Articles | Volume 18, issue 13
Research article
14 Jul 2021
Research article |  | 14 Jul 2021

Wildfire history of the boreal forest of south-western Yakutia (Siberia) over the last two millennia documented by a lake-sediment charcoal record

Ramesh Glückler, Ulrike Herzschuh, Stefan Kruse, Andrei Andreev, Stuart Andrew Vyse, Bettina Winkler, Boris K. Biskaborn, Luidmila Pestryakova, and Elisabeth Dietze

Wildfires, as a key disturbance in forest ecosystems, are shaping the world's boreal landscapes. Changes in fire regimes are closely linked to a wide array of environmental factors, such as vegetation composition, climate change, and human activity. Arctic and boreal regions and, in particular, Siberian boreal forests are experiencing rising air and ground temperatures with the subsequent degradation of permafrost soils leading to shifts in tree cover and species composition. Compared to the boreal zones of North America or Europe, little is known about how such environmental changes might influence long-term fire regimes in Russia. The larch-dominated eastern Siberian deciduous boreal forests differ markedly from the composition of other boreal forests, yet data about past fire regimes remain sparse. Here, we present a high-resolution macroscopic charcoal record from lacustrine sediments of Lake Khamra (south-west Yakutia, Siberia) spanning the last ca. 2200 years, including information about charcoal particle sizes and morphotypes. Our results reveal a phase of increased charcoal accumulation between 600 and 900 CE, indicative of relatively high amounts of burnt biomass and high fire frequencies. This is followed by an almost 900-year-long period of low charcoal accumulation without significant peaks likely corresponding to cooler climate conditions. After 1750 CE fire frequencies and the relative amount of biomass burnt start to increase again, coinciding with a warming climate and increased anthropogenic land development after Russian colonization. In the 20th century, total charcoal accumulation decreases again to very low levels despite higher fire frequency, potentially reflecting a change in fire management strategies and/or a shift of the fire regime towards more frequent but smaller fires. A similar pattern for different charcoal morphotypes and comparison to a pollen and non-pollen palynomorph (NPP) record from the same sediment core indicate that broad-scale changes in vegetation composition were probably not a major driver of recorded fire regime changes. Instead, the fire regime of the last two millennia at Lake Khamra seems to be controlled mainly by a combination of short-term climate variability and anthropogenic fire ignition and suppression.

1 Introduction

Wildfires in Siberia have become larger and more frequent in recent years (Walker et al., 2019), drawing attention from both scientists and the wider public. The occurrence of wildfires in high latitudes is closely linked to a dry climate and to rising temperatures, which have been increasing at more than twice the rate of elsewhere (Jansen et al., 2020; Lenton, 2012). With ongoing warming, the established fire regimes will likely be subject to future change. Fire has long been the most important ecological disturbance in boreal forests, shaping the appearance and composition of the world's largest terrestrial biome (Goldammer and Furyaev, 1996). It acts as a main driver behind the boreal forest's carbon pool, which comprises a third of all globally stored terrestrial carbon (Kuuluvainen and Gauthier, 2018), and strikes a fine balance between emitting carbon during combustion and sequestering it during the following regrowth (Alexander et al., 2012; Ito, 2005; Köster et al., 2018). Wildfires are thought to turn this balance towards acting as a net source of emissions (Kelly et al., 2016; Walker et al., 2019) while posing increasing risks to human livelihoods and infrastructure (Flannigan et al., 2009). Additionally, fire in boreal forests may trigger tipping points in tree mortality and tree density, as well as shifts in vegetation composition with continued global warming (Herzschuh, 2020; Lenton et al., 2008; Scheffer et al., 2012; Wang and Hausfather, 2020). For these reasons, a prediction of potential future changes in boreal fire regimes is imperative to inform and prepare adapted fire management strategies in a warming arctic and subarctic boreal environment. However, data about long-term changes in fire regimes, a prerequisite for model validation, are still very sparse for large parts of the Russian boreal forest despite it comprising more than half of the world's coniferous tree stocks (Nilsson and Shvidenko, 1998).

Macroscopic charcoal particles in sediment archives, derived from biomass burning, are commonly used as a proxy for fire activity (e.g. Clark, 1988; Conedera et al., 2009; Remy et al., 2018; Whitlock and Larsen, 2001). Charcoal records are an important tool for tracking past changes in fire regimes and searching for underlying causes and effects, from local to global scales (Marlon et al., 2008, 2013; Power et al., 2008). In recent years the evaluation of charcoal particle distributions was expanded to infer past fire frequencies (Higuera et al., 2007) and estimate source areas (Duffin et al., 2008; Leys et al., 2015), the type of vegetation burning with charcoal morphotypes (Enache and Cumming, 2007; Feurdean et al., 2020; Mustaphi and Pisaric, 2014), and fire intensity with charcoal reflectance (Hudspith et al., 2015). Yet, as the distribution of lake sediment studies in Russia is generally sparse (Subetto et al., 2017), the Global Paleofire Database (; Power et al., 2010) lists only a few continuously sampled macroscopic charcoal records across the Siberian boreal forest. Only recently have charcoal records of sufficient temporal resolution allowed for the assessment of fire return intervals (FRIs) in boreal European Russian and western Siberian evergreen forests, revealing close fire–vegetation relationships (Barhoumi et al., 2019; Feurdean et al., 2020). More studies have been conducted in North America (e.g. Frégeau et al., 2015; Hély et al., 2010; Hoecker et al., 2020; Waito et al., 2018) and boreal Europe (e.g. Aakala et al., 2018; Feurdean et al., 2017; Molinari et al., 2020; Wallenius, 2011). However, comparisons across boreal study sites are complicated by the differing predominant fire regimes in North America (high-intensity crown fires) and Eurasia (lower-intensity surface fires) (de Groot et al., 2013; Rogers et al., 2015).

The main fire regimes in the European and western Siberian evergreen boreal forest also differ markedly from those of its larch-dominated, deciduous counterpart in eastern Siberia. Many prevalent evergreen conifers (Pinus sibirica, Picea obovata, Abies sibirica) are commonly seen as fire avoiders and are more susceptible to crown fires (Dietze et al., 2020; Isaev et al., 2010; Rogers et al., 2015). The predominant eastern Siberian larches (Larix gmelinii, L. cajanderi, L. sibirica), on the other hand, can resist fires with an insulating bark protecting the cambium from heat, while their deciduous and self-pruning nature restricts fires from reaching the crown (Wirth, 2005). Moreover, larches are thought to benefit from occasional surface fires, leading to more saplings by clearing the lower vegetation layers of plant litter and mosses, although this might become a risk for young trees if fire frequencies increase above a certain threshold (Sofronov and Volokitina, 2010). Surface fires in larch forests might also play a role in the long-term preservation of permafrost (Dietze et al., 2020; Herzschuh, 2020; Herzschuh et al., 2016). These different fire strategies within the Siberian boreal forest reinforce the need for fire reconstructions towards the eastern part to evaluate changes in fire regimes depending on the prevalent tree species and to obtain a biome-specific overview of fire regimes throughout time.

The main goal of this study is to start filling a pronounced gap in the global distribution of macroscopic charcoal records by providing the first continuously sampled, high-resolution macroscopic charcoal record from eastern Siberia, using charcoal size classes and morphotypes. We specifically aim to answer the following research questions. (i) How did the fire regime in south-west Yakutia change throughout the last two millennia? (ii) How might a reconstructed fire history relate to common drivers behind changes in fire regimes (climate, vegetation, humans)?

2 Study site and methods

2.1 Location

Lake Khamra (59.97 N, 112.96 E) is located ca. 30 km north-west of the Lena River in south-west Yakutia (Republic of Sakha, Lensky District), at an elevation of 340 m a.s.l. (above sea level; Fig. 1a). Covering an area of 4.6 km2 the lake has a maximum water depth of 22.3 m. The catchment area of 107.4 km2 stretches around the gentle slopes surrounding the lake and extends in a south-western direction along the main inflow stream. It is embedded in a region of Cambrian bedrock made up of dolomite and limestone with transitions to silty Ordovician sandstone and locally restricted areas of clayey Silurian limestone (Chelnokova et al., 1988).

The surroundings of the lake were unglaciated during the Last Glacial Maximum (LGM; Ehlers and Gibbard, 2007), and discontinuous or sporadic permafrost can be present within the study area (Fedorov et al., 2018). Morphological features suggest that Lake Khamra is an intermontane basin lake that is not of thermokarst origin because of a lack of steep slopes (Katamura et al., 2009a), a maximum lake depth of 22.3 m compared to thermokarst lakes reaching mostly <10 m (Bouchard et al., 2016; West and Plug, 2008), and the absence of cryogenic features such as ice wedges indicating thermokarst processes (Katamura et al., 2009b; Séjourné et al., 2015). Accordingly, the regional soil volumetric ice content as broadly defined by Fedorov et al. (2018) is probably below the 30 % required for the development of thermokarst lakes (Grosse et al., 2013).

The region lies within the transition zone of evergreen to deciduous boreal forest. The lake catchment is covered by dense mixed-coniferous forest consisting predominantly of Larix gmelinii, together with Pinus sylvestris, P. sibirica, Picea obovata, Abies sibirica, Betula pendula, and Salix sp. (Kruse et al., 2019).

The continental climate of southern Yakutia is characterized by a short, mild growing season and extremely cold winters. The mean temperature of July at ca. 18 C is contrasted by the low mean temperature in January of about 40 C (Fedorov et al., 2018). Mean annual precipitation from 2000 to 2016 at the closest weather station in Vitim (ca. 60 km south-west of Lake Khamra) was ca. 480 mm, with the highest values in July to September (Russian Institute of Hydrometeorological Information, 2020). Wildfires most frequently occur in the central to southern regions of Yakutia, with varying stand-specific mean fire return intervals not exceeding 20–50 years at around 60 N, 120 E (Ivanova, 1996) and 15 years near Yakutsk (Takahashi, 2006). Longer estimates, increasing with higher latitude, were found by Ponomarev et al. (2016) for central Siberia of 80 years at 62 N, 200 years at 66 N, and 300 years at 71 N, with similar ranges reported by Kharuk et al. (2016).

The regional climate is mainly controlled by the large-scale Arctic Oscillation (AO) pattern emerging between cold high-latitude and mild mid-latitude air masses (Wu and Wang, 2002). It has been suggested that positive phases of AO, with low atmospheric pressure over the Arctic directing warm air from the south, potentially increase Siberian fire activity (Balzter et al., 2005; Kim et al., 2020).

Overall population density of the Lensky District is rather low at 0.5 inhabitants km−2 (Administrative Center Lensk, 2015). Lake Khamra is ca. 40 km north of one of the larger district settlements, Peleduy, with a population of ca. 5000 inhabitants. Traces of logging activity are visible in satellite imagery ca. 10 km from the lake, while in its direct vicinity only winter forest tracks are kept open. We discovered traces of recent wildfire disturbance at vegetation plots around the lake, i.e. trees with fire scars on their trunk or burnt bark. These traces likely correspond to recent fires within the catchment in 2006 and 2014 as captured by the remotely sensed forest loss data of Hansen et al. (2013) (Fig. 1b, c) and MODIS fire products (Giglio et al., 2018).

