Ecological Regime Shifts in Lake Kälksjön, Sweden, in Response to Abrupt Climate Change Around the 8.2 ka Cooling Event

A detailed diatom record from Lake Kälksjön, west-central Sweden, reveals two periods of abrupt ecological change correlative with the 8.2 ka cooling event. Using a combination of abrupt step changes and piece-wise linear regressions, the diatom data were analyzed for change points over time, and two sudden and large events that are described as regime shifts were detected. During the first event at c. 8040 cal. y BP, a doubling in diatom biomass took place over 5–10 years. This increase in primary productivity can be connected to an erosion event in the catchment that resulted in an abrupt increase in nutrient supply to the lake. The second event was characterized by a substantial shift within the planktonic diatom community from taxa indicative of colder conditions to those indicating warm over 5–10 years at c. 7850 cal. y BP. This event was superimposed on a successive change from periphytic to planktonic diatom dominance over a 250-year period and a gradual diversification of the periphytic community that spanned c. 150 years. Rapid climate warming following the 8.2 ka event likely caused these changes and both regime shifts are examples of externally driven abrupt ecological change. This study demonstrates that it is possible to detect, quantify and test for regime shifts in paleoecological data, and it highlights the need for high sampling resolution and precise chronological control. High-resolution paleoecological reconstructions of ecological regime shifts in response to climate change can provide useful analogues of future changes in ecosystem structure and functioning.


INTRODUCTION
Abrupt climate change and associated or independent abrupt ecological change has been observed in various ecosystems, including marine environments (for example , Beaugrand 2004;Collie and others 2004;Spencer and others 2011), Arctic lakes (Smol and others 2005;Ru ¨hland and others 2008), brackish lagoons (Jeppesen and others 2007), terrestrial ecosystems (Dearing 2008), coastal areas (Conley and others 2009;Duarte and others 2009) and shallow lakes (for example, Scheffer and van Nes 2007;Scheffer and Jeppesen 2007;Zimmer and others 2009).Abrupt changes, whether climatic or ecological, are often described as ''regime shifts'': sudden events that lead to rapid climate or ecosystem reconfiguration between alternative stable states (Lees and others 2006;Andersen and others 2009).When detected, past regime shifts offer rich opportunities to study the patterns and drivers of abrupt ecological change and to understand ecosystem responses to climate changes that may take place in the near future.Abrupt ecological changes have been a common phenomenon during the late Quaternary and have often been driven by abrupt climate changes, also frequent and widespread during this time (Williams and others 2011).
One manifestation of abrupt climate change is the ''8.2 ka cooling event'' which is assumed to be the highest magnitude climate change of the Holocene in the northern mid to high latitudes (Daley and others 2011).The event is believed to have been caused by the injection of glacial melt water from ice-dammed Lake Agassiz-Ojibway into the Labrador Sea (for example, Clarke and others 2004;LeGrande and others 2006;Wiersma and Renssen 2006;Daley and others 2011), slowing down the ocean thermohaline circulation (THC) and its associated heat transport (for example, Alley and others 1997; Barber and others 1999;Wiersma and Renssen 2006) and leading to a pronounced transient cooling in the North Atlantic region (Wiersma and Renssen 2006;Daley and others 2011).Analysis of Greenland ice-core records indicates that an abrupt cooling event began around 8175 ± 30 cal. y BP and lasted for c. 150 years, with maximum cooling of c. 3°lasting c. 70 years (Kobashi and others 2007).The event and its impact, however, remain to be understood in detail (Snowball and others 2010).Rohling and Pa ¨like (2005) argue that the short but intense 8.2 ka event was superimposed on a more general cooling in the Northern Hemisphere, which began at c. 8600 cal.y BP and lasted for several centuries, and it has been suggested that the climate in Fennoscandia (Norway, Sweden, Finland and parts of north-western Russia) was characterized by two or three colder episodes during a c. 400-year period centred on 8000 cal.y BP, rather than a single event (Ojala and others 2008;Zillen and Snowball 2009).A cooling episode around 8200 cal.y BP has been described in several proxy records in the Northern Hemisphere (for example, Alley and others 1997;Tinner and Lotter 2001;Veski and others 2004;Hammarlund and others 2005;Ojala and others 2008;Panizzo and others 2008;Snowball and others 2010) but due to the natural smoothing of palaeoclimate records, relatively low sampling resolution, chronological uncertainties, possible delays in ocean-atmospheric coupling and inconsistent responses of different proxy records to the assumed climate excursion, the records are difficult to compare and the total impact of the event has been difficult to study (Ojala and others 2008;Snowball and others 2010).
Annually laminated (varved) sediments enable precise and accurate dating which enables high temporal resolution studies to be carried out (Stanton and others 2010).Thus, varved sediments are a valuable tool in reconstructing rates of climatic and environmental change.Stanton and others (2010) undertook a detailed chronological study of a varved sediment sequence from Lake Ka ¨lksjo ¨n in west-central Sweden, a lake first studied by Zille ´n and others (2003).Snowball and others (2010) secured the Ka ¨lksjo ¨n varve chronology across the 8.2 ka cooling event to the dendrochronologically derived calendar-year timescale using radiocarbon wiggle-matching.Snowball and others (2010) also transferred the Greenland icecore records onto the same timescale by means of 10 Be analysis following Muscheler and others (2004).They reconstructed a period of increased erosion, interpreted as caused by increased winter precipitation, between 8066 ± 25 and 7920 ± 25 cal.y BP and drew the conclusion that aspects of climate related to precipitation changed abruptly in central Sweden at least five decades after the onset of the 8.2 ka event over central Greenland (Snowball and others 2010).
In this study, we use the precise and accurate chronology of Ka ¨lksjo ¨n to look for evidence for abrupt ecological change in the lake in response to the abrupt 8.2 ka cooling; using detailed diatom analysis, we examine how the lake ecosystem responded to the supposed cooling and subsequent warming.
High-resolution palaeoecological reconstructions, such as the one presented in this article, can help constrain the timing and trajectory of local population changes to abrupt climate change (Williams and others 2011) and provide analogues of future changes in aquatic ecosystem structure and functioning caused by climate change.The records can also serve as a testing ground for critical-threshold detection tools of interest to global-change ecologists and ecosystem managers seeking reliable and Ecological Regime Shifts and Climate Change advanced warning signs of impending regime shifts (Contamin and Ellison 2009;Williams and others 2011).

