Holocene Soil Evolution in South Siberia Based on Phytolith Records and Genetic Soil Analysis ( Russia )

Dark, coniferous hemiboreal forests in the south of West Siberia are located in the Holocene forest-steppe ecotone, where natural environmental conditions have been quite dynamic. This dynamic environment resulted in the contrasting evolution of regional soil cover and the development of unique soil profiles with the second humus horizon. The second humus horizon is assumed to be a relic from the dark-humus soil formation stage in the mid-Holocene. This article draws conclusions about changes in regional environmental conditions by analysing data from a geochemically interrelated coevolutionary soil series, obtained by using a combination of conventional soil studies, phytolith analyses, and accelerator mass spectrometry (AMS) dating of phytolith-occluded carbon (PhytOC) and humic acids. The results showed that, in general, phytocenoses changed from mire-meadow vegetation towards forest vegetation via the meadow stage. However, these stages had different durations, depending on the soil catenary position. The topographical divergence of soil phytolith profiles reflects the relief effect on the development of specific soil type combinations, accounting for the major elements of the regional mid-Holocene soil cover. The leading elementary soil-forming processes were humus accumulation and hydrogenic accumulation of calcium carbonates. In the hilltop site of Endocalcic Stagnic Albic Luvisols, the evolutionary changes were shown by the shift from the meadow phytocenosis (Calcic Stagnic Chernozem) to the forest phytocenosis. In the midslope site, the environment was more humid from the start, favouring a phytocenosis with features of the meadow-mire type. The shift from the meadow-mire environment (with Spodic Chernic Gleysols) to the forest type environment with leading profile-forming processes, acid hydrolysis and lessivage, was gradual, occurring via the meadow stage with Calcic Stagnic Chernozem. In the toeslope site (Calcic Stagnic Greyzemic Epidystric Umbrisols), the meadow-mire stage (with Spodic Chernic Gleysols) was succeeded by the forest stage of soil formation. The AMS-dating of PhytOC estimated that the dark-humus stage of soil formation began 6.5–5.7 years calBC. Despite the observed slight translocation of phytoliths down soil profiles and phytolith solubilisation, phytolith analysis can be used to reconstruct shifts in the soil formation environment for surface Holocene palaeosols.


