Temporal development and collapse of an Arctic plant-pollinator network
© Pradal et al; licensee BioMed Central Ltd. 2009
Received: 22 May 2009
Accepted: 4 December 2009
Published: 4 December 2009
The temporal dynamics and formation of plant-pollinator networks are difficult to study as it requires detailed observations of how the networks change over time. Understanding the temporal dynamics might provide insight into sustainability and robustness of the networks and how they react to environmental changes, such as global warming. Here we study an Arctic plant-pollinator network in two consecutive years using a simple mathematical model and describe the temporal dynamics (daily assembly and disassembly of links) by random mechanisms.
We develop a mathematical model with parameters governed by the probabilities for entering, leaving and making connections in the network and demonstrate that A. The dynamics is described by very similar parameters in both years despite a strong turnover in the composition of the pollinator community and different climate conditions, B. There is a drastic change in the temporal behaviour a few days before the end of the season in both years. This change leads to the collapse of the network and does not correlate with weather parameters, C. We estimate that the number of available pollinator species is about 80 species of which 75-80% are observed in each year, D. The network does not reach an equilibrium state (as defined by our model) before the collapse set in and the season is over.
We have shown that the temporal dynamics of an Arctic plant-pollinator network can be described by a simple mathematical model and that the model allows us to draw biologically interesting conclusions. Our model makes it possible to investigate how the network topology changes with changes in parameter values and might provide means to study the effect of climate on plant-pollinator networks.
The structure of plant-pollinator networks and other ecosystems has been described through features reflecting their topology and complex organization [1–3]. Also, aspects of the network dynamics have been studied, mainly from the perspective of understanding the principles that rule the number and choice of connections made by species in the network. One prominent example is the model of preferential attachment [4, 5], where new species entering the network tend to link to species already well-connected. Such models have been fitted to empirical pollination networks, see e.g. .
One key aspect of plant-pollinator networks remains little explored, namely the temporal dynamics. The temporal dynamics describes how a network is formed and modified over time, its sustainability and robustness and the mutual activities of plants and pollinators over time. Plant-pollinator networks are susceptible to environmental changes such as global warming [7–9]; hence insight into the temporal dynamics might provide valuable information about the impact of environmental changes. This is, in particular, of relevance in the Arctic where global warming is expected to have the most severe effects . In the present study, we define temporal dynamics as the development of the number of species and links in the network for a given period of time, typically a season. Here we use season in the sense of the network's activity period, i.e. the time when both flowering plants and pollinators are present. A season begins with the observation of a first link between an animal and a plant. From this starting point the network typically experiences a phase of growth where many species enter the network and with an increase in the number of connections until a slackening of the expansion appears. A collapse of the network follows with a reduction in the numbers of the links and species previously involved in the structure. In regions where climate conditions are favourable, the network can exist all year round  with fluctuations in its size throughout the year and without a complete cessation of the activity. In regions with stringent climate conditions the length of the season can be reduced to a few weeks or months because of the presence of snow and low temperatures.
A few studies explore the temporal dynamics of pollination networks [3, 6, 11–14]. Common to most of these studies is that they highlight the strong temporal dynamics of several variables, e.g. species number, species linkage level (number of links of a species to other species), total number of links in the network, network connectance and nestedness. However, these studies are based on bi-weekly, monthly or annual collection of data and comparisons, and do not propose models to describe and explain the observations, let alone the dynamics, but mainly base analysis on descriptive statistics.
Only  discusses the day-to-day dynamics and compares the fit of a (truncated) power-law and an exponential distribution to the plant and pollinator linkage levels. They use arctic plant-pollinator data collected over two consecutive seasons near Zackenberg Meteorological Station (ZMS) in Greenland. However, neither distribution is able to explain the daily development (assembly and disassembly) of links and species in the network. Here we develop a mathematical model to account for the day-to-day dynamics of the network data in . The completeness of the data enables us to follow and describe the network development using simple mathematical tools and to provide a model that reproduces aspects of the dynamics. By comparing the two seasons we assess the stability of the model parameters over the years, the influence of weather on the development of the network and draw biological consequences of the model.
Results and Discussion
The start and end of the season
Summary of weather parameters
The collapse of the network
In contrast, we found that the disappearance of pollinators from the network is highly correlated with the disappearance of plants from the network (linear regression: r2 = 0.8 for 1996, and r2 = 0.7 for 1997). Since we only know the time span a plant is visited and not for how long it is flowering, we cannot say whether links disappear A) because plants stop flowering or B) because insects disappear for other reasons. Some of the plants that stay in the network longest also enter the network very early in the season (see Additional files 2 and 3). Based on the available data, it is therefore difficult to distinguish between A and B. However, we note that the collapse of the network follows immediately after the temperature has peaked: In 1996, the temperature rises above 15°C for several days and in 1997, the temperature rises above 20°C in a single day (Additional file 4). It is plausible that flowering of plants adapted to cold climate could be affected by high temperatures .
When interpreting the results it should be kept in mind that weather conditions to some extent could be local and heterogeneous; e.g. variations in the orientation of the ground level towards the sun could impose restrictions on the hours of sun, wind and the snow coverage of plants compared to ZMS. Heterogeneous snow distribution has previsouly been reported to influence e.g. alpine ecosystems [10, 16].
Simple distributions describe the temporal development
Our main focus is to describe the temporal dynamics of the network with simple probability distributions and compare the results between the two years. The composition of the plant community is the same for the two years; in total 31 flowering plants comprising the same species each year; see  and Additional file 2. Therefore we consider the plants as fixed and study the network from the view of the pollinators. In particular, we are interested in the following features:
• How many pollinators enter the network daily?
• How many pollinators leave the network daily?
• How many links do pollinators get when they first enter the network?
• How many links do pollinators gain daily while they are in the network?
• How many links do pollinators loose daily while they are in the network?
Arrival of insects Poisson
Links for new insects Modified geometric
Death of insects Sigmoid binomial
Addition of links Geometric
Removal of links Sigmoid binomial
When pollinators first enter the network, many have few links while fewer have many links; 57% (resp. 63%) of the pollinators enter the network with one or two links and 12% (resp. 13%) have six or more links in 1996 (resp. 1997). We used a modified geometric distribution that allows the ratio of insects with one link to insects with two or more links to be higher compared to a pure geometric distribution (see details in Methods, Analysis and modelling).
For the arrival of new pollinators and new links the fit is not improved significantly if a seasonal change in the parameters are allowed (results not shown) and there is not a visual indication of seasonal change. For the death of insects and removal of links the change in parameter happens around the same time (the parameters T and t are very similar, Table 2) and a model assuming T = t fits the data (results not shown). We also note that the change in parameter values over the season is very sharp and appears to happen within a day or two (Figure 2).
Stability in the parameters over the two years
An important observation is that the estimated model parameters are similar for the two years, which indicates stability in the network dynamics. This hypothesis was validated for the arrival and leave of pollinators and for the number of links when they enter the network (Table 2). However the evolution of the number of links (addition or removal of links for pollinators while in the network) cannot be described by joined parameters. Despite of this, the parameters are still of the same order of magnitude. (In the remaining of the paper, the 'joined model' refers to a model which uses joined parameters whenever possible.)
While the plant community is the same for the two years, one fifth of the pollinator species and two thirds of all links are only observed in one of the two years; see  and Additional file 3. It is therefore interesting to see that the fitted distributions for the evolution of the number of pollinators and links through the season are identical for the two years with similar parameters. For other networks it has already been reported that overall network properties such as connectance or nestedness are conserved over the years despite the turnover in pollinator species and links [3, 18]. Our findings reinforce this observation at the more detailed level, the level of the dynamical assembly and disassembly.
Based on the above results we implemented a computer program reproducing the arrival/leaving of pollinators and the evolution of their number of links with plants during their stay in the network (see Methods, Simulation). Addition of links is assumed to be independent of the pollinator's present linkage level. Also each link is removed with a probability independent of the linkage level. When pollinators are removed from the network they are chosen with probability inversely proportional to their linkage level. This provides a better fit to the observed networks than removing pollinators randomly and might reflect sampling properties; i.e. a pollinator with many links might have higher chance of being observed over consequetive days than a pollinator with few links.
Validation of the model
Maximum number of links
Total number of insects
Total number of interactions
Maximum number of insects
Pollinator average linkage level
Number of links when at maximum
Plants per pollinator
Consequences of the model
The model assumes a Poisson number of arriving insects per day. This is compatible with a scenario where there are N available species and each insect has the same probability to visit and be observed at the study site; λ = Np i /G i where p i is the probability in year i and G i the number of good days (λ is assumed to be the same in both years). Under this scenario we expect N1N2/N = N12 insects to be observed both years. Here N i is the number of insects in year i and N12 the number of insects observed in both years. We find that N = N1N2/N12 = 61·64/49 = 79.7 pollinators (see Methods, Data sets), hence p1 = N1/N = 76.6% of the available pollinator species are observed in 1996, while p2 = 80.3% are observed in 1997. Solving λ = Np i /G i , gives p1 = 72.9% and p2 = 79.0% which are close to 76.6% and 80.3%, respectively, but derived through the Poisson model. The per day probability is independent of the year (since λ is the same), p i /G i = λ/N = 0.03.
Prolonging length of the season
Length of the season
Maximum number of links
Total number of insects
Total number of links
Maximum number of insects
Pollinator average linkage level
The model stipulates randomness in the development of the network. The linkage level of a pollinator is described by a sum of geometric (new links) and binomial (removal of links) variables. As shown in  the linkage level is far from a power-law and closer to an exponential (or geometric) distribution. Also a truncated power-law distribution (a power law restricted to a certain range defined by a cut-off) is fitted to the data in , but the cut-off is here difficult to reconcile with biological interpretation (see  for a discussion of power-laws in biology and  for a discussion and review of cut-offs/characteristic scales in ecology).
Our study highlights various interesting points. Based on our model, we have demonstrated that the plant-pollinator network shows strong dynamic stability over the two years; i.e. the dynamic features of the network are highly conserved from one year to the next. The length of the season, temperature and other weather parameters differ (Table 1) and also the visiting pollinators are not the same. Despite this we find that the development of the network could be described by very similar parameters in 1996 and 1997. We described the number of links attached to a new insect by a modified geometric distribution. This distribution does not have the characteristic power-law shape that has been reported for other types of network and is in concordance with a previous analysis of the same data . In addition, our model has distinctive random features; e.g. the number of new species per day is independent of the number already present and the number of insects being removed from the network daily is binomial such that each insect has the same probability of being removed.
Studies about the impact of habitat loss on pollination networks highlight the existence of a habitat destruction threshold at which both plants and pollinators disappear suddenly, leading to the collapse of the network . The transition from maintenance to destruction of the network is very sharp leading to the conclusion that the network's fate might change by a slight modification of the parameter controlling the transition - almost like a phase transition. We observe a similar collapse of the network with parameters changing many folds within a few days. In our case, it could be that a sudden short rise in temperature causes plants to stop flowering . Alternatively, it could be that night frost towards the end of the season creates north-facing snow patches that persist during day time, thereby reducing the resources available for the insects and making it more difficult for them to recover and regain activity. However, to test this hypothesis we need data on frost (and snow coverage) in the study site, which are not available. Finally, we demonstrated by simulation that the arrival and departure of insects in network has not reached equilibrium when the collapse of the network appears at the end of the season and we estimated that the number of available species in the area is about 80, of which 75%-80% are observed.
Climate change might have an impact on plant-pollinator networks [7, 8, 22–25]. For example an increase in the average temperature is likely to increase the length of the season and change the conditions for the existence and maintenance of the network. The effect of temperature rise is not well-understood, though some evidence is available. In , it is argued that experimental warming does not alter the length of the flowering season, whereas  and  (a study on butterflies) provide evidence that the adult life cycle of insects is unchanged with increasing temperature. However, the availability of pollinator species or plant species might change as the length of the season changes, as well as when the species are present . To study consequences of climate changes, observations over several years would help us to relate the parameters to climatic variables; however we might still need to impose more assumptions on the model, e.g. for how long are plants flowering, to what extent do pollinators overlap with their plants etc.
Remarks on and limitations of the model
We consider our model a first step in modelling the temporal dynamics of plant-pollinator networks, a complex process of network formation. As discussed in the sections above, species-specific information about plant and pollinators is not included in the model. This information includes species identity, species abundance, known species-specific interactions and weather-related parameters, such as temperature thresholds for when species are active and when they rest. In our case, species identity (Additional files 2 and 3) and the specific plant-pollinator interactions  are available. All this information potentially influence the dynamics of the network, e.g. when and which connections are made in the network [10, 16, 29]. However, it is not straightforward how to include such elaborate information in the model and the analysis we have proposed. Our model is based on network data gathered over two years and we have demonstrated that many aspects of the dynamics can be accounted for using simple mathematical tools. With detailed observations over several years we might be able to provide models that can account for further aspects of the temporal dynamics and include more detailed describtions of plants, pollinators and their characteristics.
Our study is based on the plant-pollinator network of Zackenberg in Greenland . Data were collected in 1996 and 1997 and include two full seasons from the first observation of insects visiting flowering plants to the last observation. In 1996 the season is 43 days from June 21 to August 2, while in 1997 it is 69 days from June 17 to August 24. However, bad weather reduced the number of days of observation to 24 and 26, respectively. Days of bad weather conditions were determined on site by the scientists involved in collecting the data. In the following 'good days' refers to days of observation and 'bad days' to days where observation was not possible. The study site is 500 m × 500 m.
For each day of observation we have a matrix describing the existing links between plants and pollinators present in the network at this day (Additional files 2 and 3 show summaries of the data: plant and pollinator phenology, respectively). An animal is linked to a plant species if the animal visits the flower of the plant. A link between a plant and a pollinator is supposed to exist from the first day when the link is observed to the last day when it is observed, irrespectively whether it is observed in all the intermediate days . In total, the 1996 (1997) network includes 31 (31) plant species, 61 (64) pollinators and 286 (268) different links. The plant species are the same both years. In contrast 49 pollinators are observed both years with 79 different pollinators in total. The sampling effort was assessed in  and found to be good. Climatic data were collected from ZMS including temperature, amount of radiations (in, out and net radiation), and wind speed and direction http://www.zackenberg.dk; see Table 1. We obtained full information for the 1996 season, while in the 1997 season two days (June 24 and 25) were missing.
Analysis and modelling
We used the following:
• Each day a Poisson number Po(λ) of insects enter the network
Early in the season p(x) is α, and late β. At time T, p(T) = (α + β)/2; H controls how rapid the shift from α to β is. We assume at least one insect remains in the network each day
• When an insect enters the network, it is assigned a number of links k according to a modified geometric distribution:
– The probability of k = 1 link is r1
– The probability of k >1 links is (1 - r1)r(1 - r)k-2
• Each day, each pollinator in the network receives new links according to a geometric distribution, Geo(q)
At time t, q(t) = (γ + δ)/2, where γ is the value of q(x) early in the season and δ the value late. We assume at least one link remains in the network each day.
Parameters were estimated using maximum likelihood, separately for each year and for the two years jointly. To test the validity of our models chi-square goodness-of-fit tests were used with simulated p-values because of small expected numbers. We used the Likelihood Ratio Test (LTR) to test whether parameters for the two seasons could be assumed identical. The analysis and tests were performed in R http://www.R-project.org. Logistic regression analysis on meteorological data was done using the package nnet in R.
To simulate a network from our model we wrote a program in R based on the distributions described in the previous section. We used estimated parameters, but the program allows arbitrary parameters as input, including season length (number of good days, i.e. days of observation). A simulation goes through the following steps:
• The length of the season is fixed from the beginning
• For each day, the number of new pollinators is defined according to a Poisson law. The sum of new pollinators over the season defines the pollinator list. The date of arrival is attributed to each pollinator
• For each day of the season, a number of links (modified geometric distribution) is attributed to the pollinators entering the network on this day. Links can be added to pollinators already in the network (geometric distribution) or deleted (binomial distribution)
• For each day of the season, the number of pollinators leaving the network is determined according to a binomial law. The leaving pollinators are then chosen with a probability inversely proportional to their linkage level (i.e. probability 1/k, where k is the number of links of the pollinator).
We simulated 50 networks and compared various features to the real networks. For each feature, we calculated the Z-score Z = (X - μ)/σ, where X is the value in the real network, μ the mean of the simulations and s the standard deviation of the simulations; the distribution of simulated values are approximately normal (results not shown). To compare distributions, the function ks.boot from the package Matching in R which performs a Kolmogorov-Smirnov test with bootstrapping, was used.
Julien Dutheil is thanked for helping out with R, general assistance and discussions. CW and JMO are supported by the Danish Research Council. We are also grateful to H. and B. Elberling for use of data. This work was done while CP was visiting BiRC from August 2008 to January 2009.
- Bascompte J, Jordano P, Melian CJ, Olesen JM: The nested assembly of plant-animal mutualistic networks. PNAS. 2003, 100: 9383-9387. 10.1073/pnas.1633576100.PubMed CentralView ArticlePubMedGoogle Scholar
- Basilio AM, Medan D, Torretta JP, Bartoloni NJ: A year-long plant-pollinator network. Austral Ecol. 2006, 31: 975-983. 10.1111/j.1442-9993.2006.01666.x.View ArticleGoogle Scholar
- Petanidou T, Kallimanis AS, Tzanopoulos J, Sgardelis SP, Pantis JD: Long-term observation of a pollination network: fluctuation in species and interactions, relative invariance of network structure and implications for estimates of specialization. Ecol Lett. 2008, 11: 564-575. 10.1111/j.1461-0248.2008.01170.x.View ArticlePubMedGoogle Scholar
- Barabasi AL, Albert R: Emergence of scaling in random networks. Science. 1999, 286: 509-512. 10.1126/science.286.5439.509.View ArticlePubMedGoogle Scholar
- Williams RJ, Martinez ND: Simple rules yield complex food webs. Nature. 2000, 404: 180-183. 10.1038/35006555.View ArticlePubMedGoogle Scholar
- Olesen JM, Bascompte J, Elberling H, Jordano P: Temporal dynamics in pollination network. Ecology. 2008, 89: 1573-1582. 10.1890/07-0451.1.View ArticlePubMedGoogle Scholar
- Memmott J, Craze PG, Waser NM, Price MV: Global warming and the disruption of plant-pollinator interactions. Ecol Lett. 2007, 10: 710-717. 10.1111/j.1461-0248.2007.01061.x.View ArticlePubMedGoogle Scholar
- Devoto M, Zimmermann M, Medan D: Robustness of plant-flower visitor webs to simulated climate change. Ecol Austral. 2007, 17: 37-50.Google Scholar
- Ellebjerg SM, Tamstorf MP, Illeris L, Michelsen A, Hansen BU: Inter-annual variability and controls of plant phenology and productivity at Zackenberg. Adv Ecol Res. 2008, 40: 249-273. full_text.View ArticleGoogle Scholar
- Kudo G, Hirao AS: Habitat-specific responses in the flowering phenology and seed set of alpine plants to climate variation: implications for global-change impacts. Popul Ecol. 2006, 48: 49-58. 10.1007/s10144-005-0242-z.View ArticleGoogle Scholar
- Petanidou T: Pollination fauna of a phryganic ecosystem. Versl Tech Geg Amst. 1991, 59: 1-12.Google Scholar
- Medan D, Montaldo NH, Devoto M, Mantese A, Vasellati V, Roitman GG, Bartoloni NH: Plant-pollinator relationships at two altitudes in the Andes of Mendoza, Argentina. Arctic Antarctic Alp Res. 2002, 34: 233-243. 10.2307/1552480.View ArticleGoogle Scholar
- Medan D, Basilio AM, Devoto M, Bartoloni NJ, Torretta JP, Petanidou T: Measuring generalizations and connectance in temperate, long-lasting systems. Plant-pollinator Interactions: from Specialization to Generalization. Edited by: Waser N, Ollerton J. 2006, Chicago: University of Chicago Press, 245-259.Google Scholar
- Alarcon R, Waser NM, Ollerton J: Year-to-year variation in the topology of a plant-pollinator interaction network. Oikos. 2008, 17: 1796-1807. 10.1111/j.0030-1299.2008.16987.x.View ArticleGoogle Scholar
- Hansen BU, Sigsgaard C, Rasmussen L, Cappelen J, Mernild SH, Petersen D, Tamstorf MP, Rasch M, Hasholt B: Present-day climate at Zackenberg. Adv Ecol Res. 2008, 40: 111-149. full_text.View ArticleGoogle Scholar
- Bergman P, Molau U, Holmgren B: Micrometeorological impacts on insect activity and plant reproductive success in an alpine environment, Swedish Lapland. Arctic and Alpine Res. 1996, 28: 196-202. 10.2307/1551760.View ArticleGoogle Scholar
- Arroyo MTK, Armesto JJ, Villagran C: Plant phenological patterns in the high Andean Cordillera of central Chile. J Ecol. 1981, 69: 205-223. 10.2307/2259826.View ArticleGoogle Scholar
- Jordano P, Bascompte J, Olesen JM: Invariant properties in coevolutionary networks of plant-animal interactions. Ecol Lett. 2003, 6: 69-81. 10.1046/j.1461-0248.2003.00403.x.View ArticleGoogle Scholar
- Keller EF: Revisiting "scale-free" networks. Bio Essays. 2005, 27: 1060-1068.Google Scholar
- Montoya JM, Pimm SL, Sole RV: Ecological networks and their fragility. Nature. 2006, 442: 259-264. 10.1038/nature04927.View ArticlePubMedGoogle Scholar
- Fortuna MA, Bascompte J: Habitat loss and the structure of plant-animal mutualistic networks. Ecol Lett. 2006, 9: 281-286. 10.1111/j.1461-0248.2005.00868.x.View ArticlePubMedGoogle Scholar
- Kudo G, Ito K: Plant distribution in relation to the length of the growing season in a snow-bed in the Taisetsu Mountains, northern Japan. Vegetatio. 1992, 98: 165-174. 10.1007/BF00045554.View ArticleGoogle Scholar
- Høye TT, Post E, Meltofte H, Schmidt NM, Forchhammer MC: Rapid advancement of spring in the High Arctic. Curr Biol. 2007, 17: R450-10.1016/j.cub.2007.04.047.View ArticleGoogle Scholar
- Berg TB, Schmidt NM, Høye TT, Aastrup P, Hendrichsen DK, Forchhammer MC, Klein DR: High-arctic plant-herbivore interactions under climate influence. Adv Ecol Res. 2008, 40: 275-298. full_text.View ArticleGoogle Scholar
- Klein DR, Bruun HH, Lundgren R, Philipp M: Climate change influences on species interrelationships and distributions in high-artic Greenland. Adv Ecol Res. 2008, 40: 81-100. full_text.View ArticleGoogle Scholar
- Price MV, Waser NM: Effects of experimental warming on plant reproductive phenology in a subalpine meadow. Ecology. 1998, 74: 1262-1271.Google Scholar
- Bosch J, Kemp WP, Peterson SS: Management of Osmia Lignaria populations for almond pollination: methods to advance bee emergence. Environ Entomol. 2000, 29: 874-883. 10.1603/0046-225X-29.5.874.View ArticleGoogle Scholar
- Karlsson B, Wiklund C: Butterfly life history and temperature adaptations; dry open habitats select for increased fecundity and longetivity. J Anim Ecol. 2005, 74: 99-104. 10.1111/j.1365-2656.2004.00902.x.View ArticleGoogle Scholar
- Kudo G, Suzuki S: Relationship between flowering phenology and fruit-set of dwarf shrubs in Alpine fellfields in Nothern Japan: A comparison with a subarctic heathland in Northern Sweden. Arctic, Antarctic and Alpine Res. 2002, 34: 185-190. 10.2307/1552470.View ArticleGoogle 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.