The effects of oil spills on marine fish: Implications of spatial variation in natural mortality

The effects of oil spills on marine biological systems are of great concern, especially in regions with high biological production of harvested resources such as in the Northeastern Atlantic. The scientific studies of the impact of oil spills on fish stocks tend to ignore that spatial patterns of natural mortality may influence the magnitude of the impact over time. Here, we first illustrate how spatial variation in natural mortality may affect the population impact by considering a thought experiment. Second, we consider an empirically based example of Northeast Arctic cod to extend the concept to a realistic setting. Finally, we present a scenario-based investigation of how the degree of spatial variation in natural mortality affects the impact over a gradient of oil spill sizes. Including the effects of spatial variations in natural mortality tends to widen the impact distribution, hence increasing the probability of both high and low impact events.

Spatial ecology of coastal Atlantic cod Gadus morhua associated with parasite load

Acoustic tags and receivers were used to investigate the spatial ecology of coastal Atlantic cod Gadus morhua (n = 32, mean fork length: 50 cm, range: 33-80 cm) on the Norwegian Skagerrak coast in 2012. Monthly home ranges (HR), swimming activity and depth use varied considerably among individuals and through the months of June, July and August. HR sizes for the period ranged from 0.25 to 5.20 km2 (mean = 2.30 km2. Two thirds of the tagged G. morhua were infected with black spot disease Cryptocotyle lingua parasites; these fish had larger HRs and occupied deeper water compared with non-infected fish. The infected fish also tended to be more active in terms of horizontal swimming. From an ecological and evolutionary perspective, any environmental change that modifies G. morhua behaviour may therefore also alter the parasite load of the population, and its conservation and fishery status.

A combination of hydrodynamical and statistical modelling reveals non-stationary climate effects on fish larvae distributions

Biological processes and physical oceanography are often integrated in numerical modelling of marine fish larvae, but rarely in statistical analyses of spatio-temporal observation data. Here, we examine the relative contribution of inter-annual variability in spawner distribution, advection by ocean currents, hydrography and climate in modifying observed distribution patterns of cod larvae in the Lofoten-Barents Sea. By integrating predictions from a particle-tracking model into a spatially explicit statistical analysis, the effects of advection and the timing and locations of spawning are accounted for. The analysis also includes other environmental factors: temperature, salinity, a convergence index and a climate threshold determined by the North Atlantic Oscillation (NAO). We found that the spatial pattern of larvae changed over the two climate periods, being more upstream in low NAO years. We also demonstrate that spawning distribution and ocean circulation are the main factors shaping this distribution, while temperature effects are different between climate periods, probably due to a different spatial overlap of the fish larvae and their prey, and the consequent effect on the spatial pattern of larval survival. Our new methodological approach combines numerical and statistical modelling to draw robust inferences from observed distributions and will be of general interest for studies of many marine fish species.

Within and between species competition in a seabird community: statistical exploration and modeling of time-series data

In a changing environment, the maintenance of communities is subject to many constraints (phenology, resources, climate, etc.). One such constraint is the relationship between conspecifics and competitors. In mixed colonies, seabirds may have to cope with interspecific and intraspecific competition for both space and food resources. We applied competitive interaction models to data on three seabird breeding populations: black-legged kittiwake (Rissa tridactyla), common guillemot (Uria aalge) and Brünnich's guillemot (Uria lomvia) collected over 27-years at Kharlov Island in the Barents Sea. We found a competitive effect only for the kittiwake breeding population size on the common guillemot breeding population size when kittiwakes were abundant. The timing of kittiwake breeding negatively affected the number of breeding Brünnich's guillemots. The timing of breeding was negatively correlated to biomass of the main pelagic fish in the Barents Sea, the capelin (Mallotus villosus), which suggests an indirect action. The community matrix shows that the community was not stable. The kittiwake population did not decrease as seen in north Norwegian populations. Likewise, the common guillemot population, after a crash in 1985, was recovering at Kharlov while Norwegian populations were decreasing. Only the Brünnich's guillemot showed a decrease at Kharlov until 1999. We suggest that the stability of the kittiwake and common guillemot populations at Kharlov is due to better feeding conditions than in colonies of the Norwegian coast, linked to a possible eastward shift of the capelin population with the temperature increase of the Barents Sea.

Are low but statistically significant levels of genetic differentiation in marine fishes 'biologically meaningful'? A case study of coastal Atlantic cod

A key question in many genetic studies on marine organisms is how to interpret a low but statistically significant level of genetic differentiation. Do such observations reflect a real phenomenon, or are they caused by confounding factors such as unrepresentative sampling or selective forces acting on the marker loci? Further, are low levels of differentiation biologically trivial, or can they represent a meaningful and perhaps important finding? We explored these issues in an empirical study on coastal Atlantic cod, combining temporally replicated genetic samples over a 10-year period with an extensive capture-mark-recapture study of individual mobility and population size. The genetic analyses revealed a pattern of differentiation between the inner part of the fjord and the open skerries area at the fjord entrance. Overall, genetic differentiation was weak (average F(ST)  = 0.0037), but nevertheless highly statistical significant and did not depend on particular loci that could be subject to selection. This spatial component dominated over temporal change, and temporal replicates clustered together throughout the 10-year period. Consistent with genetic results, the majority of the recaptured fish were found close to the point of release, with <1% of recaptured individuals dispersing between the inner fjord and outer skerries. We conclude that low levels of genetic differentiation in this marine fish can indeed be biologically meaningful, corresponding to separate, temporally persistent, local populations. We estimated the genetically effective sizes (N(e) ) of the two coastal cod populations to 198 and 542 and found a N(e) /N (spawner) ratio of 0.14.

Detecting genetic structure in migrating bowhead whales off the coast of Barrow, Alaska

We develop a general framework for analysing and testing genetic structure within a migratory assemblage that is based on measures of genetic differences between individuals. We demonstrate this method using microsatellite DNA data from the Bering-Chukchi-Beaufort stock of bowhead whales (Balaena mysticetus), sampled via Inuit hunting during the spring and autumn migration off Barrow, Alaska. This study includes a number of covariates such as whale ages and the time separation between captures. Applying the method to a sample of 117 bowhead whales, we use permutation methods to test for temporal trends in genetic differences that can be ascribed to age-related effects or to timing of catches during the seasons. The results reveal a pattern with elevated genetic differences among whales caught about a week apart, and are statistically significant for the autumn migration. In contrast, we find no effects of time of birth or age-difference on genetic differences. We discuss possible explanations for the results, including population substructuring, demographic consequences of historical overexploitation, and social structuring during migration.

Cryptic population structuring in Scandinavian lynx: reply to Pamilo

In a recent Commentary in this journal, Pamilo (2004) criticized our analysis of the spatial genetic structure of the Eurasian lynx in Scandinavia (Rueness et al. 2003). The analyses uncovered a marked geographical differentiation along the Scandinavian peninsula with an apparent linear gradient in the north-south direction. We used computer simulations to check on the proposition that the observed geographical structure could have arisen by genetic drift and isolation by distance in the approximate 25 generations that have passed since the last bottleneck. Pamilo disapproved of our choice of population model and also how we compared the outcome of the simulations with data. As these issues should be of interest to a wider audience we discuss them in some detail.

Cryptic population structure in a large, mobile mammalian predator: the Scandinavian lynx

The Eurasian lynx (Lynx lynx) is an example of a species that has gone through a severe bottleneck, leading to near extinction in Scandinavia around 1930-- a pattern shared with several other large carnivorous mammals. Here we extend previous genetic analyses of northern European lynx, confirming that lynx from the Scandinavian Peninsula represent a distinct clade differing clearly from European conspecifics. Furthermore, and despite a recent bottleneck and subsequent range expansion, we detect marked genetic differentiation within Scandinavia. This differentiation is largely manifested as a north-south gradient, with a linear increase in the quantity FST/(1 - FST). Aided by computer simulations we find that this pattern is unlikely to have arisen by random genetic drift in the short time since lynx started to expand in the 1950s, suggesting that the spatial structure may predate the bottleneck. Individual-based analyses indicate that, instead of a continuous gradient, Scandinavian lynx may be structured into three more or less distinct groups, possibly corresponding to northern, central and southern subpopulations. The presence of such structuring was unknown previously and was unexpected from general considerations on the mobility of the species, historical data and the absence of geographical barriers. Our study demonstrates how molecular markers may be used to detect cryptic population structure, invisible using traditional methods.

Parasites and supernormal manipulation

Social parasites may exploit their hosts by mimicking other organisms that the hosts normally benefit from investing in or responding to in some other way. Some parasites exaggerate key characters of the organisms they mimic, possibly in order to increase the response from the hosts. The huge gape and extreme begging intensity of the parasitic common cuckoo chick (Cuculus canorus) may be an example. In this paper, the evolutionary stability of manipulating hosts through exaggerated signals is analysed using game theory. Our model indicates that a parasite's signal intensity must be below a certain threshold in order to ensure acceptance and that this threshold depends directly on the rate of parasitism. The only evolutionarily stable strategy (ESS) combination is when hosts accept all signallers and parasites signal at their optimal signal intensity, which must be below the threshold. Supernormal manipulation by parasites is only evolutionarily stable under sufficiently low rates of parasitism. If the conditions for the ESS combination are not satisfied, rejector hosts can invade using signal intensity as a cue for identifying parasites. These qualitative predictions are discussed with respect to empirical evidence from parasitic mimicry systems that have been suggested to involve supernormal signalling, including evicting avian brood parasites and insect-mimicking Ophrys orchids.

Phase coupling and synchrony in the spatiotemporal dynamics of muskrat and mink populations across Canada

Population ecologists have traditionally focused on the patterns and causes of population variation in the temporal domain for which a substantial body of practical analytic techniques have been developed. More recently, numerous studies have documented how populations may fluctuate synchronously over large spatial areas; analyses of such spatially extended time-series have started to provide additional clues regarding the causes of these population fluctuations and explanations for their synchronous occurrence. Here, we report on the development of a phase-based method for identifying coupling between temporally coincident but spatially distributed cyclic time-series, which we apply to the numbers of muskrat and mink recorded at 81 locations across Canada. The analysis reveals remarkable parallel clines in the strength of coupling between proximate populations of both species--declining from west to east--together with a corresponding increase in observed synchrony between these populations the further east they are located.

Effects of age, sex and density on body weight of Norwegian red deer: evidence of density-dependent senescence

There are only a few recent studies that have demonstrated senescence in ungulates and nothing is known regarding how patterns of senescence may vary as a function of density Senescence in general is linked to the cost of reproduction, which probably increases with density in ungulates and may differ between the sexes. Further, senescence in ungulates is also regarded to be a function of tooth wear rates. As density dependence and sexual differences in food choice have been well documented, this may lead to different tooth wear rates and, thus, possibly density-dependent and sex-specific patterns of senescence. We therefore investigated the effects of age, sex, density and their possible interactions on the variability of body weight in 29,047 red deer harvested during 1965-1998 from Norway, out of which 380 males and 1452 females were eight years or older. There was clear evidence that spatio-temporal variation in density correlated negatively with body weights. In addition, there was evidence of senescence in both male and female red deer. Age at onset of senescence in females was after 20 years of age and independent of population density. In males, the onset and rate of senescence increased with increasing population density. The onset of senescence for males was at ca. 12 years of age at low density, but decreased to approximately ten years of age at high density. The pattern of density-dependent senescence in males, but not that in females, can be explained if (i) the cost of reproduction increases with density more strongly in male than in female red deer, and/or (ii) tooth wear rates are density dependent in males, but not in females. We discuss the ability of these two different, not mutually exclusive hypotheses in explaining the observed pattern of senescence.

Demography and population dynamics of the mouse opossum (Thylamys elegans) in semi-arid Chile: seasonality, feedback structure and climate

Here we present, to the authors' knowledge for the very first time for a small marsupial, a thorough analysis of the demography and population dynamics of the mouse opossum (Thylamys elegans) in western South America. We test the relative importance of feedback structure and climatic factors (rainfall and the Southern Oscillation Index) in explaining the temporal variation in the demography of the mouse opossum. The demographic information was incorporated into a stage-structured population dynamics model and the model's predictions were compared with observed patterns. The mouse opossum's capture rates showed seasonal (within-year) and between-year variability, with individuals having higher capture rates during late summer and autumn and lower capture rates during winter and spring. There was also a strong between-year effect on capture probabilities. The reproductive (the fraction of reproductively active individuals) and recruitment rates showed a clear seasonal and a between-year pattern of variation with the peak of reproductive activity occuring during winter and early spring. In addition, the fraction of reproductive individuals was positively related to annual rainfall, while population density and annual rainfall positively influenced the recruitment rate. The survival rates were negatively related to annual rainfall. The average finite population growth rate during the study period was estimated to be 1.011 +/- 0.0019 from capture-recapture estimates. While the annual growth rate estimated from the seasonal linear matrix models was 1.026, the subadult and adult survival and maturation rates represent between 54% (winter) and 81% (summer) of the impact on the annual growth rate.

A new geographical gradient in vole population dynamics

A new geographical gradient in the dynamics of small rodents is demonstrated by analysing 29 time series of density indices of the common vole (Microtus arvalis) from Poland, the Czech Republic and the Slovak Republic. This gradient extends from more stable northerly populations in coastal Poland to more variable and cyclic populations in the southernmost parts of the Slovak Republic, and is hence a reversal of the Fennoscandian gradient. All studied variables (such as mean density, cycle amplitude, density variability and the coefficients in a second-order autoregressive model) exhibit consistent latitudinal variation. Possible underlying factors are discussed. In particular, we suggest that seasonality may be a key element in explaining the observed new gradient.

Life-history traits of voles in a fluctuating population respond to the immediate environment

Life-history traits relating to growth and reproduction vary greatly among species and populations and among individuals within populations. In vole populations, body size and age at maturation may vary considerably among locations and among years within the same location. Individuals in increasing populations are typically larger and start reproduction earlier in the spring than those in declining populations. The cause of such life-history variation within populations has been subject of much discussion. Much of the controversy concerns whether the memory of past conditions, leading to delayed effects on life-history traits, resides in the environment (for example, predators, pathogens or food) or intrinsically within populations or individuals (age distribution, physiological state, genetic or maternal effects). Here we report from an extensive field transplant experiment in which voles were moved before the breeding season between sites that differed in average overwintering body mass. Transplanted voles did not retain the characteristics of their source population, and we demonstrate an over-riding role of the immediate environment in shaping life-history traits of small rodents.

Genetic structure of Siberian lemmings (Lemmus sibiricus) in a continuous habitat: large patches rather than isolation by distance

In a continuous habitat, restricted dispersal and local genetic drift are likely to create a pattern of increasing genetic differentiation with distance. Here, we describe the genetic structure of Siberian lemming (Lemmus sibiricus) populations in a continuous tundra habitat on the western coast of the Taimyr Peninsula, in order to determine the spatial scale at which genetic differentiation and isolation by distance occur. Sampling was carried out at three different geographical scales: (1) a continuous 11 km transect; (2) localities 10-30 km apart; and (3) two localities at 300 and 600 km from the main study area. Two types of genetic markers were used: partial sequences of the mitochondrial DNA control region and four microsatellite loci. On this basis the study populations were genetically quite homogeneous within patches extending over 8 km or more. Contrary to theoretical predictions, no pattern of isolation by distance among patches could be identified. This observation was interpreted as representing populations in migration-drift disequilibrium after a recent major mixing event. The lack of concordance between mtDNA haplotype phylogeny and the geographical distribution of haplotypes supported this interpretation. Spatial autocorrelation among individual genotypes on a local scale was weak and observed only in females, indicating a considerable amount of mostly male-mediated gene flow. Average gene flow per generation was estimated to be in the range of several hundred metres.

Nonlinear effects of large-scale climatic variability on wild and domestic herbivores

Large-scale climatic fluctuations, such as the North Atlantic Oscillation (NAO), have been shown to affect many ecological processes. Such effects have been typically assumed to be linear. Only one study has reported a nonlinear relation; however, that nonlinear relation was monotonic (that is, no reversal). Here we show that there is a strong nonlinear and non-monotonic (that is, reversed) effect of the NAO on body weight during the subsequent autumn for 23,838 individual wild red deer (Cervus elaphus) and 139,485 individual domestic sheep (Ovis aries) sampled over several decades on the west coast of Norway. These relationships are, at least in part, explained by comparable nonlinear and non-monotonic relations between the NAO and local climatic variables (temperature, precipitation and snow depth). The similar patterns observed for red deer and sheep, the latter of which live indoors during winter and so experience a stable energy supply in winter, suggest that the (winter) climatic variability (for which the index is a proxy) must influence the summer foraging conditions directly or indirectly.

Glacial survival of the Norwegian lemming (Lemmus lemmus) in Scandinavia: inference from mitochondrial DNA variation

In order to evaluate the biogeographical hypothesis that the Norwegian lemming (Lemmus lemmus) survived the last glacial period in some Scandinavian refugia, we examined variation in the nucleotide sequence of the mitochondrial control region (402 base pairs (bp)) and the cytochrome b (cyt b) region (633 bp) in Norwegian and Siberian (Lemmus sibiricus) lemmings. The phylogenetic distinction and cyt b divergence estimate of 1.8% between the Norwegian and Siberian lemmings suggest that their separation pre-dated the last glaciation and imply that the Norwegian lemming is probably a relic of the Pleistocene populations from Western Europe. The star-like control region phylogeny and low mitochondrial DNA diversity in the Norwegian lemming indicate a reduction in its historical effective size followed by population expansion. The average estimate of post-bottleneck time (19-21 kyr) is close to the last glacial maximum (18-22 kyr BP). Taking these findings and the fossil records into consideration, it seems likely that, after colonization of Scandinavia in the Late Pleistocene, the Norwegian lemming suffered a reduction in its population effective size and survived the last glacial maximum in some local Scandinavian refugia, as suggested by early biogeographical work.

Environmental change and rates of evolution: the phylogeographic pattern within the hartebeest complex as related to climatic variation

Global climate fluctuated considerably throughout the Pliocene-Pleistocene period, influencing the evolutionary history of a wide array of species. Using the phylogeographic patterns within the hartebeest (Alcelaphus buselaphus (Pallas, 1766)) complex, we evaluated the evolutionary consequences of such environmental change for a typical large mammal ranging on the African savannah. Our results, as generated from two mitochondrial DNA markers (the D-loop and cytochrome b), suggest an origin of the hartebeest in eastern Africa from where the species has colonized other parts of the continent. Phylogenetic analyses revealed an early diversification into southern and northern hartebeest lineages, an event that may be related to the formation of the Rift Valley lakes. The northern lineage has further diverged into eastern and western lineages, most probably as a result of the expanding central African rainforest belt and subsequent contraction of savannah habitats during a period of global warming. The diversification events appear to have coincided with major climatic changes and are highly correlated in time. These observations strongly suggest that large-scale climatic fluctuations have been a major determinant for the species' evolutionary history and that hartebeest evolution has mainly taken place in isolated yet environmentally favourable refugia during periods of global warming. Indications of sudden population expansion for two putative ancestral hartebeest populations provide further support for a refugia-based explanation of the diversification events. Reciprocal monophyly between southern and northern lineages may suggest that reproductive barriers exist and that the hartebeest complex comprises two different species.

Reproductive asynchrony increases with environmental disturbance

While it is widely recognized that the manner in which organisms adjust their timing of reproduction reflects evolutionary strategies aimed at minimizing offspring mortality or maximizing reproductive output, the conditions under which the evolutionarily stable strategy involves synchronous or asynchronous reproduction is a matter of considerable discord. A recent theoretical model predicts that whether a population displays reproductive synchrony or asynchrony will depend on the relative scales of intrinsic regulation and environmental disturbance experienced by reproducing individuals. This model predicts that, under conditions of negligible competition and large-scale environmental perturbation, evolution of a single mixed strategy will result in asynchronous reproduction. We tested this prediction using empirical data on large-scale climatic fluctuation and the annual timing of reproduction by three species of flowering plants covering 1300-population-years and four degrees of latitude in Norway. In agreement with model predictions, within populations of all three species reproductive asynchrony increased with the magnitude of large-scale climatic perturbation, but bore no relation to the strength of local density dependence. These results suggest that mixed evolutionarily stable strategies can arise from the interplay of combinations of agents of selection and the scale at which they operate; hence it is fruitless to associate synchronous versus asynchronous timing with particular single factors like climate, competition, or predation.

The impact of specialized enemies on the dimensionality of host dynamics

Although individual species persist within a web of interactions with other species, data are usually gathered only from the focal species itself. We ask whether evidence of a species' interactions be detected and understood from patterns in the dynamics of that species alone. Theory predicts that strong coupling between a prey and a specialist predator/parasite should lead to an increase in the dimensionality of the prey's dynamics, whereas weak coupling should not. Here we describe a rare test of this prediction. Two natural enemies were added separately to replicate populations of a moth. For biological reasons that we identify here, the prediction of increased dimensionality was confirmed when a parasitoid wasp was added (although this increase had subtleties not previously appreciated), but the prediction failed for an added virus. Thus, an imprint of the interactions may be discerned within time-series data from component species of a system.

Spatial structure of lemming populations (Dicrostonyx groenlandicus) fluctuating in density

The pattern and scale of the genetic structure of populations provides valuable information for the understanding of the spatial ecology of populations, including the spatial aspects of density fluctuations. In the present paper, the genetic structure of periodically fluctuating lemmings (Dicrostonyx groenlandicus) in the Canadian Arctic was analysed using mitochondrial DNA (mtDNA) control region sequences and four nuclear microsatellite loci. Low genetic variability was found in mtDNA, while microsatellite loci were highly variable in all localities, including localities on isolated small islands. For both genetic markers the genetic differentiation was clear among geographical regions but weaker among localities within regions. Such a pattern implies gene flow within regions. Based on theoretical calculations and population census data from a snap-trapping survey, we argue that the observed genetic variability on small islands and the low level of differentiation among these islands cannot be explained without invoking long distance dispersal of lemmings over the sea ice. Such dispersal is unlikely to occur only during population density peaks.

The timing of life-history events in a changing climate

Although empirical and theoretical studies suggest that climate influences the timing of life-history events in animals and plants, correlations between climate and the timing of events such as egg-laying, migration or flowering do not reveal the mechanisms by which natural selection operates on life-history events. We present a general autoregressive model of the timing of life-history events in relation to variation in global climate that, like autoregressive models of population dynamics, allows for a more mechanistic understanding of the roles of climate, resources and competition. We applied the model to data on 50 years of annual dates of first flowering by three species of plants in 26 populations covering 4 degrees of latitude in Norway. In agreement with earlier studies, plants in most populations and all three species bloomed earlier following warmer winters. Moreover, our model revealed that earlier blooming reflected increasing influences of resources and density-dependent population limitation under climatic warming. The insights available from the application of this model to phenological data in other taxa will contribute to our understanding of the roles of endogenous versus exogenous processes in the evolution of the timing of life-history events in a changing climate.

Common structure in panels of short ecological time-series

Typically, in many studies in ecology, epidemiology, biomedicine and others, we are confronted with panels of short time-series of which we are interested in obtaining a biologically meaningful grouping. Here, we propose a bootstrap approach to test whether the regression functions or the variances of the error terms in a family of stochastic regression models are the same. Our general setting includes panels of time-series models as a special case. We rigorously justify the use of the test by investigating its asymptotic properties, both theoretically and through simulations. The latter confirm that for finite sample size, bootstrap provides a better approximation than classical asymptotic theory. We then apply the proposed tests to the mink-muskrat data across 81 trapping regions in Canada. Ecologically interpretable groupings are obtained, which serve as a necessary first step before a fuller biological and statistical analysis of the food chain interaction.

Phylogeography and mitochondrial DNA (mtDNA) diversity in North American collared lemmings (Dicrostonyx groenlandicus)

Variation in the nucleotide sequence of the mitochondrial control region (250 bp) and the cytochrome b region (870 bp) was examined in collared lemmings (Dicrostonyx groenlandicus) from 19 localities in northern Alaska and the Canadian Arctic. The division of D. groenlandicus in two phylogeographical groups with limited divergence across the Mackenzie River is consistent with the separation of this species in more than one refugial area located to the northwest of the Laurentide ice sheet during the last glaciation. Populations of D.groenlandicus from formerly glaciated areas are no less variable than those in nonglaciated areas. Instead, the low intrapopulation and intraregional diversity estimates in D. groenlandicus are probably a result of regional bottleneck events due to range contractions during Holocene warming events. These results are consistent with findings previously reported on collared lemmings (D. torquatus) from the Eurasian Arctic.

Common dynamic structure of canada lynx populations within three climatic regions

Across the boreal forest of Canada, lynx populations undergo regular density cycles. Analysis of 21 time series from 1821 onward demonstrated structural similarity in these cycles within large regions of Canada. The observed population dynamics are consistent with a regional structure caused by climatic features, resulting in a grouping of lynx population dynamics into three types (corresponding to three climatic-based geographic regions): Pacific-maritime, Continental, and Atlantic-maritime. A possible link with the North Atlantic Oscillation is suggested.

Cycles and trends in cod populations

Year-to-year fluctuations in fish stocks are usually attributed to variability in recruitment, competition, predation, and changes in catchability. Trends in abundance, in contrast, are usually ascribed to human exploitation and large-scale environmental changes. In this study, we demonstrate, through statistical modeling of survey data (1921-1994) of cod from the Norwegian Skagerrak coast, that both short- and long-term variability may arise from the same set of age-structured interactions. Asymmetric competition and cannibalism between cohorts generate alternating years of high and low abundance. Intercohort interactions also resonate the recruitment variability so that long-term trends are induced. The coupling of age-structure and variable recruitment should, therefore, be considered when explaining both the short- and long-term fluctuations displayed by the coastal cod populations. Resonant effects may occur in many marine populations that exhibit this combination of traits.

Environmental variation shapes sexual dimorphism in red deer

Sexual dimorphism results from dichotomous selection on male and female strategies of growth in relation to reproduction. In polygynous mammals, these strategies reflect sexual selection on males for access to females and competitive selection on females for access to food. Consequently, in such species, males display rapid early growth to large adult size, whereas females invest in condition and early sexual maturity at the expense of size. Hence, the magnitude of adult size dimorphism should be susceptible to divergence of the sexes in response to environmental factors differentially influencing their growth to reproduction. We show that divergent growth of male and female red deer after 32 years of winter warming and 15 years of contemporaneously earlier plant phenology support this prediction. In response to warmer climate during their early development, males grew more rapidly and increased in size, while female size declined. Conversely, females, but not males, responded to earlier plant phenology with increased investment in condition and earlier reproduction. Accordingly, adult size dimorphism increased in relation to warmer climate, whereas it declined in relation to forage quality. Thus, the evolutionary trajectories of growth related to reproduction in the sexes (i) originate from sexual and competitive selection, (ii) produce sexual size dimorphism, and (iii) are molded by environmental variation.

Interspecific and intraspecific competition as causes of direct and delayed density dependence in a fluctuating vole population

A 3- to 5-year cycle of vole abundances is a characteristic phenomenon in the ecology of northern regions, and their explanation stands as a central theoretical challenge in population ecology. Although many species of voles usually coexist and are in severe competition for food and breeding space, the role of interspecific competition in vole cycles has never been evaluated statistically. After studying community effects on the population dynamics of the gray-sided vole (Clethrionomys rufocanus) in the subarctic birch forest at Kilpisjärvi, Finland, we report statistical results showing that both interspecific and intraspecific effects are important in the direct year-to-year density dependence. However, interspecific effects are not detectable in the 2-year delayed density dependence that is crucial for generating the characteristic cycles. Furthermore, we show that most of the competition takes place during the winter. The results are evaluated against two models of community dynamics. One assumes that the delayed effects are caused by an interaction with a specialist predator, and the other assumes that they are caused by overgrazing food plants. These statistical results show that vole cycles may be generated by a species-specific trophic interaction. The results also suggest that the gray-sided vole may be the focal species in the birch-forest community, as field voles may be in the taiga and as lemmings may be on the tundra.

From patterns to processes: phase and density dependencies in the Canadian lynx cycle

Across the boreal forest of North America, lynx populations undergo 10-year cycles. Analysis of 21 time series from 1821 to the present demonstrates that these fluctuations are generated by nonlinear processes with regulatory delays. Trophic interactions between lynx and hares cause delayed density-dependent regulation of lynx population growth. The nonlinearity, in contrast, appears to arise from phase dependencies in hunting success by lynx through the cycle. Using a combined approach of empirical, statistical, and mathematical modeling, we highlight how shifts in trophic interactions between the lynx and the hare generate the nonlinear process primarily by shifting functional response curves during the increase and the decrease phases.

Phase- and density-dependent population dynamics in Norwegian lemmings: interaction between deterministic and stochastic processes

We analysed two 26-year long (1970-1995) time-series on annual population growth rates of Norwegian lemmings (Lemmus lemmus) from Finse, south Norway, using a threshold autoregressive (TAR) approach. We demonstrate that the population dynamics is both phase- and density-dependent. The phase-dependence accounts for the observed nonlinearity. We used the deduced stochastic model structure as a basis for evaluating the dynamic properties of this system. The dynamics is characterized either by limit cycles or chaos (the latter with a strong semi-periodic component). Stochasticity is seen to play an important role in the determination of the periodicity. The ecological implications of these statistical and mathematical results are discussed.

Population dynamics of Norwegian red deer: density-dependence and climatic variation

We present a model on plant-deer climate interactions developed for improving our understanding of the temporal dynamics of deer abundance and, in particular, how intrinsic (density-dependent) and extrinsic (plants, climate) factors influence these dynamics. The model was tested statistically by analysing the dynamics of five Norwegian red deer populations between 1964 and 1993. Direct and delayed density-dependence significantly influenced the development of the populations: delayed density-dependence primarily operated through female density, whereas direct density-dependence acted through both female and male densities. Furthermore, population dynamics of Norwegian red deer were significantly affected by climate (as measured by the global weather phenomenon, the North Atlantic Oscillation: NAO). Warm, snowy winters (high NAO) were associated with decreased deer abundance, whereas the delayed (two-year) effect of warm, snowy winters had a positive effect on deer abundance. Our analyses are argued to have profound implications for the general understanding of climate change and terrestrial ecosystem functioning.

Global climate change and phenotypic variation among red deer cohorts

The variability of two fitness-related phenotypic traits (body weight and a mandibular skeletal ratio) was analysed among cohorts and age-classes of red deer in Norway. Phenotypic variation among cohorts was pronounced for calves, yearlings and reproductively mature adults. Fluctuations in cohort-specific mean body weights and skeletal ratios of adults correlated with global climatic variation in winter conditions influenced by the North Atlantic Oscillation while cohorts were in utero. Red deer born following warm winters were smaller than those born after cold winters, and this inter-cohort variability persisted into adulthood. Phenotypic variation among cohorts of red deer influenced by climate change may pose consequences for fitness of cohorts since body size and condition contribute to reproductive success and survival in male and female red deer. In particular, the recent trend of increasingly warm winters in northern Europe and Scandinavia may lead to reduced body size and fecundity of red deer, and perhaps other ungulates, in those areas.

Stochastic seasonality and nonlinear density-dependent factors regulate population size in an African rodent

Ecology has long been troubled by the controversy over how populations are regulated. Some ecologists focus on the role of environmental effects, whereas others argue that density-dependent feedback mechanisms are central. The relative importance of both processes is still hotly debated, but clear examples of both processes acting in the same population are rare. Key-factor analysis (regression of population changes on possible causal factors) and time-series analysis are often used to investigate the presence of density dependence, but such approaches may be biased and provide no information on actual demographic rates. Here we report on both density-dependent and density-independent effects in a murid rodent pest species, the multimammate rat Mastomys natalensis (Smith, 1834), using statistical capture-recapture models. Both effects occur simultaneously, but we also demonstrate that they do not affect all demographic rates in the same way. We have incorporated the obtained estimates of demographic rates in a population dynamics model and show that the observed dynamics are affected by stabilizing nonlinear density-dependent components coupled with strong deterministic and stochastic seasonal components.

Population regulation in snowshoe hare and Canadian lynx: asymmetric food web configurations between hare and lynx

The snowshoe hare and the Canadian lynx in the boreal forests of North America show 9- to 11-year density cycles. These are generally assumed to be linked to each other because lynx are specialist predators on hares. Based on time series data for hare and lynx, we show that the dominant dimensional structure of the hare series appears to be three whereas that of the lynx is two. The three-dimensional structure of the hare time series is hypothesized to be due to a three-trophic level model in which the hare may be seen as simultaneously regulated from below and above. The plant species in the hare diet appear compensatory to one another, and the predator species may, likewise, be seen as an internally compensatory guild. The lynx time series are, in contrast, consistent with a model of donor control in which their populations are regulated from below by prey availability. Thus our analysis suggests that the classic view of a symmetric hare-lynx interaction is too simplistic. Specifically, we argue that the classic food chain structure is inappropriate: the hare is influenced by many predators other than the lynx, and the lynx is primarily influenced by the snowshoe hare.

Is spacing behaviour coupled with predation causing the microtine density cycle? A synthesis of current process-oriented and pattern-oriented studies

Current ecological information on periodically fluctuating microtine populations are demonstrated to support a hypothesis involving both predation and intrinsic self-regulation as necessary and sufficient factors for explaining the "microtine density cycle'. The structure of the cyclic time series is largely two dimensional with strong delayed density dependence. Together with recent field studies on rodent demography, our modelling suggests that trophic interaction is a likely candidate to generate the dimensionality observed for northern microtine rodent dynamics. It is shown that the trophic interaction must be fairly strong. This suggests that specialist predation is the most likely one among the classes of trophic interactions. We also argue that some - but not too strong - self-regulation must occur to generate the structure of the available time series on northern European microtines.

A gradient from stable to cyclic populations of Clethrionomys rufocanus in Hokkaido, Japan

A total of 31 years of abundance data from 90 populations of the grey-sided vole (Clethrionomys rufocanus) in northern Hokkaido (Japan) were analysed with respect to population dynamic characteristics. Both non-periodic and multi-annually periodic fluctuations occur among the studied populations. The length of the period varies from 2 to 5 years. The dynamics appear approximately linear on a logarithmic scale, and a log-linear stochastic difference model with one time-lag is found to recreate the periodograms of the time series. The deterministic (log-linear) component of these models determines the periodicity while stochasticity sustains the cycles. The dynamics of the Hokkaidian vole populations are found to vary clinally from the western coast eastwards and towards the interior of the study area. This gradient corresponds superficially to the latitudinal gradient seen in microtine populations in Fennoscandia. However, under close scrutiny, these gradients differ greatly: the Hokkaidian gradient is caused by a cline in delayed density dependence. Statistical delayed density dependence is more negative towards the east and the interior.

A geographic gradient in small rodent density fluctuations: a statistical modelling approach

The patterns of density dependence in Fennoscandian rodents are investigated statistically using a linear autoregressive scheme. Nineteen time series of microtine abundances along a latitudinal gradient in Fennoscandia from 60 degrees N to 69 degrees N are analysed. We provide statistical evidence that there exists a latitudinal gradient in density dependence in Fennoscandian microtines. Southern populations experience significantly stronger direct density dependence than northern populations. Delayed density dependence was significantly negative throughout the region and appeared constant across the latitudinal gradient. The populations consistently exhibit dynamics of second order throughout the region. Together, the clinal direct density dependence and constant delayed density dependence give rise to a cline in cycle period from 3 to 4.5 years. The statistical results are compared to assumptions and predictions made in previous studies on the geographic gradient in the population dynamics of these rodents. The results are in agreement with the predictions of the 'generalist predator hypothesis'.

Bootstrap estimated uncertainty of the dominant Lyapunov exponent for Holarctic microtine rodents

The dominant Lyapunov exponent, as estimated from time series using the Jacobian-based method, is often used for indicating whether the underlying dynamic system is chaotic or not. The Jacobian-based method together with Response Surface Methodology has been suggested as a method for detecting chaotic dynamics in ecological time series. Besides pointing out that this may not be an appropriate method, we report on estimates of the uncertainty in the estimates of the dominant Lyapunov exponent. For this purpose, we have used time series data on Holarctic microtines. On the basis of our analyses, we are unable to find general evidence for chaotic dynamics in northern microtine populations (north of 60 degrees N) as recently suggested in the ecological literature. It seems, however, that the dynamic properties of the northern and southern populations are different. These patterns are supported by testing for nonlinearity.

Fingerprinting of diverse species with DNA probes generated from immobilized single-stranded DNA templates

Images

Age-specific optimal diets and optimal foraging tactics: a life-historic approach

By linking optimal foraging theory and optimal life history theory, we demonstrate that optimal diets, in general, may depend on the individual's age even when everything else remains the same. Older individuals (i.e., individuals with lower reproductive values) are predicted to have diets composed of highly nutritious food types that are possibly dangerous to pursue.