Macrocharcoal Signals in Histosols Reveal Wildfire History of Vast Western Siberian Forest-Peatland Complexes

Fires are a naturally cyclical factor regulating ecosystems’ function and forming new postfire ecosystems. Peat soils are unique archives that store information about ecological and climatic changes and the history of past fires during the Holocene. The paper presents a reconstruction of the dynamics of fires in the subzone of the middle taiga of Western Siberia in the Holocene. Data on fires were obtained based on the results of a study of the content of macroscopic coal particles and radiocarbon dating. The effect of fires on soil organic matter (SOM) was estimated using 13C NMR spectroscopy and the content of polyaromatic hydrocarbons (PAHs). It is shown that throughout the Holocene, the peatlands studied were prone to fires. The conducted analyses show that the maximum content of charcoal particles is observed in the Atlantic (~9100–5800 cal. B.P.) and Subatlantic (~3100 cal. B.P. to the present) periods. The high correlation dependence of the content of coals with the content of PAHs (r = 0.56, p < 0.05) and aromatic structures of SOM (r = 0.61, p < 0.05) in peat horizons is shown, which can characterize these parameters as a reliable marker of pyrogenesis.


Introduction
Since the end of the last glacial period, high-latitude landscapes of the Northern Hemisphere have accumulated vast organic carbon pools in the peat soils. In the last decades, scientists have often identified variations in climate and the extreme rate of wildfires in the widespread territory of boreal forests [1,2]. Modern climate change, together with anthropogenic influence, is dangerous because it occurs quickly and in combination with many other irreversible changes in the Earth's biosphere. However, in the geological time scale, climate change is not a new phenomenon. If we want to reveal how much climatic conditions change, it is worth paying attention to the natural climate variability during the Holocene (the last 11,800 cal. years BP). The most suitable is the Holocene climatic optimum, as it was the nearest long-term period (~10,000-5000 cal. years BP) when the average temperature exceeded present-day values [3].
Peatlands are important for the global carbon cycle, as they contain large reserves of carbon (C) in the biosphere [4]. Despite peatlands covering only 2-3% of the Earth's land surface, they store about 25% of the carbon stored in global soils [5]. Deep and dense layers of organic material accumulate in peat bog ecosystems. In addition, peat soils are unique archives that store information about climatic and environmental changes; therefore, they are suitable for studying the history of wildfires in the past [6]. The study of fire history A vertical growth rate of peat was constructed due to the model shown in Figure 1. The peat growth rate of the FH-1 site in the Boreal period of the Holocene (9350-7100 cal. years BP) was 0.5 mm/year. For a long period of 7100-2600 cal. years BP, peat growth was 0.3 mm/year. The maximum peat growth rate (up to 1 mm/year) was found for the period 2600-2160. In the periods of 2160-1400 and 1400-500 cal. years BP, the growth rate was 0.5 and 0.7 mm/year. However, in the last 500 years, the growth rate of peat has increased again to 1 mm/year. A vertical growth rate of peat was constructed due to the model shown in Figure 1. The peat growth rate of the FH-1 site in the Boreal period of the Holocene (9350-7100 cal. years BP) was 0.5 mm/year. For a long period of 7100-2600 cal. years BP, peat growth was 0.3 mm/year. The maximum peat growth rate (up to 1 mm/year) was found for the period 2600-2160. In the periods of 2160-1400 and 1400-500 cal. years BP, the growth rate was 0.5 and 0.7 mm/year. However, in the last 500 years, the growth rate of peat has increased again to 1 mm/year. Peat accumulation in FH-2 demonstrated the lowest rate of growth (0.3 mm/year) at an early stage of development (8380-5500 cal. years BP). In the period from 5500 to 4700 cal. years BP, the growth rate increased to 0.5 mm/year. In the period of ~4700-3100 cal. BP, the rate drops to 0.1 mm/year and then begins to increase sequentially. In the period of 3100-1100 cal. years BP, the average growth rate was characterized by values up to 0.21 mm/year. The maximum growth rate was detected in the period from 1100 cal. years BP to the present-0.6 mm/year.

The Botanical Composition of Peat
The analysis of the botanical composition of peat is presented in the diagram in Figure 2. The peat development in the studied area FH-1 included two stages. The initial stage of peat accumulation falls on a segment of 9300 cal. years BP. Presumably, the formation of the swamp began with the overgrowth of shallow ponds with hydrophilic plants. The next mesotrophic phase was represented by transitional sphagnum peat (180-200 cm). Oligotrophic plant associations were established at the border (transition) from the Subboreal to Subatlantic periods (~2500-3000 cal. BP) and formed a significant mass of peat (180-0 cm). Peat accumulation in FH-2 demonstrated the lowest rate of growth (0.3 mm/year) at an early stage of development (8380-5500 cal. years BP). In the period from 5500 to 4700 cal. years BP, the growth rate increased to 0.5 mm/year. In the period of~4700-3100 cal. BP, the rate drops to 0.1 mm/year and then begins to increase sequentially. In the period of 3100-1100 cal. years BP, the average growth rate was characterized by values up to 0.21 mm/year. The maximum growth rate was detected in the period from 1100 cal. years BP to the present-0.6 mm/year.

The Botanical Composition of Peat
The analysis of the botanical composition of peat is presented in the diagram in Figure 2. The peat development in the studied area FH-1 included two stages. The initial stage of peat accumulation falls on a segment of 9300 cal. years BP. Presumably, the formation of the swamp began with the overgrowth of shallow ponds with hydrophilic plants. The next mesotrophic phase was represented by transitional sphagnum peat (180-200 cm). Oligotrophic plant associations were established at the border (transition) from the Subboreal to Subatlantic periods (~2500-3000 cal. BP) and formed a significant mass of peat (180-0 cm).
The studied FH-2 peatland has also passed through several stages in its development. Its formation began with the mesotrophic phase around 8000 years ago at the transition of the Boreal and Atlantic periods of the Holocene with the waterlogged forested areas. This is evidenced by the presence of wood species residues (30%) in the peat bottom. This type of peatland is characterized by a small area of pine-shrub-sphagnum peat. The studied FH-2 peatland has also passed through several stages in its development. Its formation began with the mesotrophic phase around 8000 years ago at the transition of the Boreal and Atlantic periods of the Holocene with the waterlogged forested areas. This is evidenced by the presence of wood species residues (30%) in the peat bottom. This type of peatland is characterized by a small area of pine-shrub-sphagnum peat.

