Soil Organic Matter Accumulation and Carbon Fractions along a Moisture Gradient of Forest Soils

The aim of the study was to present effects of soil properties, especially moisture, on the quantity and quality of soil organic matter. The investigation was performed in the Czarna Rózga Reserve in Central Poland. Forty circular test areas were located in a regular grid of points (100 × 300 m). Each plot was represented by one soil profile located at the plot’s center. Sample plots were located in the area with Gleysols, Cambisols and Podzols with the water table from 0 to 100 cm. In each soil sample, particle size, total carbon and nitrogen content, acidity, base cations content and fractions of soil organic matter were determined. The organic carbon stock (SOCs) was calculated based on its total content at particular genetic soil horizons. A Carbon Distribution Index (CDI) was calculated from the ratio of the carbon accumulation in organic horizons and the amount of organic carbon accumulation in the mineral horizons, up to 60 cm. In the soils under study, in the temperate zone, moisture is an important factor in the accumulation of organic carbon in the soil. The highest accumulation of carbon was observed in soils of swampy variant, while the lowest was in the soils of moist variant. Large accumulation of C in the soils with water table 80–100 cm results from the thick organic horizons that are characterized by lower organic matter decomposition and higher acidity. The proportion of carbon accumulation in the organic horizons to the total accumulation in the mineral horizons expresses the distribution of carbon accumulated in the soil profile, and is a measure of quality of the organic matter accumulated. Studies have confirmed the importance of moisture content in the formation of the fractional organic matter. With greater soil moisture, the ratio of humic to fulvic acids (HA/FA) decreases, which may suggest an increase in carbon mobility in soils.


Introduction
Accumulation of organic carbon in soils is currently discussed for many reasons.With the increased pool of carbon stored in terrestrial ecosystems, a reduction of carbon dioxide emissions into the atmosphere and prevention from climate changes are observed [1,2].In the carbon cycle on Earth, soil is its most important reservoir.It is estimated that carbon contained in soil constitutes 75% of the total pool of organic carbon, exceeding twice the resources of carbon in the atmosphere [3].The accumulation of organic carbon in forest soils depends on many factors, but predominantly on climatic conditions [4], soil properties [5], soil moisture [6,7], the type of plant cover and type of forest management [8][9][10][11].The site conditions and plant coverage are strongly related to each other.Here, the site conditions play the main role via driving formation of certain species composition and the structure of plant cover, thus determining the type and direction of soil formation processes [12].Among the physical properties of the soil, clay fraction content acts in the formation of soil organic matter (SOM) accumulation.Many authors have studied the influence of clay and silt content on the protection of SOM, and have found a relationship between the content of SOM and soil texture, mainly with the content of fine fractions (silt and clay) [13,14].According to Wosten et al. [15], the accumulation and stabilization of soil carbon are associated with the content of clay; however, it may also depend on the soil moisture.Bauer et al. [16] and Meersmaus et al. [17] reported a strong influence of soil moisture on the decomposition of SOM, and consequently, on the carbon accumulation in the soil.SOM decomposition dynamics are influenced by pH making it an important variable, as an increase in SOM content is connected with changes of other soil parameters, namely pH [18].
By specifying the total content of soil carbon (C) [19] or by estimating the available fractions of SOM [20], the quantity and quality of SOM can be determined.To evaluate the changes in soil C dynamics due to agricultural use, the use of SOM chemical fractions was more effective than the determination of total SOM in Guimarães et al.'s [21] study.The SOM fraction can be used as an indicator of changes in the soil.The amount of SOM fractions reflects the potential mobility of C in the soil and the rate of SOM decomposition [22].
Until now, a limited number of studies reporting the relationship between the forest site condition, the accumulation of SOM, and fractions of humic substances has been published.There is a lack of studies on the effect of different moisture levels, soil properties, and natural forest plant communities on the quantity and quality of SOM in temperate climates.We hypothesize that: (1) soil moisture is an important factor in the accumulation of organic carbon in the temperate zone, and the highest accumulation of carbon was connected with a swampy variant of soils; (2) moisture has an effect on the proportion of carbon accumulation in the organic horizon to the accumulation in the mineral horizons; and (3) fractional composition of the organic matter reflects the rate of SOM decomposition in forest soils and depends on soil moisture.

