Can snow depth be used to predict the distribution of the high Arctic aphid Acyrthosiphon svalbardicum (Hemiptera: Aphididae) on Spitsbergen?
© Ávila-Jiménez and Coulson; licensee BioMed Central Ltd. 2011
Received: 13 April 2011
Accepted: 13 October 2011
Published: 13 October 2011
The Svalbard endemic aphid Acyrthosiphon svalbardicum (Heikinheimo, 1968) is host specific to Dryas octopetala L. ssp octopetala (Rosaceae). It has been hypothesized that the aphid is present on those areas with a thin winter snow cover and which therefore clear of snow earlier in the season. This early snow clearance results in a longer growing period and allows the aphid to experience at least the minimum number of degree days required to complete its life cycle. However, this hypothesis lacked a detailed field validation. We aimed to test the relationship between the aphid distribution and time of snow clearance at landscape scale, mapping snow depth at peak of snow accumulation for the two succeeding years 2009 and 2010 and examining site occupancy and plant phenology the following summers. Additionally, the distribution range mapped by Strathdee & Bale (1995) was revisited to address possible changes in range along the coast of the fjord.
A linear relation between snow depth and timing of snow melt was found but with strong inter-annual and landscape variation. Both snow depth and plant phenology were found to affect patch occupancy. In August, the aphid, at the three life stages scored (viviparae, oviparae/males and eggs), was present most frequently in those D. octopetala patches with the most advanced plant phenology and which showed shallower snow depths in spring. However, many patches predicted to contain aphids were empty. The aphid distribution range has expanded 4.7 km towards the fjord mouth from 1995.
Snow depth alone, and hence date of snow clearance, cannot precisely define species distribution at landscape scale, as this cannot explain why are they unoccupied patches under shallow snow depths with advanced plant phenology. We nonetheless present a model Arctic system that could form the basis for long term monitoring for climate- driven species shifts.
The climate envelope paradigm (the range of climatic requirements or tolerances)  is often criticised as a simplistic approach in modelling species distributions as it systematically excludes the influence of biotic factors on species distributions [2, 3]. The inherent difficulty in precisely defining the factors determining the distribution range of a species can be particularly misleading when predicting shifts in the species range as a response to climate change . Nonetheless, the fine resolution of the factors affecting successful reproduction and dispersal abilities is essential for predicting the effects of environmental change on species distribution . For Arctic species strongly dependant on the number of degree days to successfully complete their life cycle, factors such as length of summer, temperature, timing of snow melt, snow thickness and overwintering conditions  are all likely to influence species presence or absence. Surface ice layer thickness has also been shown to have a deleterious effect on soil invertebrate populations .
Depth of snow, and hence date of snow clearance and length of the summer period, is intuitively important in determining species occurrence. For many species, especially holocyclic species, the date of snow clearance may be critical. However, few studies have investigated this point with the exception of the research by Strathdee & Bale  on the Svalbard endemic aphid Acyrthosiphon svalbardicum (Heikinheimo, 1968), and studies regarding snow bed plant communities [7–10]. Climate predictions point to an increase in winter snowfall combined with increased winter snow evaporation, leading to as yet uncertain changes in the timing of the snow melt . In the case of A. svalbardicum, knowledge of the environmental factors controlling its biology and distribution is key to better understand the abilities of the species to track its ecological niche in changing times . Predicting species distribution response to environmental changes can stand as a major challenge if complex relations are to be modelled. For instance, plant phenology can also be considered among the same environmental factors that affect aphid distribution as temperature, snow thickness [7, 13] summer length, or timing of snow melt . Hence identifying the direct effect of a single isolated environmental variable can be a challenging task. Specialist phytophagous species often show a strong phenological synchronization with host plant [15–17]. Furthermore, host plant distribution, growth and phenology can be directly influenced by temperature. For example low summer temperature limit seed setting in the high Arctic , vegetative growth and germination are positively affected by warming, and peak flowering can be brought forward by increasing temperatures .
Acyrthosiphon svalbardicum is host specific to Dryas octopetala L. ssp octopetala (Rosaceae). However, the distribution of the aphid is more restricted than the distribution of the host plant, that is, not all patches of D. octopetala harbour aphids. Notable studies regarding the environmental factors effecting aphid distribution in Kongsfjord (Svalbard archipelago) were carried out during the early 1990's [6, 20], concluding that the aphid was restricted to those areas which clear up from snow earlier in the season and where the aphid will experience the minimum number of degree days required to complete its annual life cycle . This simplistic hypothesis, the direct relation between snow depth and the aphid distribution, was rapidly entrenched in polar science, becoming a modern axiom in Arctic ecology, even though the field study was constrained to six D. octopetala patches in a 2.7 by 3.5 m south facing plot in a single location within the settlement of Ny-Ålesund . Verification regarding the representativeness of that patch on the species range is lacking, even though the likelihood of misestimating the effect of unusual events is recognized to increase with the reduction on the scale of the study . In fact, later studies in the southern location of Endalen (close to Isfjord) have showed the aphid to be more abundant than expected under deep snow profiles . Consequently, whether snow depth, and hence timing of snow clearance, has a role in determining the local limit of the distribution of the aphid yet remains unclear.
Acyrthosiphon svalbardicum is locally abundant on coastal ridges along the south coast of Kongsfjord, where its occurrence has been described as decreasing towards the entrance of the fjord and with distance from the shore , even though the host plant is ubiquitous in the area. This species is thought to have low dispersal abilities since the production of alate forms appears to be rare . This has led to the hypothesis that this species may not be able to rapidly track suitable environmental conditions but that it is restrained to those few consistently microclimatically favourable sites . Acyrthosiphon svalbardicum is a holocyclic species that overwinters as an egg, and the complex life cycle includes four morphs (fundatrix, vivipara, male and female). The fundatrix gives rise directly to both sexual morphs and a viviparous form. The viviparae produces an additional generation of sexual morphs which often fails to give rise successful offspring, but which should it succeed, will increase the number of overwintering eggs by 11 fold . The production of alate forms of the aphid has been rarely reported, and this event is thought to be related to unusually long or warm summers [6, 24]. All the life stages of the aphid occur on D. octopetala, where it is more common on shoots with seed heads or flower buds than on vegetative shoots. The eggs are often found on the underside of the leaves .
Set against this background, we aimed to: 1) assess the relationship between snow depth and the distribution of A. svalbardicum; 2) survey for changes in the local range since 1995; and 3) combine field data on ice and snow depth, timing of snow melt, patch occupancy, aphid phenology and host plant phenology to determine possible relationships between these environmental variables and the distribution of the aphid at landscape scale. Moreover, we aimed to generate a baseline distribution and environmental map available for future studies both regarding snow depth and species distribution.
Snow depth and the timing of snow clearance
The minimum soil temperature recorded reached -17.6°C (recorded on the 16th March 2010 by the logger placed in a south facing site under 30 cm of snow).
Patches occupied and aphid phenology in relation to plant phenology.
Total D. octopetala patches with the given phenology
Patches free of aphids
% Patches occupied
% Oviparae scored
% Viviparae scored
% Eggs scored
Host plant phenology
Host plant phenology was examined in August 2009 in 282 patches of D. octopetala, 52 patches still bore flowers, while a total of 230 patches showed seed heads in other developmental stages. In August 2009, 11 patches had no visible seed or flower bud. The most common phenological stage at that time of the year was stage 5 (seed head twisting) with 146 patches of D. octopetala at this stage (Table 1). Host plant phenology showed no relationship with orientation of the patches (north, south slope or top of the ridge) (H = 44.23; df = 2; p > 0.1). However, a significant relationship between snow depth and plant phenology was found (H = 359.8; df = 1; p < 0.001), and the earlier phenological stages were found on areas with deeper snow profiles in winter.
Patch occupancy and aphid distribution range
From 594 points analysed in 2009, 339 represented a patch of the D. octopetala, and 98 of those were occupied by the A. svalbardicum.
From the 24 D. octopetala patches re-visited in 2010, 2 patches occupied in 2009 were now empty and 2 of the empty ones were re-colonized in 2010. The occurrence of the aphid was confirmed at the end points surveyed by Strathdee et al. .
Distribution of the aphid on the host with relation to snow depth and host phenology
Species distribution modelling is a fast advancing field , with applications that reach from population monitoring , assessing climate driven species range shift [27, 28], conservation biology , to testing biogeographical, ecological and evolutionary hypotheses . It is often possible to predict a link between climate and species distributions , although the application of climate envelope paradigms in the modelling of species distributions have repeatedly raised the need to include a range of factors/tolerances beyond only climatic variables [2, 3]. In the case of the Svalbard endemic aphid A. svalbardicum, a high Arctic species strongly dependant on number of degree days to successfully complete its life cycle, summer temperature and length of the growing season are thought to be limiting factors governing the distribution of the species [6, 20]. Previous studies indicated snow depth as an accurate predictor of A. svalbardicum distribution , indicating that the aphid succeeds in reproduction, and thus is commonly present, in those patches with snow clearance early in the season and where the aphid would experience at least the minimum number of degree days required to complete its life cycle. In our study we compared multiple transects over a series of ridges and did observe at a larger scale the predictions made by Strathdee & Bale . Eggs were found in those patches that showed shallower snow depths in winter, although the aphid was more widely distributed across areas that showed a range of snow depths. Our study also shows a statistically significant relationship between plant phenology and site occupancy that was previously overlooked, and also between plant phenology and snow depth leading to a much more complex picture. Due to strong correlation between all measured variables, neither this study, nor the previous, can infer causality between snow depth and the aphid distribution. Although the prediction by Strathdee & Bale  fits almost perfectly along a particular transect (transects 9, 10 and 12) (Figure 4), the relationship cannot be generalized for the species range across the coastal ridges in Kongsfjord. This was also the case for the distribution of the aphid in Endalen (close to Longyearbyen in Adventdalen) along a natural snow depth transect, where Dollery et al.  found a relationship between snow depth and the aphid distribution, demonstrating that shallow snow depths favour high aphid densities but while a thick winter snow pack does not preclude aphid presence. Our survey also reveals a number of cases where the species would be predicted to occur based on the Strathdee & Bale  hypothesis but where it does not actually occur (Figure 3, 4 and additional file 2). Our results show a greater number of suitable (shallow snow) patches unoccupied by the aphid than occupied. They also show site occupancy at a range of snow depths and a strong relationship between site occupancy and the phenology of the host plant.
Specialist phytophagous species sometimes rely on plant secondary metabolites as sign of stimuli during host selection, setting and parturition . Synchronization with the host plant is often essential for neonate success after hatching [15, 17]. Phenological synchronization between phytophagous species and their host plants is a common phenomenon [15, 17, 31], as in the case of the Arctic aphid Sitobion calvulus Ossiannilson, 1958, whose life cycle appears to be closely synchronized (genetically programmed) with its host plant Salix polaris. The phenology of D. octopetala nonetheless can suffer shifts even within one season due to changes in snow accumulation . In our study the aphid was found mostly in patches showing an advanced phenology, independently from the snow depth they experienced the previous winter. However, even a combination of snow depth (distribution predictor by Strathdee & Bale ) and plant phenology (additional species predictor described in the present study) cannot predict real site occupancy at a landscape scale (Figure 3, 4 and additional file 2). Strathdee & Bale  used a probabilistic approach, which although gave an insight on environmental factors affecting species distribution, cannot be used to model accurately the actual distribution range of the species. The significance of the relationship between aphid presence and plant phenology however should be interpreted with some caution due to the skew of the plant phenotype data resulting from low sample sizes from early phenological stages.
How environmental factors affect species biology, from physiology to biogeography, is also often simplified by a pragmatically reduced knowledge on functioning and dynamic of the physical environment. Even though within the same year and aspect (Figure 1A and 1B) there is evidence of a direct relationship between snow depth and timing of snow melt, this relationship weakens when the inter-annual data are considered and disappeared completely for north facing slopes. Snow depth cannot therefore be considered as a distribution predictor on north facing slopes because of the lack of any relationship between snow depth and the timing of snow melt. Variation in the timing of snow melt can be as blunt as the difference between the 100 cm of snow that melted 12th of June 2009 and the 20 cm of snow that melted with just 2 days difference (14th June) a year later; emphasizing the need for fully replicated multi-year studies and the difficulty involved in resolving cause and effect relationships.
Our data cannot yet explain why there are un-occupied patches even at shallow snow depths and possessing a suitable phenological stage. Even if aphid presence was overlooked in very low density patches that would also imply very low density patches are occurring at putative ideal patches. Inter-annual variations in snow depth, ice thickness, length of summer, plant phenology and predation pressure, together with the limited dispersal ability of the aphid might be to a certain extent driving the dynamics of a metapopulation, in which not always all suitable patches are colonized, and colonization of suboptimal habitat can occur . Acyrthosiphon svalbardicum overwinters as an egg, which has been demonstrated to survive down to -30°C . The lowest temperature recorded in the soil (recorded under 30 cm of snow in a south facing site), is greater than the recorded lower lethal temperature for the aphid in its overwintering stage. The effects of winter ground ice cover on egg hatching success has never been measured, although it has been demonstrated to reduce soil microarthropod community populations by up to 50%, and it has been suggested to have significant implications for the usually patchy distribution of these animals . Between November 2008 to April 2009 Svalbard airport weather station recorded 49.77 mm (as snow) of precipitation and 8 winter warm spells (when air temperature rises above 0°C, which combined with subsequent sub-zero temperatures leads to the creation of winter surface-ice layers in permafrost areas) comprising a total of 22 days when the winter air temperatures rose above 0°C. In the same period the following year, 83.53 mm of precipitation were registered (half of them as rain), and 15 warm spells comprising 42 days in total . This could have had an effect on eggs survival during 2009/2010 winter. Moreover it should be recognised that the snow pack is not a homogeneous column. Ice layers within the snow profile created by winter melt events, often including rain, are frequently present. Hence the passage of melt water through the snow pack is not uniform  and a combination of snow thickness and route of percolation of liquid phase water through the snow pack will influence the development of potential surface ice layers. Such surface ice layers are known to reduce the overwintering survival of oribatid mites and Collembola  and may have an effect on the hatching success of overwintering eggs. Nonetheless, these eggs appear to be very cold hardy  although duration of exposure to low temperature or the effect of icing was not assessed. Hence it is important to be aware that snow depth may have many additional effects that are correlated. Such variations in the characteristics of the snow pack may explain why snow depth in itself fails to predict accurately the distribution of the aphid.
Freeze-thaw cycles in early spring may challenge the soil invertebrate fauna, especially in alpine or sub-arctic regions . While it is known that surface freezing events can occur throughout the year , freeze-thaw cycles following snow clearance are less likely to occur in the high Arctic. Date of snow clearance is typically mid-June. At this date there has been midnight sun for over six weeks reducing soil diurnal temperature variation. As an example, in Adventfjord close to Longyearbyen the mean air temperature in June 2011 was 4.8°C  but the soil temperature at 5 mm depth only first became positive from 7th June. By June 11th, soil temperature had attained 5°C with maximum temperatures over 13°C occurring during the month. No after melting soil freeze-thaw events were observed (Coulson SJ: Soil temperature records from Adventdalen, unpublished). Whilst we cannot exclude the possible effect of freeze-thaw cycles on egg and fundatrix survival it would appear that such events are probably not a major factor.
Snow depth can be used as a proxy for timing of snow melt in Kongsfjord for south facing slopes (Figure 1A and 1B). Moreover, a significant relationship between the presence of A. svalbardicum and winter snow depth was found, and consequently snow depth can be used as a predictor of the local distribution of the aphid. This is however just a coarse approximation. Strong inter-annual variation in snow depth (Figure 2) and correlation of snow depth with other variables such as plant phenology prevents finely predicting which patches will be occupied by the aphid. The relationships described here cannot define species distribution at landscape scale (Figure 3, 4, and additional file 2), given that they do not explain un-occupancy of patches under shallow snow depths with advance phenology (Figure 5 and 6). The dispersal ability of the species might have been underestimated in the past, since the distribution range of the species has shown an expansion of 4.7 km towards the fjord mouth from 1995. Variables able to accurately predict species distribution are difficult to define yet are important to understand, especially if projections of responses to climate change are to be modelled.
The coastal ridges of Gluudneset and Gåsebu are located in an approximately east-west orientation on the south shore of Kongsfjord, Spitsbergen (78.91°N 12.05°E) (Figure 7), a location where the aphid presence has been previously recorded . Their main ridge areas face southwards, with a shorter and steeper north facing slope and only small east and west exposures. The top of the ridges show a typical polar semi- desert community dominated by D. octopetala whereas vegetation communities on the lower sections of the ridges are composed of a typical heath community dominated by Cassiope tetragona. Mean daily air temperatures greater than 0°C are recorded from June to August and the average annual precipitation is around 371 mm yr-1 mostly in form of snow in winter . Soil summer temperatures have been showed to fluctuate between 0.6-16.1°C . Dryas octopetala, the host plant of A. svalbardicum, is a circumpolar evergreen mat forming dwarf shrub species typical from, although not restricted to, dry calcareous soils which clear up of snow early in the season [19, 45]. Dryas octopetala often does not grow higher than few centimetres over the ground, and usually occurs as individuals that can vary in size from a few centimetres to over a metre in diameter.
Snow depth and ice thickness
Snow depth was measured using a 2.65 m long avalanche probe and was determined between 2nd and 6th April 2009 along a total of 16 transects (transects 1, 01, 2, 4, 5 and 5 a, b, c and d, 9, 10, and 12 on Gluudneset ridge, and transects 03, 6, 7, 8 on the Gåsebu ridge) (Figure 8), giving 594 data points with an approximate distance of 2 m between measurement points. Surface ice thickness was also measured (along each ridge by digging through the snow and the underlying ice layer if present until reaching the soil surface. The ice layer was measured with a ruler from the soil surface). The position of measurement locations were determined using a Differential GPS (DGPS) (Leica Geosystems SR20) and DGPS-position adjusted using the differential mode of the Leica Geosystems software. The position measured by the rover was adjusted to an accuracy of +/- 1 cm. In March 2010, two transects (6 and 5, one at each ridge) were re-visited and snow depth measured at each point.
Timing of snow melt
Temperature fluctuations experienced at the soil surface are less severe under the snow pack than after snow melt due to the temperature buffering effect of the snow pack . Such a buffering effect disappears as the snow melts and the sensors begin to measure a combination of air temperature and individual response to direct incoming short wave radiation. Date of snow clearance can be defined accordingly as the date the temperature data show an increase in the amplitude of daily fluctuations (See additional file 1). Five Tinytag temperature loggers (Gemini Data Loggers (UK) Ltd.) were placed on the ground surface under different snow depths (0, 10, 30, 100 and 130 cm) in south facing locations in April 2009 and collected the following August as proxy measurement of snow melt date. In March 2010, 10 additional Tinytag temperature loggers were placed on the ground surface in five south facing and five north facing locations under 10, 20, 30, 40 and 50 cm of snow, and recovered in July the same year. In all cases the loggers were placed on the soil surface under the snowpack and recorded temperature at intervals of 15 minutes.
Aphid phenology was assessed using a binocular microscope according to the system by Strathdee et al. , fundatrix, apterous viviparae, oviparae, male and egg. When present, the shiny black eggs are clearly visible on the underside of the leave and therefore it is not considered likely that they were overlooked.
Host plant phenology
Plant phenology was recorded in a total of 282 patches in accordance with Strathdee at al. : 1 (Stamens visible), 2 (Petals senescing), 3 (Seed head visible above petals), 4 (Seed head untwisting), 5 (Seed head twisting). If no bud was visible it was scored as 0. Sampling in August, although late in the season in relation to the aphid life cycle, was chosen so as to maximize the chances of detecting a relationship between delayed host plant phenology, aphid phenology and patch occupancy.
Patch occupancy by A. svalbardicum (which patches of D. octopetala were colonized by the aphid) was recorded in early August at those points where D. octopetala was found at the exact same DGPS locations where snow depth was measured the previous spring.
Binomial linear model was carried out using R version 2.10.1  with aphid presence/absence data as dependant variable and snow depth and plant phenology as environmental variables.
Mann-Whitney Rank-sum, Kruskall-Wallis tests and linear regression were performed using the statistical package included in SigmaPlot 11.0 (Systat Software Inc). Mann-Whitney Rank-sum test for non-parametric data was applied to test for relationships between snow depth and host plant phenology and between snow depth, plant phenology, and aphid life stage. Kruskall-Wallis one way analysis of variance on ranks for non- parametric data was applied 1) on the differences of snow depths covering different plant phenological stages; 2) to test for differences in the distribution of host plant phenological stages in relation to the orientation (south, north or top of the ridge); and 3) to the snow depth data grouped by plant phenological categories for occupied and non-occupied patches, with a post-hoc pairwise differences test performed implementing Dunn's method. Relationship between snow depth and timing of snow melt were tested using linear regression.
ARCFAC V Program partially funded this research through the project Snow depth and aphid distribution in Kongsfjord: Monitoring the effects of climate change on Acyrthosiphon svalbardicum (Insecta: Hemiptera) ARCFAC-026129-2008-49. UNIS internally funded the project. Jago Wallenschus, Peter Convey, Chester Sands, Lorna Little and Hanna- Kaisa Lakka for assisting in the field. Three anonymous reviewers for their enriching comments which improved the manuscript.
- Berry PM, Dawson TP, Harrison PA, Pearson RG: Modelling potential impacts of climate change on the bioclimatic envelope of species in Britain and Ireland. Global Ecology and Biogeography. 2002, 11 (6): 453-462. 10.1111/j.1466-8238.2002.00304.x.View ArticleGoogle Scholar
- Davies AJ, Jenkinson LS, Lawton JH, Shorrocks B, Wood S: Making mistakes when predicting shifts in species range in response to global warming Nature. Nature. 1998, 391 (6669): 783-786. 10.1038/35842.View ArticleGoogle Scholar
- Pearson RG, Dawson TP: Predicting the impacts of climate change on the distribution of species: are bioclimate envelope models useful?. Global Ecology and Biogeography. 2003, 12 (5): 361-371. 10.1046/j.1466-822X.2003.00042.x.View ArticleGoogle Scholar
- Strathdee AT, Bale JS: Life on the edge: Insect ecology in Arctic environments. Ann Rev Entomol. 1998, 43: 85-106. 10.1146/annurev.ento.43.1.85.View ArticleGoogle Scholar
- Coulson SJ, Leinaas HP, Ims RA, Søvik G: Experimental manipulation of the winter surface ice layer: the effects on a High Arctic soil microarthropod community. Ecography. 2000, 23: 299-306. 10.1111/j.1600-0587.2000.tb00285.x.View ArticleGoogle Scholar
- Strathdee AT, Bale JS: Factors limiting the distribution of Acyrthosiphon svalbardicum (Hemiptera: Aphididae) on Spitsbergen. Polar Biology. 1995, 15: 375-380.View ArticleGoogle Scholar
- Björk RG, Molau U: Ecology of Alpine Snowbeds and the Impact of Global Change. Arctic, Antarctic, and Alpine Research. 2007, 39 (1): 34-43. 10.1657/1523-0430(2007)39[34:EOASAT]2.0.CO;2.View ArticleGoogle Scholar
- Galen C, Stanton ML: Consequences of Emergence Phenology for Reproductive Success in Ranunculus adoneus (Ranunculaceae). American Journal of Botany. 1991, 78 (7): 978-988. 10.2307/2445177.View ArticleGoogle Scholar
- Galen C, Stanton ML: Responses of Snowbed Plant Species to Changes in Growing-Season Length. Ecology. 1995, 76 (5): 1546-1557. 10.2307/1938156.View ArticleGoogle Scholar
- Huelber KM, Gottfried M, Pauli H, Reiter K, Winkler M, Grabherr G: Phenological Responses of Snowbed Species to Snow Removal Dates in the Central Alps: Implications for Climate Warming. Arctic, Antarctic, and Alpine Research. 2006, 38 (1): 99-103. 10.1657/1523-0430(2006)038[0099:PROSST]2.0.CO;2.View ArticleGoogle Scholar
- Christensen JH, Hewitson B, Busuioc A, Chen A, Gao X, Held I, Jones R, Kolli RK, Kwon W-T, Laprise R, Magaña V, Mearns L, Menédez CG, Räisänen J, Rinke A, Sarr A, Whetton P: Regional climate projections. S. Climate change 2007. The physical science basis. Contribution of Working Group I to the fourth assessment report of the Intergovernmental Panel on Climate Change. Edited by: Solomon et al. 2007, Cambridge: Cambridge University Press, 847-940.Google Scholar
- Hullé M, Bonhomme J, Maurice D, Simon JC: Is the life cycle of high arctic aphids adapted to climate change?. Polar Biology. 2008, 31 (9): 1037-1042. 10.1007/s00300-008-0442-z.View ArticleGoogle Scholar
- Beck PSA, Kalmbach E, Joly D, Stien A, Nielsen L: Modelling local distribution of an Arctic dwarf shrub indicates an important role for remote sensing of snow cover. Remote Sensing of Environment. 2005, 98: 110-121. 10.1016/j.rse.2005.07.002.View ArticleGoogle Scholar
- Thórhallsdóttir TE: Flowering phenology in the central highland of Iceland and implications for climatic warming in the Arctic. Oecologia. 1998, 114 (1): 43-49. 10.1007/s004420050418.View ArticleGoogle Scholar
- Yukawa J: Synchronization of gallers with host plant phenology. Population Ecology. 2000, 42 (2): 105-113. 10.1007/PL00011989.View ArticleGoogle Scholar
- Guillespie M, Hodkinson ID, Cooper EJ, Bird JM, Jónsdóttir IS: Life history and host- plant relationships of the rare endemic Arctic aphid Acyrtosiphon calvulus in a changing environment. Entomologia Experimentalis et Applicata. 2007, 123 (3): 229-237. 10.1111/j.1570-7458.2007.00547.x.View ArticleGoogle Scholar
- Miñarro M, Dapena E: Resistance of Apple Cultivars to Dysaphis plantaginea (Hemiptera: Aphididae): Role of Tree Phenology in Infestation Avoidance. Environmental Entomology. 2009, 36 (5): 1206-1211.View ArticleGoogle Scholar
- Wookey PA, Parsons AN, Welker JM, Potter JA, Callaghan TV, Lee JA, Press MC: Comparative responses of phenology and reproductive development to simulated environmental change in sub- arctic and high arctic plants. Oikos. 1993, 67: 490-502. 10.2307/3545361.View ArticleGoogle Scholar
- Wookey PA, Robinson CH, Parsons AN, Welker JM, Press MC, Callaghan TV, Lee JA: Environmental constraints on the growth, photosynthesis and reproductive development of Dryas octopetala at a high Arctic polar semi- desert, Svalbard. Oecologia. 1995, 102: 478-489. 10.1007/BF00341360.View ArticleGoogle Scholar
- Strathdee AT, Bale JS, Hodkinson ID, Block WC, Webb NR, Coulson SJ: Identification of three previously unknown morphs of Acyrtosiphon svalbardicum Heikinheimo (Hemiptera: Aphididae) on Spitsbergen. Ent scand. 1993, 24: 43-47.View ArticleGoogle Scholar
- Blackburn TM, Gaston KJ: Scale in macroecology. Global Ecology and Biogeography. 2002, 11 (3): 185-189. 10.1046/j.1466-822X.2002.00290.x.View ArticleGoogle Scholar
- Dollery R, Hodkinson ID, Jónsdóttir IS: Impact of warming and timing of snow melt on soil microarthropod assemblages associated with Dryas- dominated plant communities on Svalbard. Ecography. 2006, 29: 111-119. 10.1111/j.2006.0906-7590.04366.x.View ArticleGoogle Scholar
- Hodkinson ID, Coulson SJ, Bird J, Webb NR: Discovery of the rare alate morph of Acyrthosiphon svalbardicum Heikinheimo (Homoptera; Aphididae): description and implications for species ecology. Norw J Entomol. 2002, 49: 77-80.Google Scholar
- Simon JC, Bonhomme J, Blackman RL, Hullé M: Winged morph of the high Arctic aphid Acyrtosiphon svalbardicum (Hemiptera: aphididae): abundance, reproductive status, and ecological significance. Can Ent. 2008, 140: 385-387. 10.4039/n08-002.View ArticleGoogle Scholar
- Guisan A, Thuiller W: Predicting species distribution: offering more than simple habitat models. Ecology Letters. 2005, 8 (9): 993-1009. 10.1111/j.1461-0248.2005.00792.x.View ArticleGoogle Scholar
- Shaffer HB, Fisher RN, Davidson C: The role of natural history collections in documenting species declines. Trends in Ecology & Evolution. 1998, 13 (1): 27-30. 10.1016/S0169-5347(97)01177-4.View ArticleGoogle Scholar
- Walther G.-R, Post E, Convey P, Annette M, Parmesan C, Beebee TJC, Fromentin J-M, Hoegh-Guldberg O, Bairlein F: Ecological responses to recent climate change. Nature. 2002, 416 (6879): 389-395. 10.1038/416389a.View ArticlePubMedGoogle Scholar
- Thomas CD, Cameron A, Green RE, Bakkenes M, Beaumont LJ, Collingham YC, Barend FNE, Ferreira de Siqueira M, Grainger A, Hannah L, Hughes L, Huntley B, van Jaarsveld AS, Midgley GF, Miles L, Ortega-Huerta MA, Townsted Peterson A, Philips OL, Williams SE: Extinction risk from climate change. Nature. 2004, 427: 145-147. 10.1038/nature02121.View ArticlePubMedGoogle Scholar
- Rushton SP, Ormerod SJ, Kerby G: New Paradigms for Modelling Species Distributions?. Journal of Applied Ecology. 2004, 41 (2): 193-200. 10.1111/j.0021-8901.2004.00903.x.View ArticleGoogle Scholar
- Leathwick JR: Are New Zealand's Nothofagus species in equilibrium with their environment?. J Veg Sci. 1998, 9: 719-732. 10.2307/3237290.View ArticleGoogle Scholar
- Hodkinson ID: Species response to global environmental change or why ecophysiological models are important: a reply to Davis et al. Journal of Animal Ecology. 1999, 68: 1259-1262. 10.1046/j.1365-2656.1999.00372.x.View ArticleGoogle Scholar
- Powell G, Tosh CR, Hardie J: Host plant selection by aphids: Behavioural, Evolutionary and Applied perspectives. Annual Review of Entomology. 1996, 51: 309-330.View ArticleGoogle Scholar
- Walker MD, Walker DA, Welker JM, Arft AM, Bardsley T, Brooks PD, Fahnestock JT, Jones MH, Losleben M, Parsons AN, Seastedt TR, Turner PL: Long-term experimental manipulation of winter snow regime and summer temperature in arctic and alpine tundra. Hydrological Processes. 1999, 13 (14-15): 2315-2330. 10.1002/(SICI)1099-1085(199910)13:14/15<2315::AID-HYP888>3.0.CO;2-A.View ArticleGoogle Scholar
- Hanski I: Metapopulation dynamics. Nature. 1998, 396 (6706): 41-49. 10.1038/23876.View ArticleGoogle Scholar
- Strathdee AT, Howling GG, Bale JS: Cold hardiness of overwintering aphid eggs. Journal of Insect Physiology. 1995, 41 (8): 653-657. 10.1016/0022-1910(95)00029-T.View ArticleGoogle Scholar
- Norwegian Meteorological Institute: Accessed 4th May 2011, [http://www.met.no]
- Bøggild CE, Forsberg R, Reeh N: Meltwater retention in a transect across the Greenland ice sheet. Annals of Glaciology. 2005, 40 (1): 169-173. 10.3189/172756405781813546.View ArticleGoogle Scholar
- Konestabo HS, Michelsen A, Holmstrup M: Response of springtails and mite populations to prolonged periods of soil freeze- thaw cycles in a sub- arctic ecosystems. Applied Soil Ecology. 2007, 36: 136-146. 10.1016/j.apsoil.2007.01.003.View ArticleGoogle Scholar
- Coulson SJ, Hodkinson ID, Block W, Webb NR, Worland MR: Low summer temperatures: a potential mortality factor for high Arctic soil micro-arthropods?. Journal of Insect Physiology. 1995, 41: 783-792. 10.1016/0022-1910(95)00027-R.View ArticleGoogle Scholar
- eKlima. Norwegian Meteorological Institute. Accessed 4th September 2011, [http://eKlima.met.no]
- Hodkinson ID, Bird J: Host-specific insect herbivores as sensors of climate change in arctic and Alpine environments. Arctic and Alpine Research. 1998, 30 (1): 78-83. 10.2307/1551747.View ArticleGoogle Scholar
- Hodkinson ID, Webb NR, Bale JS, Block W, Coulson SJ, Strathdee AT: Global change and Arctic ecosystems: conclusions and predictions from experiments with terrestrial invertebrates on Spitsbergen. Arctic and Alpine Research. 1998, 30 (3): 306-313. 10.2307/1551978.View ArticleGoogle Scholar
- Strathdee AT, Bale JS, Strathdee FC, Block WC, Coulson SJ, Webb NR, Hodkinson ID: Climatic severity and the response to temperature elevation of Arctic aphids. Global Change Biology. 1995, 1: 23-28. 10.1111/j.1365-2486.1995.tb00003.x.View ArticleGoogle Scholar
- Hodkinson ID, Coulson SJ, Webb NR: Invertebrate community assembly along proglacial chronosequences in the high Arctic. Journal of Animal Ecology. 2004, 73: 556-568. 10.1111/j.0021-8790.2004.00829.x.View ArticleGoogle Scholar
- Rønning OI: Features of the Ecology of Some Arctic Svalbard (Spitsbergen) Plant Communities. Arctic and Alpine Research. 1969, 1 (1): 29-44. 10.2307/1550358.View ArticleGoogle Scholar
- Strathdee AT, Bale JS, Block WC, Coulson SJ, Hodkinson ID, Webb NR: Effects of temperature elevation on a field population of Acyrtosiphon svalbardicum (Hemiptera: Aphididae) on Spitsbergen. Oecologia. 1993, 96: 457-465. 10.1007/BF00320502.View ArticleGoogle Scholar
- R Development Core Team: R: A language and environment for statistical computing. R Foundation for Statistical Computing. 2009, Vienna, AustriaGoogle 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.