Figure 1(a) Position of Lake Khamra and existing Global Paleofire Database entries (last access: 6 October 2020) within the deciduous and evergreen boreal forests of Russia (land cover classification based on © ESA Climate Change Initiative Land cover project, provided via the Centre for Environmental Data Analysis (CEDA). (b) Infrastructure and burned areas 2001–2019 CE (Giglio et al., 2018; Hansen et al., 2013) in the vicinity of the lake (centre). (c) Lake Khamra with bathymetry and catchment area (black line) (service layer credits: Esri, DigitalGlobe, GeoEye, i-cubed, USDA, USGS, AEX, Getmapping, Aerogrid, IGN, IGP, swisstopo, and the GIS User Community).

2.2 Fieldwork and subsampling

Fieldwork at Lake Khamra was conducted in August 2018. A 242 cm long sediment core, EN18232-3, was retrieved using a hammer-modified UWITEC gravity corer at the deepest part of the lake (22.3 m), based on water depth measurements using a surveying rope and a handheld HONDEX PS-7 LCD digital sounder. From the same location, a parallel short sediment core (EN18232-2, 39 cm length) was retrieved and subsampled in the field at increments of 0.5 cm (upper 20 cm) and 1 cm (lower 19 cm) for lead-210 and caesium-137 (Pb/Cs) age dating. The sediment core, stored in a plastic tube, and all sediment samples, kept in Whirl-Pak® bags, were shipped to the Alfred Wegener Institute, Helmholtz Centre for Polar and Marine Research (AWI), in Potsdam and placed in storage at 4 C. In October 2018, sediment core EN18232-3 was opened and cut in half in a cooled room under sterile conditions. While one half was archived, the other half was subsampled for proxy analysis and radiocarbon (14C) dating. The sampling scheme included one 2 mL sample from the midpoint of every 2 cm increment (n=120), used for the combined palynological analysis and charcoal extraction, and two ca. 1 mL samples between each pair of larger samples (n=234), used specifically for charcoal extraction to ensure a continuous record of charcoal concentration. Bulk sediment samples for 14C age dating were extracted every 20 cm (n=12). At 85.5 cm core depth, a ca. 2 cm long piece of wood was removed from the sediment core for 14C dating. Due to a lack of any other larger organic structures, more macrofossil samples were picked while wet-sieving bulk sediment samples (n=15). Except for one case, a clear determination of their origin was not possible because of the small size of these samples.

2.3 Laboratory analyses

2.3.1 Lithology and age dating

Sediment core EN18232-3 was visually described before sampling. Water content and bulk density were determined in subsamples from 120 sampling increments every 2 cm.

To establish a chronology, all bulk and macrofossil samples were sent to AWI Bremerhaven for accelerator mass spectrometry 14C age dating at the MICADAS (Mini Carbon Dating System) laboratory. Subsamples of the parallel short core EN18232-2 (n=19) were sent to the University of Liverpool Environmental Radioactivity Laboratory for Pb/Cs age dating, analysing 210Pb, 226Ra, 137Cs, and 241Am by direct gamma assay with ORTEC HPGe GWL series well-type coaxial low-background intrinsic germanium detectors (Appleby et al., 1986). After careful evaluation of age dating results, an age-depth model was computed using Bacon v.2.4.3 (Blaauw and Christen, 2011; package “rbacon”; Blaauw et al., 2020), combining Pb/Cs and adjusted 14C ages of all bulk samples calibrated with the IntCal20 14C calibration curve (Reimer et al., 2020).

2.3.2 Macroscopic charcoal

We developed a sample preparation protocol that allows for the extraction of both macroscopic charcoal and the smaller pollen fraction including non-pollen palynomorphs (NPPs) from the same sediment sample. Lycopodium tablets (Department of Geology, Lund University), used as marker grains in the palynological analysis, were dissolved in 10 % HCl and added to the sediment samples. These were subsequently wet-sieved at 150 µm mesh width for separation of the macroscopic charcoal from smaller fractions (e.g. Conedera et al., 2009; Dietze et al., 2019; Hawthorne et al., 2018). The suspension with the <150µm fraction was collected in a bowl below the sieve. This “pollen subsample” was then iteratively added to a falcon tube, centrifuged, and decanted prior to further preparation for palynological analysis. The >150µm fraction in the sieve was rinsed together under a gentle stream of tap water before being transferred into another falcon tube. This macroscopic “charcoal subsample” was then left to soak overnight in ca. 15 mL of bleach (<5 % NaClO) to minimize the potential counting error from darker, non-charcoal organic particles (Halsall et al., 2018; Hawthorne et al., 2018).

The counting of 304 charcoal samples was done under a reflected-light stereomicroscope at 10–40× magnification. All particles that appeared opaque and mostly jet-black with charred structures were counted in every given sample (see Brunelle and Anderson, 2003; Hawthorne et al., 2018). In addition, counted particles were grouped into three size classes (150–300, 300–500, and >500µm measured along a particle's longest axis; Dietze et al., 2019) and based on similarities in shape (charcoal morphotypes; Enache and Cumming, 2007). For size reference, preparatory needles with known diameters of 300 and 500 µm were used that could be placed next to a charcoal particle. These needles also allowed for the careful and non-destructive evaluation of the flexibility of particles of unknown origin since charcoal fragments are described as fragile and non-bendable (Whitlock and Larsen, 2001).

The grouping of particles according to their shape was based on the morphotype classification scheme by Enache and Cumming (2007) and extended by three additional types to represent the variety of charcoal particles found at the study site. The original scheme differentiates between irregular (types M, P), angular (types S, B, C), and elongated (types D, F) shapes and further divides those depending on whether they show a visible structure or ramifications. The three additional types appear as highly irregular particles (type X), elongated, fibrous particles (type E), and slightly charred, partially transparent particles (type Z), the latter of which are not included in the total charcoal sum. For correlations and visualizations, the morphotypes were grouped into their respective main categories (irregular, angular, or elongated). Within the topmost ca. 50 cm of the sediment core, relative morphotype distributions were retrospectively derived from counts of 67 subsamples. At that time, 11 samples distributed within the top 40 cm of the sediment core had already been used for other purposes and thus lack information on morphotype classification (total charcoal concentrations and size classes are available for all samples).

Eight randomly selected samples were counted a second time to obtain an estimate of counting uncertainty and ambiguity in charcoal identification.

2.3.3 Palynological samples

Established protocols following Andreev et al. (2012) were applied to the pollen subsample (n=35). Of these samples, 11 were chosen specifically from intervals with high charcoal concentrations, whereas the others were spread across the sediment core. Samples were treated with boiling potassium hydroxide for 10 min, sieved, and left to soak in 18 mL hydrofluoric acid (40 % HF) overnight. After two additional treatments with hot HF (1.5 h each), acetolysis was performed using acetic acid and in a second step a mixture of acetic anhydrite and sulfuric acid. After being fine-sieved in an ultrasonic bath, the samples were suspended in glycerol.

Pollen and NPPs were counted together with the added Lycopodium spores on pollen slides under a transmitting light microscope with a minimum count sum of 300 particles. For subsequent statistical analyses, relative frequencies of individual pollen taxa were calculated from the sum of terrestrial pollen. Spore, algal, and non-pollen palynomorph percentages are based on the sum of pollen plus either spores, algae, or non-pollen palynomorphs, respectively (Andreev et al., 2020).

2.4 Statistical methods

Statistical analysis was carried out in R v.4.0.2 (R Core Team, 2020). To assess fire history, two different approaches were applied that decompose the charcoal record into a background component, representing long-term variations in charcoal accumulation and particle taphonomy, and a peak component, representing predominantly local charcoal accumulation during fire episodes (Higuera et al., 2007, 2009; Kelly et al., 2011; Whitlock and Anderson, 2003).

First, we ran a set of analyses referred to as “classic CHAR” (R script by Dietze et al., 2019), similar to the charcoal record decomposition in the well-established “CharAnalysis” (Higuera et al., 2009). We interpolated the charcoal record to equally spaced time intervals according to its median resolution and calculated the charcoal accumulation rate (CHAR, particles cm−2 yr−1; package “paleofire” v.1.2.4, function “pretreatment”; Blarquez et al., 2014). A background component was determined by computing a locally estimated scatterplot smoothing (LOESS) at a window width of 25 % of the total record length (package “locfit” v.1.5-9.4, function “locfit”; Loader et al., 2020). This window width was found to result in an efficient distribution of a signal-to-noise index (SNI) >3 based on Kelly et al. (2011), which indicates a high degree of separation between signal and noise (Barhoumi et al., 2019; Kelly et al., 2011). A peak component was created by subtracting the background component from the time series. By fitting two Gaussian distributions into the histogram of the peak component (package “mixtools” v.1.2.0, function “normalmixEM”; Benaglia et al., 2009), a global threshold was defined at the 99th percentile of the noise distribution (Gavin et al., 2006; Higuera et al., 2011). All peak component values exceeding the threshold were subsequently identified as signal (representing fire episodes) and marked when they overlapped with periods of SNI >3, indicating a clear distinction from surrounding noise. Usually in instances of multiple consecutive samples, only the highest peak above threshold is recognized, but we included all of them as fire episodes in this case. Recent fires within the lake's catchment that were just 8 years apart (Fig. 1b), as well as a predominantly atmospheric input of charcoal and a quick recovery of the filtering wetland vegetation around the large lake, led to the suspicion that the ability of a single fire to create high charcoal counts in multiple samples is limited at this site. However, an absolute minimum estimate of the number of fire episodes was obtained by considering only one fire episode from a given peak in CHAR and only those that are clearly separated from noise (i.e. overlapping with phases of SNI >3). A global threshold was chosen to fit the steady sedimentation rate and vegetation composition around the lake. We also tested a local threshold to see how robust our approach is when compared to alternatives. For this, Gaussian distributions were fitted to samples of the peak component within a moving window of the same width as used for background component determination (package “zoo”; function “rollapply”; Zeileis and Grothendieck, 2005). Subsequent threshold values were obtained the same way as before but additionally smoothed with a default LOESS to minimize the impact of outliers.

Fire return intervals (FRIs) were determined as the temporal difference between subsequent fire episodes. An illustration of fire frequency, as the rate of identified fire episodes, was derived by counting all fire episodes within a moving window spanning 200 years before applying a LOESS of the same window width to provide a clearer visualization.

To account for accumulated uncertainties from both the chronology and the counting procedure, we used a Monte-Carlo-based (MC) approach (for detailed description and R script, see Dietze et al., 2019), referred to as “robust CHAR”. In short, it describes both the age and proxy values of each sample as Gaussian distributions, creating a pool from which random values are sampled. As inputs we used the 2σ range of Pb/Cs ages and 1σ range of calibrated and adjusted 14C ages to scale the general magnitudes of uncertainty between the two dating methods to comparable dimensions. For proxy uncertainty, the average deviation between the repeatedly counted samples was used (ca. 20 %). Following Dietze et al. (2019), we aggregated three consecutive samples to scale the high resolution of the charcoal record to the relatively high age and counting uncertainties. A total of 5000 MC runs were performed, and output data resolution was set at 3 times the record's median temporal resolution (18 years). Robust CHAR was then divided into background and peak components, similar to classic CHAR but by computing 1000 randomly sampled LOESS fits at varying window widths ranging from 5 % to 25 % of the record length, thereby not relying on an individual user-input value (Blarquez et al., 2013).

The statistical approach outlined above (classic and robust CHAR including the determination of SNI, FRIs, and fire frequency) was also applied to the charcoal size classes and morphotype categories to assess their individual contribution to the sum of all particles, whether different charcoal groups represent different types of fires, the potential relationships between charcoal particle size and source area, and/or varying fuel types over time (see Supplement).

A principal component analysis (PCA; package “stats”; R Core Team, 2020) was used to assess relationships among the various centred-log-ratio-transformed (clr) (package “compositions”; van den Boogaart et al., 2020) relative distributions of charcoal particle classes. Assuming that the various charcoal morphotypes were formed by different types of vegetation burning, we would expect an increase in a specific morphotype to coincide with changes in the distribution of some plant types in the vegetation composition, also represented by the pollen spectra. To explore this hypothesis, we applied a correlation test using Kendall's τ (package “psych”; Revelle, 2020) to clr-transformed relative distributions of pollen groups that were expected to be impacted by wildfires (pollen sums of arboreal, non-arboreal, deciduous, and evergreen taxa) and charcoal classes following Dietze et al. (2020).

3 Results

3.1 Lithological sediment properties and chronology

Sediment core EN18232-3 shows no lamination or visible changes in its brown colour or the texture of the sediment matrix. Its homogenous appearance is underlined by both uniform mean dry bulk density (184 ± 4 mg cm−3, mean ± 1σ) and water content (83 ± 4 %).

Bulk sediment 14C ages indicate a rather linear chronology, with only the deepest two samples returning similar ages. In contrast to this, the 14C ages of the macrofossils do not form a clear chronological pattern and are not in good agreement with the bulk sediment samples. Only 9 of the 15 macrofossil samples were large enough to be used for 14C dating and most were only just above the minimum amount of carbon (Table 1a). As Pb/Cs dating from the parallel short core EN18232-2 reveals no disturbances in its uniform sedimentation rate, the ages are expected to be applicable to the main core EN18232-3 from the same location within the lake and thus can be compared to the respective 14C ages. Noticeably, the topmost 14C bulk sample dates to 1415 ± 27 14C years BP (before present, i.e. before 1950 CE), whereas the Pb/Cs method from the parallel core confirms the expected recent surface age (Table 1b). This 14C age offset can have a variety of causes, such as sediment mixing processes (Biskaborn et al., 2012), the presence of old organic carbon (Colman et al., 1996; Vyse et al., 2020), or dissolved carbonate rock (hard-water effect; Keaveney and Reimer, 2012; Philippsen, 2013).

In the case of Lake Khamra, located in a zone of discontinuous permafrost and bedrock containing early Palaeozoic carbonates, the observed 14C age offset is a likely consequence of input of both old organic and inorganic carbon through the south-western inflow stream. The magnitude of this offset is documented by the difference between the 14C bulk surface sample and the corresponding Pb/Cs age. There are indicators to support the assumption that accumulation of old carbon not only happened recently but rather constitutes an ongoing process in this lake system. First, treating the surface 14C age as an outlier without also adjusting the other 14C dates would necessarily lead to a strong shift in sedimentation rates, which is reflected neither by the homogenous appearance and density nor by the uniform sedimentation rate as implied by the Pb/Cs method. This is underlined by reports of stable lake conditions and thermokarst processes in central Yakutia during the late Holocene (e.g. Pestryakova et al., 2012; Ulrich et al., 2019). Second, macrofossil 14C ages older than the surrounding sediment matrix provide direct evidence for the potential influence of old carbon on bulk sediment samples at various depths (e.g. macrofossil age of 9902 ± 97 14C years BP in sediment that dates back to only ca. 100 years BP according to the parallel core's Pb/Cs age). For these reasons, the documented age offset (mean ± 1σ) in the topmost 14C sample was used to adjust all other 14C bulk samples (Colman et al., 1996; Vyse et al., 2020). Although macrofossil samples are usually thought to be superior in precision to bulk sediment for age-dating purposes (Hajdas et al., 1995; Wohlfarth et al., 1998) and have also been used to quantify a 14C age offset throughout a sediment core (Gaglioti et al., 2014), within the present environment we cannot exclude a potential permafrost origin, as well as measurement uncertainties due to their small size. For this reason, they were not used for constructing the chronology.

The age-depth model created according to these observations (Fig. 2) shows a smooth transition from Pb/Cs dates towards adjusted bulk 14C ages. Its uniform sedimentation rate mirrors the sediment core's homogenous composition and supports the underlying assumption of a rather constant magnitude of old carbon influence on bulk 14C ages. Based on this chronology, the sediment core spans ca. 2350 years across its 242 cm length. The continuously sampled charcoal record spans ca. 2160 years and thus reaches back from 2018 CE until ca. 140 BCE. Mean (± 1σ) sampling resolution of the whole record is 7.1 ± 4.1 years (max: 27; min: 3), including the pollen samples, which cover on average 8.9 ± 3.8 years (max: 25; min: 4).

Table 1(a) 14C age results of bulk and macrofossil samples from core EN18232-3. All samples used for the chronology are marked with a * in the first column. Lab-IDs 3732.1.1, 3740.1.1, 3743.1.1, 3743.2.1, 3745.1.1, and 3753.1.1 (not shown) belong to the six macrofossil samples that are too small to be successfully dated. (b) Pb/Cs age results from parallel short core EN18232-2 (analysed by Peter Appleby, personal communication, 2019).

Download Print Version | Download XLSX

Figure 2(a) Bayesian age-depth model applying Bacon (Blaauw and Christen, 2011) to several types of determined sediment ages (red =Pb/Cs; blue =14C bulk sediment; green =14C macrofossils; grey lines = 2σ range; red line = median). (b) Model iteration log. (c, d) Prior (green line) and posterior (grey area) distributions for accumulation rate and memory, respectively. (e) Dry bulk density of sediment core EN18232-3 (grey line = mean).


3.2 Reconstructed fire regime

Samples of the charcoal record contain on average 8.1 ± 5.1 charcoal particles per cubic centimetre of sediment (min: 0; max: 38.8). Interpreting the classic peak component as temporally restricted increases in fire activity, 50 such fire episodes within the continuously sampled core segment spanning the last ca. 2200 years were identified. This results in a record-wide mean FRI of 43 years (min: 6; max: 594). A maximum estimate of this FRI, based on the reduced number of fire episodes by limiting one episode per peak with SNI >3, is 95 years (min: 12; max: 876). Results from the alternative local threshold variant slightly differ from the global threshold in expected ways (fewer fire episodes identified during periods of high CHAR and more during periods of low CHAR) but do not produce majorly contrasting results (see Appendix A).

Classic and robust CHAR analysis distinguished four phases, representing different states of the fire regime (see Fig. 3). Phase 1 (ca. 140 BCE to 600 CE) is characterized by relatively high CHAR of 0.83 ± 0.43 particles cm−2 yr−1 and numerous (n=23) peaks (i.e. fire episodes, Fig. 3c), with an SNI that is slightly above 3 for the most part but slowly decreasing towards the following phase 2 (Fig. 3b). The mean FRI is 31 years (min: 6; max: 144). Robust CHAR shows a steadily increasing background component (Fig. 3d), whereas its peak component (Fig. 3e) remains at low levels. More fire episodes (n=16) occur in the following, shorter phase 2 (600–900 CE), with a lower mean FRI of ca. 14 years. It incorporates some of the highest peaks of the whole record using either the classic or robust approach, resulting in a high CHAR of 1.3 ± 0.54 particles cm−2 yr−1 (including the maximum of 2.8 particles cm−2 yr−1), while SNI falls below 3. In the transition to phase 3 at around 900 CE a pronounced decrease in charcoal input leads up to a long period of few to no fire episodes (n=6) and longer mean FRI of >60 years (min: 6; max: 594). Hence, CHAR of phase 3 (900–1750 CE) is low with 0.53 ± 0.29 particles cm−2 yr−1, while robust CHAR background remains below average. The SNI decreases to a record-wide minimum due to the lack of any CHAR peaks above the global threshold in the second half of phase 3. Finally, phase 4 (1750–2018 CE) has higher CHAR (0.61 ± 0.47 particles cm−2 yr−1) and more frequent occurrences of fire episodes (n=5), with mean FRI sharply decreasing to 40 years (min: 6; max: 78). An outstanding peak around 1880 CE leads to a maximum SNI >6 during this phase. Although phase 4 sees increasing CHAR and fire frequency after the low CHAR of phase 3, charcoal input decreases to a minimum within the last century (mean CHAR of the last 100 years before core extraction = 0.56 ± 0.29 particles cm−2 yr−1). Similarly, the robust CHAR sum and its components show increases within phase 4 (Fig. 3d, e), with two maxima in the robust peak component around the early 1800s and 1950 CE and a following decrease in CHAR (Fig. 3d–e). In general, the older half of the record (ca. 140 BCE to 1000 CE) has higher mean CHAR and a higher variability (0.97 ± 0.5 particles cm−2 yr−1) compared to the younger half (1000–2018 CE; CHAR = 0.51 ± 0.32 particles cm−2 yr−1). Even with added uncertainties from counting and the chronology, maxima of robust CHAR mostly replicate periods of increased classic CHAR.

Over the entire record, 43.7 % of charcoal particles belong to the smallest size class (150–300 µm), while 28.8 % and 27.5 % are part of the medium (300–500 µm) and large size classes (>500µm), respectively. When assessed individually, more fire episodes are identified for smaller particles than for larger particles (Table 2). This is expected, as smaller particles tend to have a larger source area, potentially incorporating more fire events into the signal (Conedera et al., 2009).

The most prevalent morphotypes present in the sediment are types F (elongated, 31.7 %), M (irregular, 28.4 %), S (angular and black, 20.6 %), and B (angular and brownish-black, 7.2 %), with all others (X, C, D, E, P) ranking at or below 3 % each. Noticeably, two of the highest peaks of the record at ca. 650 and 1880 CE are composed primarily of large (>500µm), elongated (type F) particles. The total relative amount of type F particles seems to correlate with the largest particle size class, whereas type M is more closely associated with smaller particles (see Fig. 4 and Appendix B). Furthermore, the PCA indicates that there are rather weak grouping patterns of morphotype or size-class distributions in samples of increasing core depth and age, potentially reflecting a stable vegetation composition around the lake. The three charcoal morphotype groups show a similar temporal pattern for their background and peak component distributions (Fig. 6c, d), mostly mirroring the decreased variability in the second half of the record as described for the sum of all particles. However, when assessed individually, large particles have a generally lower variability than the other size classes, whereas the variability in irregular morphotypes is higher than that in elongated or angular particles (see Supplement).

Figure 3Overview of the charcoal record using classic and robust analysis approaches. Vertical dashed lines mark the different phases of the fire regime. (a) Classic CHAR peak component (dark-grey bars = signal, light-grey bars = noise, dashed horizontal line = threshold). (b) Signal-to-noise index (SNI) of the classic CHAR peak component based on Kelly et al. (2011) (red horizontal line = SNI cut-off value of 3). (c) Classic CHAR sum (black line = interpolated CHAR; blue line = LOESS representing the CHAR background component; red vertical lines = fire episodes with SNI >3; grey vertical lines = fire episodes with SNI <3). (d) Robust CHAR background component. (e) Robust CHAR peak component (red areas = above-average values). (f) Robust CHAR sum. For (d)(f), black line = median and grey area = interquartile range.


Table 2Reconstructed fire episodes, fire return interval (FRI) in years, and distributions of signal-to-noise index (SNI) >3 for charcoal morphotype and size classes, as well as the sum of charcoal particles over the last two millennia.

Download Print Version | Download XLSX

Figure 4Principal component analysis (PCA) of charcoal particles (size classes and individual morphotypes with a share of >5 % of the charcoal record), with coloured dots representing potential grouping patterns of charcoal assemblages with increasing age.


3.3 Vegetation history

The pollen and NPP record, covering the whole sediment core and reaching back ca. 2350 years, generally indicates a relatively stable vegetation composition (Fig. 5). The dominant arboreal pollen (AP) types comprise most of the pollen spectra (average ratio of AP : NAP =8.3:1) and include the trees and shrubs recorded around the lake. In descending order regarding their share of the pollen sum these are Pinus, Betula, Picea, Abies, Alnus, and Larix, with smaller amounts of Salix, Juniperus, and Populus. Non-arboreal pollen (NAP) types are predominantly represented by Cyperaceae, followed by less abundant Poaceae, Ericales, and Artemisia. Despite similar general palynomorph distributions, pollen assemblages can be divided into two subzones, with the upper subzone (Ib) seeing intervals of increased variability in the shares of some tree pollen and Cyperaceae (around 10 and 120 cm depth, corresponding to ca. 1950 and 700 CE, respectively). The lower subzone (Ia) demonstrates generally lower shares of Abies and Cyperaceae pollen. Charcoal particles and pollen types have mostly weak correlations without statistical significance at p>0.05, although they hint at weak associations between irregular morphotypes with AP being contrasted by those of angular morphotypes and NAP (see Appendix B).

Figure 5Pollen and non-pollen palynomorph percentage diagram from sediment core EN18232-3 at Lake Khamra (dots represent pollen taxa <1 %; A= algal remains; Z= invertebrate remains; horizontal dashed line = separation of subzones Ia and Ib).


4 Discussion

4.1 Fire regime history of the last two millennia at Lake Khamra

We use the term “fire episode” instead of “fire event” when referring to identified peaks in the record. This is to highlight that multiple fires could have contributed to any peak in the charcoal record. Consequentially, the FRIs of this study should be regarded rather as “fire episode return intervals”, marking the time span between periods of increased fire occurrence within the charcoal source area. This is because the relatively large water-surface area of Lake Khamra likely captures charcoal from a larger source area than smaller lakes. A larger source area of charcoal is directly related to an increased number of fires that were able to contribute charcoal to the present record. The gentle and densely vegetated slopes framing the catchment limit secondary charcoal input (Whitlock and Larsen, 2001), thus emphasizing a direct fire signal with predominantly primary input through the air. However, a higher number of captured fires from a larger source area also means that the comparability of FRIs reconstructed in this study to those obtained from smaller lakes or tree-ring chronologies may be limited since those usually convey direct fire impact and are more locally constrained (Remy et al., 2018).

Although it has been shown that larger charcoal particles originate generally within a few hundred metres of a lake archive (Clark et al., 1998; Higuera et al., 2007; Ohlson and Tryterud, 2000), they have also been observed to travel further depending on vegetation and fire conditions (Peters and Higuera, 2007; Pisaric, 2002; Tinner et al., 2006; Woodward and Haines, 2020). As wildfires in the Siberian boreal forest are predominantly considered low-intensity surface fires (de Groot et al., 2013; Rogers et al., 2015), the potential of the resulting convection to transport large charcoal particles is probably limited compared to high-intensity crown fires. We therefore assume a charcoal source area between a few hundred metres directly around the lake for low-intensity fires (Conedera et al., 2009) and increasing distance of up to several kilometres for more intense fires producing stronger convection, resulting in a total source area estimate of up to ca. 100 km2. Even though some extreme fires may well surpass this estimate and, occasionally, small fires within might fail to contribute sufficient amounts of charcoal, identified fire episodes in the charcoal record should still be biased towards fires closer to the lake, especially when they consist of predominantly large charcoal particles (Conedera et al., 2009; Whitlock and Larsen, 2001). The uncertainty regarding any source area estimate highlights the need for further spatial calibration studies. Also, it remains unknown whether the charcoal record is dominated more by close-proximity low-intensity fires or by high-intensity fires from a greater distance. An estimation of fire intensity from charcoal particles (e.g. measuring charcoal reflectance, Hudspith et al., 2015, or the charcoal's oxygen to carbon ratio, Reza et al., 2020) could potentially clarify the respective contributions and thus help with better constraining the source area.

The macroscopic charcoal record at Lake Khamra (Fig. 3) reveals gradually increasing but relatively stable fire activity from 140 BCE to 600 CE (phase 1). A period of high fire activity takes place between 600 and 900 CE, expressed as higher CHAR and shorter FRIs (phase 2). It then transitions into an almost 600-year period without any identified fire episodes and low CHAR (phase 3). From around 1750 CE the modern fire regime begins to take shape (phase 4), with regularly identified fire episodes marking increasing fire frequency. However, the most recent levels of CHAR are still lower than those of the maximum in phase 2 and reach a minimum in the 20th century, meaning that the amount of modern charcoal accumulation is not unprecedented within the last ca. 2200 years. The mostly rather low SNI, which is below 3 for ca. one third of the record, might result from the lower intensity surface fire regime found around Lake Khamra. Such fires probably create less distinct peaks than the high intensity crown fires of other regions, especially considering the large lake size.

Robust CHAR (Fig. 3d–f), incorporating uncertainties from the age-depth model and charcoal counting, necessarily loses the original charcoal record's short-term variability. It needs to be noted that any uncertainty potentially arising from the assumed constant rate of old carbon input to the lake, underlying the sediment core's chronology, is not included here, as any changes in magnitude of this reservoir-like effect are impossible to quantify within the scope of this study. This issue is common in studies in permafrost regions that use 14C age dating (Biskaborn et al., 2012; Colman et al., 1996; Nazarova et al., 2013; Vyse et al., 2020). In such instances, applying Pb/Cs age dating adds valuable non-carbon-related estimates of sediment accumulation rates for the upper part of a sediment core (Whitlock and Larsen, 2001). One way of potentially quantifying old carbon age offsets throughout a sediment core might be to use cleaned charcoal particles, if available in sufficient number, and evaluate their age difference to bulk 14C samples (similar to Gaglioti et al., 2014). Since charcoal is assumed to be delivered to the lake directly via the air, it might provide better estimates of the real sediment age, even when other plant macrofossils fail to do so.

The general trend of fire regime changes in south-west Yakutia over the last ca. 2200 years captured by the CHAR background component (Fig. 6c) and described above shares many similarities with a charcoal record from the evergreen forest in the Tomsk region (ca. 1500 km west of Lake Khamra; Feurdean et al., 2020). There, one period of exceptionally high CHAR was observed around 700 CE and then again starting around 1700 CE towards the present, in parallel with Lake Khamra's fire history. The onset of increased biomass burning around 1700 CE was also reconstructed using aromatic acids from an ice core on the Severnaya Zemlya archipelago (ca. 2200 km north-west of Lake Khamra), and it indicates a sudden decrease at the beginning of the 20th century (Grieman et al., 2017). However, the same study finds another maximum of biomass burning around 1500 CE, in contrast to a clear period without fire episodes from ca. 1250 to 1750 CE at Lake Khamra. Although comparisons across such long distances and between different climate zones, archives, and fire proxies are likely to show differing results, some of the described trends seem to be recorded at several sites worldwide. A global charcoal record synthesis for the last two millennia by Marlon et al. (2008) indicates decreasing biomass burning from ca. 1 CE towards the industrial era, when, after a maximum around 1850 CE, it decreases with the onset of the 20th century. A potential explanation for this similar trend during the most recent centuries, observed across many study sites from different regions of the world (e.g. Dietze et al., 2019), could be the onset of industrialization with an accompanying change in land use and subsequent fire management. However, charcoal records from Siberia are underrepresented in such synthesis studies (Marlon et al., 2016), which, together with a lack of comparable records closer to south-west Yakutia, underline the issue of sparse data in the region.

Even with the constraints imposed by the large source area at this study site, the reconstructed record-wide mean FRI of 43 years, incorporating both the exceptionally short mean FRI of phase 2 (14 years) and the long mean FRI of phase 3 (60 years, excluding the ca. 600-year-long period without identified fires), lies within the range of the few comparable studies in boreal European Russia or western Siberia. Even if the record-wide mean FRI would strongly overestimate the number of included fire episodes, we would not expect the true value to exceed the maximum estimate of 95 years. Barhoumi et al. (2019) found the shortest FRIs of the Holocene ranging from 40 to 100 years between 1500 CE and the present day in macroscopic charcoal records from the northern Ural region. A mean FRI of 45 years during recent centuries was inferred by Feurdean et al. (2020). Other studies using tree-stand ages and fire scars in tree-ring chronologies suggest mean FRIs of 80–90 years for mixed larch forests between the Yenisei and Tunguska rivers ca. 1000 km north-west of Lake Khamra since ca. 1800 CE, although the FRI of individual study sites could be as short as ca. 50 years (Kharuk et al., 2008; Sofronov et al., 1998; Vaganov and Arbatskaya, 1996). A mean FRI of 52 years was reported for the northern Irkutsk region ca. 300 km west of Lake Khamra in the 18th century (Wallenius et al., 2011) and 50–80 years for some sites in the north-eastern larch-dominated forests (Kharuk et al., 2011; Schepaschenko et al., 2008). In general, fire frequency tends to increase with decreasing latitude due to higher solar radiation, longer fire seasons, and higher flammability of dry biomass (Ivanova, 1996; Kharuk et al., 2016, 2011), which, in addition to the large source area, likely contributes to a relatively short mean FRI at Lake Khamra when compared to studies set further north. Furthermore, studies based on tree ring chronologies or fire scars are only able to detect fires where trees survived or for which deadwood within the range of existing chronologies remains, whereas sedimentary charcoal can potentially capture all types of fires.

Enache and Cumming (2007) explain how a large catchment to lake-area ratio might favour secondary deposition of compact and stable morphotypes, whereas fragile morphotypes are more prone to fragmentation during surface runoff and thus rather represent primary input through the air. However, the catchment to lake-area ratio at Lake Khamra (23:1) and the share of fragile charcoal particles (types F, M, and S alone make up >80 %) are both comparably large. This might indicate that morphotype distribution within the record is not controlled by potential filtering effects of secondary charcoal transport but rather by the type of biomass burning. This is also implied by the predominantly primary charcoal input through the air due to the densely vegetated surrounding slopes and mirrors the stable vegetation composition seen in the pollen record.

Experimental charring studies have shown how different types of vegetation produce varied charcoal appearances (Jensen et al., 2007; Mustaphi and Pisaric, 2014; Pereboom et al., 2020). Pereboom et al. (2020) found elongated charcoal particles after experimentally burning tundra graminoids, potentially hinting at the origin of the many elongated type F particles at Lake Khamra. However, these type F particles quite closely match the appearance of charred Picea needles reported by Mustaphi and Pisaric (2014). Together with the potential of more intense fires producing larger charcoal particles (Ward and Hardy, 1991), this could mean the two previously noted peaks of CHAR (ca. 1880 CE at 19.5–20.5 cm depth and 650 CE at 124.5–125 cm; both dominated by high shares of type F particles >500µm) are evidence of higher-intensity fires burning conifer trees more severely and within a few kilometres of the lake shore. Charring experiments with local vegetation and a regionally adapted morphotype classification scheme would potentially benefit future studies by providing clear ground-truthing for links between morphotypes and vegetation.

4.2 Drivers of fire regime variations

4.2.1 Vegetation

The overall stable vegetation composition during the time covered by the charcoal record, as implied by the pollen and NPP record (Fig. 5), indicates that vegetation changes were unlikely to be the main driver behind changes in the fire regime and/or that changes in fire regime did not lead to large-scale shifts in vegetation composition. Similarly, no prominent shift in charcoal morphotype composition, and hence in the type of biomass burned over time, can be inferred (Figs. 4 and 6c, d). Although some studies draw a similar conclusion (e.g. Carcaillet et al., 2001, in eastern Canada), this result contrasts with many other studies from the Eurasian and North American boreal zones, where vegetation changes were found to be closely connected to changes in fire regimes (Barhoumi et al., 2019, 2020; Feurdean et al., 2020; Gavin et al., 2007; Kelly et al., 2013). A reason for this might be that studies on longer timescales capture long-term concurrent trends in both fire and vegetation that are not observable in a record that only captures the last ca. 2200 years (e.g. glacial to interglacial changes in vegetation distribution and temperature; see Marlon et al., 2013). However, on a shorter, multi-decadal timescale, phases with more Cyperaceae pollen (sedges) in the Lake Khamra record and a higher ratio of evergreen to deciduous arboreal pollen types coincide with periods of high fire activity in phases 2 and 4 (see Fig. 6b). This could be due to either the ability of sedges to quickly settle on freshly disturbed and cleared-out forest areas, and/or sedges growing in wetter areas, possibly right at the lake shore, which are spared from fires (Angelstam and Kuuluvainen, 2004; Isaev et al., 2010; Ivanova et al., 2014). An increased proportion of evergreen trees might enable more intense crown fires. Indirectly, dry periods could lead to receding lake levels and thus an increase in both shoreland sedges and fire ignitions. However, such clear links are difficult to infer without hydrological data. In addition to differences in proxy source area and taphonomy between macroscopic charcoal and pollen grains, a variety of factors likely obscure traces of potential fire impacts: surface fires in the deciduous forests in central Yakutia mostly result in the elimination of only a share of the affected tree population depending on fire intensity (Matveev and Usoltzev, 1996). This might not be enough to leave behind a clear mark in the pollen record, which also represents a source area that is probably a lot larger than the area affected by fire. Herbs or shrubs, on the other hand, may recover too quickly for changes to be detected in our record with a median temporal sampling resolution of 6 years. Any potential mixing processes and the residence time of pollen grains before settling in the lake sediment may further diminish visibility of a fire impact (Campbell, 1999; Fægri et al., 1989). Furthermore, reconstructions of Larix dynamics based on fossil pollen are affected by the limited pollen dispersal distance of larch trees, as well as poor preservation of their pollen grains (Müller et al., 2010). This can lead to an underestimation of Larix in fossil pollen records (Edwards et al., 2000) and thus complicate the evaluation of larch tree dynamics. The relatively low proportion of fossil Larix pollen at Lake Khamra, despite Larix gmelinii being a predominant tree taxon within the study area, may well reflect that issue. This indicates how future studies aiming at specifically comparing past fire regime changes with Larix population dynamics may benefit from including plant macrofossil analysis, if possible (e.g. Birks, 2003; Stähli et al., 2006). Due to a lack of macrofossils in the Lake Khamra sediment core and their suggested prolonged terrestrial residence time, as implied by the chronology, this was not an option in the present study. These factors, together with a remaining ambiguity in morphotype classification, likely explain the rather weak correlation of the pollen and charcoal records.

4.2.2 Climate

Although it has been demonstrated that the timing and extent of supposedly ubiquitous warmer or cooler climatic phases are in fact heterogeneous (Guiot et al., 2010), evidence for their occurrence in Siberia is seen in several proxy studies (Churakova Sidorova et al., 2020; Feurdean et al., 2019; Kharuk et al., 2010; Osborn and Briffa, 2006), albeit less pronounced when it comes to vegetation response in the West Siberian Lowland (Philben et al., 2014). Neukom et al. (2019) show how such climatic periods arising from averaged reconstructions at many individual study sites are not spatially or temporally coherent on the global scale and conclude that environmental reconstructions “should not be forced to fit into global narratives or epochs”. This might be especially true for studies of a single site, using chronologies that have 14C reservoir effects. However, the low fire activity in the latter half of phase 3 (900–1750 CE) strikingly coincides with the Little Ice Age (LIA), when in many regions of the Northern Hemisphere a cooler climate prevailed from ca. 1400 to 1700 CE. In contrast to this, high fire activity during phase 2 not matching the proposed timing of the warmer Medieval Climate Anomaly (MCA; ca. 950 to 1250 CE) demonstrates the limitations of such comparisons based solely upon the 14C-dated segment of the charcoal record (estimates of LIA and MCA durations from Mann et al., 2009). Due to a lack of regional studies, the PAGES Arctic 2k temperature reconstruction (McKay and Kaufman, 2014) was used to provide a comparison between the reconstructed fire activity and large-scale changes in Arctic climate north of 60 N (Fig. 6a). This synthesis of circumpolar temperature reconstructions incorporates many records, although datasets from Siberia are sparse, and thus it is underrepresented compared to Greenland, North America, and Europe. However, climate at Lake Khamra is likely to be strongly influenced by the conditions further north as Arctic temperatures affect the strength of the AO, which has been indirectly linked to fire activity further south (Balzter et al., 2005; Kim et al., 2020). The reconstructed PAGES Arctic 2k temperature provides evidence for a colder climate around ca. 1600 CE, coinciding with the LIA and low fire activity at Lake Khamra. The following onset of increased fire frequency in phase 4 (ca. 1750 CE onwards) is concurrent with a gradual increase in Arctic temperatures during the last two centuries (Fig. 6a, c), although it does not exceed the maximum in fire frequency of phase 2. The older half of the charcoal record (before ca. 1000 CE) does not match the temperature reconstruction as clearly. This might be due to a more regionally constrained climate or a consequence of the lack of comparable data about humidity, which also affects fire activity (Brown and Giesecke, 2014; Power et al., 2008). A larch tree-ring-based, ca. 1500-year-long reconstruction of summer vapour pressure deficit in north-eastern Yakutia (ca. 2000 km from Lake Khamra) indicates that the modern, increasing level of drought stress is not yet unprecedented, being surpassed by a high vapour pressure deficit during the MCA (Churakova Sidorova et al., 2020). This is similar to trends observed in the fire reconstruction at Lake Khamra. Another possibility for a less clear relationship between Arctic temperature and fire regime changes in the older half of the record is that the assumed constant impact of old carbon on the 14C age dating might be less pronounced in this older millennium. Yet, an impact of colder Arctic temperatures on the reconstructed low fire activity in phase 2 at Lake Khamra seems probable. Since we can currently rely only on this one record, more palaeoclimatic reconstructions and high-resolution charcoal records from the region could greatly improve the validity of such inferred links between climate and the fire regime.

Figure 6Comparison of the charcoal record with climate, vegetation, and general human settlement phases. Vertical dashed lines mark the different phases of the fire regime. (a) PAGES Arctic 2k temperature (grey curve = original annual data; yellow line = LOESS). (b) Selected vegetation proxies. (c) Fire frequency: compiled, smoothed frequencies of fire episodes for the sum of particles and individual size classes and morphotypes. (d) Compiled classic CHAR background components for the sum of particles and individual size classes and morphotypes, scaled to equal extents for comparison. (e) General human settlement phases, referenced in the text. For (c)(d), lines represent size and morphotype classes, as well as the sum of charcoal particles.


4.2.3 Human activity

The discrepancy of last century's low CHAR just after a phase of increasing fire frequency and rising Arctic temperature could potentially be a sign of direct and indirect consequences of human activity around Lake Khamra. The anthropogenic influence on fire regimes may be the most difficult to quantify due to missing information about the kind and extent of human fire use and management throughout time and the complex disentanglement of other drivers such as climate and vegetation (Marlon et al., 2013). Although Lake Khamra is located in a sparsely populated region, humans have historically been shaping the surrounding landscapes by building winter forest tracks and roads (ca. 0–30 km distance), logging (ca. 10 km distance), and building villages and towns (ca. 30 and 40 km distance, respectively). It has also been shown that even very remote lake systems in Yakutia have been affected by human activity from industrialization (Biskaborn et al., 2021). Yakutia has been populated by humans since at least ca. 28 000 BCE (Pitulko et al., 2004), although the population first noticeably increased at the end of the LGM around 17 000 BCE (Fiedel and Kuzmin, 2007), eventually forming the indigenous hunting and reindeer herding tribes of Evens, Evenks, and Yukaghirs (Keyser et al., 2015; Pakendorf et al., 2006). Between 1100 and 1300 CE (in the first half of phase 3 in the charcoal record) the Sakha people moved in from the south, pressured by an expanding Mongol empire (Fedorova et al., 2013). They established a new and distinct form of semi-nomadic livelihood based on horse and cattle breeding (Pakendorf et al., 1999). Fire was mainly used at hearths to provide warmth and light but also to sustain grasslands for grazing (Kisilyakov, 2009; Pyne, 1996). Population density likely started to increase rapidly after Russian colonization in the early 17th century (Crubézy et al., 2010) and even more drastically with the onset of new industries such as wide-scale logging and mining in the 20th century (Pyne, 1996). When compared to the pastoralist societies that existed up to that point, anthropogenic influence on the fire regime likely increased at the end of phase 3 (ca. 1700 CE onwards) and throughout phase 4 after colonization and industrialization (Drobyshev et al., 2004). It has been shown that human livelihoods and the mentality towards fire use can often better explain shifts in fire regimes than population density alone (Bowman et al., 2011; Dietze et al., 2018). For example, a formerly smaller population could have relied on practices like slash-and-burn agriculture until there was a transition towards more industrialized, urban livelihoods and a new focus on active fire suppression to protect forestry resources despite an increasing population (Dietze et al., 2019; Marlon et al., 2013). This is also thought to explain a pronounced decrease in boreal fire activity within the last century in tree-ring studies from central Siberia (300 km west of Lake Khamra) and Fennoscandia (Wallenius, 2011; Wallenius et al., 2011). Forest roads and clearings could have acted as fire breaks, while the emergence of fire suppression in Russia was conceived as early as 1893 CE and later led to the founding of the first aerial firefighting unit in 1931 CE (Pyne, 1996). Adding to this, slash-and-burn agriculture was officially banned at the end of the 18th century but was likely still practised frequently up until the early 20th century (Drobyshev et al., 2004; Konakov, 1999; Kozubov and Taskaev, 1999). Higher fire activity after ca. 1750 CE, marking the onset of phase 4, therefore coincides with both a rapid increase in anthropogenic activity and land development, as well as a warming Arctic climate. Low CHAR within the last century on the other hand might be a consequence of a cultural shift towards seeing fire as a hazard to ban, control, and suppress. Whereas high fire frequency in the past corresponded with high amounts of biomass burned, the recent century has also seen increasing fire frequency with decreasing total biomass burned (Fig. 6c, d). This might indicate that the current fire regime differs from that experienced by indigenous and Sakha people a few hundred years ago, now potentially consisting of more frequent but smaller fires. A better understanding of fire use and management of the various Yakutian societies throughout history is needed to judge the extent of the human influence on fire regimes in relatively remote regions especially since it has the potential to obscure the effects of recent global warming in fire reconstructions.

5 Conclusions

With its continuous sampling scheme and high median temporal resolution, the macroscopic charcoal record at Lake Khamra provides first insights into changes in the boreal fire regime of the last ca. 2200 years in eastern Siberia where comparable data are still lacking. Recent levels of charcoal accumulation at Lake Khamra are not unprecedented within the last two millennia. The reconstructed fire regime changes do not coincide with large-scale shifts in vegetation composition, although short-term increases in evergreen trees and sedges broadly coincide with periods of increased biomass burning around 700 and 1850 CE (phases 2 and 4). Also, low fire activity from ca. 900 to 1750 CE (phase 3), expressed as long FRIs and low charcoal accumulation, corresponds to a colder Arctic climate during the LIA. Despite the generally low population density, increased anthropogenic forcing after the colonization of Yakutia by the Russians in the early 17th century might have contributed to an increase in fire frequency, together with rising temperatures. Although northern regions have been warming rapidly in recent decades, charcoal input to the lake has been minimal during the last century, coinciding with new fire management strategies and a ban on fire-related agricultural practices. The mean FRI of 43 years, with a maximum estimate of 95 years, is within the range of published literature for the wider region and incorporates a range of individual values of up to almost 600 years. The large lake size may be an important factor behind these generally shorter FRIs since it is associated with a large charcoal source area of several kilometres from the lake and thus captures more fires compared to more locally constrained studies. Overall charcoal accumulation (classic CHAR background component) and the frequency of identified fire episodes seem to be directly related to each other for the majority of the record.

Although this new charcoal record improves data availability from eastern Siberia, more reconstructions, especially from distinctly deciduous regions, are needed to form a detailed analysis of past fire regimes in the Siberian boreal forest. An improved understanding of both fire activity and its drivers throughout history will eventually enable a meaningful assessment of the presence and future of Siberian wildfires and their consequences.

Appendix A

Figure A1Classic CHAR comparing the alternative local threshold (red) to the global threshold (orange). Vertical dashed lines mark the different phases of the fire regime. (a) Classic CHAR peak component (dark-grey bars = signal; light-grey bars = noise; coloured lines = threshold versions). (b) SNI based on Kelly et al. (2011) (black horizontal line = SNI cut-off value of 3). (c) Classic CHAR sum (black line = interpolated CHAR; blue line = LOESS representing the CHAR background component; red and orange marks = fire episodes for local and global threshold versions, respectively; with colour = fire episodes with SNI >3; in grey = fire episodes with SNI <3).


Appendix B

Figure B1Correlations (Kendall's τ) of charcoal morphotype classes with selected pollen groups.


Code and data availability

The R script used to analyse the charcoal record presented in this study can be accessed in the Zenodo database: (Glückler and Dietze, 2021). The data presented in this study are available in the PANGAEA database (; Glückler et al., 2020) and will be uploaded to the Global Paleofire Database (, International Paleofire Network, 2021).


The supplement related to this article is available online at:

Author contributions

ED and UH conceived and designed the study. LP, UH, and SK organized the expedition to Yakutia. SAV and BKB collected the samples and conducted field measurements. RG and SAV sampled the sediment core and performed the lab analysis. RG conducted all the charcoal proxy analysis, supported by ED and SK. AA conducted the pollen and non-pollen palynomorph analysis. BW and ED analysed remote-sensing data and created maps of the study location. RG wrote the paper with inputs from ED. All the authors reviewed the final manuscript.

Competing interests

The authors declare that they have no conflict of interest.


Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Special issue statement

This article is part of the special issue “The role of fire in the Earth system: understanding interactions with the land, atmosphere, and society (ESD/ACP/BG/GMD/NHESS inter-journal SI)”. It is a result of the EGU General Assembly 2020, 3–8 May 2020.


We thank Nadine Bernhardt for her help working with the sediment core and Cathy Jenks for providing English proofreading, as well as the participants of the joint German–Russian expedition Yakutia 2018 for their support. We would further like to thank Philip Higuera, Angelica Feurdean, and one anonymous referee for providing valuable remarks for improvements to this paper. This study is a contribution to the PAGES-endorsed International Paleofire Network.

Financial support

This research has been supported by the European Research Council (grant no. Glacial Legacy: 772852) and the Deutsche Forschungsgemeinschaft (grant no. DI 2544/1-1: #419058007). Ramesh Glückler is funded by AWI INSPIRES (INternational Science Program for Integrative Research). Stuart Andrew Vyse is financially supported by the Earth Systems Knowledge Platform (ESKP) of the Helmholtz Foundation.

The article processing charges for this open-access publication were covered by the Alfred Wegener Institute, Helmholtz Centre for Polar and Marine Research (AWI).

Review statement

This paper was edited by Sandy Harrison and reviewed by Angelica Feurdean, Philip Higuera, and one anonymous referee.


Aakala, T., Pasanen, L., Helama, S., Vakkari, V., Drobyshev, I., Seppä, H., Kuuluvainen, T., Stivrins, N., Wallenius, T., Vasander, H., and Holmström, L.: Multiscale variation in drought controlled historical forest fire activity in the boreal forests of eastern Fennoscandia, Ecol. Monogr., 88, 74–91,, 2018. 

Administrative Center Lensk: Lensky District Map, available at: (last access: 21 October 2020), 2015 

Alexander, H. D., Mack, M. C., Goetz, S., Beck, P. S. A., and Belshe, E. F.: Implications of increased deciduous cover on stand structure and aboveground carbon pools of Alaskan boreal forests, Ecosphere, 3, 1–21,, 2012. 

Andreev, A. A., Morozova, E., Fedorov, G., Schirrmeister, L., Bobrov, A. A., Kienast, F., and Schwamborn, G.: Vegetation history of central Chukotka deduced from permafrost paleoenvironmental records of the El'gygytgyn Impact Crater, Clim. Past, 8, 1287–1300,, 2012. 

Andreev, A. A., Tarasov, P. E., Wennrich, V., and Melles, M.: Millennial-scale vegetation history of the north-eastern Russian Arctic during the mid-Pliocene inferred from the Lake El'gygytgyn pollen record, Glob. Planet. Change, 186, 103111,, 2020. 

Angelstam, P. and Kuuluvainen, T.: Boreal Forest Disturbance Regimes, Successional Dynamics and Landscape Structures: A European Perspective, Ecol. Bull., 51, 117–136, 2004. 

Appleby, P. G., Nolan, P. J., Gifford, D. W., Godfrey, M. J., Oldfield, F., Anderson, N. J., and Battarbee, R. W.: 210Pb dating by low background gamma counting, Hydrobiologia, 143, 21–27,, 1986. 

Balzter, H., Gerard, F. F., George, C. T., Rowland, C. S., Jupp, T. E., McCallum, I., Shvidenko, A., Nilsson, S., Sukhinin, A., Onuchin, A., and Schmullius, C.: Impact of the Arctic Oscillation pattern on interannual forest fire variability in Central Siberia, Geophys. Res. Lett., 32, 1–4,, 2005. 

Barhoumi, C., Peyron, O., Joannin, S., Subetto, D., Kryshen, A., Drobyshev, I., Girardin, M. P., Brossier, B., Paradis, L., Pastor, T., Alleaume, S., and Ali, A. A.: Gradually increasing forest fire activity during the Holocene in the northern Ural region (Komi Republic, Russia), The Holocene, 29, 1906–1920,, 2019. 

Barhoumi, C., Ali, A. A., Peyron, O., Dugerdil, L., Borisova, O., Golubeva, Y., Subetto, D., Kryshen, A., Drobyshev, I., Ryzhkova, N., and Joannin, S.: Did long-term fire control the coniferous boreal forest composition of the northern Ural region (Komi Republic, Russia)?, J. Biogeogr., 1–16,, 2020. 

Benaglia, T., Chauveau, D., Hunter, D. R., and Young, D.: mixtools: An R Package for analyzing finite mixture models, J. Stat. Softw., 32, 1–29, 2009. 

Birks, H. H.: The importance of plant macrofossils in the reconstruction of Lateglacial vegetation and climate: examples from Scotland, western Norway, and Minnesota, USA, Quaternary Sci. Rev., 22, 453–473,, 2003. 

Biskaborn, B. K., Herzschuh, U., Bolshiyanov, D., Savelieva, L., and Diekmann, B.: Environmental variability in northeastern Siberia during the last  13,300 yr inferred from lake diatoms and sediment–geochemical parameters, Palaeogeogr. Palaeocl., 329/330, 22–36,, 2012. 

Biskaborn, B. K., Narancic, B., Stoof-Leichsenring, K. R., Pestryakova, L. A., Appleby, P. G., Piliposian, G. T., and Diekmann, B.: Effects of climate change and industrialization on Lake Bolshoe Toko, eastern Siberia, J. Paleolimnol., 65, 335–352,, 2021. 

Blaauw, M. and Christen, J. A.: 2011 Flexible paleoclimate age-depth models using an autoregressive gamma process, Bayesian Anal., 6, 457–474,, 2011 

Blaauw, M., Christen, A., and Aquino L., M. A.: rbacon: Age-Depth Modelling using Bayesian Statistics. R package version 2.4.3, available at:, last access: 21 October 2020 

Blarquez, O., Girardin, M. P., Leys, B., Ali, A. A., Aleman, J. C., Bergeron, Y., and Carcaillet, C.: Paleofire reconstruction based on an ensemble-member strategy applied to sedimentary charcoal, Geophys. Res. Lett., 40, 2667–2672,, 2013. 

Blarquez, O., Vannière, B., Marlon, J. R., Daniau, A.-L., Power, M. J., Brewer, S., and Bartlein, P. J.: paleofire: An R package to analyse sedimentary charcoal records from the Global Charcoal Database to reconstruct past biomass burning, Comput. Geosci., 72, 255–261,, 2014. 

Bouchard, F., MacDonald, L. A., Turner, K. W., Thienpont, J. R., Medeiros, A. S., Biskaborn, B. K., Korosi, J., Hall, R. I., Pienitz, R., and Wolfe, B. B.: Paleolimnology of thermokarst lakes: a window into permafrost landscape evolution, Arct. Sci., 3, 91–117,, 2016. 

Bowman, D. M. J. S., Balch, J., Artaxo, P., Bond, W. J., Cochrane, M. A., D'Antonio, C. M., DeFries, R., Johnston, F. H., Keeley, J. E., Krawchuk, M. A., Kull, C. A., Mack, M., Moritz, M. A., Pyne, S., Roos, C. I., Scott, A. C., Sodhi, N. S., and Swetnam, T. W.: The human dimension of fire regimes on Earth, J. Biogeogr., 38, 2223–2236,, 2011. 

Brown, K. J. and Giesecke, T.: Holocene fire disturbance in the boreal forest of central Sweden, Boreas, 43, 639–651,, 2014. 

Brunelle, A. and Anderson, R. S.: Sedimentary charcoal as an indicator of late-Holocene drought in the Sierra Nevada, California, and its relevance to the future, The Holocene, 13, 21–28,, 2003. 

Campbell, I. D.: Quaternary pollen taphonomy: examples of differential redeposition and differential preservation, Palaeogeogr. Palaeocl., 149, 245–256,, 1999. 

Carcaillet, C., Bergeron, Y., Richard, P. J. H., Fréchette, B., Gauthier, S., and Prairie, Y. T.: Change of fire frequency in the eastern Canadian boreal forests during the Holocene: does vegetation composition or climate trigger the fire regime?, J. Ecol., 89, 930–946,, 2001. 

Chelnokova, S. M., Chikina, I. D., and Radchenko, S. A.: Geologic map of Yakutia P-48,49, 1: 000 000, VSEGEI, Leningrad, available at: (last access: 21 October 2020),1988 

Churakova Sidorova, O. V., Corona, C., Fonti, M. V., Guillet, S., Saurer, M., Siegwolf, R. T. W., Stoffel, M., and Vaganov, E. A.: Recent atmospheric drying in Siberia is not unprecedented over the last 1,500 years, Sci. Rep., 10, 15024,, 2020. 

Clark, J. S.: Particle motion and the theory of charcoal analysis: Source area, transport, deposition, and sampling, Quateranry Res., 30, 67–80,, 1988. 

Clark, J. S., Lynch, J., Stocks, B. J., and Goldammer, J. G.: Relationships between charcoal particles in air and sediments in west-central Siberia, The Holocene, 8, 19–29,, 1998. 

Colman, S. M., Jones, G. A., Rubin, M., King, J. W., Peck, J. A., and Orem, W. H.: AMS radiocarbon analyses from Lake Baikal, Siberia: Challenges of dating sediments from a large oligotrophic lake, Quaternary Sci. Rev., 15, 669–684,, 1996. 

Conedera, M., Tinner, W., Neff, C., Meurer, M., Dickens, A. F., and Krebs, P.: Reconstructing past fire regimes: methods, applications, and relevance to fire management and conservation, Quaternary Sci. Rev., 28, 555–576,, 2009. 

Crubézy, E., Amory, S., Keyser, C., Bouakaze, C., Bodner, M., Gibert, M., Röck, A., Parson, W., Alexeev, A., and Ludes, B.: Human evolution in Siberia: from frozen bodies to ancient DNA, BMC Evol. Biol., 10, 1–16,, 2010. 

de Groot, W. J., Cantin, A. S., Flannigan, M. D., Soja, A. J., Gowman, L. M., and Newbery, A.: A comparison of Canadian and Russian boreal forest fire regimes, Forest Ecol. Manag., 294, 23–34,, 2013. 

Dietze, E., Theuerkauf, M., Bloom, K., Brauer, A., Dörfler, W., Feeser, I., Feurdean, A., Gedminienė, L., Giesecke, T., Jahns, S., Karpińska-Kołaczek, M., Kołaczek, P., Lamentowicz, M., Latałowa, M., Marcisz, K., Obremska, M., Pędziszewska, A., Poska, A., Rehfeld, K., Stančikaitė, M., Stivrins, N., Święta-Musznicka, J., Szal, M., Vassiljev, J., Veski, S., Wacnik, A., Weisbrodt, D., Wiethold, J., Vannière, B., and Słowiński, M.: Holocene fire activity during low-natural flammability periods reveals scale-dependent cultural human-fire relationships in Europe, Quaternary Sci. Rev., 201, 44–56,, 2018. 

Dietze, E., Brykała, D., Schreuder, L. T., Jażdżewski, K., Blarquez, O., Brauer, A., Dietze, M., Obremska, M., Ott, F., Pieńczewska, A., Schouten, S., Hopmans, E. C., and Słowiński, M.: Human-induced fire regime shifts during 19th century industrialization: A robust fire regime reconstruction using northern Polish lake sediments, PLOS ONE, 14, 1–20,, 2019. 

Dietze, E., Mangelsdorf, K., Andreev, A., Karger, C., Schreuder, L. T., Hopmans, E. C., Rach, O., Sachse, D., Wennrich, V., and Herzschuh, U.: Relationships between low-temperature fires, climate and vegetation during three late glacials and interglacials of the last 430 kyr in northeastern Siberia reconstructed from monosaccharide anhydrides in Lake El'gygytgyn sediments, Clim. Past, 16, 799–818,, 2020. 

Drobyshev, I., Niklasson, M., Angelstam, P., and Majewski, P.: Testing for anthropogenic influence on fire regime for a 600-year period in the Jaksha area, Komi Republic, East European Russia, Can. J. Forest Res., 34, 2027–2036,, 2004. 

Duffin, K. I., Gillson, L., and Willis, K. J.: Testing the sensitivity of charcoal as an indicator of fire events in savanna environments: quantitative predictions of fire proximity, area and intensity, The Holocene, 18, 279–291,, 2008. 

Edwards, M. E., Anderson, P. M., Brubaker, L. B., Ager, T. A., Andreev, A. A., Bigelow, N. H., Cwynar, L. C., Eisner, W. R., Harrison, S. P., Hu, F.-S., Jolly, D., Lozhkin, A. V., MacDonald, G. M., Mock, C. J., Ritchie, J. C., Sher, A. V., Spear, R. W., Williams, J. W., and Yu, G.: Pollen-based biomes for Beringia 18,000, 6000 and 0 14C yr BP, J. Biogeogr., 27, 521–554,, 2000. 

Ehlers, J. and Gibbard, P. L.: The extent and chronology of Cenozoic Global Glaciation, Quaternary Int., 164/165, 6–20,, 2007. 

Enache, M. D. and Cumming, B. F.: Charcoal morphotypes in lake sediments from British Columbia (Canada): an assessment of their utility for the reconstruction of past fire and precipitation, J. Paleolimnol., 38, 347–363,, 2007. 

Fægri, K., Iversen, J., Kaland, P. E., and Krzywinski, K.: Textbook of pollen analysis, Edn. 4th, John Wiley & Sons Ltd., Chichester, UK, 1989. 

Fedorov, A. N., Vasilyev, N. F., Torgovkin, Y. I., Shestakova, A. A., Varlamov, S. P., Zheleznyak, M. N., Shepelev, V. V., Konstantinov, P. Y., Kalinicheva, S. S., Basharin, N. I., Makarov, V. S., Ugarov, I. S., Efremov, P. V., Argunov, R. N., Egorova, L. S., Samsonova, V. V., Shepelev, A. G., Vasiliev, A. I., Ivanova, R. N., Galanin, A. A., Lytkin, V. M., Kuzmin, G. P., and Kunitsky, V. V.: Permafrost-Landscape Map of the Republic of Sakha (Yakutia) on a Scale 1:1 500 000, Geosciences, 8, 1–17,, 2018. 

Fedorova, S. A., Reidla, M., Metspalu, E., Metspalu, M., Rootsi, S., Tambets, K., Trofimova, N., Zhadanov, S. I., Kashani, B. H., Olivieri, A., Voevoda, M. I., Osipova, L. P., Platonov, F. A., Tomsky, M. I., Khusnutdinova, E. K., Torroni, A., and Villems, R.: Autosomal and uniparental portraits of the native populations of Sakha (Yakutia): implications for the peopling of Northeast Eurasia, BMC Evol. Biol., 13, 1–18,, 2013. 

Feurdean, A., Veski, S., Florescu, G., Vannière, B., Pfeiffer, M., O'Hara, R. B., Stivrins, N., Amon, L., Heinsalu, A., Vassiljev, J., and Hickler, T.: Broadleaf deciduous forest counterbalanced the direct effect of climate on Holocene fire regime in hemiboreal/boreal region (NE Europe), Quaternary Sci. Rev., 169, 378–390,, 2017. 

Feurdean, A., Gałka, M., Florescu, G., Diaconu, A.-C., Tanţău, I., Kirpotin, S., and Hutchinson, S. M.: 2000 years of variability in hydroclimate and carbon accumulation in western Siberia and the relationship with large-scale atmospheric circulation: A multi-proxy peat record, Quaternary Sci. Rev., 226, 105948,, 2019. 

Feurdean, A., Florescu, G., Tanţău, I., Vannière, B., Diaconu, A.-C., Pfeiffer, M., Warren, D., Hutchinson, S. M., Gorina, N., Gałka, M., and Kirpotin, S.: Recent fire regime in the southern boreal forests of western Siberia is unprecedented in the last five millennia, Quaternary Sci. Rev., 244, 1–16,, 2020. 

Fiedel, S. J. and Kuzmin, Y. V.: Radiocarbon date frequency as an index of intensity of paleolithic occupation of Siberia: did humans react predictably to climate oscillations?, Radiocarbon, 49, 741–756,, 2007. 

Flannigan, M., Stocks, B., Turetsky, M., and Wotton, M.: Impacts of climate change on fire activity and fire management in the circumboreal forest, Glob. Change Biol., 15, 549–560,, 2009. 

Frégeau, M., Payette, S., and Grondin, P.: Fire history of the central boreal forest in eastern North America reveals stability since the mid-Holocene, The Holocene, 25, 1912–1922,, 2015. 

Gaglioti, B. V., Mann, D. H., Jones, B. M., Pohlman, J. W., Kunz, M. L., and Wooller, M. J.: Radiocarbon age-offsets in an arctic lake reveal the long-term response of permafrost carbon to climate change, J. Geophys. Res.-Biogeo., 119, 1630–1651,, 2014. 

Gavin, D. G., Hu, F. S., Lertzman, K., and Corbett, P.: Weak climatic control of stand-scale fire history during the Late Holocene, Ecology, 87, 1722–1732,[1722:WCCOSF]2.0.CO;2, 2006. 

Gavin, D. G., Hallett, D. J., Hu, F. S., Lertzman, K. P., Prichard, S. J., Brown, K. J., Lynch, J. A., Bartlein, P., and Peterson, D. L.: Forest fire and climate change in western North America: insights from sediment charcoal records, Front. Ecol. Environ., 5, 499–506,, 2007. 

Giglio, L., Boschetti, L., Roy, D. P., Humber, M. L., and Justice, C. O.: The Collection 6 MODIS burned area mapping algorithm and product, Remote Sens. Environ., 217, 72–85,, 2018. 

Glückler, R. and Dietze, E.: CharcoalFireReconstruction in R (Version 1.0), Zenodo [Dataset],, 2021. 

Glückler, R., Herzschuh, U., Vyse, S. A., and Dietze, E.: Macroscopic charcoal record from Lake Khamra, Yakutia, Russia, PANGAEA [Dataset],, 2020. 

Goldammer, J. G. and Furyaev, V. V.: Fire in Ecosystems of Boreal Eurasia: Ecological Impacts and Links to the Global System, in: Fire in Ecosystems of Boreal Eurasia, edited by: Goldammer, J. G. and Furyaev, V. V., Springer Netherlands, Dordrecht, 1–20, 1996. 

Grieman, M. M., Aydin, M., Fritzsche, D., McConnell, J. R., Opel, T., Sigl, M., and Saltzman, E. S.: Aromatic acids in a Eurasian Arctic ice core: a 2600-year proxy record of biomass burning, Clim. Past, 16, 395–410,, 2017. 

Grosse, G., Jones, B., and Arp, C.: 8.21 Thermokarst Lakes, Drainage, and Drained Basins, in: Treatise on Geomorphology, edited by: Shroder, J. F., Academic Press, San Diego, 325–353, 2013. 

Guiot, J., Corona, C., and ESCARSEL members: Growing season temperatures in Europe and Climate forcings over the past 1400 years, PLOS ONE, 5, 1–15,, 2010. 

Hajdas, I., Zolitschka, B., Bonani, G., Leroy, S. A. G., Negendank, J. W., Ramratht, M., and Suter, M.: AMS radiocarbon dating of annually laminated sediments from lake Holzmaar, Germany, Quaternary Sci. Rev., 14, 137–143, 1995. 

Halsall, K. M., Ellingsen, V. M., Asplund, J., Bradshaw, R. H., and Ohlson, M.: Fossil charcoal quantification using manual and image analysis approaches, The Holocene, 28, 1345–1353,, 2018. 

Hansen, M. C., Potapov, P. V., Moore, R., Hancher, M., Turubanova, S. A., Tyukavina, A., Thau, D., Stehman, S. V., Goetz, S. J., Loveland, T. R., Kommareddy, A., Egorov, A., Chini, L., Justice, C. O., and Townshend, J. R. G.: High-resolution global maps of 21st-century forest cover change, Science, 342, 850–853,, 2013. 

Hawthorne, D., Courtney Mustaphi, C. J., Aleman, J. C., Blarquez, O., Colombaroli, D., Daniau, A.-L., Marlon, J. R., Power, M., Vannière, B., Han, Y., Hantson, S., Kehrwald, N., Magi, B., Yue, X., Carcaillet, C., Marchant, R., Ogunkoya, A., Githumbi, E. N., and Muriuki, R. M.: Global Modern Charcoal Dataset (GMCD): A tool for exploring proxy-fire linkages and spatial patterns of biomass burning, Quaternary Int., 488, 3–17,, 2018. 

Hély, C., Girardin, M. P., Ali, A. A., Carcaillet, C., Brewer, S., and Bergeron, Y.: Eastern boreal North American wildfire risk of the past 7000 years: A model-data comparison, Geophys. Res. Lett., 37, 1–6,, 2010. 

Herzschuh, U.: Legacy of the Last Glacial on the present-day distribution of deciduous versus evergreen boreal forests, Glob. Ecol. Biogeogr., 29, 198–206,, 2020. 

Herzschuh, U., Birks, H. J. B., Laepple, T., Andreev, A., Melles, M., and Brigham-Grette, J.: Glacial legacies on interglacial vegetation at the Pliocene-Pleistocene transition in NE Asia, Nat. Commun., 7, 1–11,, 2016. 

Higuera, P. E., Peters, M. E., Brubaker, L. B., and Gavin, D. G.: Understanding the origin and analysis of sediment-charcoal records with a simulation model, Quaternary Sci. Rev., 26, 1790–1809,, 2007. 

Higuera, P. E., Brubaker, L. B., Anderson, P. M., Hu, F. S., and Brown, T. A.: Vegetation mediated the impacts of postglacial climate change on fire regimes in the south-central Brooks Range, Alaska, Ecol. Monogr., 79, 201–219,, 2009. 

Higuera, P. E., Gavin, D. G., Bartlein, P. J., and Hallett, D. J.: Peak detection in sediment–charcoal records: impacts of alternative data analysis methods on fire-history interpretations, Int. J. Wildland Fire, 19, 996–1014,, 2011. 

Hoecker, T. J., Higuera, P. E., Kelly, R., and Hu, F. S.: Arctic and boreal paleofire records reveal drivers of fire activity and departures from Holocene variability, Ecology, 101, 1–17,, 2020. 

Hudspith, V. A., Belcher, C. M., Kelly, R., and Hu, F. S.: Charcoal reflectance reveals early Holocene boreal deciduous forests burned at high intensities, PLOS ONE, 10, e0120835,, 2015. 

International Paleofire Network: Global Paleofire Database, available at:, last access: 12 July 2021. 

Ito, A.: Modelling of carbon cycle and fire regime in an east Siberian larch forest, Ecol. Model., 187, 121–139,, 2005. 

Isaev, A. P., Protopopov, A. V., Protopopova, V. V., Egorova, A. A., Timofeyev, P. A., Nikolaev, A. N., Shurduk, I. F., Lytkina, L. P., Ermakov, N. B., Nikitina, N. V., E?mova, A. P., Zakharova, V. I., Cherosov, M. M., Nikolin, E. G., Sosina, N. K., Troeva, E. I., Gogoleva, P. A., Kuznetsova, L. V., Pestryakov, B. N., Mironova, S. I., and Sleptsova, N. P.: Vegetation of Yakutia: Elements of Ecology and Plant Sociology, in: The Far North: Plant Biodiversity and Ecology of Yakutia, edited by: Troeva, E. I., Isaev, A. P., Cherosov, M. M., and Karpov, N. S., Springer Science+Business Media B.V., 143–260,, 2010 

Ivanova, A. A., Kopylova-Guskova, E. O., Shipunov, A. B., and Volkova, P. A.: Post-fire succession in the northern pine forest in Russia: a case study, Wulfenia, 21, 119–128, 2014. 

Ivanova, G. A.: The extreme fire season in the central taiga forests of Yakutia, in: Fire in Ecosystems of Boreal Eurasia, edited by: Goldammer, J. G. and Furyaev, V. V., Springer Netherlands, Dordrecht, 260–270, 1996. 

Jansen, E., Christensen, J. H., Dokken, T., Nisancioglu, K. H., Vinther, B. M., Capron, E., Guo, C., Jensen, M. F., Langen, P. L., Pedersen, R. A., Yang, S., Bentsen, M., Kjær, H. A., Sadatzki, H., Sessford, E., and Stendel, M.: Past perspectives on the present era of abrupt Arctic climate change, Nat. Clim. Change, 10, 714–721,, 2020. 

Jensen, K., Lynch, E. A., Calcote, R., and Hotchkiss, S. C.: Interpretation of charcoal morphotypes in sediments from Ferry Lake, Wisconsin, USA: do different plant fuel sources produce distinctive charcoal morphotypes?, The Holocene, 17, 907–915,, 2007. 

Katamura, F., Fukuda, M., Bosikov, N. P., and Desyatkin, R. V.: Charcoal records from thermokarst deposits in central Yakutia, eastern Siberia: Implications for forest fire history and thermokarst development, Quaternary Res., 71, 36–40,, 2009a. 

Katamura, F., Fukuda, M., Bosikov, N. P., and Desyatkin, R. V.: Forest fires and vegetation during the Holocene in central Yakutia, eastern Siberia, J. For. Res., 14, 30–36,, 2009b. 

Keaveney, E. M. and Reimer, P. J.: Understanding the variability in freshwater radiocarbon reservoir offsets: a cautionary tale, J. Archaeol. Sci., 39, 1306–1316,, 2012. 

Kelly, R. F., Higuera, P. E., Barrett, C. M., and Hu, F. S.: A signal-to-noise index to quantify the potential for peak detection in sediment–charcoal records, Quaternary Res., 75, 11–17,, 2011. 

Kelly, R., Chipman, M. L., Higuera, P. E., Stefanova, I., Brubaker, L. B., and Hu, F. S.: Recent burning of boreal forests exceeds fire regime limits of the past 10,000 years, P. Natl. Acad. Sci. USA, 110, 13055–13060,, 2013. 

Kelly, R., Genet, H., McGuire, A. D., and Hu, F. S.: Palaeodata-informed modelling of large carbon losses from recent burning of boreal forests, Nat. Clim. Change, 6, 79–82,, 2016. 

Keyser, C., Hollard, C., Gonzalez, A., Fausser, J.-L., Rivals, E., Alexeev, A. N., Riberon, A., Crubézy, E., and Ludes, B.: The ancient Yakuts: a population genetic enigma, Philos. T. R. Soc., 370, 20130385,, 2015. 

Kharuk, V. I., Ranson, K. J., and Dvinskaya, M. L.: Wildfires dynamic in the larch dominance zone, Geophys. Res. Lett., 35, L01402,, 2008. 

Kharuk, V. I., Im, S. T., Dvinskaya, M. L., and Ranson, K. J.: Climate-induced mountain tree-line evolution in southern Siberia, Scand. J. Forest Res., 25, 446–454,, 2010. 

Kharuk, V. I., Ranson, K. J., Dvinskaya, M. L., and Im, S. T.: Wildfires in northern Siberian larch dominated communities, Environ. Res. Lett., 6, 045208,, 2011. 

Kharuk, V. I., Dvinskaya, M. L., Petrov, I. A., Im, S. T., and Ranson, K. J.: Larch forests of Middle Siberia: long-term trends in fire return intervals, Reg. Environ. Change, 16, 2389–2397,, 2016. 

Kim, J.-S., Kug, J.-S., Jeong, S.-J., Park, H., and Schaepman-Strub, G.: Extensive fires in southeastern Siberian permafrost linked to preceding Arctic Oscillation, Sci. Adv., 6, eaax3308,, 2020. 

Kisilyakov, Y.: Prescribed Fire Experiments in Krasnoyarsk Region, Russia, IFFN, 38, 51–62, 2009 

Konakov, N. D.: Agriculture. In Traditional culture of European North-East of Russia, available at: (last access 21 October 2020), 1999. 

Köster, E., Köster, K., Berninger, F., Prokushkin, A., Aaltonen, H., Zhou, X., and Pumpanen, J.: Changes in fluxes of carbon dioxide and methane caused by fire in Siberian boreal forest with continuous permafrost, J. Environ. Manage., 228, 405–415,, 2018. 

Kozubov, G. M. and Taskaev, A. I.: Forestry and forest resources of Komi republic, DiK Publishing House, Moscow, 1999 

Kruse, S., Bolshiyanov, D., Grigoriev, M. N., Morgenstern, A., Pestryakova, L., Tsibizov, L., and Udke, A.: Russian-German Cooperation: Expeditions to Siberia in 2018, Rep. Polar Mar. Res., 734, 1–257,, 2019. 

Kuuluvainen, T. and Gauthier, S.: Young and old forest in the boreal: critical stages of ecosystem dynamics and management under global change, For. Ecosyst., 5, 1–15,, 2018. 

Lenton, T. M.: Arctic climate tipping points, Ambio, 41, 10–22,, 2012. 

Lenton, T. M., Held, H., Kriegler, E., Hall, J. W., Lucht, W., Rahmstorf, S., and Schellnhuber, H. J.: Tipping elements in the Earth's climate system, P. Natl. Acad. Sci. USA, 105, 1786–1793,, 2008. 

Leys, B., Brewer, S. C., McConaghy, S., Mueller, J., and McLauchlan, K. K.: Fire history reconstruction in grassland ecosystems: amount of charcoal reflects local area burned, Environ. Res. Lett., 10, 114009,, 2015. 

Loader, C.: locfit: Local Regression, Likelihood and Density Estimation. R package version 1.5-9.4, available at:, last access: 21 October 2020. 

Mann, M. E., Zhang, Z., Rutherford, S., Bradley, R. S., Hughes, M. K., Shindell, D., Ammann, C., Faluvegi, G., and Ni, F.: Global signatures and dynamical origins of the Little Ice Age and Medieval Climate Anomaly, Science, 326, 1256–1260,, 2009. 

Marlon, J. R., Bartlein, P. J., Carcaillet, C., Gavin, D. G., Harrison, S. P., Higuera, P. E., Joos, F., Power, M. J., and Prentice, I. C.: Climate and human influences on global biomass burning over the past two millennia, Nat. Geosci., 1, 697–702,, 2008. 

Marlon, J. R., Bartlein, P. J., Daniau, A.-L., Harrison, S. P., Maezumi, S. Y., Power, M. J., Tinner, W., and Vanniére, B.: Global biomass burning: a synthesis and review of Holocene paleofire records and their controls, Quaternary Sci. Rev., 65, 5–25,, 2013. 

Marlon, J. R., Kelly, R., Daniau, A.-L., Vannière, B., Power, M. J., Bartlein, P., Higuera, P., Blarquez, O., Brewer, S., Brücher, T., Feurdean, A., Romera, G. G., Iglesias, V., Maezumi, S. Y., Magi, B., Mustaphi, C. J. C., and Zhihai, T.: Reconstructions of biomass burning from sediment-charcoal records to improve data–model comparisons, Biogeosciences, 13, 3225–3244,, 2016. 

Matveev, P. M. and Usoltzev, V. A.: Post-fire mortality and regeneration of Larix sibirica and Larix dahurica in conditions of long-term permafrost, in: Fire in Ecosystems of Boreal Eurasia, edited by: Goldammer, J. G. and Furyaev, V. V., Springer Netherlands, Dordrecht, 366–371, 1996. 

McKay, N. P. and Kaufman, D. S.: An extended Arctic proxy temperature database for the past 2,000 years, Sci. Data, 1, 140026,, 2014. 

Molinari, C., Carcaillet, C., Bradshaw, R. H. W., Hannon, G. E., and Lehsten, V.: Fire-vegetation interactions during the last 11,000 years in boreal and cold temperate forests of Fennoscandia, Quaternary Sci. Rev., 241, 106408,, 2020. 

Müller, S., Tarasov, P. E., Andreev, A. A., Tütken, T., Gartz, S., and Diekmann, B.: Late Quaternary vegetation and environments in the Verkhoyansk Mountains region (NE Asia) reconstructed from a 50-kyr fossil pollen record from Lake Billyakh, Quaternary Sci. Rev., 29, 2071–2086,, 2010. 

Mustaphi, C. J. C. and Pisaric, M. F. J.: A classification for macroscopic charcoal morphologies found in Holocene lacustrine sediments, Prog. Phys. Geogr. Earth Environ., 38, 734–754,, 2014. 

Nazarova, L., Lüpfert, H., Subetto, D., Pestryakova, L., and Diekmann, B.: Holocene climate conditions in central Yakutia (Eastern Siberia) inferred from sediment composition and fossil chironomids of Lake Temje, Quaternary Int., 290/291, 264–274,, 2013. 

Neukom, R., Steiger, N., Gómez-Navarro, J. J., Wang, J., and Werner, J. P.: No evidence for globally coherent warm and cold periods over the preindustrial Common Era, Nature, 571, 550–554,, 2019. 

Nilsson, S. and Shvidenko, A.: Is sustainable development of the Russian forest sector possible?, IUFRO, 11, Vienna, Austria, 1–60, 1996. 

Ohlson, M. and Tryterud, E.: Interpretation of the charcoal record in forest soils: forest fires and their production and deposition of macroscopic charcoal, The Holocene, 10, 519–525,, 2000. 

Osborn, T. J. and Briffa, K. R.: The spatial extent of 20th-century warmth in the context of the past 1200 years, Science, 311, 841–844,, 2006. 

Pakendorf, B., Spitsyn, V. A., and Rodewald, A.: Genetic structure of a Sakha population from Siberia and ethnic affinities, Hum. Biol., 71, 231–244, 1999. 

Pakendorf, B., Novgorodov, I. N., Osakovskij, V. L., Danilova, A. P., Protodjakonov, A. P., and Stoneking, M.: Investigating the effects of prehistoric migrations in Siberia: genetic variation and the origins of Yakuts, Hum Genet, 20, 334–353,, 2006. 

Pereboom, E. M., Vachula, R. S., Huang, Y., and Russell, J.: The morphology of experimentally produced charcoal distinguishes fuel types in the Arctic tundra, The Holocene, 30, 1091–1096,, 2020. 

Pestryakova, L. A., Herzschuh, U., Wetterich, S., and Ulrich, M.: Present-day variability and Holocene dynamics of permafrost-affected lakes in central Yakutia (Eastern Siberia) inferred from diatom records, Quaternary Sci. Rev., 51, 56–70,, 2012. 

Peters, M. E. and Higuera, P. E.: Quantifying the source area of macroscopic charcoal with a particle dispersal model, Quaternary Res., 67, 304–310,, 2007. 

Philben, M., Kaiser, K., and Benner, R.: Biochemical evidence for minimal vegetation change in peatlands of the West Siberian Lowland during the Medieval Climate Anomaly and Little Ice Age, J. Geophys. Res.-Biogeo., 119, 808–825,, 2014. 

Philippsen, B.: The freshwater reservoir effect in radiocarbon dating, Herit. Sci., 1, 1–24,, 2013. 

Pisaric, M. F. J.: Long-distance transport of terrestrial plant material by convection resulting from forest fires, J. Paleolimnol., 28, 349–354,, 2002. 

Pitulko, V. V., Nikolsky, P. A., Girya, E. Y., Basilyan, A. E., Tumskoy, V. E., Koulakov, S. A., Astakhov, S. N., Pavlova, E. Y., and Anisimov, M. A.: The Yana RHS site: humans in the Arctic before the Last Glacial Maximum, Science, 303, 52–56,, 2004. 

Ponomarev, E., Kharuk, V., and Ranson, K.: Wildfires Dynamics in Siberian Larch Forests, Forests, 7, 125,, 2016. 

Power, M. J., Marlon, J., Ortiz, N., Bartlein, P. J., Harrison, S. P., Mayle, F. E., Ballouche, A., Bradshaw, R. H. W., Carcaillet, C., Cordova, C., Mooney, S., Moreno, P. I., Prentice, I. C., Thonicke, K., Tinner, W., Whitlock, C., Zhang, Y., Zhao, Y., Ali, A. A., Anderson, R. S., Beer, R., Behling, H., Briles, C., Brown, K. J., Brunelle, A., Bush, M., Camill, P., Chu, G. Q., Clark, J., Colombaroli, D., Connor, S., Finsinger, W., Foster, D., Frechette, J., Hallett, D. J., Higuera, P., Hope, G., Horn, S., Inoue, J., Kaltenrieder, P., Kennedy, L., Kong, Z. C., Larsen, C., Long, C. J., Lynch, J., Lynch, E. A., McGlone, M., Meeks, S., Mensing, S., Meyer, G., Minckley, T., Mohr, J., Nelson, D. M., New, J., Newnham, R., Noti, R., Oswald, W., Pierce, J., Richard, P. J. H., Rowe, C., Turney, C., Urrego-Sanchez, D. H., Umbanhowar, C., Vandergoes, M., Vanniere, B., Vescovi, E., Walsh, M., Wang, X., Williams, N., Wilmshurst, J., and Zhang, J. H.: Changes in ?re regimes since the Last Glacial Maximum: an assessment based on a global synthesis and analysis of charcoal data, Clim Dynam, 30, 887–907,, 2008. 

Power, M. J., Marlon, J. R., Bartlein, P. J., and Harrison, S. P.: Fire history and the Global Charcoal Database: A new tool for hypothesis testing and data exploration, Palaeogeogr. Palaeocl., 291, 52–59,, 2010. 

Pyne, S. J.: Wild hearth. A prolegomenon to the cultural fire history of Northern Eurasia, in: Fire in Ecosystems of Boreal Eurasia, edited by: Goldammer, J. G. and Furyaev, V. V., Springer Netherlands, Dordrecht, 21–44, 1996. 

R Core Team: R: A language and environment for statistical computing, R Foundation for Statistical Computing, Vienna, Austria, available at:, last access: 21 October 2020. 

Reimer, P. J., Austin, W. E. N., Bard, E., Bayliss, A., Blackwell, P. G., Ramsey, C. B., Butzin, M., Cheng, H., Edwards, R. L., Friedrich, M., Grootes, P. M., Guilderson, T. P., Hajdas, I., Heaton, T. J., Hogg, A. G., Hughen, K. A., Kromer, B., Manning, S. W., Muscheler, R., Palmer, J. G., Pearson, C., Plicht, J. van der, Reimer, R. W., Richards, D. A., Scott, E. M., Southon, J. R., Turney, C. S. M., Wacker, L., Adolphi, F., Büntgen, U., Capano, M., Fahrni, S. M., Fogtmann-Schulz, A., Friedrich, R., Köhler, P., Kudsk, S., Miyake, F., Olsen, J., Reinig, F., Sakamoto, M., Sookdeo, A., and Talamo, S.: The IntCal20 Northern Hemisphere radiocarbon age calibration curve (0–55 cal kBP), Radiocarbon, 62, 725–757,, 2020. 

Remy, C. C., Fouquemberg, C., Asselin, H., Andrieux, B., Magnan, G., Brossier, B., Grondin, P., Bergeron, Y., Talon, B., Girardin, M. P., Blarquez, O., Bajolle, L., and Ali, A. A.: Guidelines for the use and interpretation of palaeofire reconstructions based on various archives and proxies, Quaternary Sci. Rev., 193, 312–322,, 2018. 

Revelle, W.: psych: Procedures for Personality and Psychological Research, Northwestern University, Evanston, Illinois, USA, available at: Version=2.0.7, last access: 21 October 2020. 

Reza, M. S., Afroze, S., Bakar, M. S. A., Saidur, R., Aslfattahi, N., Taweekun, J., and Azad, A. K.: Biochar characterization of invasive Pennisetum purpureum grass: effect of pyrolysis temperature, Biochar, 2, 239–251,, 2020. 

Rogers, B. M., Soja, A. J., Goulden, M. L., and Randerson, J. T.: Influence of tree species on continental differences in boreal fires and climate feedbacks, Nat. Geosci., 8, 228–234,, 2015. 

Russian Institute of Hydrometeorological Information: World Data Center, available at:, last access: 21 October 2020. 

Scheffer, M., Hirota, M., Holmgren, M., Nes, E. H. V., and Chapin, F. S.: Thresholds for boreal biome transitions, P. Natl. Acad. Sci. USA, 109, 21384–21389,, 2012. 

Schepaschenko, D. G., Shvidenko, A. Z., and Shalaev, V. S.: Biological Productivity and Carbon Budget of Larch Forests of Northern-East Russia, Moscow State Forest University, Moscow, 2008. 

Séjourné, A., Costard, F., Fedorov, A., Gargani, J., Skorve, J., Massé, M., and Mège, D.: Evolution of the banks of thermokarst lakes in Central Yakutia (Central Siberia) due to retrogressive thaw slump activity controlled by insolation, Geomorphology, 241, 31–40,, 2015. 

Sofronov, M. A. and Volokitina, A. V.: Wildfire ecology in continuous permafrost zone, in: Permafrost Ecosystems: Siberian Larch Forests, edited by: Osawa, A., Zyryanova, O. A., Matsuura, Y., Kajimoto, T., and Wein, R. W.,, Springer Science+Business Media B.V., Dordrecht, 59–82, 2010. 

Sofronov, M. A., Volokitina, A. V., and Shvidenko, A. Z.: Wildland fires in the north of Central Siberia, Commonw. Forest. Rev., 77, 211–218, 1998. 

Stähli, M., Finsinger, W., Tinner, W., and Allgower, B.: Wildfire history and fire ecology of the Swiss National Park (Central Alps): new evidence from charcoal, pollen and plant macrofossils, The Holocene, 16, 805–817,, 2006. 

Subetto, D. A., Nazarova, L. B., Pestryakova, L. A., Syrykh, L. S., Andronikov, A. V., Biskaborn, B., Diekmann, B., Kuznetsov, D. D., Sapelko, T. V., and Grekov, I. M.: Paleolimnological studies in Russian northern Eurasia: A review, Contemp. Probl. Ecol., 10, 327–335,, 2017. 

Takahashi, K.: Future perspectives of forest management in a Siberian permafrost area, in: Symptom of Environmental Change in Siberian Permafrost Region, edited by: Hatano, R. and Guggenberger, G., Hokkaido University Press, Sapporo, 163–170, 2006. 

Tinner, W., Hofstetter, S., Zeugin, F., Conedera, M., Wohlgemuth, T., Zimmermann, L., and Zweifel, R.: Long-distance transport of macroscopic charcoal by an intensive crown fire in the Swiss Alps – implications for fire history reconstruction, The Holocene, 16, 287–292,, 2006. 

Ulrich, M., Matthes, H., Schmidt, J., Fedorov, A. N., Schirrmeister, L., Siegert, C., Schneider, B., Strauss, J., and Zielhofer, C.: Holocene thermokarst dynamics in Central Yakutia – A multi-core and robust grain-size endmember modeling approach, Quaternary Sci. Rev., 218, 10–33,, 2019. 

Vaganov, E. A. and Arbatskaya, M. K.: Climate history and wildfire frequency in the central part of Krasnoyarsky Krai. Climatic conditions in growing season and seasonal wildfires distribution, Siberian J. Ecol., 3, 9–18, 1996. 

van den Boogaart, K. G., Tolosana-Delgado, R., and Bren, M.: compositions: Compositional Data Analysis. R package version 2.0-0, available at:, last access: 21 October 2020. 

Vyse, S. A., Herzschuh, U., Andreev, A. A., Pestryakova, L. A., Diekmann, B., Armitage, S. J., and Biskaborn, B. K.: Geochemical and sedimentological responses of arctic glacial Lake Ilirney, Chukotka (Far East Russia) to palaeoenvironmental change since ∼51.8 ka BP, Quaternary Sci. Rev., 247, 106607,, 2020. 

Waito, J., Girardin, M. P., Tardif, J. C., Conciatori, F., Bergeron, Y., and Ali, A. A.: Recent fire activity in the boreal eastern interior of North America is below that of the past 2000 yr, Ecosphere, 9, e02287,, 2018. 

Walker, X. J., Baltzer, J. L., Cumming, S. G., Day, N. J., Ebert, C., Goetz, S., Johnstone, J. F., Potter, S., Rogers, B. M., Schuur, E. A. G., Turetsky, M. R., and Mack, M. C.: Increasing wildfires threaten historic carbon sink of boreal forest soils, Nature, 572, 520–523,, 2019. 

Wallenius, T.: Major decline in fires in coniferous forests – Reconstructing the phenomenon and seeking for the cause, Silva Fennica, 45, 139–155,, 2011. 

Wallenius, T., Larjavaara, M., Heikkinen, J., and Shibistova, O.: Declining fires in Larix-dominated forests in northern Irkutsk district, Int. J. Wildland Fire, 20, 248–254,, 2011. 

Wang, S. and Hausfather, Z.: ESD Reviews: mechanisms, evidence, and impacts of climate tipping elements, Earth Syst. Dynam. Discuss. [preprint],, 2020. 

Ward, D. E. and Hardy, C. C.: Smoke emissions from wildland fires, Environ. Int., 17, 117–134, 1991. 

West, J. J. and Plug, L. J.: Time-dependent morphology of thaw lakes and taliks in deep and shallow ground ice, J. Geophys. Res.-Earth, 113, F01009,, 2008. 

Whitlock, C. and Anderson, R. S.: Fire history reconstructions based on sediment records from lakes and wetlands, in: Fire and Climatic Change in Temperate Ecosystems of the Western Americas, edited by: Veblen, T. T., Baker, W. L., Montenegro, G., and Swetnam, T. W., Springer, New York, NY, 3–31, 2003. 

Whitlock, C. and Larsen, C.: Charcoal as a fire proxy, in: Tracking Environmental Change Using Lake Sediments, Vol. 3: Terrestrial, Algal, and Siliceous Indicators, edited by: Smol, J. P., Birks, H. J. B., and Last, W. M., Springer Netherlands, Dordrecht, 75–97, 2001. 

Wirth, C.: Fire regime and tree diversity in boreal forests: implications for the carbon cycle, in: Forest Diversity and Function: Temperate and Boreal Systems, edited by: Scherer-Lorenzen, M., Körner, C., and Schulze, E.-D., Springer, Berlin, Heidelberg, 309–344, 2005. 

Wohlfarth, B., Skog, G., Possnert, G., and Holmquist, B.: Pitfalls in the AMS radiocarbon-dating of terrestrial macrofossils, J. Quaternary Sci., 13, 137–145, 1998. 

Woodward, C. and Haines, H. A.: Unprecedented long-distance transport of macroscopic charcoal from a large, intense forest fire in eastern Australia: Implications for fire history reconstruction, The Holocene, 30, 947–952,, 2020.  

Wu, B. and Wang, J.: Winter Arctic Oscillation, Siberian High and East Asian Winter Monsoon, Geophys. Res. Lett., 29, 1897,, 2002.  

Zeileis, A. and Grothendieck, G.: zoo: S3 Infrastructure for Regular and Irregular Time Series, J. Stat. Softw., 14, 1–27,, 2005. 

Short summary
Data about past fire activity are very sparse in Siberia. This study presents a first high-resolution record of charcoal particles from lake sediments in boreal eastern Siberia. It indicates that current levels of charcoal accumulation are not unprecedented. While a recent increase in reconstructed fire frequency coincides with rising temperatures and increasing human activity, vegetation composition does not seem to be a major driver behind changes in the fire regime in the past two millennia.
Final-revised paper