Abstract
Fire is a major disturbance agent in the boreal forest, influencing many current and future ecosystem conditions and services. Surprisingly few studies have attempted to improve the accuracy of fire-event reconstructions even though the estimates of the occurrence of past fires may be biased, influencing the reliability of the models employing those data (e.g. C stock, cycle). This study aimed to demonstrate how three types of fire proxies – fire scars from tree rings, sedimentary charcoal and, for the first time in this context, fungal spores of Neurospora – can be integrated to achieve a better understanding of past fire dynamics. By studying charcoal and Neurospora from sediment cores from forest hollows, and the fire scars from tree rings in their surroundings in the southern Fennoscandian and western Russian boreal forest, we produced composite fire-event data sets and fire-event frequencies, and estimated fire return intervals. Our estimates show that the fire return interval varied between 126 and 237 years during the last 11,000 years. The highest fire frequency during the 18th–19th century can be associated with the anthropogenic influence. Importantly, statistical tests revealed a positive relationship between other fire event indicators and Neurospora occurrence allowing us to pinpoint past fire events at times when the sedimentary charcoal was absent, but Neurospora were abundant. We demonstrated how fire proxies with different temporal resolution can be linked, providing potential improvements in the reliability of fire history reconstructions from multiple proxies.
Introduction
Forest fires are a key disturbance in boreal forests, and characteristics of fire regimes are among the most important factors explaining the variation in forest and landscape structure and species composition (Aakala et al., 2018). Climate change is predicted to strongly influence fire occurrence in forest ecosystems in the future (Girardin et al., 2009; Khabarov et al., 2016). Boreal forests contain approximately one-third of the global forest area and one-third of terrestrial carbon stocks, and hence changes in fire activity will have a substantial impact on global carbon emissions (Flannigan, 2015). Fire regime changes can have additional impacts because high-intensity crown and low-intensity surface fires result in different net effects on climate as a consequence of their contrasting impacts on terrestrial albedo (Rogers et al., 2015). Given the current and future societal importance of forest fires and their long-term influence on many ecological processes, it is necessary to improve our understanding of fire occurrence over long time scales.
Documentary records can provide detailed information on fires and their occurrence over large areas, but their temporal coverage in the boreal zone is limited to 100–150 years at best (Wallenius, 2011). However, there are several ways of acquiring information about past fire activity from biological archives, such as tree rings or organic sediments. These archives differ in their resolution and their temporal coverage. Tree-ring-based reconstructions of past fire activity employ either information on the age structures of trees or stands that initiated after fire (Bergeron et al., 2004), the presence of fire scars on trees that were damaged but survived the fires (Aakala, 2018; Wallenius et al., 2010) or both (Lankia et al., 2012). While of high temporal resolution, tree-ring-based fire history reconstructions in the boreal forest rarely cover more than a few centuries (Wallenius et al., 2010). Low-intensity fires may not always leave fire scars, potentially causing bias in fire regime reconstructions concerning surface fires (Dieterich and Swetnam, 1984).
Charcoal particles that are formed during incomplete combustion of biomass in forest fires and become deposited within organic sediments are another widely applied proxy in fire reconstructions (e.g. Dietze et al., 2018). Every forest fire has a unique combination of fuels and temperatures that define the intensity and severity of the fire and produce a range of different charcoal forms (Feurdean et al., 2017; Keeley, 2009; Marcisz et al., 2017; Ohlson, 2012; Zackrisson, 1977). Although there is no consensus on how resistant soil charcoal particles are to fragmentation over millennia in different biogeographical settings, and the dispersal of airborne charcoal particles is somewhat difficult to assess, it is clear that charcoal provides evidence on fire occurrence in the past (De Lafontaine and Asselin, 2011; Ohlson, 2012; Oris et al., 2014; Patterson et al., 1987). Compared with fire scars, the temporal resolution of sedimentary charcoal is relatively poor, due to the uncertainty associated with the dating methods and natural sedimentation processes. In addition, evidence indicates that charcoal can be absent from the sediments during times of known forest fires (Ohlson, 2012; Ohlson and Tryterud, 2000), which leads to problems in evaluating the fire occurrence or characteristics based on the charcoal record alone.
Several authors have attempted to correlate different fire proxies such as fire scars and sedimentary charcoal for obtaining a more complete picture of the long-term fire occurrence (Brossier et al., 2014; Higuera et al., 2005; Remy et al., 2018), but only a few of them have been successful. Higuera et al. (2005) studied the calibration of fire scars and sedimentary charcoal from a small forest hollow (<0.1 ha wet depression in the forest) and concluded that the fire detection from sedimentary charcoal depends strongly on fire severity and that the fine-scale spatial patterns of lower-severity burns play an essential role in determining the charcoal signal of these events. An additional differences between these two proxies stems from differences in their origins: sedimentary charcoal is produced during medium- to high-intensity fires at temperatures between 280°C and 500°C (Whitlock and Larsen, 2001), whereas the formation of fire scars requires that some of the trees survive and are hence usually formed after low-intensity (surface) fires (Gill, 1974). Characterization of the fire regime and its changes is thus dependent on the fire intensity, and the methods applied, and integration of different proxies to a more reliable composite data set would significantly increase our understanding of the long-term fire history of the boreal forest.
Traces for the presence of pyrophilous fungal species (i.e. species whose growth depends on fire) could potentially serve as another proxy for wildfire events. The formation of ‘blossoms’ on charred wood after a fire is a well-known characteristic of some species of moulds. These include Neurospora spp., which is a mould known to grow on charred tree and shrub bark after a fire and thus could potentially serve as another fire proxy (Jacobson et al., 2004, 2006; Kuo et al., 2014; Luque et al., 2012). Ascospores of Neurospora remain dormant until exposed to moist heat above 65°C, that is, temperatures that in the boreal forests only occur during fires (Perkins et al., 1975). After a fire, ascospores of the Neurospora germinate and form characteristic ‘orange bloom’ on the black, charred tree trunks.
Here, we use composite data sets to estimate fire-event frequencies and return intervals over the Holocene in the boreal forest of northern Europe. In comparison with earlier studies that have employed two proxies (fire scars and sedimentary charcoal; for example, Higuera et al., 2005; Remy et al., 2018), we take one step further and use three independent fire event indicating proxies: fire scars, sedimentary charcoal and Neurospora from small forest hollows. We specifically select small forest hollows under a closed tree canopy to reveal local fire events. Studied time includes both periods of natural conditions and anthropogenic impacts on fire-event dynamics in the European boreal forest. We propose an approach for linking fire-scar and sedimentary fire proxies that is applicable to other similar studies.
Material and methods
Study area
The study sites are located in the boreal forest of Fennoscandia and the taiga forest in north-western Russia (Figure 1). All study sites are small forest hollows within closed-canopy forest stands. Polttiais hollow (unofficial name; 1 × 1 m; N62°16′33″, E36°45′3″) is located in the Vodlozersky National Park in western Russia and is surrounded mainly by stands of Scots pine (Pinus sylvestris) and birch (Betula sp.). Mean annual temperature is 3.4°C (17.6°C in July and –9.9°C in February), and the annual precipitation sum 660 mm (mean data 1981–2010, from the nearest city Vytegra; Peterson and Vose, 1997). Naava hollow (unofficial name; 15 × 20 m; N61°14′29″, E25°3′22″) and Kämmekkä hollow (unofficial name; 10 × 20 m; N61°14′33″, E25°3′19″) are located 90 m from each other in southern Finland. Both Kämmekkä and Naava hollow are surrounded by spruce dominated forest, situated within a protected area of old-growth forest. The forest cover consists mainly of old-growth mixed stands of Norway spruce (Picea abies) and deciduous species such as birch, aspen (Populus tremula) and alder (Alnus glutinosa and Alnus incana), and occasional large individuals of Scots pine. The mean annual temperature is around 4°C (16.5°C in July and –7°C in February), and the mean annual precipitation is 645 mm (Pirinen et al., 2012).

