Galina A. Ivanova · Elena A. Kukavskaya ·Valery A. Ivanov · Susan G. Conard ·Douglas J. McRae

Abstract Forest fuel investigations in central and southern Siberian taiga of Scots pine forest stands dominated by lichen and feather moss ground vegetation cover revealed that total aboveground biomass varied from 13.1 to 21.0 kg/m 2 . Stand biomass was higher in plots in the southern taiga,while ground fuel loads were higher in the central taiga.We developed equations for fuel biomass (both aerial and ground) that could be applicable to similar pine forest sites of Central Siberia. Fuel loading variability found among plots is related to the impact and recovery time since the last wildf ire and the mosaic distribution of living vegetation.Fuel consumption due to surface f ires of low to high-intensities ranged from 0.95 to 3.08 kg/m 2 , that is, 18—74% from pref ire values. The total amount of fuels available to burn in case of f ire was up to 4.5—6.5 kg/m 2 . Moisture content of fuels (litter, lichen, feather moss) was related to weather conditions characterized by the Russian Fire Danger Index(PV-1) and FWI code of the Canadian Forest Fire Weather Index System. The data obtained provide a strong foundation for understanding and modeling f ire behavior, emissions,and f ire effects on ecosystem processes and carbon stocks and could be used to improve existing global and regional models that incorporate biomass and fuel characteristics.

Keywords Pinus sylvestris · Tree bioma ss · Surface a nd ground f uel l oads · Fuel m oisture · Fuel c onsumption

Introduction

Boreal forest ecosystems contain large stores of carbon and occur in northern regions where climate is warming most rapidly and f ires are becoming more frequent and more severe (Groisman et al. 2017). Because Russia contains some two-thirds of the global boreal forest, accurate data on fuel characteristics and impacts of f ire in these systems are particularly important for improving our understanding of fuel consumption, emissions, and effects of f ire on regional and global carbon stocks. Even within Russia, forests vary widely in biomass and fuel characteristics and in f ire regimes—differences that are often not captured in global models or appreciated by forest managers or decisionmakers.Among about 1200 million ha of boreal forests in Russia, 115.2 million ha are dominated by Scots pine (Pinus sylvestris). About 18% of all forest f ires in Russia occur in Scots pine stands (Bartalev et al. 2015). About 65% of the Scots pine forests are in Siberia, east of the Ural Mountains(Forest Fund of Russia 1999). The ground covers most frequently burned in these stands are grass, lichen and moss(Korovin 1996). Although a number of studies have characterized specif ic fuel layers in these types of forests, few have described the loading of all components, compared fuels across or within sites, or used fuel data in conjunction with data on f ire behavior or fuel consumption from subsequent f ires. The data on fuel consumption in Siberia are very limited (FIRESCAN 1996; McRae et al. 2006; Ivanova et al.2011; Kukavskaya et al. 2013b).

The biomass in forest ecosystems is distributed vertically and includes one or more forest canopy layers (e.g.,tree foliage, living and dead tree branches), surface layers(e.g., grasses, shrubs and small trees, and dead down woody pieces), and ground layers (e.g., lichen, moss, humus, peat,and roots) (Byram 1959). Knowledge of the forest fuel complex is important for predicting forest f ire danger, potential forest f ire behavior, and the effects of f ire on carbon stocks and the atmosphere that are determined by fuel consumption. Fuel characteristics that affect f ire behavior include fuel load, continuity, structure, moisture content, and chemical composition (Davis 1959; Kurbatsky 1962; Albini 1976;Anderson 1982; Ottmar et al. 2007). A number of systems have been developed for classifying fuels and determining their inf luence on f ire characteristics or effects. These systems often combine fuel information with weather-based indices of f ire danger. In Russia, the PV-1 Index is used to estimate f ire danger (Vonsky et al. 1975). PV-1 is calculated from a combination of current and antecedent weather (temperature, humidity def icit, and precipitation), but does not include a fuel component. The Canadian Forest Fire Weather Index (FWI) System (Canadian Forestry Service 1987) is a weather-based system that divides the forest f loor (litter,moss, lichen, and duff) fuels into three layers (Van Wagner 1987) to account for different drying rates that are related to increasing compaction and decomposition with depth. This system was developed for Canada which has extensive boreal forests, but it has been applied successfully in many other regions and vegetation types (Ven?l?inen et al. 2014.). The National Fire-Danger Rating System (NFDRS) developed in the United States (Burgan and Rothermel 1984) divides fuels into two categories: dead and live. The dead fuel category consists of four diameter classes categorized by drying time(timelag), and live fuels include two classes (live herbaceous and live woody). Kurbatsky ( 1962, 1970) divided forest fuels into seven categories based on their position in the ecosystem: Group I, mosses, lichens, and small-sized plant residues; Group II, forest f loor, humus, and peat horizons;Group III, grasses and small shrubs; Group IV, large dead woody elements (fallen branches and trees, snags, stumps,logging slash); Group V, tree regeneration and tall shrubs;Group VI, living tree foliage and branches up to 7 mm in diameter; and Group VII, living tree stems and branches＞ 7 mm in diameter. This system is not linked directly to f ire behavior or fuel consumption models. Kurbatsky ( 1970)also assigned forest fuels to one of three categories based on their contribution to f ire behavior: f ire carriers, f ire supporters, or f ire retarders. Fire carriers include surface fuels that can independently carry a surface f ire, such as mosses,lichens, litter, cured grass, and logging slash. This category also includes organic layers such as peat and rootstocks of perennial grasses (sod), which can carry smoldering ground f ires. Fire supporters are fuels that can facilitate f ire spread through carrier fuels. This category includes dead and down woody branches and stems, some grasses and small shrubs(e.g.,Ledumsp.,Callunasp.,Vaccinium vitis-idaea), young conifer regeneration, standing snags, tall shrubs, and highly resinous living tree stems. Fire retarders include live plants in the grass/small shrub layer that have high moisture content at the time of f ire (e.g., bearberry [Arctostaphylos uvaursi]) as well as green foliage of deciduous woody species.