Introduction
Dark, coniferous hemiboreal forests (Milio effusi-Abietetea sibiricae) [1,2] are common in the southern taiga subzone within the West Siberian Plain, where they grow on automorphic, texturally differentiated, and semihydromorphic organic matter accumulating soils that often have residual humus horizons [3][4][5].Such forests are spread throughout the Holocene forest-steppe ecotone, and dynamic changes in the geographic environment of the ecotone in the past have resulted in the complex and contrasting evolution of its soils and soil cover.In the mid-Holocene, the soils of the forest-steppe ecotone underwent the dark-humus stage of soil formation, which was later succeeded by the active textural differentiation of automorphic soils and the hydrogenic transformation of soils on the margins of the increasing mires.Due to all these processes, mineral soils of the region display polygenetic features, recorded in the complex structure of their organic matter profiles.These include high humus content, differences in humus quality between the contemporary and the second humus horizons, occurrence of molecasts, residual solonetz, carbonates, and solodic features.
Diversity of the soil morphology, resulting from combined modern and relic features, has led to much discussion on the origin of the soil.The most heated discussions concern the second humus horizon and its residual high humus content.Presently, there exist two conflicting hypotheses: (1) The actualistic hypothesis, which considers the complexity of the organic matter profile to be the result of the illuviation of humic acids at the border between eluvial and illuvial layers [5,6]; and (2) the historico-evolutionary hypothesis, which considers the second humus horizon and its residual high humus content to be relic features of the earlier dark-humus soil formation stage of the Holocene ( [3,4,[7][8][9][10][11][12], and other).The latter hypothesis is currently the most widely accepted and has been confirmed by a series of humic acid radiocarbon dating that was obtained from the second humus horizon with modal dates of 4000-5000 year with the oldest being 6000-7000 year [13,14].
Soils with complex organic matter profiles may occupy different meso-and microrelief positions, each of which influences the morphology of the second humus horizon in a number of ways, i.e., their depth, thickness, colour, and patchiness.In addition, the topogeny in the mid-Holocene also determined the specific soil formation environment, a fact which was noted already in the second half of the last century [3,11]; however, further consideration of this aspect of the second humus horizon genesis has not obtained proper scientific consideration.
Thus, the diversity of soil polygenetic profiles and the differing views on their evolution necessitate additional studies with new techniques that will either confirm or contradict these hypotheses.Phytolith analysis is one such technique that has been increasingly employed in soil research for morphological and genetic analyses [15][16][17][18][19][20][21][22][23][24][25][26].The relatively weak migration capability of phytoliths in a landscape [27][28][29], and allows correlational study between phytolith assemblages and soil properties that reflect specific soil formation environments.Additionally, the combined use of phytolith and soil analyses provides the opportunity to study soil development due to catena positions.Employing only early methods to study complex organic matter profiles, i.e., humus structure and composition, soil micromorphology, and carbonate profiles, does not always provide unequivocal answers about the environmental conditions that favoured the development of their complex profiles, as second humus horizons display different degrees of degradation.
Therefore, to study the evolutionary peculiarities of soils with complex organic matter profiles within a geochemically dependent soil series of the Vasyugan sloping plain in the south of West Siberia, we used a set of morphological, genetic, physical, and chemical methods, as well as phytolith analysis.
In addition, we aimed to: (1) Assess the possibility of using phytoliths to determine the time of the organic matter accumulation stage of soil formation in comparison with humic acids; (2) estimate the temporal sustainability of soil phytolith profiles compared to the profile distribution of other characteristics, such as humus content, carbonate content, etc.; and (3) reconstruct soil cover during the regional Holocene optimum using polygenetic soils situated on the same catena.

Soil Sampling
A series of pits in texturally differentiated and organic matter accumulating soils were dug along a microcatena from the hilltop down to the toeslope sites on the sloping hilly transition between the watershed and the Iksa riverbed (58 • 44 04" N, 68 • 50 52" E, Bakchar district, Tomsk region, Russia).The microcatena was 80 m long with a 0.2% incline (Figure 1). between the watershed and the Iksa riverbed (58°44′04″ N, 68°50′52″ E, Bakchar district, Tomsk region, Russia).The microcatena was 80 m long with a 0.2% incline (Figure 1).Pits displaying the eluvial part and the upper portion of the illuvial layer (below the lower border of the second humus horizon, with no explicit features of windfall phenomena) were dug on both catena positions.
Soil for physical and chemical analyses was sampled every 5-10 cm by continuous coring through soil genetic horizons, while phytolith analysis samples were taken from the lower 1-2 cm of each 5 cm soil horizons [31].
For radiocarbon dating, soil samples were collected from a 1 cm thick layer at the upper and lower borders of E, Eh, and A' in each soil pit.

Hillltop and Midslope Catena Sites
In the hilltop and midslope catena sites, soils are represented by Endocalcic Stagniс Albic Luvisols (Siltic, Abruptic, Cutanic, Epiloamic) within the second humus horizon [32].The sites are covered by fir and pine (Abies sibirica L. and Pinus sibirica L.) forest with an understorey of Sorbus sibirica L. and Padus avium Mill., and a ground cover of herbs and green mosses.Shrubs are represented by Spirea media L. and Rubus idaeus L. Projective cover of the shrubs and herbs story was ca.60%, being dominated by Pteridium aquilinum L., Aegopodium podagraria L., and Vaccinium myrtillus L. The mosses were dominated by the green mosses, Hylocomium splendens (Hedw.)B.S.G. and Pleurozium schreberi (Brid.)Mitt.Pits displaying the eluvial part and the upper portion of the illuvial layer (below the lower border of the second humus horizon, with no explicit features of windfall phenomena) were dug on both catena positions.

Toeslope Catena Site
Soil for physical and chemical analyses was sampled every 5-10 cm by continuous coring through soil genetic horizons, while phytolith analysis samples were taken from the lower 1-2 cm of each 5 cm soil horizons [31].
For radiocarbon dating, soil samples were collected from a 1 cm thick layer at the upper and lower borders of E, Eh, and A' in each soil pit.

Hillltop and Midslope Catena Sites
In the hilltop and midslope catena sites, soils are represented by Endocalcic Stagnic Albic Luvisols (Siltic, Abruptic, Cutanic, Epiloamic) within the second humus horizon [32].The sites are covered by fir and pine (Abies sibirica L. and Pinus sibirica L.) forest with an understorey of Sorbus sibirica L. and Padus avium Mill., and a ground cover of herbs and green mosses.Shrubs are represented by Spirea media L. and Rubus idaeus L. Projective cover of the shrubs and herbs story was ca.60%, being dominated by Pteridium aquilinum L., Aegopodium podagraria L., and Vaccinium myrtillus L. The mosses were dominated by the green mosses, Hylocomium splendens (Hedw.)B.S.G. and Pleurozium schreberi (Brid.)Mitt.

Toeslope Catena Site
In the toeslope site, the soil is represented by Calcic Stagnic Greyzemic Epidystric Umbrisols (Siltic) under pine (Pinus sibirica L.) and birch (Betula pubescens Ehrh.) forest with some aspen (Populus tremula L.) and no understory.The shrubs are well developed, consisting of Rubus idaeus L. and Rosa acicularis Lindl.Projective cover of the herbaceous storey was ca.50%, with Calamagrostis arundinaceae L. Roth., Filipendula ulmaria L. Maxim., and Equisetum sylvaticum L. dominating.The story of lichen and mosses is poorly developed, being represented by Pleurozium schreberi (Brid.)Mitt.

Morphologic and Genetic Analyses
During the field work, the morphology of the studied soils was described following the Guidelines for soil description [33], and soils were identified according to the World Reference Base (WRB) classification [32].

Physical and Chemical Analyses
Particle-size analysis of soil samples was determined by the pipette method after preliminary treatment of soil with a 4% Na 4 P 2 O 7 solution [34] (The fraction limits were established following Kachinskiy [34].These limits are rather similar to those of WRB [32], the only discrepancy in the limits of the clay fraction, i.e., 0.001 mm according to Kachinskiy and 0.002 mm according to the WRB, does not seem to be very significant).
The studied soils had heavy or light loamy particle size composition with sharply differentiated textural profiles and a slight difference in silt differentiation coefficients.

Accelerator Mass Spectrometry 14 C Dating of Humic Acids (II) and Phytolith-Occluded Carbon
The calcium-free humic acid fraction II (HAII) was prepared following the protocol described by O. Chichagova [38].
Phytoliths for 14 C-dating were prepared according to the slightly modified method of X. Zuo et al. [39].The detailed steps were as follows: (a) Dry soil was crushed and sieved at 500 µm; (b) the sample was deflocculated with 4% Na 5 P 3 O 10 , and then washed three to four times with distilled water; (c) organic matter was first oxidised by 250 mL of H 2 O 2 (30%) for 12 h at room temperature and then heated in a water bath until the reaction stopped; d) carbonates were eliminated using 200 mL of HCl (10%) with heating for 30 min; (d) the >250 µm fraction was separated by wet sieving, and the remaining sample was disaggregated from the organic material and clay by ultrasonic treatment for 20 min; (f) clays (<5 µm) were removed by gravity sedimentation until the sample was clear (2-3 month); and (g) any remnant organic material was removed by digestion using "Schulzes Solution".50 mL KClO 3 was added to 250 mL of HNO 3 to make this highly corrosive cocktail.About 30 mL of the solution was added to the test tubes containing the sample and digestion carried out in a water bath (80 • C) in a fume cupboard.Distilled water was added to the sample, and this was centrifuged for 5 min at 5000 rpm.This procedure was repeated three times to clean the sample of any remaining solution.(h) Phytoliths were extracted three times by 200 mL of heavy liquid (CdI + KI) with a specific density of 2.35 g/cm 3 and then washed three times with distilled water; (i) extracted phytoliths were further sieved at the track membrane with pore of 5 µm to remove clay.Then, the recovered part of the remains in the sieve were treated by 20 mL of H 2 O 2 (30%) in the tube for 20 min.(j) Finally, the recovered phytoliths were dried at 60 • C for 24 h prior to testing.
The purity of the phytoliths was checked by optical microscope (X400) and SEM-EDS analysis.In this study, the extracted phytoliths were analysed using an SEM (JEOL 6610 LV, Tokyo, Japan) in association with an EDS (Energy Dispersive Spectroscopy) system (INCAWave, SEM, Tokyo, Japan).
Graphite was obtained from phytolith concentrate in which the 14 C isotope content was measured using the CAIS 0.5 MeV accelerator of the SB RAS Core Facility "Cenozoic Geochronology" (Novosibirsk, Russia) and of the Center for Applied Isotope Studies, University of Georgia (Athens, GA, USA).The radiocarbon dates were calibrated using the OxCal v4.3 software [40] and the IntCal 13 calibration curve.
Carbon content of humic acid and phytolith preparations was determined using the elemental analyser, Euro EA 3028 HT, Eurovector (Milan, Italy).