Charcoal Particles
The presence of charcoal particles in studied peat columns evidences the regular appearance of fires in the study area during the Holocene The peat bog FH-1 was exposed to 11 local fires ( Figure 3). The number of charcoal particles found during the Boreal period was 918 particles. It was found that five fires took place in the Atlantic period from 8000 to 5000 cal. years BP. During this period, 3571 pcs. charcoal were detected. The highest concentrations of charcoal particles in the peat FH-2 are observed in the period ~ 7000-6200 cal. years BP. (108-378 pcs.). The maximum coal content was detected at a depth of

Charcoal Particles
The presence of charcoal particles in studied peat columns evidences the regular appearance of fires in the study area during the Holocene The peat bog FH-1 was exposed to 11 local fires ( Figure 3). The number of charcoal particles found during the Boreal period was 918 particles. It was found that five fires took place in the Atlantic period from 8000 to 5000 cal. years BP. During this period, 3571 pcs. charcoal were detected. The highest concentrations of charcoal particles in the peat FH-2 are observed in the period~7000-6200 cal. years BP. (108-378 pcs.). The maximum coal content was detected at a depth of 327 cm (378 pcs.) and dates back to~6759 cal. years BP. In the Subboreal period, there was lower fire activity. During this period, 551 pcs. coal were found. Two local fires were recorded in the period~3500-3200 cal. years BP. (146 pcs.). High fire activity was noted in the Subatlantic period, in the time interval of~1000-2500 cal. years BP. During this  The peat bog FH-2 experienced 15 local fires during the Holocene. In total, there were 5592 pcs. charcoal found in the studied peat bog in the entire column. In the lower part of the profile, which belongs to the early Atlantic period of the Holocene (~9000-8000 cal. years BP), 598 charcoal particles were found in the studied peat soil. In addition, one local fire was detected, which is confirmed by the presence of layers with increased charcoal content (116-126 particles). Seven local fires were recorded during the Atlantic period (~8000-5000 cal. BP.) and the beginning of the Subboreal period (~5000-4000 cal. years BP). We detected 2101 pcs. charcoal in the peat column layer dated to this period. The maximum values ranged from 105 to 206 pcs. (~4275-4096 and ~6646-6364 cal. years BP.). A sharp increase in the rate of accumulation of charcoal particles at the end of the Holocene was noted, and seven local fires have been established since the end of the Subboreal period (~3000-500 cal. years BP.) and the Subatlantic period to the present (~1600-500 cal. years BP.). It was revealed that 2893 pcs. charcoal were accumulated during this period. The highest rate of coal accumulation was observed at a depth of 25 cm and dates back to the end of the Subatlantic period (~428 cal. years BP.). This layer contains the maximum content of charcoal particles (562 pcs.).

Carbon and Nitrogen Contents
The studied peat soils are characterized by high carbon content ( Table 2). The carbon content in soil FH-1 ranged from 456 to 563 g kg -1 . The nitrogen content in the peat FH-1 varied from 6.7 to 36 g kg -1 . In soil FH-2, the carbon content in peat horizons varied from 258 to 535 g kg -1 . The mineral horizon G (27 g kg -1 ) is characterized by the lowest carbon content. Soil FH-2 contains from 7.1 to 31 g kg -1 of nitrogen. The peat bog FH-2 experienced 15 local fires during the Holocene. In total, there were 5592 pcs. charcoal found in the studied peat bog in the entire column. In the lower part of the profile, which belongs to the early Atlantic period of the Holocene (~9000-8000 cal. years BP), 598 charcoal particles were found in the studied peat soil. In addition, one local fire was detected, which is confirmed by the presence of layers with increased charcoal content (116-126 particles). Seven local fires were recorded during the Atlantic period (~8000-5000 cal. BP.) and the beginning of the Subboreal period (~5000-4000 cal. years BP). We detected 2101 pcs. charcoal in the peat column layer dated to this period. The maximum values ranged from 105 to 206 pcs. (~4275-4096 and~6646-6364 cal. years BP.). A sharp increase in the rate of accumulation of charcoal particles at the end of the Holocene was noted, and seven local fires have been established since the end of the Subboreal period (~3000-500 cal. years BP.) and the Subatlantic period to the present (~1600-500 cal. years BP.). It was revealed that 2893 pcs. charcoal were accumulated during this period. The highest rate of coal accumulation was observed at a depth of 25 cm and dates back to the end of the Subatlantic period (~428 cal. years BP.). This layer contains the maximum content of charcoal particles (562 pcs.).

