Search tips
Search criteria 


Logo of plosonePLoS OneView this ArticleSubmit to PLoSGet E-mail AlertsContact UsPublic Library of Science (PLoS)
PLoS One. 2010; 5(4): e10026.
Published online 2010 April 2. doi:  10.1371/journal.pone.0010026
PMCID: PMC2848865

Quantifying Recent Ecological Changes in Remote Lakes of North America and Greenland Using Sediment Diatom Assemblages

Richard Kazimierz Frank Unsworth, Editor



Although arctic lakes have responded sensitively to 20th-century climate change, it remains uncertain how these ecological transformations compare with alpine and montane-boreal counterparts over the same interval. Furthermore, it is unclear to what degree other forcings, including atmospheric deposition of anthropogenic reactive nitrogen (Nr), have participated in recent regime shifts. Diatom-based paleolimnological syntheses offer an effective tool for retrospective assessments of past and ongoing changes in remote lake ecosystems.

Methodology/Principal Findings

We synthesized 52 dated sediment diatom records from lakes in western North America and west Greenland, spanning broad latitudinal and altitudinal gradients, and representing alpine (n = 15), arctic (n = 20), and forested boreal-montane (n = 17) ecosystems. Diatom compositional turnover (β-diversity) during the 20th century was estimated using Detrended Canonical Correspondence Analysis (DCCA) for each site and compared, for cores with sufficiently robust chronologies, to both the 19th century and the prior ~250 years (Little Ice Age). For both arctic and alpine lakes, β-diversity during the 20th century is significantly greater than the previous 350 years, and increases with both latitude and altitude. Because no correlation is apparent between 20th-century diatom β-diversity and any single physical or limnological parameter (including lake and catchment area, maximum depth, pH, conductivity, [NO3], modeled Nr deposition, ambient summer and winter air temperatures, and modeled temperature trends 1948–2008), we used Principal Components Analysis (PCA) to summarize the amplitude of recent changes in relationship to lake pH, lake:catchment area ratio, modeled Nr deposition, and recent temperature trends.


The ecological responses of remote lakes to post-industrial environmental changes are complex. However, two regions reveal concentrations of sites with elevated 20th-century diatom β-diversity: the Arctic where temperatures are increasing most rapidly, and mid-latitude alpine lakes impacted by high Nr deposition rates. We predict that remote lakes will continue to shift towards new ecological states in the Anthropocene, particularly in regions where these two forcings begin to intersect geographically.


There is mounting evidence that recent ecological and biogeochemical changes have occurred in remote lakes, defined here as those lacking any immediate, catchment-scale, anthropogenic influences. The implication is that these ecosystems can no longer be considered pristine, largely because of their high sensitivity to climate change [1], [2] which, in some regions, is compounded by significant inputs of reactive anthropogenic nitrogen (hereafter Nr; comprising all biologically-, photochemically- and radiatively- active nitrogenous compounds in the atmosphere and biosphere [3]) delivered by atmospheric deposition [4][6]. Part of the sensitivity of arctic and alpine lakes is attributable to limnological characteristics including dilute water chemistry, low primary production, and high flushing rates associated with nival hydrological regimes. The paucity of long-term climatic and environmental monitoring data in remote regions can be alleviated by the use of proxy data from high-resolution sedimentary records. Geochemical signatures and biological remains are continuously archived in the sediments accumulating at the bottoms of lakes. Diatoms (Bacillariophyceae) are unicellular aquatic photoautotrophs that respond rapidly to changes in water chemistry mediated by environmental change. Their siliceous cell walls are often preserved in lake sediments with sufficient fidelity to allow taxonomic identifications and ecological inferences, thus producing an archive of limnological history. In an evolving literature, sediment diatom records have been applied successfully to a range of global environmental issues including lake acidification [7], eutrophication [8], and climate change [2], [9]. Many of these data augment the evidence that the planet has entered the Anthropocene [10], [11], the era of human dominance over key biogeochemical cycles, with direct climatic and ecological repercussions.

In the present study we analyze the amount of diatom compositional turnover in sediments from 52 remote lakes in North America and west Greenland during three time intervals, each defined broadly by climate history: (1) the cold Little Ice Age (LIA), ~1550–1800; (2) recovery following culmination of the LIA, 1800–1900; and (3) 20th century warming. Geographically, the lakes span 38–79°N of latitude, and 12–3546 m asl of altitude, representing arctic (n = 20), alpine (n = 15), and forested boreal-montane (n = 17) ecosystems (Fig. 1, Table 1). The geographical gradients and temporal intervals captured by this sample array allow the sensitivity of remote lakes to be tested explicitly. The cornerstone of our analysis comprises estimates of diatom compositional turnover (or β-diversity) obtained by Detrended Canonical Correspondence Analysis (DCCA) constrained to time as the sole predictor variable [Materials and Methods, 12].

Figure 1
Location map of the study lakes in North America and Greenland.
Table 1
Locations, physical characteristics, pH, conductivity, and β-diversity values from lakes in this study.


Diatom Stratigraphies