Phytolith Analysis: Procedure and Interpretation
Samples for phytolith analysis were dried and subjected to standard procedures [41].After treatment with a hot 30% H 2 O 2 solution and 10% HCl, samples (approximately 50 g) were separated from sand and clay, and subjected to flotation in heavy liquid (cadmium iodide and potassium iodide solution with a gravity of approximately 2.3).After centrifugation (1500 r/min for 20 min), the floating phytoliths were collected in a tube, washed several times with distilled water, and studied under an optical microscope with 300× magnification [41].Phytoliths were counted for the cover glass (24 mm × 24 mm) area.
Phytolith analyses data were interpreted according to the methodology proposed by A.A. Golyeva [18,31,42,43].The authors identified a special type of soil profile-the phytolith profile-which is on par with other profiles (humus, carbonate, texture, etc.).
A phytolith profile develops during the entire soil life span, recording changes in the environment, soil-forming factors, or anthropogenic impacts in its qualitative and quantitative characteristics.Changes in soil-forming conditions and in phytocoenosis properties lead to the gradual replacement of one phytolith profile with a new one, but the general stratigraphic regularity of the phytolith persists, i.e., the deeper samples are older than the ones closer to the current soil surface.
To describe phytolith forms, we used the Code for Phytolith Nomenclature 1.0 [44], and then grouped the forms into diagnostic assemblages according to A. Golyeva [18,31,42,43]: The universal phytolith forms do not allow identification of plants below the class level of Dicotyledoneae; several ecological assemblages specific for certain phytocoenoses, such as taiga, meadow, steppe, or dry steppe; and phytolith forms that are specific for some families, e.g., Cyperaceae, Pinaceae, or species, e.g., Phragmites spp., called "signal forms" by some researchers [31].

