Effects of host abundance on larch budmoth outbreaks in the European Alps

相似文献   

Summer temperature dependency of larch budmoth outbreaks revealed by Alpine tree-ring isotope chronologies

Larch budmoth (LBM, Zeiraphera diniana Gn.) outbreaks cause discernable physical alteration of cell growth in tree rings of host subalpine larch (Larix decidua Mill.) in the European Alps. However, it is not clear if these outbreaks also impact isotopic signatures in tree-ring cellulose, thereby masking climatic signals. We compared LBM outbreak events in stable carbon and oxygen isotope chronologies of larch and their corresponding tree-ring widths from two high-elevation sites (1800–2200 m a.s.l.) in the Swiss Alps for the period AD 1900–2004 against isotope data obtained from non-host spruce (Picea abies). At each site, two age classes of tree individuals (150–250 and 450–550 years old) were sampled. Inclusion of the latter age class enabled one chronology to be extended back to AD 1650, and a comparison with long-term monthly resolved temperature data. Within the constraints of this local study, we found that: (1) isotopic ratios in tree rings of larch provide a strong and consistent climatic signal of temperature; (2) at all sites the isotope signatures were not disturbed by LBM outbreaks, as shown, for example, by exceptionally high significant correlations between non-host spruce and host larch chronologies; (3) below-average July to August temperatures and LBM defoliation events have been coupled for more than three centuries. Dampening of Alps-wide LBM cyclicity since the 1980s and the coincidence of recently absent cool summers in the European Alps reinforce the assumption of a strong coherence between summer temperatures and LBM defoliation events. Our results demonstrate that stable isotopes in tree-ring cellulose of larch are an excellent climate proxy enabling the analysis of climate-driven changes of LBM cycles in the long term.     

Dendroecological reconstruction and interpretation of larch budmoth (Zeiraphera diniana) outbreaks in two central alpine valleys of Switzerland from 1470 – 1990

 Outbreaks of the larch budmoth (LBM) (Zeiraphera diniana) recur cyclically approximately every 7 to 10 years in subalpine larch-cembran pine and montane to subalpine larch-Norway spruce forests of the relatively dry valleys of the European Alps. By dendroecologically analyzing increment cores from 570 host (European larch –Larix decidua) and non-host trees (cembran pine –Pinus cembra, Norway spruce –Picea abies) through the use of skeleton plots, at least 57 (59) outbreaks could be reconstructed in the optimum Upper Engadine Valley (suboptimum Goms Valley), Switzerland, during the time period 1503 (1472) to 1990. The average interval between initial years of successive outbreaks was 8.58 (8.95) years, SD 1.66 (2.13) years. Over the centuries spatial shifts of LBM activity between the two study areas occurred, probably due to climatic changes. Clear, site-specific differences in LBM attack could only be found in the suboptimum area where high-lying (>1800 m) and/or south-facing stands were infested most. LBM-afflicted trees proved to be unsuitable for climate reconstructions because the impact of the persistently recurring outbreaks on tree growth is dominant. In order to provide sufficient information for a detailed ecological interpretation of the course of an outbreak, latewood widths and/or densities have to be analyzed in addition to the ring-widths. Received: 11 February 1995 / Accepted: 19 July 1996     

Tracing the influence of larch-bud-moth insect outbreaks and weather conditions on larch tree-ring growth in Engadine (Switzerland).

Tree-ring analysis of insect-defoliated trees has so far been used for detecting past insect outbreaks only. We hypothesize that the impact of the larch-bud-moth (LBM) Zeiraphera dinian outbreaks on the growth of larch Larix decidua in the Engadine (Switzerland) is closely coupled to the spatial development of the outbreak and the ecological characteristics of the respective sites and weather conditions. We tested this hypothesis by reviewing data sets available in the literature and by analysing original data. We monitored LBM population densities and the needle phenology, growth and defoliation of larch over 28 years, i.e. over four outbreak cycles. In addition, information on defoliation patterns covering six earlier outbreaks over 50 years was matched with tree-ring information. Tree-ring chronologies of 18 larch stands were analysed with regard to abrupt growth changes and latewood events. Defoliation induces an immediate reduction in latewood, followed by a reduction in needle length and a significant decrease in radial growth in the subsequent year. We have called this tree-ring pattern the "larch-bud-moth syndrome". A careful analysis of the various parameters of the LBM syndrome for two specific population cycles enabled us to define different interaction patterns between weather conditions and tree growth. These can then be included in climate change models to help disentangle the impact of insect defoliation from that of adverse climatic conditions.     