Location of studied sites in Fennoscandia and western Russia (a), 1 – Naava hollow (b), 2 – Kämmekkä hollow (c) and 3 – Polttiais hollow (d). We compared our findings with nearby sites studied earlier: 4 – Sudenpesä (Clear et al., 2015); 5 – Vesijako (Clear et al., 2013); 6 – Kukka hollow, 7 – Larix hollow, 8 – Mosquito hollow and 9 – Olga hollow (Kuosmanen et al., 2014).
Sediment sampling
Fieldwork was conducted in August–October 2015. We extracted a 125 cm long core from Kämmekkä hollow, a 127 cm long core from Naava hollow and a 143 cm long core from Polttiais hollow, using a Russian peat sampler. We examined the cores in the field for visible charcoal layers, placed them in PVC semi-tubes and wrapped in plastic. After fieldwork, the samples were stored in a cold room (air temperature of 4–5°C) for further analysis.
Sediment dating and age–depth models
All cores were dated using accelerator mass spectrometry 14C dating (14C AMS), in Poznań, Poland (Table 1). All 14C dates were calibrated using the IntCal13 calibration data set (Reimer et al., 2013) with a two σ (95.4%) confidence level, and the age–depth models of each sequence were developed using Bacon 2.2 package (Blaauw and Christen, 2011), in the R environment (R Core Team, 2014) (Figure 2). In addition, we used a biostratigraphic correlation to define the lowermost age of Kämmekkä hollow, based on Naava hollow and Picea expansion in southern Finland (Table 1; Supplementary material 1, available online; Seppä et al., 2009b; Stivrins et al., 2017).
14C Accelerator mass spectrometry and biostratigraphic dates used in age–depth models of Polttiais, Naava and Kämmekkä hollows.