Preservation of Some Features of the Dark-Humus Stage of Soil Formation
Reconstructing soil-forming conditions and soil evolution on the basis of the properties of modern polygenetic soils presents certain methodological difficulties.Soil as an "open-closed" system is constantly exchanging energy and matter with its environment, thus compromising its relic features.Those features that are preserved represent fragments of the indigenous profile that has been transformed.Analysis of such transformation is indispensable for the studies of polygenetic soils.
The catena soils in our study are heterogenic and quite synchronic in origin, which implies their co-evolution.This circumstance allows for the assessment of the sustainability of soil features resulting from the dark-humus accumulation stage of soil formation due to the relief distribution of energy and matter flows.
Increasing textural differentiation of soils from the toeslope to the hilltop sites is the main factor limiting changes in the second humus horizons, such as decreased thickness, colouration, and spread, intrapedic mass structure, and certain other properties.

Soil Morphology and Genesis
Field examination of soil profiles and morphology revealed different thicknesses of albic material and textural differentiation of soil profiles in all catena sites (Table 1).However, the morphological traits of the processes seemed to be expressed differently depending on the relief (Figure 1), as textural differentiation is most pronounced on the hilltop and less pronounced on the toeslope.There were also signs of the gradual degradation of the upper part of the second humus horizon due to the albic material forming, which is displayed by introduction of fragments of whitish material from the E horizon into the dark-coloured material of the Eh and A' horizons.
On the hilltop, the second humus horizon (Eh) is the most degraded and could be seen as the dark-grey patches on a 10 cm thick light-grey background with a foliar structure.In the midslope and toeslope sites, the second humus horizon (A') was much better expressed morphologically, as a black-coloured continuous 25 cm thick horizon with a grainy microstructure.
The morphology of the second humus horizon on the midslope, i.e., its similar colouration and structure, allows the conclusion that humus accumulation was the leading elementary soil process, resulting in the development of the specific soil profile.A large (40-60 cm) portion of the humus horizon, currently registered in the soil profile as the humus-eluvial layer, suggests a quasi-equilibrium state between already well-developed, mature soils and their soil-forming factors at the start of textural differentiation.
Thus, the soils of the studied catena have a polygenetic structure and can be considered the relics of the main regional soil types of the optimum Holocene soil cover.

Soil Physical and Chemical Properties
The Endocalcic Stagnic Albic Luvisols had an accumulative distribution of H + cations, and a progressively eluvial distribution of Ca 2+ (Table 2).
The contribution of exchangeable Ca 2+ into the total sum of exchangeable cations in the upper soil profile was lower in the hilltop site compared to the midslope site.Exchangeable Mg 2+ was shown to have an eluvial-illuvial profile distribution.Sodium cations were also found to contribute to the cation-exchange soil complex (0.5-0.6 Cmol kg −1 ).The studied Endocalcic Stagnic Albic Luvisols were shown to have contrasting acidity, as the upper 70-80 cm of the soil profile had an acidic reaction (pH H 2 O 4.1-6.7),which became more neutral, and even alkaline, with increasing depth (pH H 2 O 7.0-8.1).Below 70-80 cm, the upper border of the high carbonate content, which reached 5%, was found.
Calcic Stagnic Greyzemic Epidystric Umbrisols with the second humus horizon occupied the toeslope site of the catena.This soil exists in both acidic and alkali conditions.The upper profile, i.e., the Aq horizon, is under the influence of acidic solutions reaching it from the acidic forest litter (pH H 2 O 4.2-6.3),whereas in the lower profile, i.e., the A' and Bgk horizons, the soil reaction is neutral or even alkaline (pH H 2 O 7.0-8.2).
Exchangeable Ca +2 was found to be 16-22 Cmol kg −1 , while the H + contribution was small, but, nevertheless, could be traced down to the B horizon.Exchangeable Na + was found to range within 0.7-0.8Cmol kg −1 .The result of particle-size analysis showed that the soil had silty clay loam/silty clay composition (Table 2).Soil organic C showed an accumulative profile distribution pattern, with a maximum (3%) in the A horizon and a gradual decrease below the A' horizon upper limit (1%) and further down to the B horizon (0.3%).The B horizon contained few carbonates (>0.4%) that were detected only below 70 cm.

Physical and Chemical Properties Preserved from the Dark-Humus Stage of Soil Formation
Soil organic C properties strongly reflected the degradation of the second humus horizon, as the SOC profile distribution was found to change depending on clay illuviation and thickness of the albic material (Table 2).Maximal and minimal second humus horizon degradation was found in the hilltop and toeslope sites studied, respectively; the SOC profile distribution had an accumulating pattern in the hilltop and the toeslope positions, reaching, in the A horizon, 1.9% and 3.0%, respectively.In the midslope site, the soil SOC profile showed a regressively accumulating pattern with relative maxima in the A and A' horizons, which suggests this soil is an intermediary variant of the degradation of the dark-humus soil formation stage.
The presence of carbonate concretia in the Bt and B horizons, in combination with a high degree of CEC saturation (Dystrict), suggests that carbonates in the soil represent a relic property, i.e., a residual from the previous soil formation stage with a predominant transpiration water regime (Tables 1 and 2).Following humus profiles, soil carbonate profiles experienced significant changes; however, despite the latter, these profiles represent solid evidence of additional water inflow of the soils in the past.The A horizon phytolith assemblage consisted of the phytoliths of dicotyledonous herbs, forest, and meadow grasses (the latter predominating), as well as of phytoliths of mosses and coniferous trees (Figures 2 and 3A).The charcoal phytoliths contributed up to 7%.The A horizon phytolith assemblage consisted of the phytoliths of dicotyledonous herbs, forest, and meadow grasses (the latter predominating), as well as of phytoliths of mosses and coniferous trees (Figures 2 and 3A).The charcoal phytoliths contributed up to 7%.The E horizon phytolith assemblage had the same phytolith composition, but with increased concentrations of some groups.The horizon upper border was found to have maximal (18%) content of charcoal phytoliths, and the ratio of forest and meadow gramineous phytoliths shifts more four-to six-fold) towards the latter.
The second humus horizon, present in the form of dark-grey patches according to the morphological description, had increased phytolith concentration, with the ratio of forest and meadow gramineous phytoliths shifted towards the former ones.Phytoliths derived from the steppe grasses were found in the lower part of the horizon (1%).The E horizon phytolith assemblage had the same phytolith composition, but with increased concentrations of some groups.The horizon upper border was found to have maximal (18%) content of charcoal phytoliths, and the ratio of forest and meadow gramineous phytoliths shifts more four-to six-fold) towards the latter.
The second humus horizon, present in the form of dark-grey patches according to the morphological description, had increased phytolith concentration, with the ratio of forest and meadow gramineous phytoliths shifted towards the former ones.Phytoliths derived from the steppe grasses were found in the lower part of the horizon (1%).
According to phytolith analysis of the soil in the eluvial position, the start of the soil formation and its most active stage were related to the meadow stage, as evidenced by the dominance of meadow grasses' phytoliths over the forest grasses' ones.The tree stands of that stage had coniferous species during the entire soil formation process, with a certain influence of fires on phytocoenosis development and structure.
Presumably, the closest analogue of the second humus horizon in the eluvial position soil can be seen in (semi)hydromorphic soils (Calcic Stagnic Chernozem) in the forest-steppe zone.
3.3.2.Endocalcic Stagnic Albic Luvisols (Siltic, Abruptic, Cutanic, Epiloamic) (Transit Position) The soil profile phytolith assemblage in the midslope site was found to be rather similar to the one in the hilltop site; that is, there were mostly forest and meadow gramineous phytoliths, the latter exceeding 3.5-3.8times the former, thus dominating in the humus-eluvial part of the soil profile (Figure 3B).
However, large parallelepiped-like morphotypes of Phragmites spp.phytoliths are found in small quantities in the E, A', and the upper A'/Bt horizons (Figures 2f and 3B), proving more periodically humid environmental conditions.Coniferous phytoliths were found throughout the entire humus-eluvial part of the soil profiles, with their contribution increasing upward.
Thus, the presence of Phragmites-type bulliform phytoliths and cubic foveolar phytoliths suggested that vegetation progressed from the intrazonal meadow-mire or mired forest phytocenoses towards the automorphic coniferous forest.The soil was identified as most likely being Calcic Stagnic Chernozem.