Study Area
The study was performed in the Czarna Rózga Reserve in Central Poland (Figure 1).The reserve area is 185.6 ha, and within the limits of the reserve stands is the dominant species common alder (Alnus glutinosa), which are of natural origin.The silver fir (Abies alba), pedunculate oak (Quercus robur), common hornbeam (Carpinus betulus) and common ash (Fraxinus excelsior) play the role of admixtures.The study area is characterized by the following climatic conditions: the average annual rainfall is 649 mm, the average daily temperature of the warmest month (July) is 17.8 • C and the coldest month (January) is −3.4 • C, and the length of the vegetation season lasts 200-210 days.Sample plots were located in the area with a predominance of fluvioglacial sand and loam with Gleysols, Cambisols and Podzols [23].The groundwater level is the basic reason for differentiation of soil subtype and vegetations communities.
mainly with the content of fine fractions (silt and clay) [13,14].According to Wosten et al. [15], the accumulation and stabilization of soil carbon are associated with the content of clay; however, it may also depend on the soil moisture.Bauer et al. [16] and Meersmaus et al. [17] reported a strong influence of soil moisture on the decomposition of SOM, and consequently, on the carbon accumulation in the soil.SOM decomposition dynamics are influenced by pH making it an important variable, as an increase in SOM content is connected with changes of other soil parameters, namely pH [18].
By specifying the total content of soil carbon (C) [19] or by estimating the available fractions of SOM [20], the quantity and quality of SOM can be determined.To evaluate the changes in soil C dynamics due to agricultural use, the use of SOM chemical fractions was more effective than the determination of total SOM in Guimarães et al.'s [21] study.The SOM fraction can be used as an indicator of changes in the soil.The amount of SOM fractions reflects the potential mobility of C in the soil and the rate of SOM decomposition [22].
Until now, a limited number of studies reporting the relationship between the forest site condition, the accumulation of SOM, and fractions of humic substances has been published.There is a lack of studies on the effect of different moisture levels, soil properties, and natural forest plant communities on the quantity and quality of SOM in temperate climates.We hypothesize that: (1) soil moisture is an important factor in the accumulation of organic carbon in the temperate zone, and the highest accumulation of carbon was connected with a swampy variant of soils; (2) moisture has an effect on the proportion of carbon accumulation in the organic horizon to the accumulation in the mineral horizons; and (3) fractional composition of the organic matter reflects the rate of SOM decomposition in forest soils and depends on soil moisture.

Study Area
The study was performed in the Czarna Rózga Reserve in Central Poland (Figure 1).The reserve area is 185.6 ha, and within the limits of the reserve stands is the dominant species common alder (Alnus glutinosa), which are of natural origin.The silver fir (Abies alba), pedunculate oak (Quercus robur), common hornbeam (Carpinus betulus) and common ash (Fraxinus excelsior) play the role of admixtures.The study area is characterized by the following climatic conditions: the average annual rainfall is 649 mm, the average daily temperature of the warmest month (July) is 17.8 °C and the coldest month (January) is −3.4 °C, and the length of the vegetation season lasts 200-210 days.Sample plots were located in the area with a predominance of fluvioglacial sand and loam with Gleysols, Cambisols and Podzols [23].The groundwater level is the basic reason for differentiation of soil subtype and vegetations communities.