Carbon and Nitrogen Contents
The studied peat soils are characterized by high carbon content ( Table 2). The carbon content in soil FH-1 ranged from 456 to 563 g kg −1 . The nitrogen content in the peat FH-1 varied from 6.7 to 36 g kg −1 . In soil FH-2, the carbon content in peat horizons varied from 258 to 535 g kg −1 . The mineral horizon G (27 g kg −1 ) is characterized by the lowest carbon content. Soil FH-2 contains from 7.1 to 31 g kg −1 of nitrogen. As a result of the analysis of water-soluble organic compounds, the maximum concentrations were found for the middle part of the peat profile (40-140 cm). The distribution of WSOC and WSON content was uneven across the FH-1 peat profile. The FH-1 peat contained 1.36-2.15 mg g −1 of water-soluble organic carbon and 0.05-0.07 mg g −1 of watersoluble organic nitrogen. The maximum content was detected at depths of 60-80 cm. The content of water-soluble carbon WSOC in the peat soil FH-2 profile varied from 0.09 to 1.63 mg g −1 , and the concentration of nitrogen WSON varied from 0.03 to 0.17 mg g −1 . There was a gradual decrease in the WSOC content from a depth of 180 cm (1.21 mg g −1 ) to the mineral horizon (0.09 mg g −1 ). The maximum content of 1.63 mg g −1 was detected at a depth of 120-140 cm.
The analysis of stable isotopes showed that the values of δ 13 C varied in the peat FH-1 from −28.21‰ to −25.54‰ (−27.08 ± 0.43) and δ 15 N varied from −2.35 to +1.63‰. The isotopic composition of peat FH-2 is characterized by higher values δ 13 C from -28.04‰ to −24.07‰ and δ 15 N from −4.6 to +1.51‰. The calculation of carbon (C stock ) and nitrogen (N stock ) stocks allowed us to establish that peat FH-1 contains 149.7 kg m -2 of carbon and 4.7 kg m -2 of nitrogen. Peat FH-2 contains 184.1 kg m -2 of carbon and 7.8 kg m -2 of nitrogen.

CPMAS 13 C NMR Data
The results of 13 C-NMR spectroscopy on the studied peats are presented in Table 3. It was found that aliphatic fragments predominate in the studied peats (67-87%). The composition of organic matter in the peat FH-1 is not homogeneous. The composition of SOM in the upper part of peat (0-200 cm) is dominated by aliphatic fragments (81-87%). Their composition is dominated by alkyl-C fragments (15.4-25.9%) and groups associated with hemicellulose C Alk-O (35.9-51.2%) and cellulose C O-Alk-O (9.9-13.3%). Fragments of aromatic structures of Aryl-C constitute 13-19%. Carboxyl groups (Carboxyl-C) account for 3-4% of the composition of the peat organic matter shifts from a depth of 200 cm. There is a gradual decrease in the proportion of fragments associated with cellulose-like substances (6.0-8.4%) and hemicellulose (23.2-32.5%) in the soil profile. The proportion of methoxyl C CH3-O (7.1-9.5%) and carboxyl (7-10%) groups and alkyl-C fragments (21.6-32.7%) increases. The proportion of aromatic fragments increased in the lower part of the soil profile (up to 27%). The ratio of alkyl/O,N-alkyl demonstrates the downward increase in profile, reflecting the stronger degree of decomposition of peat organic matter. This increase is most pronounced from 200 cm depth.
The composition of SOM at a depth of 20-230 cm is quite stable. Aliphatic fragments in peat also predominate throughout the peat profile of FH-2 peat bog (up to 84%). The largest share falls on the upper part of the profile, with a predominance of O-alkyl groups (up to 61%). The proportion of aromatic fragments in the SOM composition is about 16-23% to a depth of 215 cm and increases to 33% in the pyrogenic horizon T pyr (235-240 cm). In the composition of aliphatic fragments, the main share is accounted for by C Alk-O (12.6-41.1%) and alkyl-C fragments (16.2-37.6%). Compounds associated with cellulose-like substances C O-Alk-O (5.3-12.2%) and C CH3-O methoxyl groups (6.8-11.5%) occupy smaller shares. The degree of decomposition of organic matter, as shown by the fact that the ratio of alkyl/O,N-alkyl gradually increases along the soil profile from 0.3 to 1.5.
The upper peat layer and the pyrogenic horizon T pyr of the FH-2 profile differ significantly from the main peat mass. The horizon T (0-20 cm) is represented by the least decomposed plant residues and is characterized by a minimum ratio of alkyl/O,N-alkyl (0.3). It is dominated by O-alkyl fragments (61%) associated with cellulose-like (12.2%) and hemicellulose substances (41.4%). There is also an increased content of aromatic components due to substituted C Ar-O,N (6.2%). Pyrogenic horizon Tpyr contains an elevated portion of aromatic fragments (33%), mainly due to fragments C Ar-H(C) (26.4%). The horizon is characterized by a minimum content of fragments of hemicellulose (12.6%), cellulose (5.3%) and methoxyl groups (6.8%). The organic matter of the pyrogenic horizon T pyr is maximally decomposed in the soil profile according to the ratio alkyl/O,N-alkyl (1.5). Table 3. Percentage distribution (%) of signal intensity between selected chemical shift regions (ppm) of 13 CPMAS NMR spectra of bulk soil.

Concentration and Distribution of PAHs
The content of PAHs in the studied peats varied from 0 to 2000 ng g −1 in soil FH-1 ( Figure 4, Appendix A Table A1). Peat FH-2 contained from 0 to 1043 ng g −1 PAHs. The predominant PAHs in the peat are 2-3-4 nuclear compounds (naphthalene, phenanthrene, fluoranthrene), 5-nuclear dibenzanthracene, and 6-nuclear indenopyrene. The nature of the PAH distribution in peat profile HF-1 is uneven. The composition of PAHs of the upper 200 cm is dominated by light compounds. The content of naphthalene was 42-64 ng g −1 , phenanthrene-32-59 ng g −1 . The remaining connections were much smaller. The total PAH content in the upper part of the profile was 99-173 ng g −1 . Peat HF-1 is characterized by a sharp increase in the content of all PAHs from a depth of 200 cm. The maximum content was found for naphthalene (up to 390 ng g −1 ), phenanthrene (up to 530 ng g −1 ), dibenzanthracene (up to 590 ng g −1 ), benzoperylene (up to 510 ng g −1 ) and idenopyrene (up to 2000 ng g −1 ). The total PAH content in peat HF-1 increases with depth (99-4370 ng g −1 ). Up to a depth of 330 cm, light compounds of 2-3-4-nuclear PAHs (99-1282 ng g −1 ) predominate in the composition of PAHs, deeper there is an increase in the content of 5-6-nuclear compounds (up to 3404 ng g −1 ). Plants 2022, 11, x FOR PEER REVIEW 9 of 21 The PAH content in peat FH-2 gradually increased down the profile to a depth of 200 cm. After that, there is a gradual decrease in the PAH content. The total content varies from 71 to 1974 ng g -1 : naphthalene (33-627 ng g -1 ), phenanthrene (22-99 ng g -1 ), dibenzo[a,h]anthracene (0-126 ng g -1 ) and idenopyrene (0-1043 ng g -1 ) were characterized by a predominant content of PAHs. It was found that the concentration of high-molecular 6-nuclear dibenzo[a,h]anthracene prevailed from a depth of 100 cm, except for the lower peat horizon (0 ng g -1 ), in which light PAHs predominated. An increase in the content of naphthalene (192-545 ng g -1 ) was also detected in the pyrogenic and over-pyrogenic layers.

