Climate and spatio-temporal variation in the population dynamics of a long distance migrant, the white stork

1. A central question in ecology is to separate the relative contribution of density dependence and stochastic influences to annual fluctuations in population size. Here we estimate the deterministic and stochastic components of the dynamics of different European populations of white stork Ciconia ciconia. We then examined whether annual changes in population size was related to the climate during the breeding period (the 'tap hypothesis' sensu Saether, Sutherland & Engen (2004, Advances in Ecological Research, 35, 185 209) or during the nonbreeding period, especially in the winter areas in Africa (the 'tube hypothesis'). 2. A general characteristic of the population dynamics of this long-distance migrant is small environmental stochasticity and strong density regulation around the carrying capacity with short return times to equilibrium. 3. Annual changes in the size of the eastern European populations were correlated by rainfall in the wintering areas in Africa as well as local weather in the breeding areas just before arrival and in the later part of the breeding season and regional climate variation (North Atlantic Oscillation). This indicates that weather influences the population fluctuations of white storks through losses of sexually mature individuals as well as through an effect on the number of individuals that manages to establish themselves in the breeding population. Thus, both the tap and tube hypothesis explains climate influences on white stork population dynamics. 4. The spatial scale of environmental noise after accounting for the local dynamics was 67 km, suggesting that the strong density dependence reduces the synchronizing effects of climate variation on the population dynamics of white stork. 5. Several climate variables reduced the synchrony of the residual variation in population size after accounting for density dependence and demographic stochasticity, indicating that these climate variables had a synchronizing effect on the population fluctuations. In contrast, other climatic variables acted as desynchronizing agents. 6. Our results illustrate that evaluating the effects of common environmental variables on the spatio-temporal variation in population dynamics require estimates and modelling of their influence on the local dynamics.     

The extended Moran effect and large-scale synchronous fluctuations in the size of great tit and blue tit populations

1. Synchronous fluctuations of geographically separated populations are in general explained by the Moran effect, i.e. a common influence on the local population dynamics of environmental variables that are correlated in space. Empirical support for such a Moran effect has been difficult to provide, mainly due to problems separating out effects of local population dynamics, demographic stochasticity and dispersal that also influence the spatial scaling of population processes. Here we generalize the Moran effect by decomposing the spatial autocorrelation function for fluctuations in the size of great tit Parus major and blue tit Cyanistes caeruleus populations into components due to spatial correlations in the environmental noise, local differences in the strength of density regulation and the effects of demographic stochasticity. 2. Differences between localities in the strength of density dependence and nonlinearity in the density regulation had a small effect on population synchrony, whereas demographic stochasticity reduced the effects of the spatial correlation in environmental noise on the spatial correlations in population size by 21.7% and 23.3% in the great tit and blue tit, respectively. 3. Different environmental variables, such as beech mast and climate, induce a common environmental forcing on the dynamics of central European great and blue tit populations. This generates synchronous fluctuations in the size of populations located several hundred kilometres apart. 4. Although these environmental variables were autocorrelated over large areas, their contribution to the spatial synchrony in the population fluctuations differed, dependent on the spatial scaling of their effects on the local population dynamics. We also demonstrate that this effect can lead to the paradoxical result that a common environmental variable can induce spatial desynchronization of the population fluctuations. 5. This demonstrates that a proper understanding of the ecological consequences of environmental changes, especially those that occur simultaneously over large areas, will require information about the spatial scaling of their effects on local population dynamics.     

Predicting fluctuations of reintroduced ibex populations: the importance of density dependence, environmental stochasticity and uncertain population estimates

1. Development of population projections requires estimates of observation error, parameters characterizing expected dynamics such as the specific population growth rate and the form of density regulation, the influence of stochastic factors on population dynamics, and quantification of the uncertainty in the parameter estimates. 2. Here we construct a Population Prediction Interval (PPI) based on Bayesian state space modelling of future population growth of 28 reintroduced ibex populations in Switzerland that have been censused for up to 68 years. Our aim is to examine whether the interpopulation variation in the precision of the population projections is related to differences in the parameters characterizing the expected dynamics, in the effects of environmental stochasticity, in the magnitude of uncertainty in the population parameters, or in the observation error. 3. The error in the population censuses was small. The median coefficient of variation in the estimates across populations was 5.1%. 4. Significant density regulation was present in 53.6% of the populations, but was in general weak. 5. The width of the PPI calculated for a period of 5 years showed large variation among populations, and was explained by differences in the impact of environmental stochasticity on population dynamics. 6. In spite of the high accuracy in population estimates, the uncertainty in the parameter estimates was still large. This uncertainty affected the precision in the population predictions, but it decreased with increasing length of study period, mainly due to higher precision in the estimates of the environmental variance in the longer time-series. 7. These analyses reveal that predictions of future population fluctuations of weakly density-regulated populations such as the ibex often become uncertain. Credible population predictions require that this uncertainty is properly quantified.     

Spatial heterogeneity in climate change effects decouples the long‐term dynamics of wild reindeer populations in the high Arctic

The ‘Moran effect’ predicts that dynamics of populations of a species are synchronized over similar distances as their environmental drivers. Strong population synchrony reduces species viability, but spatial heterogeneity in density dependence, the environment, or its ecological responses may decouple dynamics in space, preventing extinctions. How such heterogeneity buffers impacts of global change on large‐scale population dynamics is not well studied. Here, we show that spatially autocorrelated fluctuations in annual winter weather synchronize wild reindeer dynamics across high‐Arctic Svalbard, while, paradoxically, spatial variation in winter climate trends contribute to diverging local population trajectories. Warmer summers have improved the carrying capacity and apparently led to increased total reindeer abundance. However, fluctuations in population size seem mainly driven by negative effects of stochastic winter rain‐on‐snow (ROS) events causing icing, with strongest effects at high densities. Count data for 10 reindeer populations 8–324 km apart suggested that density‐dependent ROS effects contributed to synchrony in population dynamics, mainly through spatially autocorrelated mortality. By comparing one coastal and one ‘continental’ reindeer population over four decades, we show that locally contrasting abundance trends can arise from spatial differences in climate change and responses to weather. The coastal population experienced a larger increase in ROS, and a stronger density‐dependent ROS effect on population growth rates, than the continental population. In contrast, the latter experienced stronger summer warming and showed the strongest positive response to summer temperatures. Accordingly, contrasting net effects of a recent climate regime shift—with increased ROS and harsher winters, yet higher summer temperatures and improved carrying capacity—led to negative and positive abundance trends in the coastal and continental population respectively. Thus, synchronized population fluctuations by climatic drivers can be buffered by spatial heterogeneity in the same drivers, as well as in the ecological responses, averaging out climate change effects at larger spatial scales.     

Climate mediated exogenous forcing and synchrony in populations of the oak aphid in the UK

Contemporary population dynamics theory suggests that animal fluctuations in nature are the result of the combined forces of intrinsic and exogenous factors. Weather is the iconic example of an exogenous force. The common approach for analyzing the relationship between population size and climatic variables is by simple correlation or using the climate as an additive covariable in statistical models. Here, we evaluated different functional forms in which climatic variables could influence population dynamics of the oak aphid Tuberculatus annulatus both in each locality and in relation to synchrony between localities. Results indicate that in at least four of eight aphid populations, climate influences population dynamics by modifying the carrying capacity of the system (lateral effect mediated by winter precipitation). Additionally, path analysis showed that synchrony in population dynamics is highly correlated with synchrony in winter precipitation regime, and the spatial scale of both processes is similar, which suggests that this is an example of the Moran effect. Our results show the key effects of precipitation on intra and inter population processes of this aphid. The methods used, mixing population dynamics modelling and test of synchrony, allowed us to connect the direct and indirect effects of exogenous variables into each population with patterns of synchrony inter populations.     

Climate impacts on Mediterranean blue tit survival: an investigation across seasons and spatial scales

In some hole nesting passerine species, long‐term monitoring data are available for several geographically independent populations. Climate forcing can then be documented and predictions made on the scale of distribution ranges. Several demographic studies of Paridae report dramatic impacts of wintertime climatic factors. However, these studies were undertaken in populations located in the northern parts of the species' ranges. Studies on the survival of Paridae in their southern ranges are necessary in order to assess potential latitudinal variation in climate forcing on survival. Based on monitoring of individual adult blue tits (Parus caeruleus), the effects of climatic factors on annual survival were assessed in three distinct Mediterranean populations. In these regions, climatic conditions in early summer might be expected to have a strong impact because they can be extremely hot and dry and because at this time of year Paridae are subjected to intrinsic constraints that stem from energetically costly postbreeding moult, recovery from reproductive costs, and from population densities inflated by the new cohort of fledglings. The impact of climatic conditions in early summer was, thus, addressed in addition to that prevailing in winter. In order to consider a large number of local climatic variables while limiting statistical power loss, integrative indices of local climate were built using multivariate techniques. In addition, the NAO and three large‐scale factors that are closely linked with atmospheric and oceanic circulation in the intertropical zone were considered as potentially influential factors in winter and early summer. Relationships between blue tit survival and indices of local temperature and precipitation in winter and in early summer were detected. Adult survival also correlated with a large‐scale tropical index in early summer: rainfall in the Sahel. This is one of the first quantitative indications that fluctuations in summer climatic conditions explain a significant part of the temporal variation in adult survival in unconnected populations of a sedentary European vertebrate. Furthermore, the results support the hypothesis that summertime local climates in Western Europe are closely linked with atmospheric and oceanic circulation in the intertropical zone.     

Metapopulation dynamics in a changing climate: Increasing spatial synchrony in weather conditions drives metapopulation synchrony of a butterfly inhabiting a fragmented landscape

Habitat fragmentation and climate change are both prominent manifestations of global change, but there is little knowledge on the specific mechanisms of how climate change may modify the effects of habitat fragmentation, for example, by altering dynamics of spatially structured populations. The long‐term viability of metapopulations is dependent on independent dynamics of local populations, because it mitigates fluctuations in the size of the metapopulation as a whole. Metapopulation viability will be compromised if climate change increases spatial synchrony in weather conditions associated with population growth rates. We studied a recently reported increase in metapopulation synchrony of the Glanville fritillary butterfly (Melitaea cinxia) in the Finnish archipelago, to see if it could be explained by an increase in synchrony of weather conditions. For this, we used 23 years of butterfly survey data together with monthly weather records for the same period. We first examined the associations between population growth rates within different regions of the metapopulation and weather conditions during different life‐history stages of the butterfly. We then examined the association between the trends in the synchrony of the weather conditions and the synchrony of the butterfly metapopulation dynamics. We found that precipitation from spring to late summer are associated with the M. cinxia per capita growth rate, with early summer conditions being most important. We further found that the increase in metapopulation synchrony is paralleled by an increase in the synchrony of weather conditions. Alternative explanations for spatial synchrony, such as increased dispersal or trophic interactions with a specialist parasitoid, did not show paralleled trends and are not supported. The climate driven increase in M. cinxia metapopulation synchrony suggests that climate change can increase extinction risk of spatially structured populations living in fragmented landscapes by altering their dynamics.     

Climate and density shape population dynamics of a marine top predator

Long-term studies have documented that climate fluctuations affect the dynamics of populations, but the relative influence of stochastic and density-dependent processes is still poorly understood and debated. Most studies have been conducted on terrestrial systems, and lacking are studies on marine systems explicitly integrating the fact that most populations live in seasonal environments and respond to regular or systematic environmental changes. We separated winter from summer mortality in a seabird population, the blue petrel Halobaena caerulea, in the southern Indian Ocean where the El Niño/Southern Oscillation effects occur with a 3-4-year lag. Seventy per cent of the mortality occurred in winter and was linked to climatic factors, being lower during anomalous warm events. The strength of density dependence was affected by climate, with population crashes occurring when poor conditions occurred at high densities. We found that an exceptionally long-lasting warming caused a ca. 40% decline of the population, suggesting that chronic climate change will strongly affect this top predator. These findings demonstrate that populations in marine systems are particularly susceptible to climate variation through complex interactions between seasonal mortality and density-dependent effects.     

Geographical gradients in the population dynamics of North American prairie ducks

1. Geographic gradients in population dynamics may occur because of spatial variation in resources that affect the deterministic components of the dynamics (i.e. carrying capacity, the specific growth rate at small densities or the strength of density regulation) or because of spatial variation in the effects of environmental stochasticity. To evaluate these, we used a hierarchical Bayesian approach to estimate parameters characterizing deterministic components and stochastic influences on population dynamics of eight species of ducks (mallard, northern pintail, blue-winged teal, gadwall, northern shoveler, American wigeon, canvasback and redhead (Anas platyrhynchos, A. acuta, A. discors, A. strepera, A. clypeata, A. americana, Aythya valisineria and Ay. americana, respectively) breeding in the North American prairies, and then tested whether these parameters varied latitudinally. 2. We also examined the influence of temporal variation in the availability of wetlands, spring temperature and winter precipitation on population dynamics to determine whether geographical gradients in population dynamics were related to large-scale variation in environmental effects. Population variability, as measured by the variance of the population fluctuations around the carrying capacity K, decreased with latitude for all species except canvasback. This decrease in population variability was caused by a combination of latitudinal gradients in the strength of density dependence, carrying capacity and process variance, for which details varied by species. 3. The effects of environmental covariates on population dynamics also varied latitudinally, particularly for mallard, northern pintail and northern shoveler. However, the proportion of the process variance explained by environmental covariates, with the exception of mallard, tended to be small. 4. Thus, geographical gradients in population dynamics of prairie ducks resulted from latitudinal gradients in both deterministic and stochastic components, and were likely influenced by spatial differences in the distribution of wetland types and shapes, agricultural practices and dispersal processes. 5. These results suggest that future management of these species could be improved by implementing harvest models that account explicitly for spatial variation in density effects and environmental stochasticity on population abundance.     

Climate variation and regional gradients in population dynamics of two hole-nesting passerines

Latitudinal gradients in population dynamics can arise through regional variation in the deterministic components of the population dynamics and the stochastic factors. Here, we demonstrate an increase with latitude in the contribution of a large-scale climate pattern, the North Atlantic Oscillation (NAO), to the fluctuations in size of populations of two European hole-nesting passerine species. However, this influence of climate induced different latitudinal gradients in the population dynamics of the two species. In the great tit the proportion of the variability in the population fluctuations explained by the NAO increased with latitude, showing a larger impact of climate on the population fluctuations of this species at higher latitudes. In contrast, no latitudinal gradient was found in the relative contribution of climate to the variability of the pied flycatcher populations because the total environmental stochasticity increased with latitude. This shows that the population ecological consequences of an expected climate change will depend on how climate affects the environmental stochasticity in the population process. In both species, the effects will be larger in those parts of Europe where large changes in climate are expected.     

Range position and climate sensitivity: The structure of among‐population demographic responses to climatic variation

Species’ distributions will respond to climate change based on the relationship between local demographic processes and climate and how this relationship varies based on range position. A rarely tested demographic prediction is that populations at the extremes of a species’ climate envelope (e.g., populations in areas with the highest mean annual temperature) will be most sensitive to local shifts in climate (i.e., warming). We tested this prediction using a dynamic species distribution model linking demographic rates to variation in temperature and precipitation for wood frogs (Lithobates sylvaticus) in North America. Using long‐term monitoring data from 746 populations in 27 study areas, we determined how climatic variation affected population growth rates and how these relationships varied with respect to long‐term climate. Some models supported the predicted pattern, with negative effects of extreme summer temperatures in hotter areas and positive effects on recruitment for summer water availability in drier areas. We also found evidence of interacting temperature and precipitation influencing population size, such as extreme heat having less of a negative effect in wetter areas. Other results were contrary to predictions, such as positive effects of summer water availability in wetter parts of the range and positive responses to winter warming especially in milder areas. In general, we found wood frogs were more sensitive to changes in temperature or temperature interacting with precipitation than to changes in precipitation alone. Our results suggest that sensitivity to changes in climate cannot be predicted simply by knowing locations within the species’ climate envelope. Many climate processes did not affect population growth rates in the predicted direction based on range position. Processes such as species‐interactions, local adaptation, and interactions with the physical landscape likely affect the responses we observed. Our work highlights the need to measure demographic responses to changing climate.     

Local and cross‐seasonal associations of climate and land use with abundance of monarch butterflies Danaus plexippus

Quantifying how climate and land use factors drive population dynamics at regional scales is complex because it depends on the extent of spatial and temporal synchrony among local populations, and the integration of population processes throughout a species’ annual cycle. We modeled weekly, site‐specific summer abundance (1994–2013) of monarch butterflies Danaus plexippus at sites across Illinois, USA to assess relative associations of monarch abundance with climate and land use variables during the winter, spring, and summer stages of their annual cycle. We developed negative binomial regression models to estimate monarch abundance during recruitment in Illinois as a function of local climate, site‐specific crop cover, and county‐level herbicide (glyphosate) application. We also incorporated cross‐seasonal covariates, including annual abundance of wintering monarchs in Mexico and climate conditions during spring migration and breeding in Texas, USA. We provide the first empirical evidence of a negative association between county‐level glyphosate application and local abundance of adult monarchs, particularly in areas of concentrated agriculture. However, this association was only evident during the initial years of the adoption of herbicide‐resistant crops (1994–2003). We also found that wetter and, to a lesser degree, cooler springs in Texas were associated with higher summer abundances in Illinois, as were relatively cool local summer temperatures in Illinois. Site‐specific abundance of monarchs averaged approximately one fewer per site from 2004–2013 than during the previous decade, suggesting a recent decline in local abundance of monarch butterflies on their summer breeding grounds in Illinois. Our results demonstrate that seasonal climate and land use are associated with trends in adult monarch abundance, and our approach highlights the value of considering fine‐resolution temporal fluctuations in population‐level responses to environmental conditions when inferring the dynamics of migratory species.     

Generalizations of the Moran effect explaining spatial synchrony in population fluctuations

The Moran effect for populations separated in space states that the autocorrelations in the population fluctuations equal the autocorrelation in environmental noise, assuming the same linear density regulation in all populations. Here we generalize the Moran effect to include also nonlinear density regulation with spatial heterogeneity in local population dynamics as well as in the effects of environmental covariates by deriving a simple expression for the correlation between the sizes of two populations, using diffusion approximation to the theta-logistic model. In general, spatial variation in parameters describing the dynamics reduces population synchrony. We also show that the contribution of a covariate to spatial synchrony depends strongly on spatial heterogeneity in the covariate or in its effect on local dynamics. These analyses show exactly how spatial environmental covariation can synchronize fluctuations of spatially segregated populations with no interchange of individuals even if the dynamics are nonlinear.     

Dynamic social system in Nubian ibex: can a second mating season develop in response to arid climate?

We studied a population of Nubian ibex Capra ibex nubiana in the eastern extreme of its range, the hyper-arid central desert of the Sultanate of Oman. Long-term data were collected from January 1983 to December 1997 by direct observation, as well as VHF telemetry on 12 animals (eight from 1987 to 1990; four from 1994 to 1996). We recorded 884 sightings: 40.4% of single animals and 59.6% of groups. Although no significant monthly variation of group size (Jarman's Typical Group Size) was found, there were distinct peaks in March (4.0 ind. group⁻¹) and September (5.1 ind. group⁻¹). Groups of males and females formed especially in March and November, and female–kid groups in February and July–August. Our data may suggest two mating periods: the first one in autumn (similar to the rut of ibex in temperate mountain areas), with kids born in spring/early summer, after winter–spring rainfall, and the second one in spring, with kids born in late summer/autumn, before winter–spring rainfalls. We suggest that the second rutting period may have evolved as a micro-evolutionary process, with the local population adapting to hyper-arid environment constraints. The spring mating season may favour only females in prime conditions, who can afford a pregnancy in the local severe summers and will deliver kids when plant greening begins, in the autumn, whereas the autumn (original) mating season may be afforded by any female, but kids will be born in an unfavourable period, before the summer drought.     

Climate variables explain neutral and adaptive variation within salmonid metapopulations: the importance of replication in landscape genetics

Understanding how environmental variation influences population genetic structure is important for conservation management because it can reveal how human stressors influence population connectivity, genetic diversity and persistence. We used riverscape genetics modelling to assess whether climatic and habitat variables were related to neutral and adaptive patterns of genetic differentiation (population‐specific and pairwise F_ST) within five metapopulations (79 populations, 4583 individuals) of steelhead trout (Oncorhynchus mykiss) in the Columbia River Basin, USA. Using 151 putatively neutral and 29 candidate adaptive SNP loci, we found that climate‐related variables (winter precipitation, summer maximum temperature, winter highest 5% flow events and summer mean flow) best explained neutral and adaptive patterns of genetic differentiation within metapopulations, suggesting that climatic variation likely influences both demography (neutral variation) and local adaptation (adaptive variation). However, we did not observe consistent relationships between climate variables and F_ST across all metapopulations, underscoring the need for replication when extrapolating results from one scale to another (e.g. basin‐wide to the metapopulation scale). Sensitivity analysis (leave‐one‐population‐out) revealed consistent relationships between climate variables and F_ST within three metapopulations; however, these patterns were not consistent in two metapopulations likely due to small sample sizes (N = 10). These results provide correlative evidence that climatic variation has shaped the genetic structure of steelhead populations and highlight the need for replication and sensitivity analyses in land and riverscape genetics.     

Density-dependent responses of fawn cohort body mass in two contrasting roe deer populations

We investigated the influence of population density on juvenile body mass in two contrasting roe deer populations, in Sweden (Bogesund) and France (Chizé), in which density was monitored for ≥15 years. We investigated the effect of population density and climatic conditions on cohort performance. We predicted that: (1) body mass of growing fawns should be sensitive to environmental changes, showing marked between-year variation (i.e., cohort effects), (2) fawns in the less productive (weakly seasonal, weakly predictable summer weather) habitat of Chizé should show stronger density-dependent responses due to more severe food competition during summer than fawns in the more productive (markedly seasonal, moderately predictable summer weather) habitat of Bogesund, and (3) fawns at Bogesund should be heavier both in absolute terms and relative to their size than their conspecifics in Chizé due to a higher degree of fat accumulation in northern environments. In both study sites we found marked cohort variation and clear effects of density, with body mass varying by as much as 29% over years. While neither summer nor winter climate influenced fawn body mass at Bogesund, fawns tended to be lighter after summers with high temperatures at Chizé. In addition, fawns were heavier after acorn mast years experienced in utero at Bogesund. As expected, the strength of the density-dependent response of fawn body mass was greater at Chizé than at Bogesund. For a given density, male fawns were consistently heavier than females in both sites. Lastly, both sexes at Bogesund had higher absolute body mass and were larger for a given body size than in Chizé. Our results clearly demonstrate that absolute density is a poor predictor of roe deer performance and supports the view that habitat quality has an overwhelming importance for determining fawn body mass in roe deer populations.     

Spatial variation in breeding phenology at small spatial scales: A stochastic effect of population size

Spatial variation in phenology can occur at small spatial scales over which individuals can disperse or forage within one generation. Previous studies have assumed that variations in phenological peaks are caused by differences in abiotic environmental characteristics. However, environments should generally be similar among local habitats over small spatial scales. When the local population size is small, the phenological peak of the local population should be strongly affected by the variation in timing expressed by individuals. If a regional population consists of small local subpopulations (e.g., a metapopulation), the stochastic processes regulated by population sizes may explain the spatial variation in phenology. In this study, we quantitatively evaluated the extent of the spatial and annual variations in the breeding phenology of the forest green tree frog, Rhacophorus arboreus habiting a small area (<10 km²). The spatial variation in phenological peaks among 25 breeding sites was large over 6 years. This spatial variation was not explained by differences in air temperature or water depth. Randomization tests revealed that a large portion of the spatial variation could be explained by differences in population size, without considering site-specific factors. Annual variations in phenological peaks tended to be greater for smaller populations. These results imply that the stochastic process might have caused the spatial and annual variations in the phenological peaks of R. arboreus observed in the study region. Understanding spatiotemporal variation in phenology determined by stochastic process would be important to better predict interspecific interactions and (meta)population dynamics at small spatial scales.     

Non-linear feedback processes and a latitudinal gradient in the climatic effects determine green spruce aphid outbreaks in the UK

The role of climatic fluctuations in determining the dynamics of insect populations has been a classical problem in population ecology. Here, we use long-term annual data on green spruce aphid populations at nine localities in the UK for determining the importance of endogenous processes, local weather and large-scale climatic factors. We rely on diagnostic and modelling tools from population dynamic theory to analyse these long-term data and to determine the role of the North Atlantic Oscillation (NAO) and local weather as exogenous factors influencing aphid dynamics. Our modelling suggests that the key elements determining population fluctuations in green spruce aphid populations in the UK are the strong non-linear feedback structure, the high potential for population growth and the effects of winter and spring weather. The results indicate that the main effect of the NAO on green spruce aphid populations is operating through the effect of winter temperatures on the maximum per capita growth rate (R_m). In particular, we can predict quite accurately the occurrence of an outbreak by using a simple logistic model with weather as a perturbation effect. However, model predictions using different climatic variables showed a clear geographical signature. The NAO and winter temperature were best for predicting observed dynamics toward the southern localities, while spring temperature was a much better predictor of aphid dynamics at northern localities. Although aphid species are characterized by complex life-cycles, we emphasize the value of simple and general population dynamic models in predicting their dynamics.     

Modelling non-additive and nonlinear signals from climatic noise in ecological time series: Soay sheep as an example

Understanding how climate can interact with other factors in determining patterns of species abundance is a persistent challenge in ecology. Recent research has suggested that the dynamics exhibited by some populations may be a non-additive function of climate, with climate affecting population growth more strongly at high density than at low density. However, we lack methodologies to adequately explain patterns in population growth generated as a result of interactions between intrinsic factors and extrinsic climatic variation in non-linear systems. We present a novel method (the Functional Coefficient Threshold Auto-Regressive (FCTAR) method) that can identify interacting influences of climate and density on population dynamics from time-series data. We demonstrate its use on count data on the size of the Soay sheep population, which is known to exhibit dynamics generated by nonlinear and non-additive interactions between density and climate, living on Hirta in the St Kilda archipelago. The FCTAR method suggests that climate fluctuations can drive the Soay sheep population between different dynamical regimes--from stable population size through limit cycles and non-periodic fluctuations.     

Effects of density-dependence and climate on the dynamics of a Svalbard reindeer population

To determine the main factors affecting the population dynamics of Svalbard reindeer, we analysed 21 yr of annual censuses, including data on population size, recruitment rate (calves per female) and mortality (number of deaths), from the Reindalen reindeer population. In accordance with previous studies on population dynamics of Svalbard reindeer, we found large inter-annual variation in population size, mortality and recruitment rates within the study area. Population size decreased in years with low recruitment rate as well as high winter mortality and vice versa. Apparently. the fluctuations were due to both direct density-dependent food limitation and variation in winter climate associated with high precipitation and icing of the feeding range. We found no delayed density-dependence or effect of climatic conditions during summer on the population dynamics. The mortality during die-off years was mainly of calves and very old individuals, indicating that the population was more vulnerable to high die oft in years following high recruitment rate. These results suggest an unstable interaction between the reindeer population and its food supply in these predator-free environments.