Site Description
Lake Ka ¨lksjo ¨n (60°09¢N, 13°03¢E) is located in the province of Va ¨rmland in west-central Sweden (Figure 1) and has an area of 0.4 km 2 .Maximum water depth is 14.2 m, and the lake is located below the highest shoreline at an altitude of 98 m a.s.l.(Zille ´n and others 2003).The lake has four inflows entering from the east and north, and it drains through a stream to the west.The lake catchment covers 4 km 2 and is dominated by managed boreal forest composed of spruce, pine and birch (Zille ´n and others 2003).Areas close to the lake margin are dominated by silty and clayey soils, and arable land occurs along the western shore of the lake.The remainder of the catchment is covered by sandy and silty till (Zille ´n and others 2003).
The study area was deglaciated around 10,000 cal.y.BP.The immediate surroundings of Lake Ka ¨lksjo ¨n, covered by fine-grained deposits, were initially part of Ancient Lake Va ¨nern (Bjo ¨rck 1995), and the lake became isolated at c. 9000 cal.y.BP due to isostatic uplift (Zille ´n and others 2003).The lake provides an ideal sedimentary environment for the study of rapid environmental change, because the seasonal cycle of sedimentation is preserved in the lake in the form of biogenic-clastic varves (Stanton and others 2010;Snowball and others 2010).By counting annual layers, such records can be precisely dated, providing highly resolved quantitative microfossil stratigraphies for studies of environmental change (Renberg and Segerstro ¨m 1981;Ojala 2001).

Core Collection and Chronology
Four parallel sediment profiles were collected from the lake ice during fieldwork in the late winters of 2002 and 2003 using a modified rod-operated Russian corer (Snowball and others 2010).Each profile consisted of a combination of five 1.5-m long overlapping cores.The sediment profiles cover a sediment thickness of 6.75 m each and include the entire post-isolation period (Snowball and others 2010).The sediments are stored at 4°C in a cold humidified room at the Department of Geology, Lund University.An initial varve chronology for Lake Ka ¨lksjo ¨n sediments was established by Stanton and others (2010).Incremental seasonal layer counting was conducted from the sedimentwater interface.Historically dated atmospheric lead pollution isochrones were detected and established ages of palaeomagnetic secular variations were combined with the annual layer counts to produce

1338
L. Randsalu-Wendrup and others a varve chronology that extended from the present day to the beginning of varve formation at 9193 ± 186 cal.y BP.The uncertainty of the varve chronology established by Stanton and others (2010) at 8000 cal.y BP was on the order of 150 years.Snowball and others (2010) reduced this uncertainty to 25 years through the radiocarbon wiggle-matching of a 400-year-long section of Stanton and others (2010) varve chronology that spanned the 8.2 ka cooling event.This method finds the best statistical fit of a series of 14 C measurements to the calendar-year timescale that forms the basis of the IntCal04 radiocarbon calibration curve (Reimer and others 2004).Snowball and others (2010) also matched the GRIP and NorthGRIP ice-cores to the same absolute timescale by comparing 10 Be data and tree-ring D 14 C.This procedure synchronized the Lake Ka ¨lksjo ¨n data and Greenland chronology ice-core data to the absolute timescale, without any dependency on changes in proxy-climate parameters.We emphasize that possible future changes to the tree-ringderived radiocarbon calibration curve will not affect the relative synchronization of the Ka ¨lksjo ¨n and Greenland ice-core data, although absolute ages might be altered.