Peat Formation in the Context of Pyrogen History
The paludification process and the accumulation of peat in our study territory began at 14,700-12,700 cal. years BP [14]. The research area is unique in terms of coverage of the processes of peat bog formation and peat accumulation. The vegetation of the studied territory is represented by communities of wetland pine forests. The average waterlogging of the study territory was about 50%; in some areas, it reached 70-75% [15]. Pine sphagnum and shrub-sphagnum forests replace each other, which are then replaced by pine shrub-sphagnum forested peat bogs [16,17]. It was revealed that the studied peatlands differ in the rate of vertical growth in certain periods of the Holocene. These differences reflect the development of the studied peatlands at different stages of the Holocene and their current state. At the same time, it is important to take into account that the rate of vertical growth in boreal ecosystems is due to a variety of natural factors.
Considering the botanical composition of the peat FH-1, it can be argued that at first, a sedge-sphagnum phytocenosis was formed with the dominance of Eriophorum on a cover of Sphagnum magellanicum with an admixture of S. angustifolium. Then a sphagnum community with the dominance of S. fuscum emerged. Fuscum peat contains the remains of hydrophilic Scheuchzeria palustris. During the upper stage of the development of the peat, which began to be established at the border Subboreal and Subatlantic periods and continues to the present, the forest areas surrounding the swamp were regularly exposed The PAH content in peat FH-2 gradually increased down the profile to a depth of 200 cm. After that, there is a gradual decrease in the PAH content. The total content varies from 71 to 1974 ng g −1 : naphthalene (33-627 ng g −1 ), phenanthrene (22-99 ng g −1 ), dibenzo[a,h]anthracene (0-126 ng g −1 ) and idenopyrene (0-1043 ng g −1 ) were characterized by a predominant content of PAHs. It was found that the concentration of high-molecular 6-nuclear dibenzo[a,h]anthracene prevailed from a depth of 100 cm, except for the lower peat horizon (0 ng g −1 ), in which light PAHs predominated. An increase in the content of naphthalene (192-545 ng g −1 ) was also detected in the pyrogenic and over-pyrogenic layers.

Peat Formation in the Context of Pyrogen History
The paludification process and the accumulation of peat in our study territory began at 14,700-12,700 cal. years BP [14]. The research area is unique in terms of coverage of the processes of peat bog formation and peat accumulation. The vegetation of the studied territory is represented by communities of wetland pine forests. The average waterlogging of the study territory was about 50%; in some areas, it reached 70-75% [15]. Pine sphagnum and shrub-sphagnum forests replace each other, which are then replaced by pine shrubsphagnum forested peat bogs [16,17]. It was revealed that the studied peatlands differ in the rate of vertical growth in certain periods of the Holocene. These differences reflect the development of the studied peatlands at different stages of the Holocene and their current state. At the same time, it is important to take into account that the rate of vertical growth in boreal ecosystems is due to a variety of natural factors.
Considering the botanical composition of the peat FH-1, it can be argued that at first, a sedge-sphagnum phytocenosis was formed with the dominance of Eriophorum on a cover of Sphagnum magellanicum with an admixture of S. angustifolium. Then a sphagnum community with the dominance of S. fuscum emerged. Fuscum peat contains the remains of hydrophilic Scheuchzeria palustris. During the upper stage of the development of the peat, which began to be established at the border Subboreal and Subatlantic periods and continues to the present, the forest areas surrounding the swamp were regularly exposed to fires [2,18]. This is indicated by the presence of charcoal particles diagnosed in all layers of peat. However, their small number and the absence of pronounced pyrogenic layers serve as indirect evidence that fires did not affect the watered part of the peat, and their intensity was not high.
At the initial stage (I) of development, peat FH-2 was affected by a fire (7535 ± 120 cal. BP). This is confirmed by the presence of burnt bark and wood residues in the botanical composition of peat (charcoals 10%). The remaining peat-forming plants have retained their anatomical structure. The peat is dominated by the remains of Betula sp., Pinus sibirica, P. sylvestris (5-15%), Equisetum, Carex caespitosa and Eriophorum (55% in total). At the next mesotrophic stage (II), covering the entire Atlantic and Subboreal period, paleophytocenoses are formed with the dominance of Eriophorum (50-90%) in the botanical composition of peat (230 to 90 cm). In addition, the samples contain the remains of woody plants (Betula, Pinus, Picea), sedges (Carex caespitosa, Carex lasiocarpa, C. limosa), Equisetum and Sphagnum mosses (Sphagnum obtusum, S. teres and S. magellanicum). Communities existed throughout the Atlantic and Subboreal periods of peat development (8050-2915 cal. BP). The presence of charcoal particles serves as evidence that the adjacent forest territories regularly experienced the burning. At the same time, the studied swamp was not affected by fires due to its high water content, as evidenced by the abundance of hydrophilic plants in the composition of phytocenoses. The absence of large fires is indirectly confirmed by the high degree of preservation of peat-forming plants. At the end of the Subboreal and the beginning of the Subatlantic period, the swamp passed to a new stage of development (III-IV). Successions occurred gradually without a sharp change in humidification conditions. Mesotrophic Carex lasiocarpa, C. caespitosa, Sphagnum obtusum, S. teres and others disappeared from the associations. They were replaced by oligotrophic Scheuchzeria palustris, Sphagnum fuscum and S. angustifolium. At the final stage, marsh shrubs (Ledum palustre, Andromeda polifolia, Chamaedaphne calyculata, etc.) began to appear. As a result, a shrubby-sphagnum swamp was formed, similar in composition and structure to the one currently existing. In the sub-Atlantic period (approximately the last 2500 years), the impact of fires and adjacent forests on the swamp was regular.
The data obtained by us are consistent with the results of L.V. Karpenko [18,19] for this region (11,800-10,200 cal. BP) and the central part of the West Siberian lowland (9000-11,500 cal. years BP) [20]. According to the authors [14], the differences can be explained by the methods of determining the age and the calibration parameters of the method. Nevertheless, in most studies, approximately the same age of peat accumulation in swamp ecosystems, for example, in Western Siberia (8500 cal. BP) [21] and the European North of Russia (11,010 cal. BP), in the Holocene [7,22] was observed. All researchers agree that fires affect the development of peat soils in the boreal zone. The pyrogenic features of peatlands reflect the modern development of bogs and their individual stages [2]. Knowledge about the rate of carbon accumulation and its changes over time is an important step in modern research. According to [23,24], peat and tundra ecosystems are able to shift from a carbon sink into its source due to increased soil respiration under warmer and drier conditions. Therefore, in order to better understand current and future ecosystem pyrogenic changes, it is necessary to study and understand their relationship with vegetation and climatic changes in the past [25].