Fires never consume the entire biomass of a living forest,because much of it is too moist to burn, or too discontinuous to sustain a f ire. The mass of surface and ground fuels combusted is determined by their biomass loading, structure,and moisture content (Kurbatsky and Ivanova 1987). Some researchers in Russia have classif ied all dead plant materials (sometimes including living mosses) above the mineral soil as forest f loor (Rodin and Bazelevich 1965; Pozdnyakov et al. 1969; Vedrova et al. 2000). From a f ire viewpoint, it is very important to distinguish between the litter layer and the decomposing forest f loor (F and H layers) below, since the litter layer is the most important fuel component in actively carrying f ire along with lichens, mosses, and dried grasses,leaves, and needles. These are the primary fuels consumed under the surface f ire regime that accounts for about 80% of the annual area burned in Russia (Korovin 1996). Litter can consist of dead and down woody pieces, which are deposited on the surface of the moss and lichen layer, as well as any cured herbs and dead foliage that have fallen onto the layer. Often dead herbs and fallen conifer needles get caught and hang-up on whortleberry (V. vitis-idaea) shrubs found in Scots pine (P. sylvestris)—mixed herbs—V. vitis-idaeaforest stands. Litter can sometimes be directly deposited on top of the forest f loor soil surface in forest stands with no ground cover (i.e., sites devoid of any lichen or moss cover).The rate of litter and forest f loor accumulation is strongly dependent on forest type, understory composition and cover, stocking, and age class (Melekhov 1957). Fuel load and fuel properties are dynamic and can vary throughout the year. Herb biomass increases in the spring, while the moisture content gradually decreases during the growing season. Cured and fallen vegetation from previous years is compressed and incorporated into the forest f loor and its fuel load then decreases as a result of decomposition. These constant changes determine the natural fuel dynamics in a forest stand, especially those with large amounts of herbaceous and shrub vegetation types that create a loosely packed fuel layer in spring (Kurbatsky and Ivanova 1987). Increases in the biomass of living herbs and shrub foliage during the growing season can alter the microclimate and the overall moisture of the fuel bed sufficiently to inhibit ecosystem f lammability and f ire spread, especially on herb- and shrubrich sites (Deeming et al. 1972; Baranov 1976).

Fuel layers in Siberian taiga forests include f ine-textured hygroscopic materials (litter, mosses, lichens), whose moisture content can range from low (10—30%) to highly saturated (up 300% and more) and may change relatively rapidly in response to rain or temperature and relative humidity.On the other hand, living shrubs and herbaceous vegetation respond more slowly to the environment and may often have moisture contents of 120—400% when the hygroscopic fractions of litter, mosses and lichen are quite dry. Forest fuel layer development and structure are impacted by a number of factors, such as vegetation type, precipitation pattern, and f ire frequency. Needle-fall at the end of the growing season leads to a relatively even distribution of fallen needles across the litter surface of central Siberian Scots pine stands by mid-August (Elagin 1976). The rate of decomposition of this litter is signif icantly reduced by the high resin content of newly fallen needles and twigs (Atkin 1984). This results in the accumulation of duffthat is characterized by high water-retention capacity (Bannikova 1967). Amosov ( 1964)reported that loadings of f ine surface fuels in pine forests with lichen and grass ground cover that are lower than about 0.2 kg/m 2 are insufficient to allow f ire to spread well. However, Sofronov ( 1988) found that an evenly distributed, cured sedge layer with a fuel load of only 0.05 kg/m 2 could carry f ire if quite dry. According to Valendik and Gevel ( 1975) f ire spreads through dry grasses with moisture content less than 26%. The dynamics, structure, and composition of surface and ground fuels are complex and site-specif ic. Fuel loading and structure are determined by vegetation type, moisture and temperature conditions, soil type, permafrost presence,and time since the last f ire (Melekhov 1947, 1957; Kurbatsky 1962; Morgan et al. 2001; Ivanova et al. 2007). Site characteristics, terrain, and weather conditions at the time of the f ire, including wind speed, are also important factors in limiting or promoting f ire spread (Valendik 1990).

Our study focused on quantifying the forest fuel structure,characteristics and consumption in typical Scots pine—lichen(Cladoniasp.)—feather moss [Pleurozium schreberi(Brid.)Mitt.] forests of the Siberian central and southern taiga.These data improve the knowledge of the biomass and fuel characteristics that are necessary for the effective f ire management and accurate modeling in these extensive and frequently burned Scots pine forest types. This work is part of a larger interdisciplinary study, the Fire Effects in the Boreal Eurasia Region (FIRE BEAR) Project that was developed to quantify the impacts of f ire on carbon cycle, emissions,and forest health and sustainability in Central Siberia. The data discussed here provided a basis for studying f ire behavior, and f ire impact on carbon stocks and various ecosystem components. In addition, fuel consumption data are essential for improving f ire emission estimates in Siberia.

Materials a nd methods

Study s ites

We investigated forest fuel loads on three sites in Scots pine stands of central Siberia (Fig. 1). For convenience, we identify the study sites by the names of the nearest villages. The Yartsevo study site is in the central taiga subzone, west of the Yenisey River on the Sym plain of the Western Siberian Lowland (60°38′N, 89°41′E) 55 km northwest of the village of Yartsevo and is an area of intermixed sand islands,bogs, and lakes typical of vast areas of western Siberia. The climate is characterized by cool and humid conditions, with an average annual air temperature ranging from - 3.2 to- 5.7 °C and total precipitation of 400—500 mm. Precipitation occurs throughout the year, with up to 60% of it falling in the summer (Zhukov et al. 1969). There is no underlying permafrost. The dry Scots pine stands with lichen and feather moss-dominated ground vegetation are characterized by sandy podzolic soils. The site quality classes are IV and V (Anuchin 1982), which indicates poor growing conditions.Pine regeneration, consisting of up to 0.5-m high Scots pine,was uniformly distributed across the site. Taller regeneration(2.0—3.5 m high) occurred in small pockets across some of the plots, where gaps in the canopy had resulted from past f ires. Sparse tall shrubs included wild rose (Rosa acicularisLindl.) and willow (Salix capreaL.).