Calcic Stagnic Greyzemic Epidystric Umbrisols (Siltic) (Toeslope Site)
The soil profile was found to have a phytolith assemblage similar to those in the hilltop and midslope catena sites, being, in this respect, closer to that in Endocalcic Stagnic Albic Luvisols in the midslope site than to the hilltop type, as meadow gramineous phytoliths dominated by three to five times over the forest gramineous phytoliths, and Phragmites-type phytoliths were detected in the Aq and A' horizons (Figures 2f and 3C).
The phytolith assemblage profile in the toeslope site was similar to the profile in the midslope site, but Phragmites-type bulliform cell phytoliths were found throughout the entire A and A' horizons, thus confirming the longer period of hydromorphic soil formation in the trans-accumulative position and the relatively later start of the lessivage and albic material formation compared to the soils in the higher positions of the catena.The second humus horizon should be considered a residual from the Spodic Chernic Gleysols.

AMS 14 C Dating of the Eluvial and Humus Horizons
The optical and electron micropscopy data and EDS analysis of phytolith surfaces (Figure 4) showed no additional organic matter that could have compromised the radiocarbon dating results.According to phytolith analysis of the soil in the eluvial position, the start of the soil formation and its most active stage were related to the meadow stage, as evidenced by the dominance of meadow grasses' phytoliths over the forest grasses' ones.The tree stands of that stage had coniferous species during the entire soil formation process, with a certain influence of fires on phytocoenosis development and structure.
Presumably, the closest analogue of the second humus horizon in the eluvial position soil can be seen in (semi)hydromorphic soils (Calcic Stagnic Chernozem) in the forest-steppe zone.