Charcoal Particles as a Marker of Pyrogenesis
Charcoal particles are present throughout both analyzed peat soil profiles. This indicates that fires were constantly occurring during all periods of the Holocene when the studied peatlands were being developed. Probably, the presence of charcoals reflects the natural causes of fires associated with climatic conditions in Western Siberia [26]. Differences in the content of coal particles were revealed between the studied peat soils. Peat bog FH-2, in comparison with peat bog FH-1, is characterized by an increased frequency of local fires. This is probably due to the fact that the peat FH-2 has a smaller area and is surrounded by forest, which has been exposed to fires more often. The FH-1 peat bog developed in an open area and at a distance from the forest and was less often subjected to pyrogenic influence. In contrast, the peat bog FH-1 is characterized by a higher concentration of pyrogenic coals in general. This may indicate a higher intensity of fires or the ability to capture more charcoals. The increased content of pyrogenic coals was found in the deepest horizons of the studied peatlands. This indicates fires (or a series of fires) could trigger the formation of peats in this area.
The longest period of time with a high rate of coal accumulation occurred during the Atlantic period. According to the literature [27][28][29][30][31][32], in the Atlantic period, the temperature exceeded the current values by 3.4-5 • C, which probably affected the increase in the frequency of fires. Similar data are presented in the works of other authors [21,33]. It was also revealed that the summer values of solar radiation in the boreal zone reached their maximum values during the Holocene optimum (~9000-5000 cal. years BP), and, probably, this led to an increase in fires. The dynamics of fires could also be influenced by the taxonomic composition of plants in the Boreal biome in the Holocene [34]. Probably, plant communities resistant to fires were replaced by fire-hazardous ones. During the Subboreal and the beginning of the Subatlantic periods, the proportion of fire-hazardous plants was less, which led to a lower rate of accumulation of pyrogenic plant residues.
Thus, with an increase in summer air temperature, a decrease in climate moisture and a change in plant communities, we can expect an increase in the number and scale of fires in the middle taiga of Western Siberia. Both peatlands are characterized by an increase in the rate of charcoal accumulation at the end of the Holocene. Similar results were published in works on the European part of Russia [8,35], where the authors associate an increase in the rate of charcoal accumulation at the end of the Holocene with an increase in anthropogenic influence. A similar increase in the frequency of fires over the past two millennia has also been recorded in southern Siberia, such as Lake Baikal [25] and in the Altai Mountains [36,37].