Sampling Scheme
Forty circular test areas with a surface of 0.1 hectare were established within the area of the reserve.Test surfaces were located in a regular grid of points (100 × 300 m) (Figure 1).Each plot was represented by one soil profile located at the plot's center.The surface horizons (organic-O, mineral-A, or organic/mineral-AM) and the mineral topsoil (AE or AG) horizon were sampled.The subsoil gleyic horizon (G) was sampled from the bottom of the topsoil horizon to 60 cm deep.The surface horizons' samples constituted a cumulative sample from five subsamples around the soil profile.
On each test surface (Table 1), site conditions including vegetation types, soil types, and moisture levels (based on the groundwater depth and using soil moisture monitoring sensors in the daily interval in May and June) were recorded.The obtained results revealed three types of vegetation communities (Abietetum albae, Tilio-Carpinetum ficarietosum and Fraxino-Alnetum) within the reserve and water table from 0 to 100 cm in soil profiles.The three variants of soil moisture content (first variant with water table 80-100 cm (WT 80-100 ), second variant with water table 40-50 cm-moist, and last variant with water table 0-30 cm-swampy) were distinguished.In further analysis, the above sequence of soil moisture content (WT 80-100 , moist and swampy) was considered.Ten areas that represent the Abietetum albae from Stagnosols and Podzols have been identified for the WT 80-100 variant.The moist variant represents the Tilio-Carpinetum ficarietosum (19 areas).Gleysols dominated this variant.The third swampy variant (11 areas) is represented by the Fraxino-Alnetum dominated by Mollic Gleysols.

Laboratory Analysis of Soil
Soil samples obtained in the field were dried and sieved through 2.0 mm mesh.Using the potentiometric method, the pH of the samples was analyzed in H 2 O and KCl.By using a laser diffraction technique (Analysette 22, Fritsch, Idar-Oberstein, Germany), the soil texture was determined.Carbon (C t ) and nitrogen (N) contents were measured with an elemental analyzer (LECO CNS TrueMac Analyzer) (Leco, St. Joseph, MI, USA).Using the method of Kononowa and Bielczikowa, fractional composition of humus was determined in which extraction is performed in a mixture of 0.1 M NaOH and 0.1 M Na 4 P 2 O 7 [24].In order to obtain a humin (Hm), humic acid (HA), and fulvic acid (FA) fraction, the chemical fraction was conducted.The relationships HA/FA and (HA+FA)/Hm were also calculated.Exchangeable Ca, Mg, K and Na were determined by inductively coupled plasma optical emission spectrometry (ICP-OES) (iCAP 6500 DUO, Thermo Fisher Scientific, Cambridge, UK) in 1 M ammonium acetate at pH 7.0 extracts.The sum of base cations (BC) and effective base saturation (BS) were calculated.Available phosphorus was determined by Bray-Kurtz's method.The hydrolytic acidity was determined in the samples after extracting 5 g of soil with 30 mL 1 mol•L −1 (CH 3 COO) 2 Ca (shaking time 1 h), followed by filtration.After extracting 5 g of soil with 30 mL KCl, exchangeable acidity was determined.Soil on the filter was washed several times by extractant solution, up to a volume of 200 mL.About 25 mL of the obtained solution was titrated by potentiometric titration (automatic titrator, Mettler Toledo Inc., Columbus, OH, USA) to a pH of 8.2.For determining the exchangeable aluminum and hydrogen, NaF was added to the KCl extract and titrated with NaOH.Total cation exchange capacity (CECt) was calculated as a sum of base cations (BC) and hydrolytic acidity (Ah); and effective cation exchange capacity (CECe) was calculated as a sum of base cations (BC) and exchangeable acidity (Aex).Using samples with intact structure collected to metal cylinders, bulk density was determined via the drying-weighing method [25].
The soil organic carbon stock (SOC S ) was calculated as a sum of its total content at particular soil horizons.In each of the analyzed soil samples, the SOC S was calculated in the soil block area of 1 m 2 and a depth from the surface to 60 cm.The calculation of SOC S for the particular depths was made by summing the SOC S at subsequent genetic soil horizons according to the formula: where C is the organic carbon content at the subsequent genetic horizon (%), D is the bulk density (g•cm −3 ), and m is the thickness of the horizons (cm).Moreover, a Carbon Distribution Index (CDI) was provided and calculated from the ratio of the carbon accumulated in organic horizons and the amount of organic carbon accumulated in the mineral horizons up to 60 cm according to the formula: (2)