The Govorkovo and Khrebtovo sites are in the southern taiga subzone of the Lower Angara region. Compared to the central taiga, stands of the southern taiga are more productive (tree volume can be as much as 50—60% higher), and living ground cover is more diverse (Zhukov et al. 1969). The Govorkovo site (58°35′ N, 98°55′ E) is south of the Angara River and 60 km southeast of the village of Govorkovo. The Khrebtovo site (58°42′ N, 98°25′ E) is north of the Angara River and 20 km north of the village of Khrebtovo (Fig. 1).The climate of this region is strongly continental, with an average annual air temperature of - 2.0° to - 2.4 °C and annual precipitation of 320—380 mm (Central Siberia 1964).The soils are illuvial-ferruginous sandy podzols, which are more fertile than soils at the Yartsevo site. The site quality class is III (Anuchin 1982). Uniformly distributed regenerated pines consisted of up to 0.5-m high Scots pine, while taller trees had regenerated only in Plot 2 in Govorkovo. The Scots pine—lichen—feather moss stand on the Govorkovo site is on a f lat bench of a gradual slope. It is a xeric forest type,similar to the Yartsevo site. The tall shrub layer contains willow (S. capreaL.), dog rose (R. acicularisLindl.), and mountain ash (Sorbus sibiricaHedl.). The Khrebtovo site is on a gentle slope in more mesic conditions than at the other two sites, as indicated by the herb—feather moss ground cover. This site is level with minimal microrelief and with nanorelief formed by down dead wood of varying ages. The tall shrub layer averages 3.5 m high and contains mainly alder (Alnus fruticosaRupr.), with scattered individuals of willow, dog rose, honeysuckle (Lonicera tataricaL.), and spiraea (Spiraea mediaFranz Schmidt).

Fig. 1 Map of location of study sites and view of the Scots pine forests in Russia

Fig. 2 Schematic representation of fuel categories (crown, surface,and ground) used and fuel types within each category

Summer drought periods, which lead to periods of high f ire danger, are common on all study sites (Valendik 1990;Valendik and Ivanova 1996). The pure pine forest stands at Yartsevo were of f ire origin and included several age groups:450, 275, 185, 130, and 45 years of age. The average landscape wildf ire interval was approximately 50 years (Ivanova 2005). The last wildf ire that occurred on the experimental site was in 1956. The Scots pine stand at Govorkovo last burned in 1946 and included the following age groups:200, 130, 80, and 50 years. The landscape wildf ire interval was 30—70 years. The last f ire on the Khrebtovo site was in 1948. The forest stand included age groups: 205, 170, 90,and 55 years. The average landscape wildf ire interval was approximately 46 years.

We analyzed fuel loads above the mineral layer, including crown, surface, and ground fuels (Fig. 2). We obtained fuel data from 13 plots measuring that were established at the experimental sites in preparation for a series of planned prescribed f ires: nine plots in Yartsevo, two plots in Govorkovo, and two plots in the Khrebtovo study area. These plots represented a range of forest stand conditions within this Scots pine forest type. Plots were approximately 200 × 200 m(4 ha). A 25 × 25 m sampling grid was laid out on each plot using numbered metal pins (normally 49 pins in a 7 × 7 grid)(Fig. 3).

We estimated mean density of trees ≥ 10 cm DBH by a point-centered quarter (PCQ) method (Morisita 1957). The centers of PCQ plots were offset by 1 m from the points on the sampling grid for each plot (Fig. 3). We measured height,diameter, and distance from the center point for four trees at every other sample point, for a total of 25 sample grid points (100 trees) per plot. Plot characteristics for both live and dead trees are summarized in Table 1. We calculated stand diameter at breast height (DBH) as a quadratic mean of the average basal area. Mean height is reported as the mean height of trees at the quadratic mean diameter, and basal area (BA) was based on cumulative BA of individual trees (Vagin et al. 1978).

Fig. 3 Diagram of basic sampling layout for each plot that illustrates the location of various types of fuel sampling relative to the grid points.Each plot was bordered by a f ire line cleared to mineral soil. Point-centered quarter (PCQ) plots were established at alternate grid points

Table 1 Stand structure characteristics of Scots pine forests on individual 4-ha study plots at study sites in the central and southern taiga in Russia

Tree fue l s ampling

We estimated aboveground fuels for Scots pine by developing biomass regressions from harvested trees across a range of size classes (e.g., Walker and Stocks 1975; Stocks 1980). Sample trees covered the range of DBH classes present at the site. Sixteen trees, with DBH from 4 to 44 cm,were felled and destructively sampled. The stump portion of each tree was cut offat the ground surface and included in the total tree stem biomass. Any tree pieces broken offduring the felling operation were retrieved, and all tree components were placed on tarps to ensure that no materials were lost. For each tree, we measured DBH, diameter at stump height, total tree height, stump height, height to living crown, height to the lowest live and dead branch, and crown width. Stem diameter was measured at 2-m intervals starting at the point the tree was cut. The 2-m increments of tree bole were weighed in the f ield. Tree age was determined by counting growth rings on each stump. All branches were removed from the main stem of each tree and separated into standard fuel diameter classes (0.0—0.49, 0.5—0.99, 1.0—2.99,3.0—4.99, and 5.0—6.99-cm size classes) using a go-no-go gauge (McRae et al. 1979). All branch diameter classes for dead and living fractions and needles on small twigs, were weighed separately. After weighing, three subsamples were taken of each branch diameter class (for both live and dead branches) and of each bole section for dry mass determination. Cones found on branches were collected, weighed,and subsampled to weigh dry mass. All subsamples were returned to the laboratory and dried at 100 °C until they reached constant mass. Crown needle biomass was calculated using conversion coefficients. For each tree, 15—20 living twigs with needles were partitioned into twig and needle components to determine the ratio of needle biomass (dry mass) to that of the live 0.0—0.49 cm twig. The average ratio for each sample tree was used to estimate overall dry mass of needles and small twigs for that tree.