Diatom β-diversity summarizes the amount of compositional turnover having occurred in a core's successive diatom assemblages over a specified interval of time. It is estimated in units of standard deviations (SD). We first illustrate examples of the raw data used in these computations by showing the relative frequencies of dominant diatoms in dated sediments from four lakes expressing variable amounts of assemblage change and representing a range of lake characteristics (Fig. 2). In the Arctic (CF-11, Baffin Island) and alpine lakes (Curator, Jasper National Park, Alberta, and Emerald, Beartooth Wilderness, Wyoming), the increased success of planktonic Cyclotella spp. and Discostella spp. in the 20th century is evident (Fig. 2A–C). Both of these genera become competitive as the water column stratifies, and thus benefit from prolongation of the summer ice-free period [13][15]. During earlier intervals of lowered planktonic diatom abundance, small benthic taxa (e.g. Pseudostaurosira, Staurosirella, Staurosira and Achnanthes spp.) predominate in these lakes. These diatoms are successful when light attenuation is low, ice-free growing seasons are short, and production is concentrated in the littoral zone, often within ice-free moats [2], [9]. In these and other examples, it is the relative success of planktonic diatoms over benthic taxa in the 20th century that drives β-diversity (Fig. 2).

Figure 2
Examples of diatom stratigraphic records.

An additional shift is observed in the sediments of Emerald Lake: after ~1950, Asterionella formosa becomes completely dominant, replacing the Discostella stelligera complex as the principal planktonic diatom. The magnitude of changes in Emerald Lake's diatom flora produced the highest 20th century β-diversity in our survey (2.00 SD). However, recent increases of A. formosa are in no way isolated, being relatively common in alpine lakes of the American Cordillera [4], where they have been linked causally to increase Nr availability from atmospheric deposition [16].

In contrast to lakes that preserve high (CF-11, Curator) to extreme (Emerald) degrees of diatom assemblage compositional turnover, a number of lakes possess relatively complacent diatom stratigraphies. Such lakes are commonly, but not exclusively, from the montane-boreal subset of sites (Table 1). For example, diatom assemblages in sediments from North Barrière Lake (British Columbia interior) have changed very little over recent centuries, resulting in low, near constant β-diversity (0.42–0.53 SD) over the length of the record (Fig. 2D).

Diatom β-Diversity

The estimated 20th century β-diversities from both arctic (p = 0.004) and alpine (p = 0.003) diatom records are significantly higher than those from forested montane-boreal lakes (Fig. 3), as tested using Wilcoxon Rank Sum tests with Bonferonni-adjusted p-values [Materials and Methods]. The alpine sites range from 0.52–2.00 SD (mean: 1.01±0.43), the arctic sites range from 0.44–1.34 SD (mean: 0.84±0.22), and the montane-boreal sites range from 0.24–1.17 SD (mean: 0.63±0.27). However, during the 19th century and the 1550–1800 intervals, β-diversities for arctic (p = 0.29 and p = 0.24, respectively) and alpine (p = 0.71 and p = 0.49, respectively) lakes do not differ significantly from montane-boreal counterparts (Fig. 3). Complete β-diversity results are reported in Table 1.

Figure 3
Boxplots of diatom β-diversity from each lake region.

Additionally, 20th century diatom β-diversities from alpine lakes are significantly greater than either interval of the prior 350 years (p = 0.01 for both time periods). These increase from a mean value of 0.52 SD during the LIA (1550–1800) to 0.61 SD during the 19th century, rising to 1.01 SD in the 20th century (Fig. 3A). A parallel increase in 20th century diatom β-diversity is evident in the arctic lakes, with a lower level of significance (p = 0.05 for both time periods). Here, the succession of mean values is: 0.75 SD (1550–1800 AD); 0.71 SD (19th century); and 0.84 SD (20th century). In contrast, none of the temporal changes in diatom assemblage turnover from montane-boreal lakes represent statistically significant trends (p = 0.76 for the 19th century and p = 0.89 for the LIA).

While we accept this result as a first approximation, we caution that the limited available dating control strongly limits inferences prior to the 20th century from these sites, relative to the arctic and alpine records (Fig. 3, Table 1). Indeed, large changes in diatom assemblages from north-boreal lakes during the 19th century have been documented elsewhere [14], [17].

Summarily, diatom β-diversity for all 52 lakes during the 20th century is significantly greater than the 19th century (p = 0.03), while there is only a small and non-significant difference in turnover between the 19th century and 1550–1800 intervals (p = 0.86). Alpine and arctic lakes reveal greater diatom assemblage turnover in the 20th century relative to forested montane-boreal sites, even though the previous 350 years reveal no significant differences between the three biomes (Fig. 3).

Diatom β-Diversity – Environment Relationships

To investigate potential influences on diatom compositional turnover, we tested relationships between β-diversity and a number of physical and limnological characteristics of the lakes, including: lake and catchment area, maximum depth, pH, conductivity, [NO3], modeled Nr deposition, gridded ambient summer and winter air temperatures, and modeled temperature trends 1948–2008 [Materials and Methods]. For each lake, we first standardized 20th century β-diversity estimates to the 19th century value, in order to account for between-lake differences in natural variability. We failed to uncover statistically significant relationships between standardized 20th century diatom β-diversity and any one of the physical or environmental variables considered. For example, 20th century diatom β-diversity is not correlated to lake-water pH (r = 0.01, p = 0.47). One reason there are no immediate correlations is the upward inflection of diatom β-diversity values at both extremes of the latitudinal gradient (i.e. arctic and alpine lakes, respectively), a pattern that is not mirrored by any single variable tested.


Diatom β-Diversity, Latitude, and Altitude