Return of the moth: rethinking the effect of climate on insect outbreaks

The sudden interruption of recurring larch budmoth (LBM; Zeiraphera diniana or griseana Gn.) outbreaks across the European Alps after 1982 was surprising, because populations had regularly oscillated every 8–9 years for the past 1200 years or more. Although ecophysiological evidence was limited and underlying processes remained uncertain, climate change has been indicated as a possible driver of this disruption. An unexpected, recent return of LBM population peaks in 2017 and 2018 provides insight into this insect’s climate sensitivity. Here, we combine meteorological and dendrochronological data to explore the influence of temperature variation and atmospheric circulation on cyclic LBM outbreaks since the early 1950s. Anomalous cold European winters, associated with a persistent negative phase of the North Atlantic Oscillation, coincide with four consecutive epidemics between 1953 and 1982, and any of three warming-induced mechanisms could explain the system’s failure thereafter: (1) high egg mortality, (2) asynchrony between egg hatch and foliage growth, and (3) upward shifts of outbreak epicentres. In demonstrating that LBM populations continued to oscillate every 8–9 years at sub-outbreak levels, this study emphasizes the relevance of winter temperatures on trophic interactions between insects and their host trees, as well as the importance of separating natural from anthropogenic climate forcing on population behaviour.     

The potential for using Larix decidua ring widths in reconstructions of larch budmoth (Zeiraphera diniana) outbreak history: dendrochronological estimates compared with insect surveys

In the Alps, larch (Larix decidua Mill.) forests show periodic discolouration due to larch budmoth (LBM) outbreaks (Zeiraphera diniana Guénée, Lepidoptera: Tortricidae). Tree defoliation causes severe reductions in radial growth, visible in tree rings. This paper aims at reconstructing LBM outbreak history, and critically examining the potential for using dendrochronological data by comparing tree-ring estimates with insect surveys. The occurrence of LBM outbreaks was investigated using 249 cores from larch growing near the timberline in three regions of the French Alps (Briançonnais, Maurienne, and Tarentaise). Years with an abrupt decrease in radial growth (-40%) were considered as negative pointer years reflecting the potential impact of LBM. The comparison with three non-host conifers (Norway spruce, stone and mountain pines) permitted us to distinguish growth reductions in larch due to climatic effects from those due to defoliation by LBM. The dendrochronological data were matched with information reporting conspicuous discolouration in old forestry reports or recovered from systematic field observations. Twenty-two outbreaks are discernible within the period 1800-1983. A peak reduction in larch growth occurred at intervals of 8.86ǃ.01 years in the following years: 1801, 1811, 1820, 1830, 1838, 1846, 1857, 1867, 1875, 1884, 1892, 1901, 1909, 1918, 1926, 1936, 1945, 1953, 1963, 1972, 1980, and 1987. Our objective method based on ring measurements and comparison with non-host trees was compared with qualitative techniques based on the visual detection of conspicuous latewood anomalies. Larch in the Briançonnais (which experiences a continental climate) are infested first, whereas the Tarentaise region exhibits a much weaker impact of LBM. Complete tree recovery was observable 3 years after outbreak peaks.     

Larch- and pine-feeding host races of the larch bud moth (Zeiraphera diniana) have cyclic and synchronous population fluctuations

Population cycles of many forest-defoliating insects often show synchronous fluctuations at both intra-specific (spatial synchrony) and inter-specific levels. However, population dynamics of different host-associated biotypes of the same species, such as those of the larch budmoth (LBM), Zeiraphera diniana (Lepidoptera: Tortricidae), have never been compared. This species causes extensive defoliation of larch forests every 8 to 9 years in the Alps, but it consists of two genetically-differentiated host races, the first one developing on European larch, Larix decidua , and the other one developing on Swiss stone pine, Pinus cembra . The dynamics of Zeiraphera populations have been extensively studied on larch, whereas little is known about larval density and possible population fluctuations on sympatric pines. A larval census on Swiss stone pine was conducted in the Swiss Alps intermittently between 1958 to 2004 and in the French Alps from 1992 to 2004. Population density of Zeiraphera on pine varied up to 5000-fold and showed periodic oscillations, with five peaks in Switzerland and one in France. Because the feeding activity of the pine race is restricted to the elongating shoot of the current year, no conspicuous defoliation of pine trees was noted during years of high larval densities. Zeiraphera populations on pine oscillated in significant synchrony with larch-associated populations, and peak densities were observed either the same year or shifted by±one year. Our results did not allow any explanation for cyclic fluctuations of LBM on pine, but the synchrony with the larch race's cycle suggests that studies on genetics as well as on parasitism should be intensified.     