Surface and gr ound fue l s ampling

The line intersect method of Van Wagner ( 1968), as adapted by McRae et al. ( 1979), was modif ied slightly for measuring the dead and down woody fuels (DDWF) before and immediately after each f ire. Five-meter transects were placed at all 49 grid sample points on each plot (Fig. 3). Due to the random natural fall of material, the transect orientations were kept similar to simplify the sampling set-up and relocation of transect for later postf ire sampling. The number of fuel pieces ＜ 7.0 cm in diameter were tallied into standard 0.0—0.49, 0.5—0.99, 1.0—2.99,3.0—4.99, and 5.0—6.99 cm diameter fuel size classes along the entire length of each 5-m transect using a metal “go-no-go”gauge that had templates for each size class. The jack pine(Pinus banksiana) multiplication factor of McRae et al. ( 1979)was used to obtain the fuel load (kg/m2), as jack pine and Scots pine have similar small branch diameters. The diameters of woody pieces ≥ 7.0 cm in diameter were determined using calipers, and decomposition level (solid vs. rotten) was recorded.Surface vegetation and ground fuel loadings were sampled using methods of Walker and Stocks ( 1975) and Alexander et al. ( 1991). Twenty-f ive 20 × 25-cm fuel quadrats were laid out per plot, one at every other grid point (Fig. 3). The sample location was offset 2 m away from each grid point so as not to interfere with other sampling. The aboveground components of small vegetation (herbs and small shrubs, e.g.,Vacciniumspp.) within the sampling frame were collected before forest f loor samples. Litter (needles, cones, and bark pieces) was collected as a separate fraction from on top of any living ground cover. The projected covers of moss and lichen were determined by ocular estimation to the nearest 5%. The outer edge of the sample quadrat was then cut using a knife and duffwas stripped offin 2-cm layers until the mineral soil was reached.Forest f loor depth (to mineral soil) for each sample quadrat was determined by averaging four measurements per quadrat(one in each corner). Because these methods do not estimate fuel loads for live mosses and lichens separately from the dufflayer, we collected an additional 10 fuel samples near each experimental plot for which lichens and mosses were separated from the duff. The ratio of lichens + mosses to duffobtained was used to determine the contribution of different fuel types to the total loads on our experimental plots. To determine fuel consumption, we placed two depth-of-burn pins 1 m on either side of the end DDWF sampling pin.

Samples for surface and ground fuel moisture content were taken immediately before experimental burns, on areas adjacent to each experimental plot and with similar vegetation composition and structure. In addition, in 2003 we collected daily samples for fuel moisture for a month (July), and during other years, we collected up to 150 samples of different fuels annually adjacent to our experimental plots in similar vegetation. Herbs and small shrubs, mosses and lichens (collected from both sunny and shaded ground locations), and litter were sampled. Six samples per fuel category were taken and averaged to obtain data on moisture content. The samples were oven-dried in the laboratory, and moisture content of each sample was calculated as percentage of dry mass. Weather characteristics (temperature, relative humidity, pressure, wind speed and direction, and precipitation) were measured at an onsite weather station that was established near the experimental plots and maintained during June and July. For the early part of the f ire season, we used data from the nearby weather stations.

The moisture content of different surface and ground fuel types (litter, lichen, and moss) was analyzed as a function of weather conditions, as characterized by the Russian Fire Danger Index (PV-1). PV-1 was calculated as:

wheretis temperature (°C) at 13:00 h DST, τ is dew point temperature (°C), and κ is the precipitation coefficient (Vonsky et al. 1975).

Burning p rocedures

Experimental plots were burned over four consecutive years in June and July, which corresponds to the main f ire season for this region. Yartsevo sites were burned in 2000—2002(Plots 13 and 14 in 2000; Plots 2, 3, 6, 19 in 2001; and Plots 4, 20, 21 in 2002); Govorkovo and Khrebtovo sites were burned in 2002—2003 (Plots 1 and 2 in 2002; Plots 3 and 5 in 2003). All plots were burned using line ignition along the windward side to mimic equilibrium f ire behavior conditions. Fires were primarily surface f ires, but they ranged widely in f ire behavior and in f ire severity (McRae et al.2006; Ivanova et al. 2011).

Byram’s ( 1959) equation was used to calculate fireline intensity as a product of rate of spread, dry mass of consumed fuel, and heat of combustion of consumed fuel(McRae et al. 2006 updated). Fires were classif ied into lowintensity f ires with f ireline intensity being up to 2000 kW/m,moderate-intensity (2001—4000 kW/m), and high-intensity f ires (more than 4001 kW/m).

Results and di scussion

Tree biomas s and c rown fue ls

Both individual fuel components of the trees and total tree biomasses showed close relationships with tree DBH, height,and age (Table 2), with the highest correlation coefficients obtained for tree DBH. We explored a number of independent variables, and concluded that regressions based on the logarithm of DBH provided the best models for representing biomass components across the range of our sampled trees.Our goal was to f ind a monotonic functional form that f it the data as well as possible, while minimizing both proportional error and the number of parameters in the equation. The general form of the resulting series of equations for estimating total Scots pine tree biomass and biomass of individual components is:

whereBrepresents biomass (kg);b0 ,b1 andb2 are coeffi-cients that were estimated for the logarithms of the various biomass components; and DBH is the diameter at breast height (cm). The f itted values of the parameters (b0 ,b1 ,b2 ) for the logarithms of biomass components are listed in Table 3.

Table 2 Correlations ( R)between Scots pine component and tree DBH, height and age

We developed the same regressions for the dry Yartsevo and Govorkovo sites because they have similar biomass and separate regressions for the higher-moisture Khrebtovo site.Figure 4 shows illustrative plots of several of these regression relationships and the data points used for developing them. For ease of interpretation, these are plotted using the nontransformed variables rather than the log 10 values that we used for the regressions.

The total tree biomass on each plot was calculated based on Eq. 2 and tree diameter size distribution determined from the PCQ plots. Aboveground stand biomass of living trees varied from 8.0 to 16.0 kg/m2among our plots (Table 4).Trees with higher DBH (32 and 36 cm DBH), which were common on our site, accounted for most of the tree biomass over the site. Nine of 13 experimental plots had dead standing trees, ranging from 12 to 45 cm in DBH and accounting for 1—15% from the total stand basal area (Table 1). Their aboveground biomass accounted for 0.12—0.18 kg/m2over experimental plots (Table 4). The largest component of the living aboveground tree biomass was the stems (up to 80%)followed by living branches (10.6—13.2%), dead branches(2.1—5.4%), and needles (3.4—4.9%). Dead trees accounted for no more than 2% of the total aboveground stand biomass.Total stand biomass on the study plots in the southern taiga was higher than in the central taiga, most likely due to more favorable climate conditions, higher stand productivity and density (Table 1). However, as the stand in the southern taiga contained younger trees with smaller diameters (DBH was 25% less), the difference in stand biomass between our southern and central taiga sites was only 12%.

As mentioned in the “ Introduction”, crown fires are relatively rare on Scots pine sites in Siberia. This rarity is attributable to the generally low stand density, lack of ladder fuels, and high crowns, which are particularly characteristic of dry Scots pine forests. The typical height to live crown of the PCQ trees was 11.1 ± 0.3 m in the central taiga and 12.5 ± 0.3 m in the southern taiga (Table 1). The difference in height to crown is attributed to the greater degree of self-pruning in the higher-density southern taiga stands.If a crown f ire were to occur in these forests, the f ire would normally only consume the foliage and living and dead branches up to 1.0 cm in diameter (Kurbatsky 1970). This category would only account for 0.60—1.35 kg/m 2 on our experimental plots (Table 4). The remaining portions of a living tree are usually too large, moist, and widely spaced to be consumed, especially given the normally short residence time of a crown f ire.

While not a major component of our stands, small pockets of Scots pine saplings (＞ 1.5 m tall), caused when gaps were created by past f ires that killed the mature trees and allowed full sunlight to reach the ground, were scattered irregularly across many of the plots. These pockets occupied from 0.01 to 0.08 ha (0.3—2.0%) in Yartsevo on seven of nine plots and 0.31 ha on Plot 2 in Govorkovo, but not at Khrebtovo. The Scots pine saplings in these pockets varied across the experimental plots from averages of 2.1—5.3 cmin DBH and 4.0—7.6 m in height. The total biomass for this regeneration varied from 0.18 to 1.59 kg/m 2 in the pockets,but was very small at the plot level (0.006—0.101 kg/m 2 on plots where sapling pockets occurred; Table 5). These conifer regeneration pockets can be an important ladder fuel to initiate torching or crown f ire in the canopy trees under high wind conditions.

Table 3 Coefficients used to determine the logarithm of biomass of individual trees and their components as related to the logarithm of DBH on Scots pine forest sites (see Eq. 2 in the text)

Surface and gr ound fue ls

The living surface ground cover on the experimental plots ref lected the microsite conditions on the plot. The total projected cover of ground vegetation depended on stand canopy closure. The surface vegetation layer had a well-pronounced structure, which was controlled by small (1—2 m) variations in relief as well as impacts from the last wildf ire. Moist microdepressions were occupied by aLedumsp.—V. vitisidaeacommunity.

The total projected cover of the herb—grass—small shrub layer at Yartsevo ranged from 15 to 40%. The height of this layer was typically 20—35 cm. The layer was dominated by small shrubs, whose proportions varied signif icantly depending on growing conditions:V. vitis-idaeadominated on better drained sites (ecotopes);Vaccinium myrtillusL.dominated mesic sites,Ledum palustreL. dominated wetter sites, and small bog shrubs andSphagnumdominated very wet sites. Mosses (60—100% of projected cover) were dominated byP. schreberi. Although lichens occurred on all plots,the projected cover varied widely (20—100%) with complete dominance on sites severely burned by previous wildf ires,where low tree density allowed more full sunlight to reach the ground leading to lower soil moisture.

In Govorkovo, the projected cover of the grass—small shrub layer, which was dominated byV. vitis-idaea,V. myrtillus, andLedum palustrevaried from 30 to 80% of the total plot area depending on stand canopy closure. Feather mosses, mainlyP. schreberi, covered up to 100% of the plot areas. Lichens, dominated byCladonia rangiferina, covered up to 40% of the plots. On the Khrebtovo sites, lichens(mainlyC. rangiferina)covered less than 1% of the plots.Well-developed mosses, mainlyPleurozeum schreberiandHylocomium splendens(Hedw.) B. S. G., formed a layer 5—7 cm high covering up to 100% of the plot areas. The grass—small shrub layer was dominated by small shrubs (V.vitis-idaea,V. myrtillus)and herbs (Calamagrostis arundinacea(L.) Roth.,Iris ruthenicaKer-Gawl.,Pyrola rotundifoliaL.,Carex macrouraMeinch.).