There are two geographic subsets within our data for which 20th century diatom β-diversity is highest: high latitude-low altitude lakes and high altitude-low latitude lakes (Table 1). It has previously been demonstrated that diatom community turnover since 1850 increases with latitude in the Arctic [2]. Our present results confirm this finding, and places it in the longer context of the last 450 years. Furthermore, we show that 20th century diatom turnover rates in high altitude lakes are of similar or greater magnitude to those recorded in the Arctic. When the data from arctic and montane-boreal lakes are collated (i.e. alpine lakes excluded), there is a significant linear relationship between 20th century β-diversity and latitude (r = 0.52; p<0.001; n = 37, Fig. 4A). Similarly, when alpine and montane-boreal lakes are combined (i.e. sites >60°N excluded), a significant linear relationship (r = 0.51; p<0.003; n = 32) emerges between 20th century β-diversity and altitude (Fig. 4B). In contrast, no relationships are evident between β-diversity and either altitude or latitude during the 19th century or the LIA, indicating that these ecosystems do not have naturally high rates of diatom assemblage turnover. Thus, the observed relationships between diatom β-diversity and both latitude and altitude (Fig. 4) confer greater ecological responses among arctic and alpine lakes to 20th century environmental changes than to natural perturbations, including the LIA cold spell.

Figure 4
Scatterplots of 20th century β-diversity against latitude and altitude.

These observed geographic patterns raise the hypothesis that 20th century diatom community dynamics were driven alternately by changes in two external factors: climate warming and enhanced Nr deposition. Summarizing the 20th century β-diversity trend over the complete latitudinal range of lakes in relation to the mean value (0.82 SD) enables a graphic delineation of the arctic and alpine regions characterized by the highest β-diversity values (Fig. 5A). When the climate change and Nr deposition gradients are added for comparison, Nr deposition rises for the alpine locations, while climate warming is greatest in the Arctic (Figs. 5B,C). Thus, the failure of 20th century diatom β-diversity to correlate statistically with any individual physical or limnological variable, including surrogates for Nr deposition and climate warming, can be explained by the sensitivity of β-diversity to more than one primary forcing factor across the 52-lake population. Given the geographical breadth of lakes considered, this is perhaps not surprising. However, the reality of multiple stressors is their potential to act independently in some regions, and synergistically in others. For example, the diatom stratigraphy from Emerald Lake (Fig. 2C), which produced the highest 20th century β-diversity of any site considered, reveals an early response consistent with warming (Discostella spp.) followed by later changes associated with Nr deposition (Asterionella formosa).

Figure 5
Latitudinal distribution of 20th century diatom β-diversity in relation to modeled Nr deposition and recent climate change.

Climate Change as a Driver of Diatom β-Diversity

A compelling range of data reveals the legacy of pronounced climate warming during the last century, which is firmly imprinted in both arctic [18], [19] and North American alpine [20], [21] environments. Regional compilations of arctic temperature reanalysis data [22] and paleoclimatic proxies [18] suggest that mean temperatures since ~1970 exceed any period within the prior 400 years (Fig. 6 A and B). The most recent arctic analysis [19] indicates that four of the five warmest decades of the past two millennia have occurred since 1950, and that recent summer temperatures are much as 1.4°C warmer relative to the long-term trend modulated by solar insolation. Furthermore, it has been established that the rate of arctic warming increases with latitude [23]. Although similar trends are evident in North American alpine regions, the rate and amplitude of warming is greater in the Arctic, primarily because positive feedbacks involving the cryosphere become accentuated with increasing latitude. Recent arctic warming corresponds closely with the regional increase of diatom β-diversity observed in high-latitude lakes over the 20th century (Fig. 5), especially given that the character of diatom assemblage shifts in these lakes is consistent with those anticipated from warmer summers and longer growing seasons: increased diatom biomass, species richness and species size, and in many cases increased representation by planktonic forms [2].

Figure 6
Temperature anomalies for arctic (A and B) and alpine (C and D) regions.

However, it is important to note that climate trends in the Arctic are spatially heterogeneous [24], [25]. Accordingly, regions that have experienced the least recent warming, including northern Québec, Labrador, and west Greenland, also contain the arctic lakes with the lowest 20th century diatom β-diversity values (e.g. Tasikutaaq Lake: 0.64 SD; Lake B: 0.76 SD, Table 1). This is concordant with the contention that climate warming is the primary driver of diatom β-diversity in the Arctic.

While climate change must be viewed as a potential driver of diatom shifts across each of the regions considered here, montane-boreal lakes tend to have considerably longer open-water seasons, so that thresholds for diatom taxonomic shifts are more gradually met. These ecosystems are potentially not as sensitive to small changes in climate or nutrient input [14], in part due to the edaphic stability of their catchments [26]. Relatively low mean 20th century diatom β-diversity in these lakes (Fig. 3A) is entirely consistent with these attributes. Nonetheless, certain montane-boreal lakes produce very high 20th century β-diversity, including Toquart Lake (Vancouver Island, 1.17 SD) and Namur Lake (northern Alberta, 1.16 SD). Such responses are, however, restricted to a minority of sites: 20th century diatom β-diversity exceeds the grand mean of 0.82 SD in only 3 of 17 montane-boreal lakes considered (Table 1).