Endocalcic Stagniс Albic Luvisols (Siltic, Abruptic, Cutanic, Epiloamic) (Transit Position)
The soil profile phytolith assemblage in the midslope site was found to be rather similar to the one in the hilltop site; that is, there were mostly forest and meadow gramineous phytoliths, the latter exceeding 3.5-3.8times the former, thus dominating in the humus-eluvial part of the soil profile (Figure 3B).
However, large parallelepiped-like morphotypes of Phragmites spp.phytoliths are found in small quantities in the E, A', and the upper A'/Bt horizons (Figures 2f and 3B), proving more periodically humid environmental conditions.Coniferous phytoliths were found throughout the entire humus-eluvial part of the soil profiles, with their contribution increasing upward.
Thus, the presence of Phragmites-type bulliform phytoliths and cubic foveolar phytoliths suggested that vegetation progressed from the intrazonal meadow-mire or mired forest phytocenoses towards the automorphic coniferous forest.The soil was identified as most likely being Calcic Stagnic Chernozem.

Calcic Stagnic Greyzemic Epidystric Umbrisols (Siltic) (Toeslope Site).
The soil profile was found to have a phytolith assemblage similar to those in the hilltop and midslope catena sites, being, in this respect, closer to that in Endocalcic Stagniс Albic Luvisols in the midslope site than to the hilltop type, as meadow gramineous phytoliths dominated by three to five times over the forest gramineous phytoliths, and Phragmites-type phytoliths were detected in the Aq and A' horizons (Figures 2f and 3C).
The phytolith assemblage profile in the toeslope site was similar to the profile in the midslope site, but Phragmites-type bulliform cell phytoliths were found throughout the entire A and A' horizons, thus confirming the longer period of hydromorphic soil formation in the trans-accumulative position and the relatively later start of the lessivage and albic material formation compared to the soils in the higher positions of the catena.The second humus horizon should be considered a residual from the Spodic Chernic Gleysols.