Fig. 4 Plots of data and regressions used to estimate biomass components of individual Scots pine trees based on diameter at breast height (DBH). We used the same regressions for the dry Yartsevo and Govorkovo sites because they have similar biomass. Note the difference in vertical scale for the f ine fuel biomass, which included small diameter (Class 1) branches and foliage

The types of surface vegetation cover were similar across the experimental plots in different sites. However,the relative importance of individual cover types differed among plots, which is clear from the proportions of structural elements of the ground vegetation layer. The distribution of a few major cover types (e.g., feather moss and lichen) ref lects differences in microsite moisture regimes.The transitional ground cover types (e.g., whortleberry,lichen, feather moss) occurred in most cases as a result of the inf luence of ecosystem (or vegetation) factors related to the influence of the previous fire and edaphic factors. In Yartsevo and Govorkovo, fewer than 15 herb and small shrub species were found, but up to 10 varieties of feather moss and lichen species were noted. A similarity coefficient (Vasilevich 1969) of 60% found for the plots suggested that all the plots are f loristically comparable(Kovaleva et al. 2011). In the more mesic Khrebtovo site,there were 28 species of herbs and shrubs, and only three species of lichens, which grew on rotten woody debris rather than directly on the forest f loor.

The total surface and ground fuel loading (Table 6),including herbs and small shrubs, litter, moss, lichen, duff(F and H layers), and DDWF varied from 4.0 to 6.2 kg/m2in the central taiga and from 3.8 to 5.0 kg/m 2 in the southern taiga. Although stands are generally more productive in the southern taiga, fuel loads in the central taiga exceed those in the southern taiga by 13% (P＜ 0.1), which might be a result of lower decomposition rates in northern sites(Johansson 1994). Variability in fuel loading among the plots within the same site location was related to the impact and recovery time since the last wildf ire, microsite variability, and the mosaic distribution (heterogeneity) of living ground vegetation that establishes after f ires (Ivanova and Perevoznikova 1996). Mosses, lichens and duffaccounted for the major proportion (58.1—89.1%) of the total surface and ground fuel loading with small shrubs and herbs, litter, and dead down woody debris making up 0.3—1.7%,1.6—6.3%, and 8.3—34.7%, respectively (Table 6). Mosses and lichens accounted for 0.7—1.7 kg/m 2 . Duffloading was 71—158% higher than that of moss and lichen. While litter,lichen and moss are commonly consumed almost completely during most f ires even of low to moderate severity in Scots pine forest, the consumption of the dufflayer depends on its moisture content at the time of burning. Kurbatsky ( 1962)stated that where the feather moss is moist, or dry only at the tips, little to no duffconsumption occurs.

Surface and ground fuel loading varied both within and between individual plots. The within-plot coefficient of variation ranged from 24 to 46%. With this in-plot variation, 25 sample points provided sufficient accuracy (standard error 5.1—9.1% of the mean). The coefficient of variation for the herb/small shrub layer, with the latter as the dominant component, was high (50—170%) due to the non-uniform distribution and low fuel loading (up to only 0.07 kg/m 2 ) of this component of the fuel—vegetation complex.

The relatively small amount of litter (＜ 0.3 kg/m 2 ) on the Yartsevo site was dominated by pine cones (43%) and needles (32%) (Fig. 5). Bark contributed 24% and cured herbs including grasses and small shrubs only 1% to the litter biomass. The dominance of cones in the litter loading at this site could be a result of more prolif ic pine cone production that can occur at low stand densities (Pobedinsky 1979) as well as lower rates of the decomposition process in the central taiga compared to the southern taiga (Vedrova 2002). In Govorkovo, litter loading was dominated by needles (47%),but in Khrebtovo, leaves of the well-developed understory layer (mainly alder) contributed the most to the total litter loads (40%). Higher litter biomass (5—15% of the total ground and surface fuel load) is often found where Scots pine occurs on more productive sites farther to the south(Ivanova 2005).

Table 4 Aboveground stand biomass and biomass of individual crown components for living trees on our sample plots

The total ground fuel load on the more mesic pine/moss forests in central Siberia was 40—55% higher than we observed. In birch—pine forests with grasses as dominant ground cover, ground fuel loads varied from 1.8 to 2.5 kg/m 2(Buryak et al. 2003), 30—55% lower than loads at our Scots pine sites. In the northern taiga spruce forest dominated by moss vegetation, ground fuels loads can reach 6.8 kg/m 2(Marchenko and Karlov 1961). Both ground and surface fuel loads in mixed larch—pine forest with a dark-coniferous subcanopy layer in the Angara region are up to 65% higher than on the sites reported on here (Zhila et al. 2011). For boreal forests of Canada, Stocks and Hartley ( 1995) reported surface and ground fuel loads in an Ontario mature jack pine forest of 7.9 kg/m 2 , almost twice as high as on our study plots.

Table 5 Characteristics and biomass of young coniferous regeneration over 1.5 m tall in pockets across the study plots

We observed a strong relationship between ground fuel depth and ground fuel loading (Fig. 6). While fuel loads up to 6 cm depth in the central and southern taiga have similar relations with fuel depth, fuel loads ＞ 6 cm in the southern taiga are typically 10—25% lower than loads in the central taiga, perhaps a result of the higher rate of the decomposition in southern sites (Vedrova 2002). Ground fuel loading and stand basal area were highly correlated, but only for the plots situated on drier soils where lichen dominated the ground cover (Fig. 7).