The impacts of recent climate warming on alpine lakes has been more closely investigated, revealing a range of hydrological and limnological changes that ultimately confer high degrees of sensitivity to climate change [21], [27]. Reanalysis data from the Rocky Mountains show similar late 20th century temperature trends in both the Canadian and American sectors of the cordillera (Fig. 6C). The two warmest intervals of the 20th century (1950–1970 and 1990-present) are the only decades that exceed mean 1901–1960 temperatures inferred from the 400 year dendroclimate record (Fig. 6D). Thus, recent warming in the Rocky Mountains is of lesser magnitude than that witnessed in the Arctic. And yet mean 20th century diatom β-diversity from the alpine lakes is slightly higher than arctic counterparts (Fig. 3A), including the highest individual values and the only ones >1.4 SD (Table 1). In the few cases where analyses have been conducted in sufficient detail, the climatic legacy, while present, is less pronounced than that attributable to Nr deposition, with respect to both diatoms and stable isotopes [28]. The widespread recent proliferation of nitrophilous diatoms in American alpine lakes suggests that this assertion is regionally valid [4], [16]. The diatom succession from Emerald Lake aptly illustrates this type of complex response: Discostella spp. increase dramatically following the LIA, only to be swamped by Asterionella formosa after 1950 (Fig. 2C). It is on the basis of these observations that we consider in greater detail Nr deposition in relation to diatom community turnover in alpine lakes.

Nr Deposition as a Driver of Diatom β-Diversity

Atmospheric emissions of Nr have increased 9-fold over the last century, now exceeding 140 Tg N yr−1 globally. Humans fix more N2 than the sum of natural processes [29]. In aquatic ecosystems, greater Nr availability from atmospheric subsidies becomes ecologically relevant by alleviating N limitation; lakes that are naturally N-limited or co-limited by N and P are thus most susceptible. It has been estimated that deposition rates in excess of 1.5 kg N ha−1 yr−1 are sufficient to induce rapid changes in the diatom communities of lakes [30]. However, a pronounced latitudinal gradient exists with respect to gridded estimates of Nr deposition increases between 1860 and 1993, mirroring patterns of human settlement and the intensity of agricultural activities [3]. With respect to the lakes under consideration here, the greatest increases in Nr deposition rates occur in the region immediately adjacent to low-latitude, high-altitude lakes of the Rocky Mountains (Fig. 5B). Here, 20th century diatom β-diversity values are typically high, including the lakes in our survey that have changed the most Table 1). Given that climate warming is of lower amplitude than that witnessed in the Arctic, and considering the ecology of diatom taxa involved in these changes [4], [6], [16], [28], we surmise that Nr deposition has contributed directly to the inflection of 20th century diatom β-diversity values observed in low-latitude alpine lakes (Fig. 5A). This is not to say that climate is unimportant in driving recent ecological changes in the alpine, but rather that the combined influences of climate warming and Nr deposition may conspire to produce higher 20th century diatom β-diversities than those attributable to climate alone. The inclusion of additional alpine lake populations [1] is required to fully assess these nuances.

Synthesis of Results by Principal Components Analysis (PCA)

PCA allows an effective integration of both the climate and Nr deposition gradients for all 52 lakes, for which the additional variables of pH and lake:catchment area ratio are included to summarize the chemical and physical environments of each site. The independently-generated 20th century β-diversity values are included in the analysis as passive variables, and subsequently contoured onto the resulting plot (Fig. 7A). The first PCA axis accounts for 38.1% of total variance explained, and conveniently distinguishes sites with high Nr deposition (positive scores) from those more influenced by climate warming (negative scores). This is the only significant principal component in the analysis, when compared to a broken-stick model based on randomly-generated matrices of identical proportions and total explanatory power [31]. The limnological variables (pH and lake:catchment area ratio) are more closely associated with the second axis, which is not statistically significant. This suggests that regional climate change and Nr deposition rates, which define the primary axis, are more important than site-specific factors in the analysis.

Figure 7
PCA of site-specific temperature and Nr deposition trends, lake pH, and lake[ratio]catchment area ratio.

When 20th century diatom β-diversity is passively superimposed on the PCA and contoured, values increase concentrically away from the origin (Fig. 7A). This implies greater diatom change at sites towards the extremities of axis 1, confirming that climate warming and enhanced Nr deposition alternately represent the dominant influence on diatom assemblage turnover within subsets of our lake population (Fig. 7A). The ordering of lake scores on the first PCA axis provides a visualization of site distribution along a gradient opposing climate change with Nr deposition (Fig. 7B). Although there are some exceptions, the arctic and montane-boreal lakes generally produce negative scores and hence are more strongly associated with temperature change, whereas alpine lakes have negative scores that confer a greater influence to Nr deposition. Montane-boreal lakes cluster more closely with the arctic population, suggesting that they have changed primarily in responses to climate warming. This is corroborated by recent increases in centric planktonic diatom populations in many of these lakes, and the absence of similar trends among taxa responsive to Nr deposition, such as Asterionella formosa [32]. The PCA thus provides a concise a posteriori verification of the pattern identified earlier: changes in diatom assemblages over the 20th century relate to more than one environmental forcing, and these are expressed in a geographically coherent manner.