The role of larch budmoth (Zeiraphera diniana Gn.) on forest succession in a larch (Larix decidua Mill.) and Swiss stone pine (Pinus cembra L.) stand in the Susa Valley (Piedmont, Italy)

In the Alps, larch (Larix decidua Mill.) is severely affected by larch budmoth (Zeiraphera diniana Guénée) (LBM) attacks. The impact of these outbreaks on the Swiss stone pine (Pinus cembra L.) and on the dynamic processes acting in subalpine forest stands are still not well known. Dendroecological methods were used in this study to reconstruct past LBM outbreaks in Susa Valley, Piedmont, Italy. The analysis was carried out on 62 cores from larch and 101 cores from stone pine. The length and severity of each outbreak was quantified for both species and for each tree by means of the programme OUTBREAK. The frequency of the outbreaks was determined using singular spectral analysis and superposed epoch analysis was used to test the significance of the associations between outbreaks and tree-ring growth. In order to verify if trees belonging to different age classes are differently affected by LBM, the reconstructed outbreaks are then classified taking into account the cambial age of the tree at the moment of the outbreak. From 1760 to 1999, 19 outbreaks were recorded in the larch chronologies, while only three outbreaks in the stone pine chronologies. The larch growth is strongly influenced by LBM and the identified outbreaks are equally distributed in all the three age classes. On the stone pine the sporadic occurrence of the identified events made difficult any interpretation of the eventual effect of LBM. Our results lead us to argue that LBM has not played an important role both in determining the stone pine growth rate and in influencing the present observed succession from the stage dominated by larch, to a stage dominated by stone pine or by a mixed stone pine-larch forest.     

Consequences of larch budmoth outbreaks on the climatic significance of ring width and stable isotopes of larch

Tree-ring widths and stable carbon and oxygen isotopes of five European larch trees from Lötschental, Switzerland were investigated for the period 1900–2004. The objective was to test the suitability of each of these parameters for high-frequency climate reconstructions. This is of special interest with regard to the problem of cyclic larch budmoth (LBM) infestations of alpine larch trees. The results clearly demonstrate that tree-ring width chronologies are not suitable for high-frequency reconstructions because infestations lead to variably reduced tree-ring increments, largely suppressing climate signals. On the other hand, the stable isotope chronologies proved less affected by larch budmoth outbreaks, independent of the strength of the infestations. The correlation of the carbon isotopes with summer temperatures was especially high (r = 0.73) and with precipitation lower but nevertheless significant (r = ?0.43). Oxygen isotopes were also correlated with summer temperature (r = 0.46); however, a certain perturbation of normal oxygen isotope signatures due to LBM outbreaks was evident. Contrary to tree-ring widths, none of the LBM outbreaks caused a significant disturbance of the current year’s isotopic climate signal and, most importantly, there were no delayed effects in the following years. Thus, stable carbon isotopes in tree-ring chronologies of the European larch provide an excellent opportunity for high-frequency temperature reconstructions.     

European larch phenology in the Alps: can we grasp the role of ecological factors by combining field observations and inverse modelling?

Vegetation phenology is strongly influenced by climatic factors. Climate changes may cause phenological variations, especially in the Alps which are considered to be extremely vulnerable to global warming. The main goal of our study is to analyze European larch (Larix decidua Mill.) phenology in alpine environments and the role of the ecological factors involved, using an integrated approach based on accurate field observations and modelling techniques. We present 2 years of field-collected larch phenological data, obtained following a specifically designed observation protocol. We observed that both spring and autumn larch phenology is strongly influenced by altitude. We propose an approach for the optimization of a spring warming model (SW) and the growing season index model (GSI) consisting of a model inversion technique, based on simulated look-up tables (LUTs), that provides robust parameter estimates. The optimized models showed excellent agreement between modelled and observed data: the SW model predicts the beginning of the growing season (B(GS)) with a mean RMSE of 4 days, while GSI gives a prediction of the growing season length (L(GS)) with a RMSE of 5 days. Moreover, we showed that the original GSI parameters led to consistent errors, while the optimized ones significantly increased model accuracy. Finally, we used GSI to investigate interactions of ecological factors during springtime development and autumn senescence. We found that temperature is the most effective factor during spring recovery while photoperiod plays an important role during autumn senescence: photoperiod shows a contrasting effect with altitude decreasing its influence with increasing altitude.     