Statistical Analysis
Basic descriptive statistics, arithmetic mean and standard deviation were used to estimate the soil properties' variability.Tukey's honest significant difference (HSD) test was used to assess differences between means for carbon accumulation in soils.To evaluate the relationships between stock of soil carbon, other soil properties, forest species, and moisture, the Principal Components Analysis (PCA) method was used.Multiple regression models were developed to describe the relationship between the estimated values of carbon content and other soil properties.The statistical significance of the results was verified at the significance level of alpha = 0.05.All statistical analyzes were performed with Statistica 12 software (2012).

Results
The investigated soils are characterized by diversification of properties (Table 2).The values of pH H 2 O and pH KCl were in the range from 3.85 to 7.30 and from 3.04 to 5.78, respectively.A significantly greater pH was noted in soils of moist and swampy variants (Table 2).For cation content, a similar dependence was recorded (Table 2).The greater Ah and Aex were observed in soils with water table 80-100 cm (Table 3).The texture of the investigated soils was dominated by sand (on average 50-72% in the deeper horizons) and silt (on average 23-42% in deeper horizons) with addition of clay (on average 5-8% in the deeper horizon) (Table 3).The sand content (only in the deepest horizons) was significantly greater in WT 80-100 variant (Table 3).The silt and clay content (only in the deepest horizons) was higher in soils of moist and swampy variants (Table 3).The analyzed soils clearly differ in their moisture content (Table 4).Statistically significant higher soil moisture content was recorded in the soils of the swampy variant (mean soil moisture of the surface humus horizon was 46.85%).The lowest soil moisture content was recorded in the soils of the WT 80-100 variant (Table 4).The differences in the carbon content of the examined soils were clearly noted.A statistically higher carbon content in the surface horizons was noted in the soils of the WT 80-100 variant.The average content of carbon in the surface horizon of this variant was 241.87 g•kg −1 , and was 101.10 g•kg −1 and 135.96 g•kg −1 in the soils of the moist and swampy variants, respectively (Table 4).The highest carbon content in the lower lying horizons was recorded for soils of the swampy variant.The mean carbon content at these horizons was 44.13 g•kg −1 and 8.51 g•kg −1 , respectively.In the case of WT 80-100 and moist variant of soils in the deeper horizons, the carbon content was 35.07-4.87g•kg −1 and 24.29-3.95g•kg −1 , respectively.Similar dependencies concerned nitrogen content (Table 4).The lowest C/N, which expresses the degree of organic matter distribution, was recorded in the soils of the swampy and moist variants (the C/N ratio ranged from 11.8 to 14.2).The highest C/N ratio, from 19.6 to 20.3 depending on the horizon, was recorded in the WT 80-100 variant of soils (Table 4).Soil organic matter (SOCs) concentration in the 0-60 cm soil horizon differed significantly among moisture gradients and followed the order: swampy > WT 80-100 > moist variants.The average SOCs up to a depth of 60 cm in the soils of moist variant is lower and is estimated at 12.61 Mg•ha −1 .The higher accumulation was found in the soils of WT 80-100 variant, where the average SOCs was 14.09 Mg•ha −1 .A significantly greater accumulation was noted in the soils of swampy variant (19.02Mg•ha −1 ).The ratio between carbon accumulation in the organic horizon and carbon accumulation in the mineral horizons, up to 60 cm of soils (CDI) of estimated variants, is notably diverse (average 0.299-0.837).The lowest values of CDI are found in the soils of moist variant in which the CDI indicator is 0.299.The soils of swampy variant are characterized by higher values of the CDI index-0.329, and the CDI index of soils of WT 80-100 variant approached 1.0 and is 0.837 (Table 4).
Concentrations of C in fulvic acid (FA), humic acid (HA) and humin fractions in the different soil horizons of different moisture variants are presented in Table 4. Mean values of HA, FA and humin varied from 11.9% to 0.79%, 4.43% to 0.63%, and 9.79% to 0.25%, respectively.In the fractions of SOM, humic acids dominate over fulvic acids and humins.Soils of the WT 80-100 variant are characterized by the highest content of humic and fulvic acid (Table 4).Significant differences between the soils in terms of humic acid were observed in two surface horizons (swampy variant = WT 80-100 variant > moist variant = swampy variant).No significant differences between soil variant in terms of fulvic acid and humin were noted (Table 4).All SOM fractions decreased significantly with depth for all variants of soil.The HA/FA ratio ranged from 4.70 to 1.05.The lowest HA/FA ratio regardless of depth was in the soil of WT 80-100 variant, and the highest in the soil of moist and swampy variant.The relationship of the sum of C humic and fulvic acids to C humin (HA+FA/Hm) assumes the values of 54.3-0.68 (Table 4).In the soil of WT 80-100 variant, the highest HA/FA ratio and HA+FA/Hm, as compared to moist and swampy variant, was recorded and was statistically significant (Table 4).
Multiple regression models explained 89% of the variance in the C content (Table 5).C content was largely dependent on hydrolytic acidity and moisture.The relationships between carbon content, hydrolytic acidity and moisture are shown in Figure 2. The increase in C accumulation in soils correlated with high moisture and high acidity.
A projection of the variables on the factor-plane clearly demonstrated a correlation between the variables and soil moisture (supplementary variables).Two main factors had a significant total impact (60.62%) on the variance of the variables in soil.Factor 1 explained 35.76% of the variance of the examined properties and Factor 2 explained 24.86% of the variance (Figure 3).The soil of moist and swampy variants had a higher pH.In contrast, the soils of WT 80-100 variant were characterized by a high C/N ratio.The soils of WT 80-100 variant contained more humic acids.The fulvic acids and humins were connected with nitrogen content.