Bacon age–depth models for (a) Polttiais hollow, (b) Naava hollow and (c) Kämmekkä hollow. The calibrated 14C AMS dates are shown in transparent blue, and the age–depth model is shown in grey (darker greys indicate more likely calendar ages; grey stippled lines show 95% confidence intervals). The red curve is the posterior mean chronology. The laboratory numbers for the 14C AMS dates are also shown (see Table 1). Biostratigraphic marker of Picea pollen values according to Naava hollow data and southern Finland Picea pollen-stratigraphical patterns (Supplementary material 1, available online).
Sediment analyses
We analysed consecutive 1-cm samples for sedimentary charcoal and non-pollen palynomorphs from Kämmekkä, Naava and Polttiais hollows. We extracted a record of sedimentary charcoal and categorized them according to size: (1) microscopic charcoal (20–150 μm), (2) macroscopic charcoal (>150 μm–1 mm) and (3) large charcoal pieces (>1 mm).
Microscopic charcoal was counted from the pollen slides. These particles are in the range from 20 to 150 μm, but commonly 20–80 μm in size due to pollen preparation method which involves centrifuge usage that can break larger charcoal fragments into smaller pieces. We prepared each sample (volume 1 cm3, thickness 1 cm), using common pollen preparation procedures (HCl 10%, KOH 10%, acetolysis 3 min, mounted in glycerol; Berglund and Ralska-Jasiewiczowa, 1986). We added known quantities of Lycopodium spores to each sample to allow calculation of charcoal concentrations (Stockmarr, 1971) and counted microscopic charcoal particles according to Finsinger et al. (2008).
Macroscopic charcoal analysis was used to estimate the concentration of charcoal pieces > 150 μm in size in each sample. For this, we treated each sediment sample (volume 1 cm3, thickness 1 cm) with dilute NaOCI to promote sediment bleaching and disaggregation before sieving at 150 µm. The sediment residue was added to 20 mL distilled water and decanted to a petri dish for charcoal analysis. Charcoal was identified as brittle, black crystalline particles with angular broken edges using a stereomicroscope at 30–60× magnification, (Swain, 1973).
Large charcoal pieces >1 mm in size were recorded during plant macrofossil analysis, which followed the procedure described by Gałka et al. (2017). Sediment samples (volume 3 cm3) were wet-sieved and analysed using a Nikon SMZ800 stereoscopic microscope at magnifications of 10–200 and a transmitted light microscope.
Finally, we identified and counted the Neurospora (HdV-55c in size 22–28 × 15–18 μm; with a grooved surface with longitudinal ribs; Miola, 2012; Shear and Dodge, 1927; Van Geel, 1978) and Gelasinospora (HdV-1, HdV-2, HdV-528 in size ca. 22–33 × 14–20 μm; with a pitted wall; Dowding, 1933; van Geel, 1972, 1978) ascospores from the pollen slides. Based on molecular methods, Gelasinospora and Neurospora species were treated as a synonym of Neurospora genus within the order Sordariales (Dettman et al., 2001; García et al., 2004). Hence, identified spores of Gelasinospora and Neurospora were merged and treated as Neurospora spp. Fungal spores were counted alongside pollen analysis and we stop counting them when at least 500 pollen grains were identified.
Fire dates from scarred trees
We sampled fire scars at the Polttiais hollow site in August 2015, at the same time when sediment samples were taken. For this, we searched the vicinity of Polttiais hollow for fire-scarred trees, at approximately 100 m radius around the hollow, but so that we remained in the same stand in which the hollow was located. When encountered, we extracted a partial stem disc that contained the tip of the fire scar for dating the year of fire. In the lab, samples were dried and sanded to fine grit (600). We visually cross-dated the year of fire under a microscope, using marker rings (Yamaguchi, 1991) from samples obtained from live trees in the same stand.
For the Naava and Kämmekkä hollow, we used the fire dates from Tuominen (1990). Based on these data, fires occurred in 1773, 1827, 1857 and 1886. Fires have been absent from the stand since 1886. These fire dates were originally obtained from tree-ring samples from 18 fire scarred pines in the study stand, by inspecting the changes in tree growth rates typically associated with surface fires as described by Mikola (1950).
Data analyses
We used CharAnalysis (Higuera, 2009), to identify fire events from the sedimentary charcoal record. This analysis is based on the charcoal peak screening process (Gavin et al., 2006), in which a threshold value is used to separate the background charcoal deposition (noise) from the occurrence of peaks that are indicative of actual fire events. CharAnalysis was originally designed for the analysis of lake sediments, and hence we adjusted the procedure for using it for small forest hollow data (see Supplementary material 2, available online). CharAnalysis was used for screening the macroscopic and microscopic charcoal, whereas the observations of Neurospora and large charcoal pieces were treated as a dichotomous fire indicator.
For assessing the usefulness of Neurospora as a fire indicator, we tested its connection with other fire indicators (fire scars, macro- and microcharcoal) using a permutation test. In the test, we used the number of samples (slices), where Neurospora and other indicators were observed, as our test statistic. As the slicing of the sediment core may cause different indicators to be in separate slices even if they indicate the same fire event, Neurospora and another indicators were considered to match if they were in the same or in the neighbouring samples. We assigned the fire-scar-based fire events into a sediment sample whose dating appears to match it the best. In the permutations, Neurospora observations were randomly assigned to sediment samples while keeping sequential Neurospora observations in sequent slices. The p-value was obtained as the proportion of random permutations where the number of matches were at least as high as observed. We conducted the test separately for each sequence and by pooling the area-wise test statistics together. Furthermore, we examined whether the existence of charcoal or fire scar in a sample increased the probability of observing Neurospora. This was estimated as the risk ratio RR:
where C is proportion of slices with Neurospora among a sample with charcoal in the same or neighbouring sample and N is the proportion of slices with Neurospora among samples without another fire indicator in the same or neighbouring sample.
Risk ratio above 1 indicates that Neurospora was observed more often when another fire indicator is present than when another fire indicator was absent. The confidence intervals (CIs) of the risk ratio were calculated using the function riskratio in package fmsb of R software.
While combining the three separate fire indicators into a composite indicator, a sample was considered as a ‘fire sample’ if at least one proxy indicated a fire. Using CharAnalysis, we then estimated fire frequencies, as well as the mean and median fire return (mFRI) intervals. Usually (and as implemented in the CharAnalysis software), the CIs for the mFRI are obtained, using bootstrapping. However, as these intervals omit the effect of timing uncertainty which is available from the Bayesian age–depth model, we augmented the computation of CIs by considering both the sampling distribution of the mFRI and the timing uncertainty. To quantify and compare these two sources of uncertainty, we also report the uncertainty intervals obtained only considering the timing uncertainty and the bootstrap intervals without the timing uncertainty.
While computing the time-varying fire frequency, CharAnalysis sums the total number of fires within a 1000-year period and smooths the obtained series with a Lowess smoother. Preliminary analyses showed that Naava, Kämmekkä and Polttiais sediment sequences all had higher sedimentation rate towards the present day. The topmost section of sediment sequence was formed mainly of moss peat matrix that was less decomposed and not suppressed as subsequent lowermost sediment layers (similarly to bogs). As a result, the sampling density was considerably higher during the last 200 years (sediment is less compacted at the top – higher temporal resolution than for lowermost sediment), and we therefore binned the samples into 120-year bins that correspond to the smallest of the 95% quantiles of the sampling resolution over the period. We show the timing uncertainty also in the fire frequency estimates in the form of highest posterior density intervals.
Results
Sediment chronology
Age–depth modelling of all three sediment sequences reveals overall superposition sedimentation. Few samples were excluded from the age–depth models as they indicated outliers (Table 1). Because plant remains for AMS 14C dating were not available throughout the sediment sequences, several bulk peat samples were dated. A bulk sample consists of a mixture of terrestrial plant remains and tend to be slightly older than the date of plant macrofossil. These differences may arise due to a composition of the bulk sample (e.g. roots, fungi, humic acids), trampling of animals or change in peat decomposition rate (see, for example, Väliranta et al., 2014). Hence, establishing a chronology of small forest hollow is somewhat challenging as possible hiatus and non-linear sedimentation may occur.
The chronology for the basal section of Kämmekkä hollow was established based on the rise of spruce (Picea) pollen values, as the AMS 14C date of wood was an outlier (too young age). According to the Naava hollow pollen data (Stivrins et al., 2017) and general knowledge of spruce migration into Finland (Seppä et al., 2009b), the appearance of spruce 4800 years ago was used as a biostratigraphic marker. Considering the proximity of Naava hollow, we adjusted the basal age of Kämmekkä hollow according to this biostratigraphic marker. An example of circumstances that result in younger-than-expected wood remains is the fall of a tree branch that penetrates into an older sediment section. Except for this younger piece of wood, we found no signs of significant disturbance in sediment or the pollen spectra (Stivrins et al., 2017).
Tree-ring and sedimentary data on fire events
The tree-ring record from the stand surrounding Polttiais hollow indicated four fires: 1654, 1734, 1914 and 1949. While the number of fires from tree rings was the same for both Finnish and Russian sites, their spread in time was very different. The record from Tuominen (1990) shows fires for Naava and Kämmekkä hollows in 1773, 1827, 1857 and 1886.
All sediment sequences contained charcoal (microscopic to large charcoal pieces) which were further analysed individually and combined into a composite data set (Figure 3). Interestingly, we observed several samples with macrocharcoal pieces in size of > 1 mm but no significant microcharcoal (20–150 μm) values and vice versa.