Fire Impact on Soil Organic Matter
Knowledge about the rate of carbon accumulation and its changes over time is an important step in modern research of soil organic matter. According to [23,24,38], peat ecosystems from a pure carbon sink are able to turn into a source due to increased soil respiration under warmer and drier conditions. The high content of carbon, nitrogen and their water-soluble fraction is typical for peatlands, which is consistent with the data on Western Siberia, Finland and Canada (420-630 g kg −1 ) [39][40][41][42]. The results of this study are comparable with those previously published for the study region (460-570 g kg −1 ) [43]. The studied peats largely depend on peat-forming plants in terms of carbon and nitrogen content. Peat from vascular plants is characterized by a higher content of C and N than from sphagnum [43][44][45]. So, for peat FH-1, the upper part (0-180 cm) is characterized by lower concentrations of carbon and nitrogen than in deeper layers. There is a sharp transition to an increase in the carbon content in the deeper layers of peat (180-420 cm). This indicates differences in the botanical composition of the peat. Similar regularities were revealed for peat FH-2, where vascular plants (Scheuchzeria, sedge, etc.) increased in the depth of the profile in the composition of peat. The nitrogen concentration also increases with depth as the degree of decomposition increases. At the same time, the C/N ratio narrows. The C/N ratio in the peat FH-1 was 14-95. In the peat FH-2, the ratio varied from 20-79. The highest values of the C/N ratio are characteristic of the upper peat horizons. There is a noticeable transition from the sphagnum type of peat with a low nitrogen content to the type of peat with more significant nitrogen concentrations, which is probably due to the predominance of vascular plants in the deep layers of peat [19]. In general, the nitrogen content in the peat of the temperate and boreal zones varies between 4.0-29 g kg −1 [44]. There is no obvious effect of pyrogenesis on the carbon and nitrogen content. However, we assume that the impact of wildfire may affect carbon and nitrogen stocks.
In terms of peat reserves, Russia ranks first in the world with 3.691 million km 2 or 21% of the country's territory. The world reserves of carbon accumulated in swamps (in an area of 6.41 million km 2 ) amount to 329.0-528.0 Gt of carbon [46]. Significant differences in carbon reserves in the studied soils can be explained by the density of peat. The peat density of the HF-2 site is significantly higher than that of the FH-1 site. The lower bulk density of the soil FH-1 is explained by the higher rate of accumulation of young sphagnum mosses [47,48]. At the FH-1 site, the sphagnum gained not only height but also carbon faster compared to the FH-2 soil. The influence of fires on soil density should not be excluded. The works [45] show an increase in density after fires due to the accumulation of ash. The impact of fires may also be reflected in an increase in pyrogenic carbon (PyC) content [2]. More frequent fires at the FH-2 site are also likely to affect the increase in carbon stocks.
The molecular composition of peat and the degree of its humification is determined by the geographical location, vegetation and environmental conditions. The most common components in the organic matter of peat soils include cellulose, hemicellulose and lignin, which are mainly derived from plants [49]. Thus, peat is very heterogeneous and rich in organic carbon material.
The molecular composition of the studied peat is characterized by the predominance of aliphatic fragments. As is known, pyrogenesis is an important factor that determines the content of aromatic fragments in the composition of organic matter [50,51]. Horizons with clear signs of fire impact at depths of 219, 331, 399 and 413 cm are distinguished in the soil of FH-1. The contents of aromatic fragments in them are 20%, 27%, 26% and 27%, respectively. The maximum values of the degree of aromaticity (AR/AL) were found in the same horizons. Also, in these horizons, as the proportion of aromatic fragments increases, an increased content of coal particles is revealed. It can be assumed that fires led to an increase in the aromaticity of organic matter. The correlation coefficient between the content of coals and CAr-H(C) is significant (r = 0.63, p < 0.05) and for CAr-O,N (r = 0.56, p < 0.05) and for the general degree of aromaticity (fa; r = 0.61, p < 0.05). Data on molecular fragments in the organic composition of soil FH-2 showed an increase in aromatic fragments down the profile. The lower peat horizons stand out especially, one of which is pyrogenic. As it was noted earlier, the studied peat bog was exposed to fire at the beginning of its formation. This explains the significant increase in Aryl C and the overall aromaticity of organic matter. Fires and pyrogenic residues (charcoal particles) are probably the reason for the high aromaticity of peat horizons of the studied soils. The preservation or retarded decomposition of these residues was facilitated by the fact that they are least actively processed by destructors [9,52,53].