Finally, we note that lake scores on the primary PCA axis form a near-continuum of values (Fig. 7B), which mandates that some lakes are concommitantly influenced by climate change and Nr-deposition, even though their 20th century β-diversity values may be modest. While we have attempted to separate these forcings, ultimately such efforts may prove futile because Nr delivery to lakes and climate warming are not mutually exclusive of each other. For example, climate warming influences precipitation, and in turn Nr deposition rates from wet deposition [33]. Furthermore, in both the Artic [34] and the alpine [26], melting of glaciers and perennial snowpacks has the potential to relegate anthropogenic Nr archived in snow and ice back to surface waters, irrespective of Nr deposition rates from precipitation. In our view, it is the potential for such synergistic processes that represents the most ominous threat for aquatic ecosystems in the 21st century, given stratigraphic evidence that such combined changes are already underway in some regions (Fig. 2C).


During the 20th century, changes in diatom assemblages from both arctic and alpine lakes have accelerated relative to the previous 350 years. Elevated 20th century β-diversity is associated primarily with climate warming and Nr deposition, with strong regional variations in the degree of influence attributable to either forcing factor. Given future scenarios for both climate change [35], [36] and Nr deposition [3], the diatom β-diversity trends are unlikely to be reversed, and in our view will only become exacerbated as the 21st century progresses. Our study contributes further proof that distinctive biological fingerprints exist for the Anthropocene [37], while extending this notion to include the microbiota of remote lakes. As with organisms such as birds [38] and higher plants [39], future trajectories of lake diatom communities include states for which no prior analogs exist. Although the full range of ecological implications remain poorly understood, changes at the base of food webs necessarily entail consequences for higher trophic levels, while modifying the biogeochemical cycling of major elements including, but not limited to, carbon, nitrogen, phosphorus and silicon. Increases of primary production and organic matter sedimentation may also influence the recruitment of metals, both natural and anthropogenic, to sediments. The prognosis for truly unperturbed lake ecosystems, if indeed any still remain, is that they are highly susceptible to marked biological reorganizations. Lakes that have already entered new biological regimes will continue to change as humankind tightens its grip on both the global climate system and key biogeochemical cycles.

Materials and Methods

Site Selection, Core Chronology, and Diatom Analysis

Fifty-two diatom stratigraphies were compiled from lakes in western and northern North America and west Greenland (Fig. 1), spanning latitudes from 37.67°N to 79.33°N, and altitudes from 12 m asl to 3546 m asl (Table 1). The range of lake-water pH is 5.9–8.4, representing a large environmental gradient that is captured by diatom assemblages ranging, accordingly, from acidophilous to circumneutral to alkaliphilous. Sites can be categorized into arctic (>60°N; n = 20), alpine (above altitudinal tree-line; n = 15), and temperate montane-boreal lakes (<60°N and forested; n = 17) (Table 1). Despite being restricted to a single cordillera (the Rocky Mountains), the alpine lake population spans 15° of latitude, from southern Colorado to Alberta. The arctic lakes include sites from the continental Northwest Territories, much of the Canadian archipelago, northern Quebec, and west Greenland, together spanning 18° of latitude and 61° of longitude. To close the geographical gap (Fig. 1), and for meaningful comparisons with arctic and alpine sites, we also included a number of low elevation lakes in forested catchments of western Canada (Alberta, central British Columbia, and Vancouver Island). None of the lakes in our analysis is affected by direct point-source anthropogenic activities (e.g. shoreline development, effluent discharge, acidification). Most of the sites have been previously published upon in the context of local and regional studies (Table 1).

Continuous lake-sediment cores with an intact sediment-water interface were collected using a gravity-type corer [40]. Once recovered, cores were sampled continuously and immediately at either 0.25 cm or 0.5 cm resolution. Chronostratigraphy of the cores was established by 210Pb dating with alpha or gamma spectroscopy [41]. Sedimentation rates were extrapolated to an estimated age of 450 years only where constant stable values were reached in the unsupported 210Pb inventory of the core. Adequate chronology was the primary criterion for inclusion in our compilation. However, for many sites only the 20th century is dated with acceptable accuracy and precision, so that the number of sites for which earlier intervals are considered is somewhat less than 52 (Table 1). The three sites with previously unpublished chronologies have unsupported 210Pb inventories that can be reliably interpreted using the constant rate-of-supply model (Fig. S1).

Preparation of sediment samples for diatom enumeration followed standard protocols [42]. Diatom enumeration was carried out on selected sediment intervals, at resolutions ranging from 0.25 cm to 3.0 cm. Two taxonomic complexes were used include small (<15µm) and ecologically comparable Discostella and Cyclotella spp., which can not be resolved by light microscopy: these are the ‘Discostella stelligera complex’ (D. stelligera and D. pseudostelligera) and the ‘Cyclotella comensis complex’ (C. comensis, C. rossii, and C. tripartita).

Numerical Analyses