AMS 14 C Dating of the Eluvial and Humus Horizons
The optical and electron micropscopy data and EDS analysis of phytolith surfaces (Figure 4) showed no additional organic matter that could have compromised the radiocarbon dating results.The radiocarbon dating results of phytolith-occluded carbon (PhytOC) and HAII are given in Table 3.The sample, collected from the lower boundary of the A' horizon within Calcic Stagnic Greyzemic Epidystric Umbrisols had PhytOC dated as 6404-6466 cal year BP (median 6435 year cal BP) and humic acid dated as 4051-4185 cal year BP (median 4118 year cal BP) (Figure 5).The difference between the dates exceeded 2300 year and most likely results from the different rates of organic carbon rejuvenation in the soil and PhytOC.The carbon contained in different humic acid fractions, which is the most resistant fraction in HAII under increased biochemical activity and lessivage, is rejuvenated, thus decreasing the radiocarbon dates [45]; conversely, the carbon occluded in phytoliths is preserved, since it is subject to a much less aggressively destructive soil-forming environment.
The radiocarbon dating results of phytolith-occluded carbon (PhytOC) and HAII are given in Table 3.
The sample, collected from the lower boundary of the A' horizon within Calcic Stagnic Greyzemic Epidystric Umbrisols had PhytOC dated as 6404-6466 cal year BP (median 6435 year cal BP) and humic acid dated as 4051-4185 cal year BP (median 4118 year cal BP) (Figure 5).The difference between the dates exceeded 2300 year and most likely results from the different rates of organic carbon rejuvenation in the soil and PhytOC.The carbon contained in different humic acid fractions, which is the most resistant fraction in HAII under increased biochemical activity and lessivage, is rejuvenated, thus decreasing the radiocarbon dates [45]; conversely, the carbon occluded in phytoliths is preserved, since it is subject to a much less aggressively destructive soil-forming environment.The samples taken at the lower border of the second humus horizons of the hilltop and midslope sites were of a close PhytOC age: 5321-5220 cal year BP (median 5270 year cal BP) and 5843-5711 cal year BP (median 5777 year cal BP), respectively.If the difference in 5 cm between sample depths is taken into consideration, then the relative difference of 500 year between the radiocarbon dates can be discarded, and the samples can therefore be considered to have been formed synchronically.The same is true for the HAII radiocarbon dates, which were determined as 4600-4462 cal year BP (median 4531 year cal BP) for Eh horizons and 4696-4548 cal year BP (median 4622 year cal BP) for the A' horizons.However, the difference between the PhytOC and HAII radiocarbon dates exceeds 1000 year.The samples taken at the lower border of the second humus horizons of the hilltop and midslope sites were of a close PhytOC age: 5321-5220 cal year BP (median 5270 year cal BP) and 5843-5711 cal year BP (median 5777 year cal BP), respectively.If the difference in 5 cm between sample depths is taken into consideration, then the relative difference of 500 year between the radiocarbon dates can be discarded, and the samples can therefore be considered to have been formed synchronically.The same is true for the HAII radiocarbon dates, which were determined as 4600-4462 cal year BP (median 4531 year cal BP) for Eh horizons and 4696-4548 cal year BP (median 4622 year cal BP) for the A' horizons.However, the difference between the PhytOC and HAII radiocarbon dates exceeds 1000 year.
Comparison of the radiocarbon date differences in samples from the lower border of the second humus horizon of soils in the different catena sites showed that the carbon rejuvenation rate in HAII was almost 2.5 times higher in hilltop and midslope sites than in the toeslope site.
In the sample taken from the eluvial horizon of the midslope site soil, the PhytOC radiocarbon date was determined as 4005-3936 cal year BP (median 3970 year cal BP), whereas the date for HAII was estimated as 2632-2544 cal year BP (median 2588 year cal BP), the difference being ca.1400 year.This confirms the process of carbon rejuvenation both in soil and in phytoliths, the rate being slower in phytolith-occluded carbon compared to soil organic carbon.
In the sample taken from the E horizon of the eluvial position soil, the PhytOC and HAII radiocarbon dates were rather close, i.e., 2151-1902 cal year BP (median 2026 year cal BP) and 2410-2031 cal year BP (median 2220 year cal BP), respectively, with the difference being approximately 200 year.
As Table 3 and Figure 5 show, soil phytolith profiles are rather resistant to transformation, as can be concluded from their close radiocarbon dates of PhytOC from the second humus horizons in the different catena positions, namely, 6.5-5.3 years cal BP.Therefore, soil phytolith concentrations and phytolith assemblages should be considered synchronic to the dark-humus Mid-Atlantic stage of soil formation.
Recent studies investigating sources of PhytOC and assessing phytolith potential for radiocarbon dating showed multiple sources for PhytOC [46,47].Some researchers established that one of the PhytOC sources may be soil carbon, which significantly decreases the age of recurrent plants' phytoliths.However, to establish the starting date of soil formation, the use of phytoliths from soil samples taken at the lower border of the humus horizon is justified, despite it being a source of carbon in phytoliths.
The PhytOC radiocarbon date would be especially close to the start of humus formation in those cases when soil organic carbon, which undergoes rejuvenation, is analysed in different HA fractions.
In addition, the radiocarbon date regression with soil depth (R 2 = 0.808, Figure 5) in all catena positions allows the assumption that soil carbon was not the leading source of PhytOC; if this had been the case, the profile distribution would have been more even with no explicit regression pattern.The main source of phytolith profile rejuvenation is provided by new phytoliths becoming mixed into the soil due to plant litter decomposition and by their subsequent translocation down towards the lower border of the eluvial portion of the soil profile [48].
Thus, different materials used for radiocarbon dating were shown to rejuvenate with depth, but the depth and the rate of the new carbon input depends on the carbon-preserving environment and the rate of soil textural differentiation, resulting in the redistribution of organic matter downwards.The largest degree of rejuvenation manifested in the soil in the hilltop site, whereas that of the toeslope showed the lowest degree of rejuvenation.The phytolith-occluded carbon was found to be more recalcitrant to rejuvenation than the soil organic carbon.