Isotopic Composition and PAH Content as a Criterion for Assessing the Impact of Fire
The stable isotope composition of soils can be considered a criterion for detecting environmental changes. This is especially true in peat soils, which are important archives containing information about the history of peat formation, climate change and past fires [6,54,55].
Peat consists of various plant residues that have undergone transformations as a result of the activity of microorganisms, and, first of all, the isotopic composition is influenced by peat-forming plants. In the botanical composition of peat FH-1, Sphagnum spp. dominated in the upper two meters. The values of δ 13 C in peat FH-1 ranged from −27.84 to −26.57 to a depth of 180 cm after δ 13 C began to increase (−26.66‰-−25.54‰). Deeper peat horizons (220-420 cm) are characterized by a decrease of δ 13 C values from -26.95‰to -28.03‰. An increase in the isotopic composition may also indicate a change in the hydrological regime towards a decrease in the humidity of peat horizons. An increase in δ 13 C with a depth in sphagnum peat in the upper 100 cm suggests wetter conditions during the transition to another stage of peat development [43]. In peat FH-2, the δ 13 C values increased to −25.31%-−24.07‰ up to a depth of 120 cm, followed by a decrease in deeper horizons (−25.14‰ . . . −28.04‰). In this peat at depths from 100 to 240 cm in the botanical composition, there was a change in the dominant peat-forming plants; namely, the Sphagnum disappeared, and Eriophorum appeared. In the mesotrophic stage of peat bog formation, fluffy phytocenoses prevailed, which deposited the bulk of peat (50-90%).
The δ 15 N values are characterized by similar patterns and are within the range of reported values for peat bogs from the Subarctic and Arctic Regions [56,57]. For the peat ecosystems considered, there is an obvious period of increase in total N and δ 15 N with depth. This period is characterized by the dominance of Eriophorum. This is consistent with the data of [58], where the predominance of Scheuchzeria palustris L. was noted (up to 50-70%). Similar results were obtained by the authors when studying different peatlands [59][60][61]. The characterization of the isotopic composition allows us to assess changes in moisture conditions, the type of peat vegetation and the degree of decomposition of organic matter. Nevertheless, we found a negative correlation coefficient between the δ 13 C and the charcoal content in the peat profile FH-1 (r= −0.47, p < 0.05). The negative correlation of δ 13 C values is shown with the aromatic structures of SOM (r= −0.53, p < 0.05). In contrast, δ 13 C has a positive relationship with aliphatic fragments of SOM (r = 0.53, p < 0.05).
It Is difficult to assess the direct impact of fires on δ 13 C and δ 15 N of peat. The revealed correlations should be discussed from the standpoint of the transformation of plant residues and the formation of peat by different types of peat-forming plants, which more clearly reflect the content of the isotopic composition. However, some authors suggest that fires lead to an increase in δ 13 C values [59]. Polyaromatic hydrocarbons are high-molecular structures, the sources of which can be both natural and man-made processes. The natural sources of PAHs in peat soils are usually attributed to the intake from rocks and oils [62], vegetation [59,63,64] and fires [2,65,66]. The content of PAHs in the studied peat soils is characterized by a tendency to increase the content of polyarenes deep into the profile. The predominance of light PAHs was revealed ( Figure 4, Table A1). Their share in peat FH-1 ranges from 1 to 52%. However, after 180 cm, there is an increase in the content of heavy PAHs, in particular, indeno[1,2,3-c,d]pyrene. The proportion of heavy PAHs after it varied from 1 to 58%. Similar patterns have been identified in peat FH-2. The prevalence of heavy PAHs after 140 cm was up to 65%.
The features of PAH distribution are the result of the influence of various factorsbotanical composition, redistribution of substances in peat, the degree of decomposition of SOM and pyrogenesis. Fires lead to an increase in the PAH content in soils [1,2]. Some authors note that the sources of PAHs are not only modern fires but also fires that have passed a long time ago (paleofires) [65,67]. In our work, we revealed an increase in the content of PAHs in horizons with pyrogenic features. High correlation values were obtained between the content of PAHs and the content of coal particles in peat horizons. This is especially clear for soil FH-1. The correlation coefficient between ∑PAH and the content of coal particles was r = 0.56, p < 0.05. Several individual PAHs with significant correlation coefficients with the charcoal contents were revealed: naphthalene (r = 0.53, p < 0.05), chrysene (r = 0.65, p < 0.05), dibenzo[a,h]anthracene (r = 0.70, p < 0.05). Presumably, coal particles formed as a result of fires (paleofires) can participate in the formation of PAHs. At the same time, significant concentrations of high-molecular PAHs are associated with their content in certain plant species. Deep peat horizons are characterized by an increase in herbaceous vegetation containing most of the lignin. Lignin is the main source of organic matter of no-pyrogenic origin [68]. In addition, the PAH content significantly correlates with aromatic fragments (markers of pyrogenesis) of organic matter (r = 0.77-0.78, p < 0.05) and negatively with aliphatic fragments (-0.78, p < 0.05). In [67], it was noted that pyrogenesis is a key factor in the formation of PAHs, and most of them are pyrogenic. The most frequently encountered PAH ratios in the literature were calculated to assess whether the isolated PAHs are pyrogenic [69][70][71][72]. The calculated coefficients are presented in Table A2.
The ratios ANT/178, ANT/ANT+PHE, FL/PYR and PYR+BaP/CHR+PHE showed borderline values of PAH origin. This indicates that the observed PAHs have a mixed origin from oil and from fires. In the studied peats, the ANT/178 and ANT/ANT+PHE ratios show that PAHs in the upper part of the soil profile (up to the "turning point") do not have pyrogenic origins. From the depth (200 cm for peat FH-1 and 120 cm for peat FH-2), boundary values (>0.1) were detected. This may indicate the involvement of pyrogenesis in their formation. The ratio PYR+BaP/CHR+PHE also characterizes the PAHs of the upper peat horizons as non-pyrogenic. However, in the lower part, the pyrogenic origin of PAHs is confirmed by this coefficient. The FL/PYR ratio has opposite indicators. PAHs in the upper part of the studied peats are pyrogenic, and in the lower horizons, they have a mixed origin. It can be assumed that despite the boundary values of the coefficients, PAHs can be formed as a result of fires. Probably, these ratios should be applied with caution because it is impossible to reliably assert the origin of PAHs.
The ratios FLA/PYR, FLA/(FLA+PYR) and IcdP/(IcdP+BghiP) also serve to assess the formation of PAHs [70,73]. The obtained coefficients characterize the studied PAHs in the majority as pyrogenic. The studied peatlands are characterized by similar patterns of distribution of these diagnostic coefficients. According to the ratios FLA/PYR (<1.4) and IcdP/(IcdP+BghiP; >0.5), the predominance of pyrogenic PAHs in the profile of peats was revealed. However, mainly from the middle of the soil profile (approximately from a depth of 100 cm). According to the literature [74][75][76], these coefficients reflect the burning of grasses, wood and coal, shrubs and other plant residues. The ratio FLA/(FLA+PYR) more accurately conveys the pyrogenic origin of PAHs. The obtained values (>0.5) can be considered evidence as the main source of PAHs of vegetation, coal and wood. Thus, the obtained features of the vertical distribution of PAH concentrations in peat make it possible to more accurately identify the role of paleofires in the formation of the studied peatlands.

Area Description and Soil Sampling
The studies were carried out in the boreal zone of the north-west of the eastern part of Russia. The research area belongs to the Sym-Dubsky Middle Taiga cedar-pine region and is located on the left bank of the Malaya Khoyba River (a tributary of the Yenisei) on the territory of the Sredneenisei station of the V.N. Sukachev Institute of Forest SD RAS (ZOTTO Observatory). The study area is characterized by a sharply continental climate. According to the data of the Bor weather station (the closest to the objects of the study), the average annual air temperature is −3.5 • C, and the average precipitation is 594 mm/g.
The objects of the study were forested peatlands which formed Fibric Histosol (FH). Site FH-1 (60 • 75 61.11 N, 89 • 40 45 E) is a vast swamp with an area of 160 hectares. Site FH-2 (60 • 81 31.39 N, 89 • 32 99.44 E) was a sphagnum partially forested swamp with an area of 0.7 hectares. Two soil sections were laid for sampling from peat horizons. In addition, peat columns were taken using an Eijkelkamp soil drill with a sampler diameter of 5 cm and a length of 50 cm.