Discussion
Results present in this study demonstrate that soil moisture is a factor that most strongly influences the carbon accumulation in forest soils.Compared to the soil with water table from 80 to 100 cm and moist soils, those with the highest moisture content (swampy variant) are characterized by the higher C accumulation.Soil organic matter (SOCs) concentration differed significantly among moisture gradient and followed the order: swampy > WT80-100 > moist variants.The greater SOM accumulation in the soils is affected not only by the high moisture, but also periodic deficiency of water.Our results show the high SOC stock on the sites with low moisture (WT80-100 variant) and in swampy variant, and that site condition promoted lower rate decomposition.The high SOC stock in the soils of WT80-100 variant is affected by greater acidity (surface horizons are not enriched due to low groundwater level).In the soils of swampy variant, there is an effect of anaerobic processes: the soils better aerated in the moist variants are characterized by faster and more efficient decomposition of SOM.Jandl et al. [26] and Stendahl et al. [27] related the accumulation of SOC with site conditions.A large stock of carbon can be associated with the accumulation of organic matter in acidic environments.A similar situation was noted for the soils of the WT80-100 variant where surface horizons were under coniferous species influence.Tree species affect pH, and the specific composition of organic matter may lead to the modification of soil acidity and nutrient pools [28][29][30].pH is an important variable that influences SOM decomposition dynamics [31].Our results confirmed the importance of moisture and acidity of soils on soil organic matter accumulation.In our study, 89% of the variance in the C content was explained by the hydrolytic acidity and moisture of the soils.Li et al.'s [32] results indicated that moisture level had an effect on carbon accumulation.At the same time, the increase of soil's acidity results in deterioration of the conditions of decomposition of dead organic matter, slowing down their humification rate and microbial activity.