To estimate compositional turnover of sediment diatom assemblages, or β-diversity, relative abundances were analyzed using detrended canonical correspondence analysis (DCCA) constrained to time [12]. DCCA models assemblage composition as a unimodal response to one or more environmental variable. By constraining the DCCA to time, as inferred from 210Pb dating, we retain the biostratigraphic integrity of each sedimentary diatom sequence. The larger the β-diversity value obtained over the interval under consideration, the greater the assemblage turnover. This strategy allows for explicit comparisons between assemblage changes during the 20th century and discrete earlier intervals. The selection of time slices for diatom β-diversity calculations is based on climate change over recent centuries. Paleoclimatic reconstructions consistently reveal the LIA (~1550–1850 AD) as one of the coldest intervals of the Holocene, when glacial advances were widespread and low summer temperatures prevailed [43]. The onset of the LIA was gradual, and maximum glacier expansions occurred asynchronously between the late 17th and 19th centuries. The LIA terminated by the end of the 19th century, and temperatures have increased since, although not monotonically. The latter part of the 20th and the first decade of the 21st century have probably witnessed the highest temperatures of at least the last millennium [44]. Accordingly, we adopted the following time slices for diatom β-diversity calculations: 1550–1900, the 19th century, and the 20th century. Owing to dating uncertainty and variable sediment recovery from site to site, the number of diatom records decreases with age (Table 1). Where dating control was adequate, we included only those lakes with a minimum of three sediment intervals within any one time period. The DCCA protocol included square-root transformation of relative frequencies to stabilize between-taxon variance, no down-weighting of rare taxa, inclusion of all taxa with >1% relative abundance, detrending by segments, and non-linear rescaling. DCCA was performed using CANOCO 4.51 [45].

We subsequently used multiple linear regression to test for correlations between selected physical and limnological variables and 20th century diatom β-diversity. Variables include lake and catchment area, lake:catchment area ratio (a surrogate for water residence time), maximum depth, winter (December–February) and summer (June–August) ambient air temperatures, surface air temperature trends for each site from 1948–2008, measured lake-water pH, specific conductance, [NO3], and modeled atmospheric reactive nitrogen (Nr) deposition between 1860 and 1993 [46]. Three gridded air temperature parameters were extracted for each lake's coordinates: (1) winter (December–February) and (2) summer (June–August) air temperatures from high resolution (10′) interpolations from station means for the period of 1961–1990 [25], and (3) surface air temperature trends 1948–2008 extracted from the NCEP/NCAR Reanalysis Project [22], summarized using the linear trend (slope) of the data period. Nr deposition estimates were extracted from gridded (5° longitude by 3.75° latitude) global data-sets using the nearest grid centroid to the corresponding lake coordinates [3], [46].

Non-parametric Wilcoxon rank sum tests were used to assess the statistical significance of β-diversity between lake categories within each time period and between time periods for each lake category. Levels of statistical significance were evaluated using Bonferroni-adjusted values of p. Finally, in order to summarize and evaluate collective environmental influences on β-diversity lake-by-lake, ordination by principal components analysis (PCA) was undertaken on the following standardized variables: 1948–2008 temperature trend, ΔNr deposition 1860–1993, pH, and lake:catchment area ratio (with 3 sites excluded due to insufficient catchment data). The latter two variables were included as conservative measures of lake chemical and physical characteristics, respectively. These analyses, as well as the extraction of gridded climate data, were conducted using the open-source software environment R [47].

Supporting Information

Figure S1

210Pb chronology from three unpublished alpine lakes included in this study. Lakes are detailed in Table 1. All three age models met our criteria for inclusion in the study, yielding reliable dates through the unsupported 210Pb section of the core. Constant sedimentation rates near the limit of 210Pb dating are shown by the linear relationships of log 210Pb activity against cumulative dry mass and were used to extrapolate dates to approximately 1550.

(1.77 MB TIF)


This compilation could not have been realized without data contributions from D. Antoniades, B. Cumming, M. Douglas, B. Keatley, K. Laird, N. Michelutti, R. Pienitz and K. Rühland. We also thank E. Doxsey-Whitfield, B. Meneken, K. Nydick and R. Vinebrooke for field assistance, M. Anderson for aid in retrieval and management of climate data, J. Smol for insightful comments on early versions of the manuscript, as well as the journal reviewers and the Academic Editor for suggestions that have led to substantial improvements.


Competing Interests: The authors have declared that no competing interests exist.

Funding: Funding was provided by the Alberta Ingenuity Centre for Water Research, the Natural Sciences and Engineering Research Council of Canada, the U.S. National Science Foundation, and the U.S. Department of the Interior. W.O.H. received a travel grant from the Faculty of Graduate Studies and Research, University of Alberta, for the visit to the University of Bergen from which this paper emerged. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.