The DDWF load within individual plots for our study area had a coefficient of variation ranging from 187 to 302%. This variation resulted from the non-uniform distribution of fuels across the study area and differences in the various diameter size class loadings (Table 7). DDWF ≥ 7 cm in diameter accounted for 90% of the DDWF load. However, these largepieces were widely scattered and present on only 22 ± 3% of the line intersects. The highest DDWF load for pieces ＜ 7 cm in diameter in the 1.0—2.99 cm diameter size class (up to 11% of the total). A signif icant relationship (R= - 0.64,P＜ 0.05) between the loading of down woody debris and stand basal area at our Yartsevo site in the central taiga(Fig. 8) showed that DDWF loads are negatively correlated with BA values, perhaps a result of old trees dying and their transition from the category “stand” to “down woody debris”which accumulates on the ground due to low decomposition rates (Vedrova 2002). We did not f ind the same relationship for plots in the southern taiga because they are representedby younger stands where the majority of dead trees have DBH of 10—15 cm compared to the average DBH of 26 cm in the central taiga plots.

Table 6 Surface and ground fuel loadings for central and southern taiga Scots pine stands in Siberia

Fig. 5 Composition of litter loading on different study sites

Fig. 6 Relationships between ground fuel depths and loadings in Scots pine stands of the central and southern taiga. Means for the central taiga are based on data from 199 sample quadrats and 88 sample quadrats for the southern taiga

Fig. 7 Relationship between ground fuel loading and stand basal area for the dry plots with lichens dominating the ground cover (Plots 2, 3, 13, and 14 in Yartsevo; Plot 2 in Govorkovo). Ground fuels include litter, live lichens, mosses and duff

Table 7 Down woody fuel loading (kg/m 2 ), by diameter size class for the central and southern taiga Scots pine stands in Siberia

Fig. 8 Relationship between the loading of down woody debris and stand basal area for the central taiga (Yartsevo) plots

Surface a nd ground f uel m oisture c ontent

Fuel moisture plays a key role in f ire ignition and behavior,since it determines fuel f lammability levels and the capability of the f ire to spread (carry), which ultimately affects rate of spread and f ireline intensity (Kurbatsky 1970). Living herbs/grasses and small shrubs are characterized by relatively high and stable moisture content and are relatively independent of short-term variability in weather conditions compared to dead forest fuels with a similar diameter.Important shrubs (Vaccinium vitis idaea,V. myrtillus,Ledum palustre) on our experimental plots had lower moisture content (110—140% during the growing season) than herbs(250—350%), primarilyC. arundinacea. The presence of large amounts of herbs and small shrubs can affect the microclimate, as leaf emergence in the spring can reduce air temperature and increase the relative humidity. Such an interaction can increase the fuel moisture in the forest stand,thus reducing f ire danger (Kurbatsky and Ivanova 1987).However, herbs and small shrubs on our plots did not contribute substantially to decreasing the natural f ire danger,because their total cover was low and distribution was spotty.

Feather moss moisture content on our plots was at least twice that of lichens, since the absorption capacity of lichens is lower compared to that of mosses (Vonsky et al. 1975;Volokitina 1985). Moisture content differences between sun and shady conditions for feather moss (60—70%) were higher than for lichen (13—44%). Lichens dried rapidly to below the moisture of extinction after precipitation in the central taiga Scots pine stands. According to Sofronova ( 1970), the moisture content variability in lichens over 24 h is about 30%higher than mosses. The moisture content of lichens and mosses in the Scots pine stands was lowest from 13:00 to 16:00 h; after 17:00 h, fuel moisture content began gradually to increase due to decreased air temperature and increased relative humidity. Moisture content in the upper lichen layer can change from between 150 and 160—40% (the ignition threshold for lichen) over an 8—12-h period (Sofronov et al.2005). One of our experimental burns on a site with predominantly lichen ground cover was conducted successfully only 2.5 h after 5 mm of rain had fallen (McRae et al. 2006),illustrating the importance of rapid changes in moisture content of lichen layers for potential f ire ignition and spread.

The correlation coefficient characterizing the relationship between PV-1 and fuel moisture content was - 0.55(P＜ 0.0005) for lichen, and - 0.54 (P＜ 0.005) for litter(Fig. 9). The moisture content of the litter layer at a given PV-1 value was 5—15% lower than that of the lichen. This result corresponds to results of Zhukovskaya ( 1970) in a laboratory after drying fuel samples and is a result of the high moisture absorption capacity of lichen compared to the litter layer. In addition, lichen contacts the wetter mineral layer, while litter is on the top of the surface layer. We found correlations between the moisture content of theP. schreberilayer at different depths with the Russian Fire Danger Index(PV-1), with correlation coefficients ranging between - 0.59 and - 0.69 (P＜ 0.0005) (Fig. 10). Our data agree with layerby-layer moss drying data obtained for Scots pine stands in the northern part of European Russia and in Buryatia (Kurbatsky 1962; Vonsky et al. 1975; Volokitina 1985).

Fig. 9 Fuel moisture differences related to the PV-1 Index (see Eq. 1 in the methods) for the litter ( a) and lichen ( b) in a Scots pine—lichen—feather moss forest stand

Of the different Canadian Forest Fire Weather Index System codes (DC, BUI, FWI, FFMC), the best correlation with fuel moisture content was found for FWI (from- 0.68 to - 0.81; Fig. 10). McRae et al. ( 2006) also found the best correlation between f ire parameters (fuel consumption,depth of burn) with FWI in Scots pine forest of Siberia.Based on these relationships between moisture content and PV-1, the moisture contents of hygroscopic litter, lichens,and mosses are clearly dependent on timing and amount of precipitation. The regressions developed could be used to forecast the moisture content of the different fuels depending on the weather conditions and the time when fuels approach the ignition threshold range. Ignition thresholds have been reported to range from 13 to 27% for litter, 25—40% for lichens, 10—40% for mosses, and 16—25% for cured herbs and grasses (Melekhov 1947; Kurbatsky 1962; Sofronov 1967; Zhukovskaya 1970; Valendik and Gevel 1975; Konev 1977).

Fuel c onsumption