Conclusions
The combined use of radiocarbon dating, phytolith analysis, and conventional soil analysis methods to study co-evolved soils with complex organic matter profiles revealed the changes in the soil formation environment and the corresponding shifts in soil types.
The topographical divergence of soil phytolith profiles reflects the relief effect on the development of specific soil type combinations, accounting for the major elements of the regional mid-Holocene soil cover.The leading elementary soil-forming processes were humus accumulation and hydrogenic accumulation of calcium carbonates.
In the hilltop site in Endocalcic Stagnic Albic Luvisols, the evolutionary changes were shown by the shift from the meadow phytocenosis (Calcic Stagnic Chernozem) to the forest phytocenosis, i.e., from the soddy type of soil formation to textural differentiation of the soil profile with albic material formation.
In the midslope site, the environment was more humid from the start, favouring a phytocenosis with features of the meadow-mire.The shift from the meadow-mire environment (with Spodic Chernic Gleysols) to the forest type with clay illuviation was gradual, occurring via the meadow stage with Calcic Stagnic Chernozem.
Comparison of the soil radiocarbon date profiles, obtained using different organic carbon containing material, revealed a higher recalcitrance to rejuvenation of phytolith-occluded carbon compared to humic acid carbon.Carbon rejuvenation of different materials depends on the soil textural differentiation along a particular catena.

Figure 1 .
Figure 1.Schematic presentation of the catena studied.Vegetation composition and polygenetic soils with the second humus horizon on the study sites on the Vasyugan sloping plain: А-Endocalcic Stagniс Albic Luvisols (Siltic, Abruptic, Cutanic, Epiloamic) on the hilltop; B-Endocalcic Stagniс Albic Luvisols (Siltic, Abruptic, Cutanic, Epiloamic) in the midslope; C-Calcic Stagnic Greyzemic Epidystric Umbrisols (Siltic) on the toeslope.A pit displaying the eluvial part and the upper portion of the illuvial layer (below the lower border of the second humus horizon, with no explicit features of windfall phenomena) was dug on both catena positions.1-Soils with the second humus horizon [30].

Figure 1 .
Figure 1.Schematic presentation of the catena studied.Vegetation composition and polygenetic soils with the second humus horizon on the study sites on the Vasyugan sloping plain: A-Endocalcic Stagnic Albic Luvisols (Siltic, Abruptic, Cutanic, Epiloamic) on the hilltop; B-Endocalcic Stagnic Albic Luvisols (Siltic, Abruptic, Cutanic, Epiloamic) in the midslope; C-Calcic Stagnic Greyzemic Epidystric Umbrisols (Siltic) on the toeslope.A pit displaying the eluvial part and the upper portion of the illuvial layer (below the lower border of the second humus horizon, with no explicit features of windfall phenomena) was dug on both catena positions.1-Soils with the second humus horizon [30].

Figure 4 .
Figure 4. Images of phytoliths extracted from soils: (a)-Optical microscopy images of phytoliths; (b)-SEM image of phytoliths with EDS analysis of phytolith surface.

Figure 4 .
Figure 4. Images of phytoliths extracted from soils: (a)-Optical microscopy images of phytoliths; (b)-SEM image of phytoliths with EDS analysis of phytolith surface.

Table 1 .
Field description of soils.

Table 2 .
Soil chemical and physical properties.

Table 3 .
Results of AMS14C dating of the eluvial and second humus horizons.