Fire proxies in Polttiais (a), Naava (b) and Kämmekkä (c) hollows, as a function of depth. Fire proxies are represented by quadrats – microscopic charcoal (20–150 μm), stars – macroscopic charcoal (150 μm–1 mm), diamonds – large charcoal pieces (>1 mm), asterisk – fire scar (assigned into a sediment sample whose age appears to match it the best) and filled circles – Neurospora.
The connection between Neurospora and other fire indicators
Neurospora abundance in samples was low and did not exceed eight spores per sample. In Kämmekkä hollow, nine out of 11, Polttiais hollow six out of seven and in Naava hollow 18 out of 19 Neurospora observations were in the same or neighbouring slice with another fire indicator. The permutation test yielded p-values 0.16 for Kämmekkä, 0.21 for Polttiais and 0.09 for Naava, whereas the p-value for pooled data was 0.01. Hence, even though the connection was not significant in separately analysed sequences, a statistically significant connection was found when all the sequences were pooled together.
The risk ratios for Kämmekkä hollow was 3.5 (95% CI 0.8–15.5), for Polttiais hollow 3.4 (CI 0.4–27.7) and Naava hollow 5.3 (CI 0.7–38.2) and when combining all the areas, the risk ratio was 4.3 (CI 1.5–11.8). All the risk ratios were well above one (3.4–5.3), but the CIs were wide and above one only for the pooled data.
Composite data for fire-event reconstructions
We used the composite of fire proxies (Figure 3) to compute the fire return intervals (FRIs) (Table 2). The median FRIs were smaller than the mean FRIs indicating a skewed distribution of FRIs. In addition, it appeared that the uncertainty related to the sampling distribution of median or mean FRIs was more prominent than the uncertainty related to the timing.
The mean and median fire return intervals with three different sources of uncertainty. Timing: posterior mean of the mean or median fire return interval along with the credibility interval computed from the Bacon chronologies. Bootstrap: the mean and median of the fire return intervals using the mean chronology along with the bootstrap-based confidence interval. Both: The posterior mean of the mean (median) fire return interval along with the uncertainty limits that considers the timing uncertainty and the uncertainty from the estimation of the sampling distribution of the mean (median).
The composite fire record from Polttiais hollow in western Russia showed that the fire frequency was lower 11–8 ka ago but increased 8–4 ka ago (Figure 4). Fire frequency slightly decreased 4–2 ka ago but increased substantially over the last 1000 years. In Finland, the FRI was lower in Kämmekkä hollow from 6.5 to 4 ka ago and increased from two to seven fires per 1000 years in a period from 4 to 1.5 ka ago. Fire frequency decreased from 1.5 ka to 600 years ago followed by an abrupt rise up to nine fires per 1000 years towards the present day. Naava hollow had on average four to five fire events per 1000 years 7.5 to 4 ka ago, with a slight decrease 3–2 ka ago and distinct increase in fire frequency afterwards, particularly over the last 600 years. A common feature in all fire reconstructions was the abrupt increase in fire frequency over the last 1000 to 500 years.