Diatom Analysis
Continuous and quantitative samples for diatom analysis were taken using an aluminium tray of known volume.The sediments in the tray were cut into 78 samples under a macroscope using an Aniol electronic tree-ring measuring system in the Laboratory for Wood Anatomy and Dendrochronology at Lund University.They were cut to varves along the stratigraphic boundaries marking the division between winter and spring layers, and based on the wiggle matched varve chronology (Snowball and others 2010), the samples span the time period from c. 8500 to 7500 cal.y BP, with a resolution of 20 varves between c. 8500-8300 cal.y BP and c. 7700-7500 cal.y BP and 10 varves in between.The diatom subsamples were freezedried and approximately 0.01 g of dry sediment was digested using the water-bath technique described by Renberg (1990) with digestion by 30% H 2 O 2 (Battarbee and others 2001).Diatom concentrations (valves mg -1 ) were determined by adding divinylbenzene spheres to the cleaned samples (Battarbee and Kneen 1982).Cleaned subsamples were dried onto cover slips and mounted onto microscope slides using Zrax (RI $ 1.7+).At least 400 diatom valves were counted in each subsample under a light microscope at 10009 using phase-contrast optics.Diatom taxonomy followed Krammer and Lange-Bertalot (1986, 1988, 1991a, b), Round and Bukhtiyarova (1996), Lange-Bertalot andMetzeltin (1996), Ha ˚kansson (2002), Houk and Klee (2004) and Siver and others (2005) and the diatoms were grouped according to their habitats into planktonic and periphytic life forms.In total, 16 planktonic and 78 periphytic taxa were identified.All diatom data are expressed as relative percentage abundances and for further numerical analyses these were reduced to 25 taxa with more than 3% relative abundance.The diatom accumulation rate (DAR), PCA of the planktonic and periphytic communities and Hill's N2 diversity index (Hill 1973) were calculated separately using all identified taxa.

Pollen Analysis
Pollen samples at intervals of 5-9 cm were prepared using standard methods (Berglund and Ralska-Jasiewiczowa 1986) with the addition of six lycopodium spike tablets for the calculation of pollen influx (Stockmarr 1971).Pollen grains were well preserved and between 800 and 1100 pollen were counted per sample.Selected pollen data are shown as influx (grains cm -2 y -1 ).

Organic Carbon Content
The total carbon content of the same samples used for diatom analysis was analyzed using a Costech Instruments ECS 4010 elemental analyzer at the Department of Geology, Lund University.The subsamples were acid-washed in 10% HCl to remove trace amounts of carbonates, rinsed with deionized water and freeze-dried prior to being homogenized and crushed with a mortar and pestle.Total inorganic carbon content is negligible (Snowball and others 2010) and the total carbon content is thus assumed to equal the total organic carbon (TOC) content of the sediments.

Numerical Analysis
The Stratigraphic diatom assemblage and production zones (DAZ) were determined through regression tree analysis (De'ath and Fabricius 2000) using the mvpart package (De'ath 2011) for the R statistical software (R Development Core Team 2011) in combination with change-point analysis (described below).The smallest tree within one standard error of the tree with the lowest cross-validation relative error was chosen, producing four zones.A detrended correspondence analysis (DCA) performed Ecological Regime Shifts and Climate Change using the vegan package (Oksanen and others 2009) in R revealed that the gradient length in the diatom data was short (<2 SD units) and, therefore, principle component analysis (PCA) was performed on Hellinger transformed diatom data using the same package.Hellinger transformation preserves the Hellinger distances, a dissimilarity coefficient shown to have good properties for ecological data, between samples in the resulting PCA (Legendre and Gallagher 2001).One PCA axis was identified as explaining significant proportions of the variance in the stratigraphic data when compared to the broken stick (null) distribution (Jackson 1993).This axis explained 47% of the variance.Principle components analysis was also performed on the planktonic and periphytic communities separately.
In both the analyses one axis was significant, and these axes explained 61% of the variance in the planktonic community and 45% of the variance in the periphytic community, respectively.
Total DAR, the proportion of planktonic to total DAR, and the first principal components of both the planktonic and the periphytic communities were analyzed for change points over time using a combination of abrupt step changes and piece-wise linear regressions (for details see Carstensen and Weydmann 2012).The abrupt step-change model has a discontinuity where the mean of the time series changes from one time step to another, whereas the piece-wise linear regression is continuous with changing slopes for different periods.Logarithmic and logistic transformations were applied to total DAR and the planktonic proportion of DAR, respectively, to achieve time series with normal distributions.The principal components were approximately normally distributed and thus were not transformed.An autoregressive term (AR(1)-process) for the error process was included in the change-point analyses to account for potential correlation in the time series, acknowledging that samples were not equidistantly distributed with time (10, 15, 20, or 25 years in-between samples) but were multiples of 5 years.This implied that if two observations were, for example, 15 years apart (that is, 3 lags), the autoregressive parameter would be raised to the power of 3. Thus, the autoregressive parameter described temporal correlation at a 5-year scale.The significance of the change-point models was assessed using likelihood ratio statistics, and the abrupt and piece-wise linear models were compared using Bayes Information Criterion (BIC).Non-significant slopes in the piece-wise linear model were set to zero.The change-point detection analyses were carried out using PROC MODEL in SAS v.9.2.