Botanical Composition
The analysis of the botanical composition and the degree of decomposition of peat was carried out in the Laboratory of Phytocenology and Forest Resource Studies of the Sukachev Institute of Forest, Siberian Branch of the Russian Academy of Sciences (Krasnoyarsk, Russia) according to GOST 28245-89 using a microscope "Leitz Wetzlar" with magnification × 20 and × 40. Stratigraphic diagrams of peat composition were constructed using the software «Korpi» [77]. Atlases were used to identify plant residues in peat [78,79]. Sampling for the analysis of the botanical composition of peat FH-1 was carried out only to a depth of 210 cm.

Charcoal Analysis and 14 C Dating
The calculation of macroscopic particles of charcoal in peat soils has been out, according to [80]. From the raw column sample, peat was selected with a volume of 1 cm 3 every 2 cm. Samples were soaked with a 100 cm 3 5% NaOCl for 24 h at room temperature to be able to identify black char particles. Then samples for charcoal analysis were gently washed with distilled water through a sieve (0.125 mm), diluted in water in Petri-dishes, and researched for macroscopic charcoal under 40 x stereoscopes. The char particle analysis results were processed in the CharAnalysis program in the R software environment.
Radiocarbon dating of peat samples was carried out at the equipment provided by Shared Research Facilities of the Tomsk Scientific Center SB AS (Tomsk, Russia) by liquid scintillation method using a spectrometer-radiometer Quantulus 1220 (PerkinElmer Life Sciences/Wallac Oy, Turku, Finland). Calibration of the radiocarbon age to calendar age was performed using the CALIB REV-7.10 programs.

Carbon and Nitrogen Analysis
The total C tot and N tot were determined by dry combustion on an EA-1100 analyzer (Carlo Erba, Milano, Italy). Carbon and nitrogen pools were calculated using the formula according to [81]: C stock (kg m -2 ) = 0.1 × Ctot × BD × h, where h is the thickness of the soil layer (cm), C tot is the organic carbon concentration of the soil layer (g kg −1 ), BD is the bulk density of the soil layer (kg m -3 ). A 47.7 cm 3 drill was used to determine the bulk density of soil horizons.
Water-soluble organic carbon (WSOC) and nitrogen (WSON) were extracted with deionized water (ELGA Lab Water, Lane End, UK) at room temperature at a ratio of 1:50 (soil: water) for mineral horizons and 1:100 for organic and pyrogenic horizons in BIOFIL test tubes with quartz filters (MN, Germany, pore size 0.4 µm). Total carbon (TC) and nitrogen (TN) were assessed using the TOC-VCPN analyzer (Shimadzu, Kyoto, Japan) with a TNM-1 module. The obtained results were recalculated on a dry soil mass basis.

13C-NMR Spectroscopy
The composition of the organic matter was determined by solid-state CP-MAS 13 C-NMR. Solid-state 13 C-NMR spectra of the two upper soil horizons were recorded on a 100.53 MHz Bruker Avance III 400WB (Bruker, Ettlingen, Germany) with a rotation frequency of 12.5 kHz, a contact time of 5 ms, and a 2-s recycle delay at the resource center of the research park "Magnetic Resonance Research Methods" of Saint-Petersburg State University, Russia. Before analysis, the samples were treated with 10% hydrofluoric acid to remove paramagnetic iron, according to [82]. Chemical shifts of fractions were determined relative to a tetramethylsilane shift (0 ppm). The contribution of main carbon forms to the total spectral intensity was determined by the integration of the corresponding chemical shift regions according to [83,84]. The total content of aryl C (AR) was calculated as the sum of the signals at 110-165 ppm fields. Signals from alkyl C (AL) were recorded in the 0-110 ppm range. The degree of aromaticity (fa) was determined as the 100% * proportion of total content of aryl C components (110-145 and 145-165 ppm) on total C (excluding the contribution in the range of 165-220 ppm). The hydrophobicity index (Hb) was calculated by the sum (%) of the signals 0-45 and 110-145 ppm.

Statistics
Correlation analysis was used to determine the relationship between the obtained data. Correlation coefficients (r-Pearson) were calculated using the STATISTICA 10.0 (Stat. Soft Inc, Tusla, OK, USA); differences were considered significant at the significance level p < 0.05.

Conclusions
In our study, a reconstruction of the fire activity of the Western Siberian forest-peatland complexes of the Holocene time period was presented. The studied peatlands have a close history of past fires in the Holocene. It has been established that for 9000 years, the studied peats have been repeatedly confirmed by fires (from 11 to 14 local fires). Two peaks with a maximum content of charcoal particles were characteristic of the Atlantic (~9100-5800 cal. years B.P.) and Subatlantic (~3100 cal. years B.P. to the present) periods of the Holocene. It is shown that the Subboreal period (~5000-2500 cal. years BP) is characterized by a significantly smaller number of fires. There are certainly many reasons for such dynamics of fires (Holocene climate optimum, botanical composition of peat and modern anthropogenic impact). The Holocene climate optimum was a warm period that occurred between about 9000 and 5000 years ago, which contributed to the higher frequency and intensity of fires in the past. A high content of coal particles is characteristic of the lower horizons. It can be assumed that fires in them were a factor that led to the initial waterlogging of the territories. Anthropogenic factors are probably related to the increase in the frequency of fires recorded in recent times (the last 500-1000 years). The effect of fires on soil organic matter is manifested in an increase in carbon reserves in soils and aromatic components in the composition of the SOM. It is suggested that pyrogenesis may affect the increase δ 13 C. There is a significant correlation between the content of PAHs with charcoal particles and with the aromatic structures of SOM in pyrogenic peat layers. Probably, the isolated PAHs are of pyrogenic origin. This is also confirmed by the calculated coefficients of the origin of PAHs. Thus, the content of charcoal particles and the increased content of aromatic fragments of SOM and PAHs in peat horizons can serve as reliable markers of pyrogenesis in peat soils.