Fire frequency for Polttiais (a), Naava (b) and Kämmekkä (c), hollows. Dots at the top of boxes indicate the sampling density, a higher row of crosses – the inferred fires and the lower row of crosses – the fires in 120-year bins. Black curve – the posterior mean of the fire frequency and grey band – the 95% highest density interval. Horizontal axis – time in years before present. Vertical axis – fire frequency.
Discussion
Proxies of fire
Prior to further exploration of the results and their interpretation, it is necessary to discuss the validity and strength of proxies used in fire reconstructions. In palaeoecology, sedimentary charcoal has been used in numerous publications, but to a lesser extent as a compilation of different size categories. Empirical studies show that larger charcoal particles (>150 μm) fall out relatively close to its emission source (<100 m), and smaller particles (20–150 μm) can be windblown from a broader region (>100 m) (Clark and Patterson, 1997; Conedera et al., 2009; Whitlock and Larsen, 2001). Recently, Adolf et al. (2018) provided the first European-scale geospatial training set relating the charcoal signal in surface lake sediments to fire parameters recorded by moderate satellite resolution imaging spectroradiometer sensors. According to these findings, the source area for both microscopic and macroscopic charcoal particles is very similar and can regionally be 40 km in diameter. However, these results were obtained from somewhat open European vegetation conditions and lakes. In the current study, only small (< 0.1 ha) forest hollows from densely forested conditions were analysed with a closed forest canopy. It is difficult to assess the source area of the smallest charcoal particles as they can spread over long distance in the air, but considering that for pollen, the source area in forest hollow usually is limited to 100 m (Overballe-Petersen and Bradshaw, 2011), it seems likely that also the sedimentary charcoal from such depository might reflect local fire events. Small-scale variability in fire occurrence can be seen from our results where two nearby hollows Naava and Kämmekkä (located 90 m apart) did not record the same fire events. Nevertheless, as noted by Remy et al. (2018) and Pitkänen et al. (2001), fire events detected from terrestrial sedimentary environments are more spatially and temporally precise and robust than those detected from lake sediments.
The spatial and temporal precision of the fire record is not an issue for fire scar-based reconstructions. Trees record fires at specific locations and, when cross-dated, at an annual resolution. However, compared with the sediments, the time span of the tree-ring-based fire record is in most cases more limited rarely dating back more than the last several centuries (but see Wallenius et al., 2010). This time limitation is generally due to the short life span of the tree species and the disappearance of fire scars with advancing decomposition following tree death. In addition, also fire-scar-based reconstructions may suffer from imperfect detection, depending on the fire regime. In high-intensity fires, all trees may be killed, leaving no fire-scar record. In such instances, fires can still be deduced from tree age structures (Dansereau and Bergeron, 1993), but this information gradually disappears as post-fire age cohorts die. On the contrary, low-intensity fires may pass through a forest without leaving a scar (Piha et al., 2013), as large fire-adapted trees can sustain low-intensity fires without being damaged. In such a case, fires could locally produce charcoal and/or raise temperatures enough to induce growth of pyrolytic fungi (such as the Neurospora spp.), producing evidence of fire that is visible only in the sediment data.
In the field of non-pollen palynomorphs, ascospores of Neurospora (HdV-55c; with longitudinal ribs) and Gelasinospora (HdV-1, HdV-2, HdV-528; pitted spores) are usually counted separately. However, genetic data suggest that these two genera of Sordariaceae family should be united (Dettman et al., 2001; García et al., 2004). In a palaeoecological context, ascospores of Gelasinospora found in sediments have been associated with fires when found in sediment layers containing charcoal (Dietre et al., 2017; Kuhry, 1985, 1997; Van Geel, 1978) or during dry conditions (Van Geel, 1972). Shumilovskikh et al. (2015) demonstrated correspondence of Gelasinospora to charred layers and dry phases of Sphagnum peat bog development. Some studies, however, indicate that Gelasinospora species can be coprophilous as well as carbonicolous and lignicolous (Krug et al., 2004; Lundqvist, 1972). Likewise, Neurospora (HdV-55c) ascospores in palaeoecology are associated with fires, as they have been found in charred sediment layers (Van Geel, 1978). The palaeoecological context of these ascospores usually refer to bog sediment sequences and thus far have not been used as a forest fire proxy.
Most field collections of Neurospora species contain isolates from tropical and subtropical regions, originating from either prescribed or natural fires. In general, natural populations of Neurospora often occur where fires are an essential part of the ecosystem (Perkins et al., 1975). The ascospores do not germinate under ambient temperatures, but they do germinate after being exposed to moist heat at 65–70°C for a few minutes. Therefore, for comparatively long time, Neurospora spp. were thought to belong only to moist tropical and subtropical regions, but more recent literature has made it clear that particular species are common primary colonizers of trees and shrubs after the forest fires also in North America and Europe (Jacobson et al., 2006; Luque et al., 2012). In particular, Kuo et al. (2014) showed that Neurospora spp. can grow as a symbiont (endophyte) within Scots pine, rapidly shifting to a parasitic or saprophytic lifestyle following a fire. Moreover, Neurospora spp. ascospores were not detectable in the soil after wildfire in samples taken after wildfire from Indonesia and Finland, rather Neurospora spp. survive the wildfire within the tree trunk (Kuo et al., 2014).
Typical ‘orange bloom’ developing on the charred tree trunks consists mainly of conidia containing conidium spores – micro and macro conidiospores with a typical size of 3 μm × 4 μm and 5 μm × 9 μm, respectively (Lee, 2012; Maheshwari, 1999). In contrast to the ascospores, increased heat is not necessary for conidiospores to develop (Maheshwari, 1999). In our samples, none of the particles detected corresponded to a size characteristic to the conidium spores, only ascospores with their characteristic surface morphology (ribbed or pitted) and size 20–30 μm. All in all, the use of Neurospora spp. as a fire proxy requires a completion of a full life cycle, starting with germination of ascospore initiated by heat, continued by the formation and development of perithecium, and culminating in the maturation of new asci that produce significant amounts of ascospores.
Compared with the use of charcoal as a fire proxy, Neurospora spp. has the additional advantage that it can potentially be used to constrain the season of fire occurrence. The minimum temperature for Neurospora spp. development from the germinated ascospores is 4°C (Dix and Webster, 1995). In laboratory settings, the development of perithecium, generation of asci and spore maturation has taken approximately 1 month. It is likely that in a forest setting this process takes an even longer time (Lee, 2012). Therefore, hypothetically, if a wildfire has occurred in the boreal forest during late autumn/winter time, it is most likely that no Neurospora spp. perithecium will develop even if the heat treatment from a wildfire was present, leaving no trace from such late season fires. Hence, findings of Neurospora spp. indicate fires occurring during spring/summer.
Fire events in the Northern European boreal forest
For the period for which we had data for fire proxies, Neurospora was linked with the two other, well-known proxies. Based on this relationship and the autecology of this fungus (i.e. its dependence on high temperatures), we suggest that integrating these different proxies complements the fire record and helps to improve our understanding of fire occurrence in the boreal forest. The fire history reconstructions developed here indicated that, on average, the fire occurred every 126–237 years (median 52–137 years) in these studied forests (Table 2).
Analysed separately, the mean and median FRIs in the Polttiais hollow in western Russia were 237 and 137 years over the last 11 ka. For the Finnish sites, the mean FRIs were 126–143. These long-term estimates of FRI are within the range for modern boreal forest fire intervals of 50–200 years (Bonan and Shugart, 1989), but shorter than many previous long-term fire history reconstructions from these regions indicate (Pitkänen et al., 2002). For example, Kuosmanen et al. (2014) used microscopic charcoal particles to study the effect of local fire on the forest in western Russia (in the same region as Polttiais hollow; Figure 1) and suggested relatively few fire events, or even an absence of local fires during the Holocene (last ca. 12 ka), depending on the site. In a similar study, Clear et al. (2013) used macroscopic charcoal data from the small forest hollow to study the fire frequency variability in Vesijako, Finland. The mean fire frequency was estimated as 430 years in semi-natural conditions (5–2 ka ago), 180 years during anthropogenic influence (2 ka to 750 years ago) and no fires over the last 750 years (Clear et al., 2013). In fire-scar-based reconstruction, Wallenius et al. (2007) estimated a 50-year fire interval for the 17th and 18th century for the region, before the cessation of fires in the mid-19th century.
Variability of the FRI can be partly explained by the substantial changes in climate over the entire Holocene, and changes in human impact especially in the more recent past. During the 11.7–8.2 ka ago, air temperatures were lower than present in western Russia and Finland (Kuosmanen et al., 2016). The warmest period appeared within the 8.2–4.2 ka ago and during the last 4.2 ka temperatures decreased gradually (Seppä et al., 2009a). However, while analyses on fire-climate relationships often rely on temperature reconstructions as the primary driver of the forest dynamics (Kuosmanen et al., 2016), the use of recently accomplished precipitation reconstructions have shown that concomitant decreases in precipitation may induce increased fire occurrence during colder time periods in the European boreal forest (Aakala et al., 2018). This makes explaining changes in fire frequency as a function of climate variability less straightforward.
A much stronger signal in the fire record was the increased fire frequency in the last ca. 1000 years that is likely linked to increased human activities. While the timing and the intensity of changes in human activities have differed between the study regions, both have been subjected to similar fire-conducive human impact in the past (Tikkanen and Chernyakova, 2014; Wallenius, 2011). These have included slash-and-burn agriculture which has been an active use of fire to alter the landscape but also probably unintentionally from other human activities such as fire spread from hunters’ campfires. Huttunen (1980) notes that the Evo region, where Naava and Kämmekkä hollows are located, comprises only 1–2% arable land and slash and burn cultivation with the clearing of forest for agricultural purposes was carried out rarely, whereas around 20 km towards the south at Lammi, the land becomes more fertile and as such more suitable for cultivation. In line with this, Stivrins et al. (2017) showed that the first crop pollen indicating agricultural practices appeared only 400 years ago at Naava and Kämmekkä hollows, while in other areas of Finland, an intensification of slash and burn activity can be seen in the charcoal records after 1000 cal. yr BP due to the expansion of cultivation and establishment of more permanent settlements (Alenius et al., 2013; Lagerås, 1996; Taavitsainen et al., 1998).
Earlier studies have demonstrated that using various charcoal size categories with a combination of fire-scar data may aid in gaining more comprehensive information about the occurrence of forest fires. However, we suggest that including Neurospora as a complementary proxy for fire reconstructions can further improve our understanding of fire occurrence. An additional advantage of Neurospora is that the germination and spread of their spores occurs over more extended time (e.g. days to weeks) compared with charcoal production during a fire event and are hence less dependent on short-term weather (particularly wind) conditions, increasing the probability of detection from sediments. Using multiple types of evidence for tracking past fires still needs more experimental studies to understand the strengths and weaknesses of these different fire proxies.
Conclusion
In this study, we presented a novel approach for past fire event reconstruction. We utilized three types of independent fire proxies – fire scar, charcoal and fungal spore of Neurospora spp. – to create a composite data set that was statistically analysed and used to reconstruct fire events in two study areas in the boreal forest. While each proxy has its strengths and weaknesses that require further research, our findings suggest that adding an additional complementary proxy (here: the spores of Neurospora spp.) complements the fire record, suggesting that fire reconstructions based on a single proxy technique may provide underestimates about past fire activity.
Supplemental Material
Supplementary_material_1 – Supplemental material for Integrating fire-scar, charcoal and fungal spore data to study fire events in the boreal forest of northern Europe
Supplemental material, Supplementary_material_1 for Integrating fire-scar, charcoal and fungal spore data to study fire events in the boreal forest of northern Europe by Normunds Stivrins, Tuomas Aakala, Liisa Ilvonen, Leena Pasanen, Timo Kuuluvainen, Harri Vasander, Mariusz Gałka, Helena R Disbrey, Janis Liepins, Lasse Holmström and Heikki Seppä in The Holocene
Supplemental Material
Supplementary_material_2 – Supplemental material for Integrating fire-scar, charcoal and fungal spore data to study fire events in the boreal forest of northern Europe
Supplemental material, Supplementary_material_2 for Integrating fire-scar, charcoal and fungal spore data to study fire events in the boreal forest of northern Europe by Normunds Stivrins, Tuomas Aakala, Liisa Ilvonen, Leena Pasanen, Timo Kuuluvainen, Harri Vasander, Mariusz Gałka, Helena R Disbrey, Janis Liepins, Lasse Holmström and Heikki Seppä in The Holocene
Footnotes
Acknowledgements
We would like to thank the referees.
Funding
The research was funded by the Academy of Finland (project nos. 276255, 252629 and 275969), by the Kone Foundation. Additional support was provided by the University of Latvia project ‘Studies of the fire impact on the bog environment and recovery’ with partners JSC ‘Latvia’s State Forests’, The Nature Conservation Agency and Latvian Peat Association, the Latvian Council of Science project no. LZP-2018/1-0171, National basic funding for science Y5-AZ03-ZF-N-110, the Estonian Research Council grants IUT1-8 and PRG323, and COST CA18135 Fire in the Earth System: Science & Society.
Supplemental material
Supplemental material for this article is available online.
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