Discussion
Results present in this study demonstrate that soil moisture is a factor that most strongly influences the carbon accumulation in forest soils.Compared to the soil with water table from 80 to 100 cm and moist soils, those with the highest moisture content (swampy variant) are characterized by the higher C accumulation.Soil organic matter (SOCs) concentration differed significantly among moisture gradient and followed the order: swampy > WT 80-100 > moist variants.The greater SOM accumulation in the soils is affected not only by the high moisture, but also periodic deficiency of water.Our results show the high SOC stock on the sites with low moisture (WT 80-100 variant) and in swampy variant, and that site condition promoted lower rate decomposition.The high SOC stock in the soils of WT 80-100 variant is affected by greater acidity (surface horizons are not enriched due to low groundwater level).In the soils of swampy variant, there is an effect of anaerobic processes: the soils better aerated in the moist variants are characterized by faster and more efficient decomposition of SOM.Jandl et al. [26] and Stendahl et al. [27] related the accumulation of SOC with site conditions.A large stock of carbon can be associated with the accumulation of organic matter in acidic environments.A similar situation was noted for the soils of the WT 80-100 variant where surface horizons were under coniferous species influence.Tree species affect pH, and the specific composition of organic matter may lead to the modification of soil acidity and nutrient pools [28][29][30].pH is an important variable that influences SOM decomposition dynamics [31].Our results confirmed the importance of moisture and acidity of soils on soil organic matter accumulation.In our study, 89% of the variance in the C content was explained by the hydrolytic acidity and moisture of the soils.Li et al.'s [32] results indicated that moisture level had an effect on carbon accumulation.At the same time, the increase of soil's acidity results in deterioration of the conditions of decomposition of dead organic matter, slowing down their humification rate and microbial activity.On the other hand, the increase in organic matter accumulation results in an increase in sorption capacity, which, under acidic forest soil condition and lower level of groundwater, includes aluminum and hydrogen ions.The increase in soil moisture was favorable for the accumulation of water soluble carbon.In our study, soils in the swampy variant are characterized by the highest accumulation of C, which can be combined with the accumulation of organic debris in the peat-forming process.These studies confirm earlier findings, that although the peat covers only 3% of the Earth's surface, it constitutes an important reservoir of carbon [33,34].Under anaerobic conditions, 2-6 times less CO 2 was produced than under aerobic conditions [35].Furthermore, the presence of water in the soil profile facilitates dissolution of minerals, and during the capillary transport, also their transportation to higher horizons where they can form stable conjunctions to the humic substances.In addition, higher soil moisture results in higher biomass production of trees and plant cover of lower layers of the stand, which results in increased accumulation of organic material into the forest floor and soil.
In our study, CDI was used to evaluate the distribution of carbon accumulated in the soil profile and as a measure of quality of the accumulated organic matter.The advantage of accumulation of carbon in the mineral horizons over accumulation in the organic horizons indicates a greater efficiency of the processes of organic matter decomposition.In numerous studies [36], the greater resistance of SOM, in combination with soil minerals, is observed.In comparison to the accumulation occurring in the organic horizons, the predominance of carbon accumulation in the mineral horizons can be observed in soils of moist and swampy variants.The soils of moist variants were occupied by a very diverse species composition; the presence of admixture species-mainly hornbeam, alder and maple-provides abundant input of easily decaying organic matter to the soil.The positive impact of hornbeam on the accumulation of SOM was also described by other authors [31,37].Stand mixing improves soil properties, especially the quality of SOM and biochemical properties [38].The dominance of coniferous species in the forest (WT 80-100 variant) changes the CDI index.In coniferous forests with fir, the C accumulation in the mineral horizons is reduced in relation to the total accumulation in the whole soil profile.The reason for the higher accumulation of C at the organic horizons is undoubtedly the slow decomposition of organic debris.Litter decomposition is particularly influenced by decomposer organisms [39].The conifer litter decomposes slower than that from hardwood species [40].The soils of coniferous forests are characterized by low pH, and the decomposition of litter is mainly conducted by fungal organisms.Fungi dominate in acidic soils, whilst bacteria in soils are characterized by a neutral or alkaline pH [41].
Forest soils in the studied moisture content gradient differ in the accumulation of carbon and the fractional composition of the soil organic matter.In the studied soils, a large proportion of humic acids were found, which prevailed in the moist variants.This is in line with Tavares and Nahas' study [42].More humic acids than fulvic acids indicate the potentially low mobility of carbon accumulated in the surface soil horizons [21].In our studies, the ratio of humic to fulvic acids (HA/FA) changes as moisture content of soils increases, suggesting an increase in carbon mobility in these soils.More humic acids than fulvic acids indicate high humification rates.The differences in soil moisture contents influences the amounts and chemical composition of the soil humic fraction [32].The results implied that a high moisture level was unfavorable for the accumulation of humic fraction in the soil.In our experiment, we used the ratio between the sum of humic and fulvic acids and humin ((HA+FA)/humin) as an indicator of loss of SOM through the soil profile.The highest ratio ((HA+FA)/humin) was determined in the soils with the lowest moisture.As the soil moisture content increased, the ratio decreased.The lower ratio ((HA+FA)/humin) indicates the sufficient humin level in the soils with the highest moisture content.