RESULTS
The studied sediment section has been subdivided into four zones, which have different diatom community structures: the bottom zone is dominated by periphytic diatoms, the next is defined by the appearance of Puncticulata radiosa (Lemmerman) Ha ˚kansson 2002, the third encompasses a shift to a higher lake productivity and the top zone records a shift in dominant taxa from periphyton to plankton.All absolute ages that refer to Ka ¨lksjo ¨n hereafter have an uncertainty of ±25 years.
The pollen record shows a domination of Betula and Pinus during this time period with an increase in total pollen influx from around 8400 cal.y BP, determined by Betula and Pinus, with Pinus expanding slightly earlier than Betula.A subsequent decrease in total pollen influx was recorded at c. 8260 cal.y BP (Figure 3).

DAZ 2 (8206-8036 cal. y BP)
Periphytic diatoms dominate this zone, as well, and make up 60-75% of the total diatom assemblage (Figure 2).The periphyton is still dominated by Planothidium frequentissimum and small Fragilariaceae, but the planktonic assemblage is characterized by the sudden appearance and rise in abundance of Puncticulata radiosa, a centric planktonic diatom that makes up 60-70% of the plankton (Figure 2).The appearance of Puncticulata radiosa defines this DAZ.Other abundant planktonic species are Cyclotella DAR values remain relatively low throughout this zone, with DAR total between 50 and 175 9 10 4 valves cm -2 y -1 , reaching a minimum of 39 9 10 4 valves cm -2 y -1 at c. 8040 cal.y BP (Figure 3).At the same time, there is a drop in the TOC content from c. 9 to 6% (Figure 3).The samples belonging to this zone plot on the left side of the ordination diagram (Figure 4) and the change on PCA axis 2 is heavily influenced by the appearance of Puncticulata radiosa, which has its highest abundance in the lower left quadrant.
The total pollen influx values are low and relatively stable, fluctuating between 24000 and 28000 grains cm -2 y -1 throughout this period, with only minor variations in Betula and Pinus percentages (not shown).Relative abundance within periphytic assemblage (%) Figure 2. Diatom stratigraphy of the major (>3%) taxa in Lake Ka ¨lksjo ¨n as relative abundance (%), diatom assemblage and production zones (DAZ), and proportion of plankton and periphyton (%).Major (>5%) periphytic taxa as relative abundance (%) within the periphytic community and Hill's N2 for all species as well as plankton and periphyton, respectively.Grey line z-varve at 8066 cal.y BP.
Ecological Regime Shifts and Climate Change DAZ 3 (8036-7856 cal.y BP) This diatom assemblage is very similar to DAZ 2, with Puncticulata radiosa, C. distinguenda and T. flocculosa as the dominant planktonic species and Planothidium frequentissimum and small Fragilariaceae dominating the periphyton (Figure 2).The proportion of periphyton is approximately equal to plankton, with plankton reaching 50-55% towards