Endurance of larch forest ecosystems in eastern Siberia under warming trends

The larch (Larix spp.) forest in eastern Siberia is the world's largest coniferous forest. Its persistence is considered to depend on near‐surface permafrost, and thus, forecast warming over the 21st century and consequent degradation of near‐surface permafrost is expected to affect the larch forest in Siberia. However, predictions of these effects vary greatly, and many uncertainties remain about land – atmosphere interactions within the ecosystem. We developed an integrated land surface model to analyze how the Siberian larch forest will react to current warming trends. This model analyzed interactions between vegetation dynamics and thermo‐hydrology, although it does not consider many processes those are considered to affect productivity response to a changing climate (e.g., nitrogen limitation, waterlogged soil, heat stress, and change in species composition). The model showed that, under climatic conditions predicted under gradual and rapid warming, the annual net primary production of larch increased about 2 and 3 times, respectively, by the end of the 21st century compared with that in the previous century. Soil water content during the larch‐growing season showed no obvious trend, even when surface permafrost was allowed to decay and result in subsurface runoff. A sensitivity test showed that the forecast temperature and precipitation trends extended larch leafing days and reduced water shortages during the growing season, thereby increasing productivity. The integrated model also satisfactorily reconstructed latitudinal gradients in permafrost presence, soil moisture, tree leaf area index, and biomass over the entire larch‐dominated area in eastern Siberia. Projected changes to ecosystem hydrology and larch productivity at this geographical scale were consistent with those from site‐level simulation. This study reduces the uncertainty surrounding the impact of current climate trends on this globally important carbon reservoir, and it demonstrates the need to consider complex ecological processes to make accurate predictions.     

Anomalous outbreaks of an invasive defoliator and native bark beetle facilitated by warm temperatures,changes in precipitation and interspecific interactions

Biotic disturbance agents such as insects can be highly responsive to climatic change and have widespread ecological and economic impacts on forests. Quantifying the responses of introduced and native insects to climate, including how dynamics of one agent may mediate those of another, is important for forecasting disturbance and associated impacts on forest structure and function. We investigated drivers of outbreaks by larch casebearer Coleophora laricella, an invasive defoliator, and eastern larch beetle Dendroctonus simplex, a native, tree‐killing bark beetle, on tamarack Larix laricina from 2000 to in Minnesota, USA. We evaluated the utility of temporal, spatial and climatic variables in predicting the presence/absence of outbreaks of each insect in cells of rasterized aerial survey data. The role of defoliation by larch casebearer in outbreaks of eastern larch beetle was also investigated. For both species, the most important predictors of outbreak occurrence were proximity of conspecific outbreaks in space and time. For larch casebearer, outbreak occurrence was positively associated with spring precipitation and warmer growing seasons. Outbreak occurrence of eastern larch beetle was positively associated with warmer and dryer years and was more likely in cells with prior defoliation by larch casebearer. Our results demonstrate that climate can drive large scale outbreaks of introduced and non‐native disturbance agents on a single host species, and that interactions at the tree level between such agents may scale up to manifest across large temporal and spatial scales.     

Changes in the population structure and sporulation behaviour of <Emphasis Type="Italic">Phytophthora ramorum</Emphasis> associated with the epidemic on <Emphasis Type="Italic">Larix</Emphasis> (larch) in Britain

During a decade of invasion, the exotic pathogen Phytophthora ramorum has undergone an unexpected change in behaviour during spread into the woodlands and forests of Great Britain. From 2002 to 2008 most outbreaks centred on nurseries and managed gardens with affected hosts almost exclusively broadleaf shrubs and trees. However 2009 saw a major shift as larch tree plantations (Larix) were affected by widespread infection and mortality incited by P. ramorum. To understand the processes underlying the host jump to larch, isolates of the EU1 lineage of P. ramorum collected from 2002 to 2012 were investigated using seven polymorphic microsatellite markers. Analysis of 347 isolates resolved 51 multilocus genotypes (MLGs) which partitioned into two distinct clusters. One comprised MLGs unique to Britain and unknown elsewhere in Europe, the other cluster was primarily of MLGs already known in other European countries but dominated by one genotype, EU1MG1. Pre-2009 isolates were predominantly of the unique British cluster with only a few typical of the European cluster. This reversed after 2009 with European MLGs, especially EU1MG1, becoming increasingly common as the larch epidemic expanded. We hypothesise that the growing dominance of EU1MLG1 has been an important driver in the emergence of the epidemic on larch, aided by its ability to sporulate more abundantly compared with the dominant unique British MLG. European MLGs appear closely associated with the distribution of larch along the west coast of Britain whereas unique British MLGs tend to be concentrated in south west England. The two population clusters suggest at least two separate introductions of the EU1 lineage into Britain with subsequent diversification.     