Conclusions
In the soils under study, in the temperate zone, moisture is an important factor in the accumulation of organic carbon in the soil.The highest accumulation of carbon was observed in soils of swampy variant, while the lowest in the soils of moist variant.In the soil of swampy variant, the anaerobic conditions affect the organic matter decomposition leading to slower processes.The increased accumulation of C in soils of WT 80-100 variant is caused by other factors.Large accumulation of C in the soils of the WT 80-100 variant results from the thick organic horizons that are characterized by lower organic matter decomposition and higher acidity.In the WT 80-100 variant, the higher acidity is simultaneously an effect of the stand species composition and lack of water neutralization.The proportion of carbon accumulation in the organic horizons to the total accumulation in the mineral horizons expresses the distribution of carbon accumulated in the soil profile, and is a measure of quality of the organic matter accumulated.The advantage of accumulation at the mineral horizons over accumulation at the organic horizons confirms a greater efficiency of decomposition processes of organic matter, which, thanks to more advanced transformations and connecting with mineral soil, can be considered stable.Studies have confirmed the importance of moisture content in the formation of the fractional organic matter.The increase in soil moisture content changes the ratio of humic to fulvic acids (HA/FA), which may suggest an increase in carbon mobility in soils.

Figure 2 .
Figure 2. Three-dimensional plot of the relationship between soil organic carbon stock (SOCs), moisture and hydrolytic acidity (Ah) of soils (circles with WT80-100-group of plots with water table 80-100 cm, circles with M-group of moisture plots, circles with S-group of swampy plots); regression plane marked by black lines.

Figure 2 .
Figure2.Three-dimensional plot of the relationship between soil organic carbon stock (SOCs), moisture and hydrolytic acidity (Ah) of soils (circles with WT 80-100 -group of plots with water table 80-100 cm, circles with M-group of moisture plots, circles with S-group of swampy plots); regression plane marked by black lines.

Figure 3 .
Figure 3. Diagram of Principal Components Analysis (PCA) with projection of variables on a plane of the first and second factor in the surface horizon (WT80-100, moist and swampy-soil moisture gradient; HA-humic acids, FA-fulvic acids, Hm-humins).

Figure 3 .
Figure 3. Diagram of Principal Components Analysis (PCA) with projection of variables on a plane of the first and second factor in the surface horizon (WT 80-100 , moist and swampy-soil moisture gradient; HA-humic acids, FA-fulvic acids, Hm-humins).

Table 1 .
Characteristic of vegetation and soils in groups of studied areas.

Moisture Gradient Plant Communities Dominant Stand Species Dominant Ground Cover Species Type of Soil Depth of Groundwater Level (in the Spring) (cm) Humus Type
WT 80-100 -variant with water table 80-100 cm.

Table 2 .
Chemical properties of study soils.
n.d.no determine, mean ± standard deviation, small letters (a, b, c) in the upper index of the mean values mean significant differences in soil horizons between the moisture variant, WT 80-100 , moist and swampy-soil moisture gradient; VCS-very coarse sand, CS-coarse sand, MS-medium sand, FS-fine sand, VFS-very fine sand, CSi-coarse silt, FSi-fine silt; BD-bulk density.

Table 4 .
Humus substances, organic carbon and nitrogen content, moisture of study soils.

Table 5 .
Multiple regression analysis for carbon content based on soil characteristic.R 2 describes the percentage of explained variance, β is the regression coefficient for the given equation parameter and p is the significance level for the equation parameter.

Table 5 .
Multiple regression analysis for carbon content based on soil characteristic.R 2 describes the percentage of explained variance, β is the regression coefficient for the given equation parameter and p is the significance level for the equation parameter.