1342
L. Randsalu-Wendrup and others the top of the zone.The main characteristics of the zone are the increases around 7990 cal.y BP in both DAR plankton and DAR periphyton up to a DAR total of 360 9 10 4 valves cm -2 y -1 at c. 7970 cal.y BP (Figures 3, 5A, B).The periphytic accumulation rate reaches a peak of 240 9 10 4 valves cm -2 y -1 at this time and then drops to the same low levels as in DAZ 1 (Figure 3).The DAR plankton reaches high levels at the same time and then remains high, albeit fluctuating throughout the sequence (Figure 3).The samples plot in the lower left quadrant of the PCA diagram (Figure 4).The samples from the bottom of the zone are difficult to distinguish from the samples of DAZ 2, whereas the samples move to the right on PCA axis 1 further up in the zone.Because changes in DAR, the defining characteristic of DAZ 3, are not represented in the PCA, it is not surprising that the two groups are quite similar in the ordination diagram.
TOC values are low at the beginning and middle of this zone, increasing to c. 8% at the top.Total pollen influx values remain low throughout the zone, reaching their lowest values of 20300 grains cm -2 y -1 at c. 7890 cal.y BP (Figure 3), although Quercus influx values are already showing some recovery slightly prior to the other taxa.
DAZ 4 (7856-7500 cal.y BP) At c. 7850 cal.y BP, there is a substantial shift in the diatom assemblage of Lake Ka ¨lksjo ¨n, with planktonic diatoms becoming completely dominant.Planktonic species make up 70-90% of the total diatom assemblage in this upper part of the record and are mainly composed of Aulacoseira taxa. A. distans (Ehr.)Simonsen 1979 is the dominant, with an abundance of up to 50% of the planktonic species (Figure 2).Other important planktonic species are T. flocculosa and small centric D. stelligera (Cleve et Grunow) Houk and Klee 2004 and D. pseudostelligera (Hustedt) Houk and Klee 2004 (Figure 2).Puncticulata radiosa, which dominated two previous assemblage zones, disappears at the beginning of this zone and is not present in any large quantities, although it returns at the top of the sequence.In the periphytic assemblage, Planothidium frequentissimum decreases in abundance along with the small Fragilariaceae, and a slightly more diverse flora develops, with Achnanthidium minutissimum (Ku ¨tzing) Czarnecki 1994 and small Eunotia species increasing in abundance.Accumulation rates are relatively high, with DAR total fluctuating around 150-630 9 10 4 valves cm -2 y -1 .DAR periphyton remains stable at the same low values as in DAZ 1, whereas DAR plankton is varying between 75 and 500 9 10 4 valves cm -2 y -1 (Figure 3).The samples plot on the right side of the PCA diagram and are clearly set apart from the earlier samples (Figure 4).
The TOC content rises at the beginning of the zone and reaches a maximum towards the top (Figure 3).At c. 7860, the pollen influx values increase, with Betula and Pinus as the most Ecological Regime Shifts and Climate Change common taxa, but the deciduous trees (Alnus, Quercus, Ulmus and Corylus) also show increasing abundances throughout the zone (Figure 3).

Change Points Over Time
The step-change and the piece-wise linear models were equally good in describing the total DAR variations (Figure 5A), but the step-change model used one parameter less, resulting in a lower BIC, and was therefore considered more suitable to describe the trends in total DAR.The step-change model suggested that DAR total increased abruptly from 107 9 10 4 to 269 9 10 4 valves cm -2 y -1 between 8041 ± 25 and 8031 ± 25 cal.y BP.Conversely, the piece-wise linear model was much better at describing the trends in the DAR proportion of planktonic diatoms (Figure 5B).The piecewise linear model suggested a constant proportion of 33% planktonic diatoms that increased to 84% over a 250-year period between 8012 ± 33 and 7762 ± 24 cal.y BP, that is, almost immediately after the abrupt step change in DAR total .
Changes were also apparent at the community level.The first principal component of the plankton community reflects the dominance of Puncticulata radiosa at positive values and small Discostella spp.and A. distans at negative.The axis was best described with a step-change model, suggesting abrupt changes in the planktonic community around both 8201 ± 25 and 7851 ± 25 cal.y BP (Figure 5C).A piece-wise linear model was selected for the first principal component of the periphyton community showing a stable community, with Fragilariaceae and Planothidium frequentissimum as the main contributors to the negative axis values, before 8237 ± 102 cal.y BP.The community changed slightly over the following 300 years and then shifted more drastically over the next 150 years to positive axis values that are associated mainly with small Eunotia spp., attaining a new and stable community after 7787 ± 12 cal.y BP (Figure 5D).For all the four models, the 5-year correlation was estimated to between 0.64 and 0.74.Consequently, the 10-year correlation was the square of the AR(1) parameter estimates, the 15-year correlation was the cube and so forth.