Landscape geometry and travelling waves in the larch budmoth

Travelling waves in cyclic populations refer to temporal shifts in peak densities moving across space in a wave‐like fashion. The epicentre hypothesis states that peak densities begin in specific geographic foci and then spread into adjoining areas. Travelling waves have been confirmed in a number of population systems, begging questions about their causes. Herein we apply a newly developed statistical technique, wavelet phase analysis, to historical data to document that the travelling waves in larch budmoth (LBM) outbreaks arise from two epicentres, both located in areas with high concentrations of favourable habitat. We propose that the spatial arrangement of the landscape mosaic is responsible for initiating the travelling waves. We use a tri‐trophic model of LBM dynamics to demonstrate that landscape heterogeneity (specifically gradients in density of favourable habitat) alone, is capable of inducing waves from epicentres. Our study provides unique evidence of how landscape features can mould travelling waves.     

Waves and synchrony in Epirrita autumnata/Operophtera brumata outbreaks. I. Lagged synchrony: regionally, locally and among species

1. In 1990-2003, during a complete 10-year outbreak cycle, the synchrony of the birch defoliating outbreaks of the geometrids Epirrita autumnata and Operophtera brumata was studied quantitatively in the northern part of the Fennoscandian mountain chain (the Scandes). Data were supplemented with similar data from 1964 to 1966 and historical information. A 30-year series of field data from one locality in southern Scandes made possible interregional comparisons. 2. In 1991, outbreaks started in north-eastern Fennoscandia and moved westward like a wave and reached the outer coast of north-western Norway in about 2000. This wave is a new observation. In the same years, a previously documented outbreak wave moved southward along the Scandes. 3. Outbreak periods have usually occurred around the middle of each decade. Seemingly unrelated population peaks at the decadal shift 2000 were reported from islands at the coast of north-western Norway. They are shown here to have been the final ripples of the east-west wave. 4. At some localities, O. brumata peaked 2 years after E. autumnata. A lag of 1 or 2 years also occurred at the locality in southern Scandes. This interspecific time lag is a new observation. In accordance with the north-south wave, a time-lag of 1-2 years occurred between the fluctuations of northern and southern E. autumnata and O. brumata populations. 5. The population peak of E. autumnata occurred 1 year earlier at one locality than at a nearby locality. This pattern and particular altitudinal shifts of the O. brumata population density at these localities repeated in two outbreak periods. This indicates that, for example, local climate may modify outbreak synchrony between nearby localities. 6. At the same localities, O. brumata peaked first at one altitude and 1 or 2 years later at another altitude. This vertical lag is a new observation. 7. E. autumnata shows fluctuation traits similar to some other cyclic animals, e.g. the larch budmoth in the European Alps, some European tetraonid birds and the Canadian snow-shoe hare. These similarities (and dissimilarities) in intra- and interspecific synchronies and causes of E. autumnata and O. brumata synchronies, regionally, locally and among the two species are discussed.     

Molecular epidemiology and population structure of the honey bee brood pathogen Melissococcus plutonius

Melissococcus plutonius is the causative agent of European foulbrood (EFB), which is a serious brood disease of the European honey bee (Apis mellifera). EFB remains a threat because of a poor understanding of disease epidemiology. We used a recently published multi-locus sequence typing method to characterise 206 M. plutonius isolates recovered from outbreaks in England and Wales over the course of 2 years. We detected 15 different sequence types (STs), which were resolved by eBURST and phylogenetic analysis into three clonal complexes (CCs) 3, 12 and 13. Single and double locus variants within CC3 were the most abundant and widespread genotypes, accounting for 85% of the cases. In contrast, CCs 12 and 13 were rarer and predominantly found in geographical regions of high sampling intensity, consistent with a more recent introduction and localised spread. K-function analysis and interpoint distance tests revealed significant geographical clustering in five common STs, but pointed to different dispersal patterns between STs. We noted that CCs appeared to vary in pathogenicity and that infection caused by the more pathogenic variants is more likely to lead to honey bee colony destruction, as opposed to treatment. The importance of these findings for improving our understanding of disease aetiology and control are discussed.     