Total fuel consumption completeness ranged from 18 to 74%compared to pref ire values (Table 8). Grasses and herbs were mainly consumed completely during the f ire. Mosses and duffwere partially consumed depending on their moisture and f ire intensity. The highest fuel consumption (3.08 kg/m 2 )was registered for Plot 14 where almost total tree mortality had recently been observed (Kukavskaya et al. 2014). The consumption was similar to the average fuel consumption value of 3.44 kg/m 2 measured during a high-intensity surface and crown f ire in a dry Scots pine forest on Bor Island just north of our Yartsevo site (FIRESCAN 1996).

Fuel consumption due to low-intensity f ires varied from 0.95 to 1.66 kg/m 2 , compared with 1.29 to 2.09 kg/m 2 for fires of moderate-intensity. Low- to moderate-intensity f ires result in slightly higher fuel consumption values in the southern taiga compared to central taiga due to more complete fuel combustion, that seems to be a result of drier fuels in southern Siberia. There are just two f ires of highintensity—not enough for a valid comparison. However, fuel consumption in the central taiga due to high-intensity f ire potentially could be greater due to the greater accumulation of fuels on the ground (Table 6).

Conclusion

The total aboveground biomass in Scots pine forest types of central and southern Siberia varied widely (from 13.1 to 21.0 kg/m 2 ). Stand biomass was somewhat higher on plots in southern taiga, while ground fuel loads were higher in central taiga. Variations in total biomass among plots within the same study area were due mainly to differences in stand density. Scattered groups of young Scots pine saplings, while not contributing substantially to stand biomass,have the potential to serve as ladder fuels that can carry surface f ires into the tree canopy. Needles and small-diameter branches that we would expect might be consumed by crown f ire, accounted for 18.4 ± 1.3% of potentially available fuel,suggesting that crown fuels will not contribute a great deal to overall fuel consumption and emissions over that from moderate to high-severity surface f ires.

Fig. 1 0 Fuel moisture changes in feather moss related to the PV-1 Index (I) and Canadian Fire Weather Index (FWI) (II) in a Scots pine—small shrub—feather moss forest stand ( a 0—2 cm layer; b 2—4 cm layer; c 4—6 cm layer)

Surface and ground fuels are the main contributor of fuel consumed during f ire in these Scots pine forest types(McRae et al. 2006). Their biomass varied nearly twofold across our central and southern taiga study sites. This layer was composed of widely varying amounts of mosses,lichens, duff, and DDWF (tree stems and branches). The high variability in fuel loads of these key surface and ground fuels indicates the importance of quantifying fuel structure and loads in different areas, rather than using broad regional averages, for input into models of emissions, carbon stocks,and carbon cycle impacts of f ire. Fuel consumption due to f ires of low to high intensity varied from 0.95 to 3.08 kg/m 2 at our experimental sites.

In high-intensity f ires, forest f loor materials (primarily surface vegetation, lichens and mosses, duffand DDWF up to 7-cm diameter), live foliage, and small branches may be completely consumed, while large downed material, live tree stems and big branches are usually not. Estimated biomass available to burn varied from 4.5 to 6.5 kg/m 2 across our experimental plots and included the loads of needles and small branches in the tree canopy, surface vegetation, moss,lichen, duff, DDW of up to 7-cm diameter and rotten stems as well as the half of the solid stem loads on the ground. This accounts for about a third of the total aboveground biomass.These fuel loads can support surface f ires over a large range of severities, as well as the occasional crown f ire in Scots pine—feather moss—lichen forest stands, depending on the burning conditions. After long droughts resulting in high f ire hazard (PV-1 ＞ 4000), surface and ground fuels in Scots pine forest can burn very severely, and where ladder fuels such aspockets of smaller trees that regenerated in openings from previous f ires are present, torching or transition to crown f ire may occur (Kurbatsky 1970).

Table 8 Fuel consumption values at the central and southern taiga Scots pine stands in Siberia

The equations we have developed for estimating separate crown fuel components from tree diameter and relating surface and ground fuel loadings to stand characteristics such as basal area and fuel depth, as well as average values on fuel loadings should be applicable to similar pine (P. sylvestris)forests in the central and southern taiga of Siberia, where specif ic f ield data are not available.

Feather moss surface fuels remain moist at much higher PV-1 than lichen surface fuels. Thus, lichen fuels are f lammable under a much wider range of weather conditions.Lichen fuels also respond more rapidly to short-term changes in weather. An understanding of relationships between fuel moisture, weather patterns, and f ire hazard indices such as PV-1 and FWI provides a foundation for modeling potential effects of seasonal weather patterns on f ire behavior for Scots pine sites in central and southern Siberia.

Data on biomass, fuel characteristics and consumption like we presented here can be used to estimate the risk of f ire ignition and spread, emissions, and f ire intensity, and to support fire suppression efforts. Site- or forest typespecif ic data are required for understanding and modeling f ire effects, emissions, and carbon stocks on individual sites and for improving existing global and regional models that require biomass or fuel characteristics. Similar data to those presented here are needed for other forest types in Siberia,including the extensive larch and spruce/f ir forests. Ongoing studies (Chernykh and Furyaev 2011; Ivanova et al. 2011;Valendik et al. 2013) are helping to f ill this gap. Postf ire data from this and other studies (Buryak et al. 2011; Zhila et al. 2011; Kukavskaya et al. 2013a) are required to assess the long-term f ire impacts on biomass fuel characteristics,carbon stocks, and vegetation dynamics.

AcknowledgementsCooperation and logistical support of the Russian Aerial Forest Protection Service (Avialesookhrana) and Russian Forest Service (Regional and Local Forestry Committees) is greatly appreciated. A special thanks to L. Bobkova, N. Koshurnikova, and E. Krasnoshchekova for their assistance in fuel sampling and to D.Randall for statistical analysis of tree data.

Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License ( http://creat iveco mmons.org/licen ses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.