The influence of weather conditions on the activity of high-arctic arthropods inferred from long-term observations
© Høye and Forchhammer; licensee BioMed Central Ltd. 2008
Received: 31 January 2008
Accepted: 02 May 2008
Published: 02 May 2008
Climate change is particularly pronounced in the High Arctic and a better understanding of the repercussions on ecological processes like herbivory, predation and pollination is needed. Arthropods play an important role in the high-arctic ecosystem and this role is determined by their density and activity. However, density and activity may be sensitive to separate components of climate. Earlier emergence due to advanced timing of snowmelt following climate change may expose adult arthropods to unchanged temperatures but higher levels of radiation. The capture rate of arthropods in passive open traps like pitfall trap integrates density and activity and, therefore, serves as a proxy of the magnitude of such arthropod-related ecological processes. We used arthropod pitfall trapping data and weather data from 10 seasons in high-arctic Greenland to identify climatic effects on the activity pattern of nine arthropod taxa.
We were able to statistically separate the variation in capture rates into a non-linear component of capture date (density) and a linear component of weather (activity). The non-linear proxy of density always accounted for more of the variation than the linear component of weather. After accounting for the seasonal phenological development, the most important weather variable influencing the capture rate of flying arthropods was temperature, while surface-dwelling species were principally influenced by solar radiation.
Consistent with previous findings, air temperature best explained variation in the activity level of flying insects. An advancement of the phenology in this group due to earlier snowmelt will make individuals appear earlier in the season, but parallel temperature increases could mean that individuals are exposed to similar temperatures. Hence, the effect of climatic changes on the activity pattern in this group may be unchanged. In contrast, we found that solar radiation is a better proxy of activity levels than air temperature in surface-dwelling arthropods. An advancement of the phenology may expose surface-dwelling arthropods to higher levels of solar radiation, which suggest that their locomotory performance is enhanced and their contribution to ecological processes is increased.
The earth is undergoing climatic changes and the model predictions for the high-arctic climate are particularly dramatic . This underscores the urgency for achieving a better understanding of how climate affects ecological processes like herbivory, predation and pollination . The major controls of such processes are the population size of the functional group (e.g. herbivores, predators or pollinators) and the activity level of the individuals. A particularly strong response to climatic conditions can be expected in high-arctic arthropods since they are poikilothermic and because the abiotic conditions of the High Arctic are generally close to the minimum requirements for locomotion . The seasonal emergence of adult arctic arthropods (i.e. their phenological development) is determined mainly by the timing of snowmelt [4–7], whereas their activity pattern is a function of body temperature. The microclimatic conditions that determine body temperature in arthropods is an integration of a suite of weather variables (e.g. radiation, precipitation and wind speed) in addition to temperature [8–11]. Hence, the magnitude of arthropod-related ecological processes in the High Arctic depends on several potentially uncorrelated components of climate. We have previously shown that high-arctic arthropods have advanced their emergence phenology in recent years in response to earlier timing of snowmelt . This shift in phenology may expose arthropods to a time of the season with different weather conditions and may therefore affect their level of activity.
Capture rates from passive traps like pitfall traps integrate arthropod density and activity levels  and therefore provides a proxy for the magnitude of the ecological processes in which the organisms are involved. For instance, the capture rate of wolf spiders is related to their level of activity as well as their density. As such it is a correlate of the predation pressure by wolf spiders. Likewise, the capture rate of pollinators provides a proxy for the magnitude of pollination. The effect of weather on activity patterns is immediate whereas the physiological processes controlling density involves emergence and mortality. Hence, the time scale of fluctuations in activity is likely to be shorter than the time scale of fluctuations in density. Previous work from seasonal environments has shown that these two factors can be separated statistically . Here, we extend this technique by applying generalized additive models . Using this modelling approach, it is possible to separate the variation in capture rate due to a component which varies over a short time scale (activity) from variation due to a component which varies over a longer time scale (density) without a priori assumptions about its exact shape. We aim at elucidating broader patterns in the relation between climate and arthropod ecology, and therefore we take a higher taxon approach. We addressed the following questions: Which meteorological variable most strongly influences the locomotory performance of different taxa of flying and surface-dwelling arthropods in a high-arctic location? What is the relative influence of weather parameters and capture date on capture numbers of these taxa?
Study area and data
Arthropods were monitored during 10 consecutive years (1996–2005) with samples from one window trap plot (plot 1) and six pitfall trap plots (plot 2–7) collected weekly during June, July and August. All plots were operated during the period 1996–2005 except for plots 6 and 7 which were operated during the periods 1996–1998 and 1999–2005, respectively. Each pitfall trap plot (10 × 20 m2) consisted of eight pitfall traps and each window trap plot consisted of two window traps . The window trap plot was located next to a pond. The pitfall traps were yellow plastic cups 10 cm in diameter and each window trap consisted of two chambers and the two traps were placed perpendicular to each other [see  for details]. Trapping started in June once the snow at each trap had melted. All specimens caught in all years were sorted to the family level except for collembolans and mites who were only counted. Here, we focus on: Chironomidae (Diptera), Muscidae (Diptera), Sciaridae (Diptera), Ichneumonidae (Hymenoptera), Nymphalidae (Lepidoptera), Lycosidae (Aranea), Linyphiidae (Aranea), collembolans and mites. The family Nymphalidae is represented by two species (Boloria chariclea, (Schneider) and B. polaris, (Boisduval)), the Lycosidae by one species (Pardosa glacialis, (Thorell)) and the Linyphiidae by five species (Collinsia thulensis, (Jackson), Hilaira vexatrix, (O. P.-Cambridge), Erigone arctica (White), Erigone psychrophila, Thorell and Mecynargus borealis, (Jackson)) . The data set included a total of 531,036 individuals, each belonging to one of the nine taxa. This corresponds to 93.6% of all arthropod specimen caught during the study period. Traps were occasionally flooded, emptied by arctic foxes (Alopex lagopus, (Linneaus)), or trampled by muskoxen (Ovibos moschatus, (Zimmermann)), so the capture numbers in each plot were transformed to individuals caught per trap per day for each trapping period .
Where Y i are the log10-transformed number of individuals per trap per day, β0 is the intercept and trapping date (DAY) is modelled by a spline function s(·). The term WEATHER refers to one of the four different weather variables (TDD, SR, PREC and WIND) and was modelled as a linear continuous predictor and ε i is the error term. The weather variables were calculated as daily averages for each trapping period. We repeated the model estimation for each species group (n = 9) in each plot (n = 7 for flying arthropods, n = 6 for surface-dwelling arthropods, because this group was excluded from the window trap plot) in each year (n = 10) with each of the weather parameters (n = 4) as linear predictors and date as a spline function (with df = 4). This resulted in a total of 500 sets of four competing models. Some taxa were only caught in small numbers in some plots in some years, and we therefore restricted our analyses to years and plots where at least 100 individuals were caught except for the Nymphalidae and the Ichneumonidae where the limit was set to 50 due to the larger size of individuals in these families. This reduced the number of sets of models to 385. We had an average of 11.5 trapping periods within each season. Although it is likely that the capture numbers were affected by several weather variables, we included only one in each model to reduce the risk of over-parameterization. Since all models with the same response variable had the same number of degrees of freedom and the same null deviance, the lowest residual deviance indicated the best fit to data. In this way, we identified for each species group for each plot for each year the weather variable that was best able to explain variation in capture numbers. A previous study  separated seasonal patterns of abundance from activity by first fitting a Gaussian curve to the data and using weather variables to explain variation in the residuals. In addition to fitting GAM models we adopted this approach. Hence, we fitted Gaussian curves to the entire set of 500 time series and identified the weather variable that best explained variation in the residuals from linear regression models.
Statistical analysis of capture data in relation to weather patterns
Weather variable in best model
Number of models
Average explained deviance in %
Our analyses demonstrated that the temporal variation in capture rates of high-arctic arthropods from pitfall and window traps consisted of two components indicating the effects of density and activity. By using generalized additive modelling to separate variation in capture rates, we found one component that was unimodally related to capture date and one that was a function of short-term weather fluctuations through the season. Air temperature was the most important weather variable followed by solar radiation. The influence of temperature was most pronounced on the taxa of flying arthropods like Lepidoptera and Diptera, which is consistent with previous findings . The surface-dwelling taxa like the spiders (Linyphiidae and Lycosidae) responded most strongly to variation in solar radiation. Solar radiation is probably a better predictor of near-surface temperatures than ambient air temperature two metres above the ground. Although, capture rates may be influenced by factors other than density and activity, such as seasonally-determined changes in behaviour related specifically to reproduction, our flexible non-linear functions would have included such effects in the function of capture date. When we reran the entire analysis using Gaussian curves to model the seasonal development the same weather variables were identified as the most important in eight of twelve cases. The deviation between the model approaches could arise if the seasonal development does not follow a Gaussian curve. Alternatively, the spline function of the GAM could have captured some of the variation due to variation in activity levels, but with the degrees of freedom in the spline (i.e. the flexibility of the curve) set at a constant and low level this is rather unlikely. The good correspondence between the two methods lends further support to the conclusion that the identified weather variables are indeed affecting the activity level of the individual species groups.
Species-specific responses may be hidden by pooling multiple species in larger taxonomic units. However, we found no bimodality in the seasonal capture of the taxa included in this study which could suggest that individual species within a taxon occurred during different parts of the season. Also, since flying and surface-dwelling arthropods responded to different components of climate, general differences in the environmental conditions experienced by these two groups may be more important than behavioural differences between species within taxa. Clearly, this higher taxon approach precludes a detailed understanding of the species-environment interactions. However, since our aim here is to elucidate broader patterns in the relation between climate and arthropod ecology, we refer to these interactions simply as arthropod-related ecological processes. It is possible that sampling with shorter time intervals would have identified a stronger influence of weather variables on activity of the arthropods. However, without increasing the number of traps this would also have reduced the sample size from each sampling and thereby reduced the statistical power of the analyses. Indeed we are confident that by using data from ten seasons we identified the most important weather variable for the activity level of each taxon.
Our non-linear statistical modelling approach also made it possible to quantify the relative importance of density and activity. It has long been acknowledged that passive traps like pitfall traps provide biased estimates of density, because differences in activity levels between trapping periods are difficult to control for . It has stimulated the development of alternative sampling techniques to provide unbiased estimates of population density [23, 24]. However, in some cases no alternative sampling technique is feasible. Indeed, some organisms are simply too small and fragile for applying mark recapture techniques or too mobile, cryptic or non-randomly distributed to allow for counting using a transect approach. In these cases, parallel data on dominating meteorological variables make it possible to control for the effect of differences in locomotory performance among trapping periods. Also, capture rates in passive traps may actually provide a proxy for the rate of ecological processes involving arthropods.
We were able to statistically separate the variation in capture rates of high-arctic arthropods into a non-linear component of capture date (density) and a linear component of weather (activity). Consistent with previous findings, air temperature best explained variation in the activity level of flying insects. An advancement of the phenology in this group due to earlier snowmelt will make individuals appear earlier in the season, but parallel temperature increases could mean that individuals are exposed to similar temperatures. Hence, the effect of climatic changes on the activity pattern in this group may be unchanged. In contrast, we found that solar radiation is a better proxy of activity levels than air temperature in surface-dwelling arthropods. An advancement of the phenology may expose surface-dwelling arthropods to higher levels of solar radiation, which suggest that their locomotory performance is enhanced and their contribution to ecological processes is increased. Therefore, the integrated effect of climatic changes on ecological processes depends on the change in both temperature and timing of snowmelt and may vary between processes involving surface-dwelling and flying arthropods.
We acknowledge the Zackenberg Basic monitoring programme at the Department of Arctic Environment, National Environmental Research Institute, University of Aarhus, Denmark, for providing access to ecosystem monitoring data and the Danish Polar Center for logistic support at the research station at Zackenberg. Niels M. Schmidt is thanked for fruitful discussions on the analysis and comments on an earlier draft by Nicholas Tyler are much appreciated.
- Kattsov VM, Källén E, Cattle H, Christensen J, Drange H, Hanssen-Bauer I, Jóhannesen T, Karol I, Räisänen J, Svensson G, Vavulin S: Future climate change: modeling and scenarios for the Arctic. Arctic Climate Impact Assessment. 2005, Cambridge: Cambridge University Press, 99-150.
- Walther GR, Post E, Convey P, Menzel A, Parmesan C, Beebee TJC, Fromentin JM, Hoegh-Guldberg O, Bairlein F: Ecological responses to recent climate change. Nature. 2002, 416: 389-395. 10.1038/416389a.View ArticlePubMed
- Strathdee AT, Bale JS: Life on the edge: insect ecology in arctic environments. Annu Rev Entomol. 1998, 43: 85-106. 10.1146/annurev.ento.43.1.85.View ArticlePubMed
- Danks HV: Some effects of photoperiod, temperature, and food on emergence in three species of Chironomidae (Diptera). Can Entomol. 1978, 110: 289-300.View Article
- MacLean SFJr, Pitelka FA: Seasonal patterns of abundance of tundra arthropods near Barrow. Arctic. 1971, 24: 19-40.View Article
- Corbet PS, Danks HV: Seasonal emergence and activity of mosquitoes (Diptera: Culicidae) in a High-Arctic locality. Can Entomol. 1973, 105: 837-872.View Article
- Høye TT, Forchhammer MC: Phenology of high-arctic arthropods: effects of climate on spatial, seasonal and inter-annual variation. Adv Ecol Res. 2008, 40: 299-324. 10.1016/S0065-2504(07)00013-X.View Article
- Peng RK, Fletcher CR, Sutton SL: The effect of microclimate on flying dipterans. Int J Biometeorol. 1992, 36: 69-76. 10.1007/BF01208916.View Article
- Edde PA, Phillips TW, Nansen C, Payton ME: Flight activity of the lesser grain borer, Rhyzopertha dominica F. (Coleoptera: Bostrichidae), in relation to weather. Environ Entomol. 2006, 35: 616-624.View Article
- Waringer JA: Phenology and the influence of meteorological parameters on the catching success of light-trapping for Trichoptera. Freshwater Biol. 1991, 25: 307-319. 10.1111/j.1365-2427.1991.tb00493.x.View Article
- Service MW: Effects of wind on the behaviour and distribution of mosquitoes and blackflies. Int J Biometeorol. 1980, 24: 347-353. 10.1007/BF02250577.View Article
- Høye TT, Post E, Meltofte H, Schmidt NM, Forchhammer MC: Rapid advancement of spring in the High Arctic. Curr Biol. 2007, 17: R449-R451. 10.1016/j.cub.2007.04.047.View ArticlePubMed
- Southwood TRE, Henderson PA: Ecological methods. 2000, Oxford: Blackwell Science
- Briers RA, Cariss HM, Gee JHR: Flight activity of adult stoneflies in relation to weather. Ecol Entomol. 2003, 28: 31-40. 10.1046/j.1365-2311.2003.00480.x.View Article
- Hastie TJ, Tibshirani RJ: Generalized additive models. 1990, Boca Raton: Chapman & Hall/CRC
- Meltofte H, Thing H: Zackenberg Ecological Research Operations 1st Annual Report 1995. 1996, Copenhagen; Danish Polar Center, Ministry of Research & Information Technology
- Rasch M, Caning K: Zackenberg Ecological Research Operations, 10th Annual Report, 2004. 2005, Copenhagen: Danish Polar Center, Ministry of Science, Technology and Innovation
- Jónsson E, Gardarsson A, Gíslason G: A new window trap used in the assessment of the flight periods of Chironomidae and Simuliidae (Diptera). Freshwater Biol. 1986, 16: 711-719. 10.1111/j.1365-2427.1986.tb01012.x.View Article
- Meltofte H, Berg TB: BioBasis – Conceptual design and sampling procedures of the biological programme of Zackenberg Basic. 2006, Roskilde: National Environmental Research Institute, Department of Arctic Environment
- Larsen S: Diversity and composition of the spider fauna in different vegetation types at Zackenberg. Zackenberg Ecological Research Operations, 6th Annual Report, 2000. Edited by: Caning K, Rasch M. 2001, Copenhagen: Danish Polar Center, Ministry of Research and Information Technology, 54-56.
- Venables WN, Ripley BD: Modern applied statistics with S. 2002, New York: SpringerView Article
- Topping CJ, Sunderland KD: Limitations to the use of pitfall traps in ecological studies exemplified by a study of spiders in a field of winter wheat. J Appl Ecol. 1992, 29: 485-491. 10.2307/2404516.View Article
- Perner J, Schueler S: Estimating the density of ground-dwelling arthropods with pitfall traps using a nested-cross array. J Anim Ecol. 2004, 73: 469-477. 10.1111/j.0021-8790.2004.00821.x.View Article
- Thomas CFG, Parkinson L, Marshall EJP: Isolating the components of activity-density for the carabid beetle Pterostichus melanarius in farmland. Oecologia. 1998, 116: 103-112. 10.1007/s004420050568.View Article
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.