Climate change and outbreaks of the geometrids Operophtera brumata and Epirrita autumnata in subarctic birch forest: evidence of a recent outbreak range expansion

1. Range expansions mediated by recent climate warming have been documented for many insect species, including some important forest pests. However, whether climate change also influences the eruptive dynamics of forest pest insects, and hence the ecological and economical consequences of outbreaks, is largely unresolved. 2. Using historical outbreak records covering more than a century, we document recent outbreak range expansions of two species of cyclic geometrid moth, Operophtera brumata Bkh. (winter moth) and Epirrita autumnata L. (autumnal moth), in subarctic birch forest of northern Fennoscandia. The two species differ with respect to cold tolerance, and show strikingly different patterns in their recent outbreak range expansion. 3. We show that, during the past 15-20 years, the less cold-tolerant species O. brumata has experienced a pronounced north-eastern expansion into areas previously dominated by E. autumnata outbreaks. Epirrita autumnata, on the other hand, has expanded the region in which it exhibits regular outbreaks into the coldest, most continental areas. Our findings support the suggestion that recent climate warming in the region is the most parsimonious explanation for the observed patterns. 4. The presence of O. brumata outbreaks in regions previously affected solely by E. autumnata outbreaks is likely to increase the effective duration of local outbreaks, and hence have profound implications for the subarctic birch forest ecosystem.     

Short-term signals of climate change in Italian summit vegetation: observations at two GLORIA sites

Short-term changes occurring in high mountain vegetation were analysed using the data from two Italian sites already part of the GLobal Observation Research Initiative in Alpine environments (GLORIA – central Apennines and southwestern Alps). The study focused on a set of floristic (endemics), structural (life forms) and ecological (thermic vegetation indicator) variables. Vegetation data were collected according to the GLORIA multi-summit standardized method during the last decade. The re-visitation revealed a moderate decrease in regional endemic flora and significant variations in structural and ecological parameters. The increase in caespitose hemicryptophytes in both sites, in suffruticose chamaephytes in the central Apennines and in rosette-forming hemicryptophytes in the southwestern Alps emerged, highlighting the rapid responses of the alpine vegetation to climate warming. The increase in perennial life forms is related with the expansion of graminoids and small woody plants. These life forms seem to be most suitable to face climate warming in Italian summits. The increase in the thermic vegetation indicator exceeds the mean European summits increment, and this is due to the expansion of thermophilic species. Short-term analyses with fine spatial and temporal resolutions are still necessary to improve our understanding concerning species behaviour in high-elevation ecosystems.     

Evolutionary history of the European whitefish Coregonus lavaretus (L.) species complex as inferred from mtDNA phylogeography and gill-raker numbers

We compared mitochondrial DNA and gill-raker number variation in populations of the European whitefish Coregonus lavaretus (L.) species complex to illuminate their evolutionary history, and discuss mechanisms behind diversification. Using single-strand conformation polymorphism (SSCP) and sequencing 528 bp of combined parts of the cytochrome oxidase b (cyt b) and NADH dehydrogenase subunit 3 (ND3) mithochondrial DNA (mtDNA) regions, we documented phylogeographic relationships among populations and phylogeny of mtDNA haplotypes. Demographic events behind geographical distribution of haplotypes were inferred using nested clade analysis (NCA) and mismatch distribution. Concordance between operational taxonomical groups, based on gill-raker numbers, and mtDNA patterns was tested. Three major mtDNA clades were resolved in Europe: a North European clade from northwest Russia to Denmark, a Siberian clade from the Arctic Sea to southwest Norway, and a South European clade from Denmark to the European Alps, reflecting occupation in different glacial refugia. Demographic events inferred from NCA were isolation by distance, range expansion, and fragmentation. Mismatch analysis suggested that clades which colonized Fennoscandia and the Alps expanded in population size 24 500-5800 years before present, with minute female effective population sizes, implying small founder populations during colonization. Gill-raker counts did not commensurate with hierarchical mtDNA clades, and poorly with haplotypes, suggesting recent origin of gill-raker variation. Whitefish designations based on gill-raker numbers were not associated with ancient clades. Lack of congruence in morphology and evolutionary lineages implies that the taxonomy of this species complex should be reconsidered.