DISCUSSION
Changes in diatom composition and abundance provide information on the prevailing conditions in Lake Ka ¨lksjo ¨n and its catchment in the centuries before and after the 8.2 ka event.The lake formed approximately 9000 years ago (Stanton and others 2010) and from c. 8500 to 8200 cal.y BP periphytic diatoms, in particular, small Fragilariaceae, dominated the diatom assemblage (Figure 6A).These species are common in lake sediments deposited in newly deglaciated terrain (Smol 1983) and are considered as adaptable and cosmopolitan (Lotter and others 2000).They have short life spans and are common in lakes with a prolonged ice cover and a short ice-free season (Lotter and Bigler 2000;Jones and Birks 2004).Rising pollen influx values indicate developing vegetation in the catchment from c. 8500 cal.y BP, consisting mainly of birch and pine forest.The presence of heavily silicified Aulacoseira species in the planktonic diatom assemblage suggests wind-induced turbulence and sustained water-column mixing, conditions that are required by the heavy Aulacoseiras to stay buoyant (Jones and others 2011).
At c. 8200 cal.y BP Puncticulata radiosa, a centric planktonic diatom, appeared (Figure 6B).This is a species indicative of cold summer conditions (Rimet and others 2009;Rioual and others 2006).It tolerates low light and low nutrient levels and is characterized by its fast growth (Rimet and others 2009).Change-point analysis of the planktonic community shows that there is an abrupt and significant change on the first principle component at c. 8200 cal.y BP (Figure 5C).PC1 can be interpreted as a climate proxy, with species such as Puncticulata radiosa, indicating cold summer conditions at positive values and species indicative of warmer stratified summer conditions, for example, D. stelligera and D. pseudostelligera (Ru ¨hland and others 2008; Enache and others 2011), at negative values.The overall DAR in the lake was still low, and the dominating periphytic community was made up of Fragilariaceae, which, in combination with the change in planktonic species composition, suggest weak stratification during the open-water period and deep water-column mixing.Low pollen influx values indicate less vegetation in the catchment and/or a stressed terrestrial environment.Comparable declines in tree pollen influx (Corylus, Alnus, Tilia, Ulmus) were recorded in Lake Ro ˜uge, Estonia (Veski and others 2004) between c. 8250 and 8150 cal.y BP but, in contrast to our results, a simultaneous increase in Betula influx was found.This was attributed to a higher tolerance of Betula to lowered temperatures during the flowering period in early spring.An increase in Betula while other tree taxa declined (excepting Pinus) was also recorded between 8200 and 8000 cal.y BP at Højby Sø, Denmark (Hede and others 2010), together with an indication of cooler summer temperatures.The changes in our diatom and pollen records can 1344 L. Randsalu-Wendrup and others be connected to the 8.2 ka event as defined by oxygen isotope records from the Greenland icecores (Kobashi and others 2007).Kobashi and others (2007) concluded that the event started at 8175 ± 30 years BP and that it took less than 20 years to reach the coldest period, with a cooling of 3.3 ± 1.1°C in central Greenland.In general, Northern Hemisphere conditions were cooler, drier and perhaps windier during the event (Alley and Agustsdottir 2005), and many palaeoclimate records from various parts of Europe suggest that it was predominantly a winter and spring event (for example, Alley and Agustsdottir 2005;Wiersma and Renssen 2006;Seppa ¨and others 2007).In Arctic and high-altitude lakes, cold and long winters with thick ice cover and delayed ice break-up in spring lead to a higher abundance of periphytic diatoms as compared to planktonic (for example, Lotter and Bigler 2000), and these observations provide a model for shifts in the diatom community at the 8.2 ka event.Many Fennoscandian lakes, for example, Lake Nautaja ¨rvi in southern central Finland, showed a shift from planktonic species dominance to more littoral species dominance at the 8.2 ka cooling (Ojala and others 2008) and in a small west Greenland lake an abrupt shift from a planktonic to a periphytic community, including the sudden disappearance of Cyclotella (Puncticulata) radiosa, took place at the same time (Andersen and others 2009).The ecological response to the 8.2 ka cooling event is somewhat different in Lake Ka ¨lksjo ¨n: here a prolongation of the early, postglacial, Fragilariaceae dominated conditions occurred, and planktonic diatom composition changed abruptly to taxa more tolerant of cold conditions, weak stratification and deep watercolumn mixing in the summer months.The differential response of lakes in different locations points to the role of local conditions, such as vegetation and morphometry, in mediating the nature of limnological response to a given environmental forcing.
Snowball and others (2010) reconstructed a period of enhanced erosion in the catchment of Lake Ka ¨lksjo ¨n from 8066 to 7920 cal.y BP.They argued that this anomalous period of erosion was caused by increased winter precipitation and higher spring meltwater discharge, increasing the input of minerogenic matter into the lake.This process caused the formation of the z-varve at 8066 cal.y BP (Figure 2): a varve with a distinctly thick mineral spring layer, which marks a transition from relatively organic-rich to relatively mineral-rich sediments (Snowball and others 2010).The timing of this coincides with a short period of increased global methane emissions and, by extension, Ecological Regime Shifts and Climate Change warming and wetting as reconstructed by Kobashi and others (2007) at c. 8070 BP.This short methane spike is also evident in the Greenland temperature record, which suggests that the short-term warmer and wetter event was a large-scale climatic excursion (Kobashi and others 2007).The Lake Ka ¨lksjo ¨n data suggest an increase in nutrient supply associated with elevated erosion and a resultant increase in lake productivity (Figure 6C).This is evident in the DARs, which show a doubling in 5-10 years time, centered at c. 8040 cal.y.BP (Figure 3 and 5A).The abrupt step-change character of the trend in total DARs, however, suggests an abrupt increase in nutrient input into the lake rather than a gradual increase, which supports the hypothesis of a major erosion event (initially recorded as the z-varve) in the catchment at c. 8060 cal.y BP.At this time, pollen influx had already been low for around 200 years, possibly indicating decreased vegetation cover and more exposed soils.The diatom community structure did not change significantly during this time period, suggesting that the overall climate was still cold.
PCA and regression tree analysis data (Figures 2  and 4) show that the diatom assemblage of DAZ 4 (7856-7500 cal.y BP) differs significantly from the rest of the sequence, which implies an important environmental shift in the lake (Figure 6D).The lake became plankton dominated, and within the planktonic community Puncticulata radiosa virtually disappeared, accompanied by an increase in lowalkalinity species such as A. distans.The planktonic community was dominated by Aulacoseira taxa, which suggests ample supplies of dissolved silica and nutrients, high turbulence keeping these heavy diatoms afloat (Ojala and others 2008;Panizzo and others 2008) and prolonged open-water seasons (Lund 1955).Also important were small Discostella species (D. stelligera and D. pseudostelligera), which are characteristic of stable lake stratification (Ru ¨hland and others 2008; Enache and others 2011), and T. flocculosa, which often increases with initial lake eutrophication (Anderson and others 1995).Taken together, the diatom assemblage suggests that the lake was stratified during summer; favoring the Discostellas, but that Aulacoseira taxa were abundant during the isothermal period in spring (for example, Rautio and others 2000).An alternative explanation, given that the samples are aggregates of 10 or 20 years, is that limnological conditions varied from year to year in this interval, with conditions more favorable for Discostella species in some years and those favoring Aulacoseira during other years.In general, rises in Cyclotella abundance indicate warming-induced lakewater changes and inferred higher summer temperatures in northern Sweden at this stage (Rosen and others 2000;Bigler and Hall 2002;Bigler and others 2006).The shift from Fragilariaceae to Cyclotella, and other planktonic taxa, in Lake Ka ¨lksjo ¨n mirrors data from Lake Vuolep Njakajaure in northern Sweden (Bigler and others 2006), where an appearance and dominance of planktonic Cyclotella taxa was recorded at c. 7800 cal.y BP.A similar pattern also has been observed in many Arctic lakes in response to recent climate warming (Smol and others 2005;Ru ¨hland and others 2008).The change within the planktonic community was abrupt at c. 7850 cal.y BP, whereas there was a more gradual change within the periphyton (Figure 5C, D).The periphytic community changed over a period of c. 150 years centered at c. 7870 cal.y BP, and a more diverse flora replaced the Fragilariaceae (Figure 2).The pollen data suggest that the vegetation in the catchment recovered, with rapidly increasing abundance of tree pollen (mainly birch and pine) after a distinct low at c. 7860 cal.y BP.The appearance of A. distans, a diatom indicative of increased humic contents of lake waters (Eloranta 1986;Turkia and others 1998), was probably favored by the development of organicrich soils in the catchment.Alnus (alder) also expanded at about the same time (c.7850) in Lake Vuolep Njakajaure (Bigler and others 2006).In this case, the authors suggested that the establishment of Alnus enhanced the high in-lake productivity because Alnus is able to fix atmospheric nitrogen, and as a result nitrogen concentrations in soils and near-shore waters probably increased (Bigler and others 2006).In newly deglaciated lakes in Alaska (Engstrom and others 2000;Fritz and others 2004), the establishment of Alnus has major impacts on lake chemistry and nutrient availability with pronounced consequences for diatom assemblages.
Lake Ka ¨lksjo ¨n displays two distinct periods of abrupt ecological change during the investigated time period: the change in accumulation rate at c. 8040 cal.y BP and the shift to planktonic dominance at c. 7850 cal.y BP.In this article, we use ''abrupt change'' in the broad phenomenological sense described by Williams and others (2011): ''abrupt change occurs when one or more variables in a system rapidly jump from one state to another, at rates that are rapid relative to background rates''.A ''regime shift'' has occurred if the new state is persistent (Folke and others 2004;Williams and others 2011).Regime shifts have been described as sudden and large events that lead to rapid ecosystem reconfiguration (for example, Lees and others 2006;Andersen and others 2009;Williams and

1346
L. Randsalu-Wendrup and others others 2011), hence we choose to describe both periods as regime shifts.
Both shifts are rapid, but the second one is more complex with different parameters involved.The shift at c. 8040 cal.y BP is an abrupt step change occurring over just 5-10 years (Figure 5A).The rapid increase in nutrient supply associated with elevated catchment erosion led to increased diatom productivity in the lake, but not a change in species composition.In lakes in the high Arctic, increased nutrient supply commonly does not lead to changes in diatom species assemblage (Douglas and Smol 2000;Smol and Douglas 2007) but has been inferred to affect the diatom productivity (Antoniades and others 2011), as seen here in Ka ¨lksjo ¨n.The change in species composition at c. 7850 cal.y BP is complex and can be seen as a successive change from periphytic to planktonic diatom dominance, accompanied by a continued increase in planktonic DAR, over a 250-year period (Figures 3,5B).Within the planktonic community, however, the shift from taxa indicative of colder conditions to those indicating warm was abrupt and took place over 5-10 years around 7850 cal.y BP (Figure 5C), whereas the change in periphytic taxa to a more diverse community took c. 150 years (Figure 5D) from c. 7950 to 7800.The change in species composition was most likely caused by climate warming leading to stable lake stratification, shallower mixing depth and increased nutrient recycling, which favored planktonic species over the periphyton.
Both regime shifts are examples of extrinsic changes where strong and rapid forcing (nutrient input and climate, respectively) is the primary causative agent of the abrupt ecological change (Williams and others 2011).They are also both examples of regime shifts resulting from smooth pressure-status relationships, where an abrupt change in external forcing leads to an equally abrupt ecosystem change: the regime shift in the driver is linearly mediated to the ecosystem state (Andersen and others 2009).However, the abrupt shift within the planktonic community can be described as a threshold-like or tipping-point regime shift: the change in driver (climate warming) does not seem to affect the planktonic community until a critical threshold is reached, at which point the system state changes rapidly (Andersen and others 2009).
These kinds of rapid ecosystem changes caused by extrinsic factors have been described before.They have, for example, been found at the Younger Dryas cooling and subsequent warming (for example, Ammann and others 2000; Birks and Ammann 2000;Yu 2007;Steffensen and others 2008;Gonzales and Grimm 2009;Williams and others 2011).The regime shifts in Lake Ka ¨lksjo ¨n are also consistent with the changes described by Smol and others (2005) and Ru ¨hland and others (2008), who have observed that recent climate warming caused an increase in planktonic diatoms, such as Cyclotella species, and a compensatory decrease in small periphytic species in several Arctic lakes.However, although these studies established that abrupt changes occurred during the time periods investigated, the shifts were only evaluated using qualitative assessments, and no inferential analyses were made.In this article, we demonstrate that using change-point models can help detect, quantify and test for regime shifts in palaeoecological data, and we propose that other regime shift studies can benefit from this approach.
A key element for rigorous quantitative methods, however, is high sampling resolution and precise chronological control.

Figure 1 .
Figure 1.A Map of Scandinavia showing the location of the study site Lake Ka ¨lksjo ¨n, located in Va ¨rmland, west central Sweden.B Map of the province of Va ¨rmland.C Close-up map of the lake showing inlets, outlets and coring site (redrawn from Stanton and others 2010).
Figure3.Diatom accumulation rates (DAR), weight % total organic carbon (TOC), total pollen influx and influx of selected tree pollen, and d 18 O records from the GRIP and NorthGRIP ice-cores(Rasmussen and others  2006).The GICC05 timescale was synchronized to the IntCal04 tree-ring timescale (seeSnowball  and others 2010).Diatom assemblage and production zones (DAZ) as in Figure2.

Figure 5 .
Figure 5.Time series of A the total diatom accumulation rate, B the planktonic proportion of the total diatom accumulation rate, C first principle component of the planktonic community, and D first principle component of the periphytic community.Log-likelihood and BIC values are given for the different models.The optimal change-point models (assessed by the lowest (BIC)) are shown with solid lines, and the alternative less suitable models with dashed lines.Grey lines in accordance with DAZ as in Figure 2.

Figure 6 .
Figure 6.Illustration of the hypothetical development of Lake Ka ¨lksjo ¨n across the 8.2 ka event.