1. Battarbee RW, Thompson R, Catalan J, Gytnes J, Birks HJB. Climate variability and ecosystem dynamics of remote alpine and arctic lakes: the MOLAR project. J Paleolimnol. 2002;28:1–6.
2. Smol JP, Wolfe AP, Birks HJB, Douglas MSV, Jones VJ, et al. Climate-driven regime shifts in the biological communities of arctic lakes. Proc Nat Acad Sci U S A. 2005;102:4397–4402. [PubMed]
3. Galloway JN, Dentener FJ, Capone DG, Boyer EW, Howarth RW, et al. Nitrogen cycles: past, present, and future. Biogeochemistry. 2004;70:153–226.
4. Wolfe AP, Van Gorp AC, Baron JS. Recent ecological and biogeochemical changes in alpine lakes of Rocky Mountain National Park (Colorado, USA): a response to anthropogenic nitrogen deposition. Geobiology. 2003;1:153–168.
5. Bergstöm A, Jansson M. Atmospheric nitrogen deposition has caused nitrogen enrichments and eutrophication of lakes in the northern hemisphere. Glob Change Biol. 2006;12:635–643.
6. Nanus L, Williams MW, Campbell DH, Elliott EM, Kendall C. Evaluating regional patterns in nitrate sources to watersheds in National Parks of the Rocky Mountains using nitrate isotopes. Environ Sci Technol. 2008;42:6487–6498. [PubMed]
7. Battarbee RW. Diatom analysis and the acidification of lakes. Phil Trans R Soc London B. 1984;305:415–477.
8. Bennion H, Juggins S, Anderson NJ. Predicting epilimnetic phosphorus concentrations using an improved diatom-based transfer function and its application to lake eutrophication management. Environ Sci Technol. 1996;30:2004–2007.
9. Smol JP, Douglas MSV. From controversy to consensus: making the case for recent climate change in the Arctic using lake sediments. Front Ecol Environ. 2007;5:466–474.
10. Crutzen PJ. Geology of mankind. Nature. 2002;415:23. [PubMed]
11. Zalasiewicz J, Williams M, Smith A, Barry TL, Coe AL, et al. Are we now living in the Anthropocene? Geol Soc Am Today. 2008;18:4–8.
12. Birks HJB. Estimating the amount of compositional change in late-Quaternary pollen-stratigraphical data. Veg Hist Archaeobot. 2007;16:197–202.
13. Bradbury JP. A climatic-limnologic model of diatom succession for paleolimnological interpretation of varved sediments at Elk Lake, Minnesota. J Paleolimnol. 1988;1:115–131.
14. Rühland K, Paterson AM, Smol JP. Hemispheric-scale patterns of climate-related shifts in planktonic diatoms from North American and European lakes. Glob Change Biol. 2008;14:1–15.
15. Winder M, Reuter JE, Schladow SG. Lake warming favours small-sized planktonic diatom species. Proc Roy Soc London B. 2009;276:427–435. [PMC free article] [PubMed]
16. Saros JE, Michel TJ, Interlandi SJ, Wolfe AP. Resource requirements of Asterionella formosa and Fragilaria crotonensis in oligotrophic alpine lakes: implications for recent phytoplankton community reorganizations. Can J Fish Aquat Sci. 2005;62:1681–1689.
17. Moser KA, Smol JP, MacDonald GM, Larsen CPS. 19th century eutrophication of a remote boreal lake: a consequence of climate warming? J Paleolimnol. 2002;28:269–281.
18. Overpeck J, Hughen K, Hardy D, Bradley R, Case R, et al. Arctic environmental change of the last four centuries. Science. 1997;278:1251–1256.
19. Kaufman DS, Schneider DP, McKay NP, Amman CM, Bradley RS, et al. Recent warming reverses long-term arctic cooling. Science. 2009;325:1236–1239. [PubMed]
20. Luckman BH, Wilson RJS. Summer temperatures in the Canadian Rockies during the last millennium: a revised record. Clim Dyn. 2005;24:131–144.
21. Parker BR, Vinebrooke RD, Schindler DW. Recent climate extremes alter alpine lake ecosystems. Proc Nat Acad Sci U S A. 2008;105:12927–12931. [PubMed]
22. Kalnay E, Kanamitsu M, Kistler R, Collins W, Deaven D, et al. The NCEP/NCAR 40-year reanalysis project. Bull Am Met Soc. 1996;77:437–471.
23. Moritz RE, Bitz CM, Steig EJ. Dynamics of recent climate change in the Arctic. Science. 2002;297:1497–1502. [PubMed]
24. Box JE. Survey of Greenland instrumental temperature records: 1873–2001. Int J Climatol. 2002;22:1829–1847.
25. New M, Hulme M, Jones P. A high-resolution data set of surface climate over global land areas. Clim Res. 2002;21:1–25.
26. Heegaard E, Lotter AF, Birks HJB. Aquatic biota and the detection of climate change: are there consistent aquatic ecotones? J Paleolimnol. 2006;35:507–518.
27. Schindler DW, Donahue WF. An impending water crisis in Canada's western prairie provinces. Proc Nat Acad Sci U S A. 2006;103:7210–7216. [PubMed]
28. Enders SK, Pagani M, Pantoja S, Baron JS, Wolfe AP, et al. Compound-specific stable isotopes of organic compounds from lake sediments track recent environmental changes in an alpine ecosystem, Rocky Mountain National Park, Colorado (U.S.A.). Limnol Oceanogr. 2008;53:1468–1478.
29. Galloway JN, Cowling EB. Reactive nitrogen and the world: 200 years of change. Ambio. 2002;31:64–71. [PubMed]
30. Baron JS. Hindcasting nitrogen deposition to determine and ecological critical load. Ecol Appl. 2006;16:433–439. [PubMed]
31. Bennett KD. Determination of the number of zones in a biostratigraphical sequence. New Phytol. 1996;132:155–170.
32. Hazewinkel RRO, Wolfe AP, Pla S, Curtis C, Hadley K. Have atmospheric emissions from the Athabasca oil sands impacted lakes in northeastern Alberta, Canada? Can J Fish Aquat Sci. 2008;65:1554–1567.
33. Williams M, Baron JS, Caine N, Sommerfeld R, Sanford R. Nitrogen saturation in the Rocky Mountains. Environ Sci Technol. 1996;30:640–646.
34. Anderson RK, Miller GH, Briner JP, Lifton NA, DeVogel SB. A millennial perspective on arctic warming from 14C in quartz and plants emerging from beneath ice caps. Geophys Res Lett. 2007;35:L01502.
35. Lean JL, Rind DH. How will Earth's surface temperature change in future decades? Geophys Res Lett. 2009;36:L15708.
36. IPCC. Climate change 2007: the physical science basis. In: Solomon S, Qin D, Manning M, Chen Z, Marquis M, et al., editors. Working group 1 to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change. New York: Cambridge University Press; 2007.
37. Root TL, Price TL, Hall KR, Schneider SH, Rosenzweig C, et al. Fingerprints of global warming on wild animals and plants. Nature. 2003;421:57–60. [PubMed]
38. Stralberg D, Jongsomjit D, Howell CA, Snyder MA, Alexander JD, et al. Re-shuffling of species with climate disruption a no-analog future for California birds. PLoS One. 2009;4:e6825. [PMC free article] [PubMed]
39. Williams JW, Jackson ST. Novel climates, no-analog commmunities, and ecological surprises. Front Ecol Environ. 2007;5:475–482.
40. Glew JR, Smol JP, Last WM. Sediment core collection and extrusion. In: Last WM, Smol JP, editors. Tracking environmental change using lake sediments vol 1. Dordrecht: Kluwer; 2001. pp. 73–106.
41. Appleby PG. Chronostratigraphic techniques in recent sediments. In: Last WM, Smol JP, editors. Tracking environmental change using lake sediments vol 1. Dordrecht: Kluwer; 2001. pp. 171–203.
42. Battarbee RW, Jones VJ, Flower RJ, Cameron NG, Bennion H, et al. Diatoms. In: Smol JP, Birks HJB, Last WM, editors. Tracking environmental change using lake sediments vol 3. Dordrecht: Kluwer; 2001. pp. 155–202.
43. Bradley RS, Briffa KR, Cole J, Hughes MK, Osborn TJ. The climate of the last millennium. In: Alverson K, Bradley RS, Pedersen TF, editors. Paleoclimate, global change and the future. Berlin: Springer; 2003. pp. 105–141.
44. Mann ME, Zhang Z, Hughes MK, Bradley RS, Miller SK, Rutherford S, Ni F. Proxy-based reconstructions of hemispheric and global surface temperature variations over the past two millennia. Proc Nat Acad Sci U S A. 2008;105:13252–13257. [PubMed]
45. ter Braak CJF, Šmilauer P. CANOCO: software for canonical community ordination (version 4.5) Ithaca: Microcomputer Power; 2002.
46. Dentener FJ. Global maps of atmospheric nitrogen deposition, 1860, 1993, and 2050. Oak Ridge National Laboratory Distributed Active Archive Center; 2006. Available online:
47. R Development Core Team. R: A language and environment for statistical computing. Vienna, Austria.: R Foundation for Statistical Computing,; 2008. Available online:
48. Cleveland WS. Robust locally weighted fitting and smoothing scatterplots. J Am Stat Assoc. 1979;74:829–836.
49. Perren BB, Bradley RS, Francus P. Rapid lacustrine response to recent High Arctic warming: A diatom record from Sawtooth Lake, Ellesmere Island, Nunavut. Arct Antarct Alp Res. 2003;35:271–278.
50. Antoniades D, Douglas MSV, Smol JP. Quantitative estimates of recent environmental changes in the Canadian High Arctic inferred from diatoms in lake and pond sediments. J Paleolimnol. 2005;33:349–360.
51. Michelutti N, Douglas MSV, Wolfe AP, Smol JP. Heightened sensitivity of a poorly buffered high arctic lake to late-Holocene climatic change. Quat Res. 2006;65:421–430.
52. Keatley BE, Douglas MSV, Smol JP. Early-20th century environmental changes inferred using subfossil diatoms from a small pond on Melville Island, NWT, Canadian high Arctic. Hydrobiologia. 2006;553:15–26.
53. Wolfe AP, Cooke CA, Hobbs WO. Are current rates of atmospheric nitrogen deposition influencing lakes in the eastern Canadian Arctic? Arct Antarct Alp Res. 2006;38:465–476.
54. Michelutti N, Wolfe AP, Briner JP, Miller GH. Climatically controlled chemical and biological development in arctic lakes. J Geophys Res. 2007;112:G03002.
55. Perren BB, Douglas MSV, Anderson NJ. Diatoms reveal complex spatial and temporal patterns of recent limnological change in West Greenland. J Paleolimnol. 2009;42:233–247.
56. Saulnier-Talbot E. Impacts de l'évolution climatique postglaciaire sur les lacs de l'extrême Nord de l'Ungava, Québec. 2007. PhD thesis, Université Laval, Québec, Canada.
57. Hobbs WO, Wolfe AP. Recent paleolimnology of three lakes in the Fraser River Basin (BC, Canada): no response to the collapse of sockeye salmon stocks following the Hells Gate landslides. J Paleolimnol. 2007;40:295–308.
58. Laird K, Cumming BF. A regional paleolimnological assessment of the impact of clear-cutting on lakes from the central interior of British Columbia. Can J Fish Aquat Sci. 2001;58:492–505.
59. Laird K, Cumming BF, Nordin R. A regional paleolimnological assessment of the impact of clear-cutting on lakes from the west coast of Vancouver Island, British Columbia. Can J Fish Aquat Sci. 2001;58:479–491.
60. Saros JE, Interlandi SJ, Wolfe AP, Engstrom DR. Recent changes in the diatom community structure of lakes in the Beartooth Mountain Range, USA. Arct Antarct Alp Res. 2003;35:18–23.

Articles from PLoS ONE are provided here courtesy of Public Library of Science