Soil Respiration in the Profiles of Forest Soils in Inland Dunes

Forest soil profiles of two dunes within the European belt of inland dunes were analysed in the laboratory. We carried out respirometric measurements of carbon dioxide production and oxygen consumption for every horizon of the studied soils while simultaneously quantifying the organic matter and humidity. Oxygen consumption and carbon dioxide excretion decreased exponentially with depth. The oxygen consumption decrease was less rapid than the decrease in carbon dioxide production. We found a statistical significant linear dependence between oxygen consumption and carbon dioxide excretion, and organic matter content and soil water capacity. Respiration processes in the profiles were divided into two strata; oxygen respiration dominated in the first and fermentation processes in the second. We estimated total respiration in the studied profiles for an area of 1 m down to around 1 m depth. We concluded that when assessing the soil’s role in carbon cycling in an ecosystem, it is necessary to consider both the respiratory and fermentation strata, as both produce large quantities of carbon dioxide. The main factor determining carbon dioxide production intensity is organic matter content; thus the distribution of organic matter in the soil profile determines carbon cycling intensity.

marily concern the surface horizon or measurements of CO 2 excretion from the surface. The review of Kuzyakow [3] discussed the diversity of methods applied in this field, which results from, amongst others, the range of factors driving CO 2 excretion from the soil. Moreover, the paper warned against making generalisations and suggested caution in drawing conclusions from published results. Field investigations of soil respiration in the soil profile are usually based on probe measurements of changes in amounts of carbon dioxide [4] [5]. On the other hand, laboratory investigations are sparse due to methodical limitations (difficulties relating to the maintenance of anaerobic conditions) and primarily concern the reactions of soil to changes in environmental factors [6].
There have been many laboratory investigations of soil (soil samples) and field investigations of the soil surface in different ecosystems, and these have shown that a large portion of the organic carbon in soil is found in the surface horizon [1] [7]. However, to understand overall carbon cycling in an ecosystem it is necessary to understand the dynamics of carbon in the deeper soil horizons [6] [8]. In reality, these deeper soil horizons excrete low amounts of carbon dioxide, but they are often deep [9] and their features depend on the type of soil in question. The excretion of CO 2 from the soil surface can tell us much, but it must be remembered that excretion is the product of two processes, carbon dioxide production in the soil and diffusion [10] [11]. Diffusion depends on the physical characteristics of soil [12]. Metabolic processes are, of course, also widely dependent on the physical and chemical characteristics of soil. However, Fierer and co-authors [6] showed that basic environmental factors like temperature, moisture, and organic substance content affect the metabolism of the soil surface differently than that of the deeper horizons. This difference is caused by differences in the microorganism populations between the horizons. Deeper soil horizons diverge from the surface in a variety of ways, but mainly in their oxygenation and organic matter content. In the present study, it was assumed that to understand respiration processes in the deeper soil horizons it is necessary to understand their reactions to basic environmental factors under anaerobic conditions, especially humidity and organic matter content [6] [13]. It was also assumed that the thermal conditions in the deeper horizons are relatively stable and that their impact on respiration is generally understood [14]. Humidity and organic matter content in both the surface and deeper horizons depend on their location within an ecosystem and the relief of the terrain. Location within the relief also to a large extent defines the soil profile structure. The research of Wiaux [15] showed that carbon dioxide excretion from a soil profile at the peak of a loess hill significantly differs from that at the foot of the hill. This was evidently due to differences in the thermal-humidity conditions, in both the surface and deeper horizons. The investigations of Wiaux [15] primarily concerned CO 2 excretion from the soil surface, but measurements carried out at different depths showed how the physical characteristics of every soil horizon can impact the intensity of CO 2 excretion.
The main objective of the present study was to investigate under laboratory conditions respiration in soil sampled from various horizons of two soil profiles. We analysed the oxygen consumption, CO 2 excretion and respiratory quotient RQ. The investigations were carried out in anaerobic conditions while maintaining the appropriate humidity level (GWC, WHC) and soil organic matter content (SOM). To minimise differences in physical characteristics between given soil horizons, two soil profiles from the peaks of inland dunes in the European dune belt were chosen for analysis. These soils were composed of aoelian sands with similar particle size compositions and were from locations with similar relief. These were forest soils-Podzol soils. In the European dune belt, aeolian sands have similar mineral compositions (68% -98% quartz) [16] [17] [18].

The Study Site
The studied soil profiles were sampled from the European dune belt (European sand belt) [17] [19] [20] [21] et al. Basic profile 1 was from Ukraine (Male Polissia region), and basic profile 2 was from the ancient Vistula River valley in Poland. Additionally, to complement the analysis, we included profile 3, which was also obtained from Male Polissia in Ukraine (2017 May). The climate of the study area is temperate. Some years it varies from Atlantic type to more continental. Average air temperature in 2017 was equal 10.3˚C, and average yearly precipitation amounts to 745.2 mm. Site 1. The dune is located in Ukraine (geographic coordinates: 50˚15'20.97''N, 24˚19'36.06''E). This is a parabolic type dune with a distinct shape. Profile 1 lies on the dune's peak. The top is covered in a coniferous forest dominated by pine-Pinus sylvestris L. with an admixture of oak-Quercus robur L. and hornbeam-Carpinus betulus L.

Soil Sampling
Soil samples were taken from every genetic horizon of the studied profiles (Table 1). Each soil sample weighed about 1 kg. Every sample was placed in a waterproof container to maintain natural humidity levels. consumption-OCD or PCO 2 ). The respective regression coefficients were all r > 0.95. They were later corrected for pressure and temperature, multiplied by the calibration factor and related per gram of dry soil per hour. The method was described in details by Fischer and Blažka [22]. Measurements were repeated 12 times for each soil sample.
2) Humidity (GWC-gravimetric water content, WHC-water holding capacity, dm-dry mass). Our procedure was essentially based on Grace et al. [23], Priha and Smolander [24] with minor modifications. GWC and other contents are given as ratios, not in %. Dry mass was determined by drying at 110˚C to constant weight.
3) Soil Organic Matter (SOM). SOM was estimated as the loss of weight in ignition to constant weight (LOI) at 500˚C on selected samples. 4) Granulomentric composition. Particle-size analysis was done with the laser particle measuring device Analysette 22 Comfort. Particle-size measurements were done with a step of 0.5 phi.

Statistics
Means (m), standard deviations (s), coefficients of variation (cv) and standard errors of the mean (se) were determined for all data series. A T-test for equality of means was used to determine the statistical significance of mean differences.
A test for regression co-linearity determination was also carried out.

Results
SOM, GWC, WHC and granulometric composition were determined for six horizons in each of profiles 1, 2 and 3. The characteristics of these profiles are presented in Table 1. Respiromentric investigations PCO2, OCD, RQ were carried out for all horizons of profiles 1 and 2.

Soil Profile Characteristics
The investigated profiles were classified as Podzol soil. Six horizons were demarcated in each of the profiles (O/A-A-Es-Bs-B/C-C). In the horizons of the investigated soil profiles, fine or medium grain, poorly or medium sorted sands dominated ( Table 2). In the profiles, the granulomentric fractions were divided into three main groups-sand, powder and silt. The percentage contents of these groups were analysed in each horizon except in the O/A-A horizon, in view of its high organic matter content. The type of relation was similar for all three soil profiles. As expected, the amount of sand increased with depth. Meanwhile, the proportion of powder and silt decreased. However, it is worth noting that in the investigated soils an increase in the proportion of the finest particles was not observed in horizon Bs. It can be assumed that this is due to the location of the profiles on the peaks of dunes, where are the best interior drainage conditions. Open Journal of Soil Science The highest organic matter contents (SOM) ( Table 1)

Humidity, Soil Organic Matter Content and Respiration in the Soil Profile
Both the oxygen consumption and carbon dioxide production decreased exponentially with profile depth (Figure 1).  An analysis of regression co-linearity was carried out to clarify this problem.
The results showed there was a closer relationship between WHC and SOM than between WHC and GWC (Table 3). This relation indicates the close link between organic matter and water holding capacity, which has been found before in numerous studies [25] [26] and others. An analysis of regression co-linearity was also carried out to analyse the correlation between the granulometric per-  given horizon. To obtain data on how respiration changes in response to changes in soil humidity, each horizon was hydrated up to 80% of its maximal water capacity. 80% of water capacity is generally considered to be optimal. The obtained results are presented in Figure 3; the calculations were based on percentage values, assuming that respiration at the moment of sampling was 100%.
These results show that oxygen consumption during moistening increases more than carbon dioxide production. This could be due to the addition of extra air with the water. However, the deeper horizons (Bs, B/C, C) reacted less to moistening than the other horizons. Using microbiological methods, Fierer and his team concluded that "The relative responses of C mineralization rates to changes in soil moisture, temperature, and nutrient levels varied with soil depth" [6].
Our results showing the effect of increased humidity on soil respiration demonstrate that the surface and sub-surface horizons responded differently. Moreover, there were further differences between the subsurface horizons. It is likely that the poor organic matter content in horizons Bs, B/C, C had some significance. These results fully confirm the conclusion of Fierer [6], that there exists a different carbon cycling mechanism in the deeper soil horizons than in the surface and subsurface horizons.
The parameter RQ can indirectly tell us about this mechanism. In research on soil metabolism this parameter is very important, but usually does not give definitive answers on what is dominating-respiration or fermentation; it only gives information about tendencies. Angert [11] showed how for soils the parameter RQ is difficult to interpret, as it is dependent on many factors and varies widely. However, by assuming a simplification it is possible to see some tendencies; namely, if the RQ parameter varies between 0.5 and 1.0 then respiration (aerobic) processes dominate, and if it is smaller than 0.5 or larger than 1.0 then fermentation processes dominate. The data on oxygen consumption and carbon dioxide production for each horizon of profile 2 was used to calculate the RQ parameter. The results were divided into 2 groups. The first included results higher than 1.0 and lower than 0.5, for which fermentation processes were assumed to dominate, and the second results from 0.5 to 1.0, for which respiration processes were assumed to dominate (Figure 4). In horizon B/C respiration was at a lower level that could be detected using our method, so it is not included in

Respiration in the Soil Profile
Although oxygen consumption and carbon dioxide excretion in the soil surface  Table 4. Bearing in mind that these results only have an approximate character, we can deduce that the analysed profiles were similar in their granulometric composition, relief and soil organic matter content but different with regards to their metabolic processes. Oxygen consumption and carbon dioxide production in profile 2 were most intensive in the surface stratum (77% and 72% Table 4).
On the other hand, in profile 1 respiration in the surface stratum (55% and 51%) was close to that of the fermentation stratum (44% and 48%). This resulted in profile 1 producing slightly more carbon dioxide and consuming much more oxygen. This is likely because although the soil organic matter content in both profiles is similar, its distribution within the profiles is different ( Table 4). The fermentation stratum of both profiles contains more than 90% of the whole soil mass; however, the respiration stratum of profile 1 had only about 22% of the organic matter, whereas that of profile 2 had 60% ( Table 4). The differing vertical distributions of the organic matter are perhaps caused by differences in the filtration abilities of the A horizons, which depend on granulometries. The average particle size ( Different soil filtration abilities led to differences in organic matter decomposition that caused the amount of carbon dioxide excretion in the respiration stratum of profile 1 to be 55%, whereas in profile 2% -77%. Davidson et al. [27] cautiously estimated that in their studied mixed forest in Massachusetts, profile O was the source of 40% -48% of annual carbon dioxide production. On the other hand, investigations of boreal forest in Finland showed that horizon A was responsible for 69.9% of carbon dioxide production [14]. However, it is clear  [4]. In our investigations, which were just a snapshot and did not look at multiple seasons, and which only described production in the studied profiles and not excretion from the soil surface, the carbon dioxide production in the surface strata was quite high. We obtained similar results for oxygen consumption. The transfer of part of the soil organic matter in profile 1 to the fermentation stratum led to this stratum having a much higher oxygen consumption and carbon dioxide excretion. The interrelations between organic matter and CO 2 excretion in a soil profile are discussed in, amongst others, Delogu et al. [1]. Finally, our results indicate that the distribution of organic matter in a soil profile significantly affects the intensity of carbon cycling. An amount of organic matter located in the surface respiration stratum leads to a noticeably higher amount of CO 2 production and oxygen consumption than an equivalent amount located in deeper strata. However, this does not mean that in investigations of carbon cycling in an ecosystem the fermentation stratum can be left out, as its respirometric processes can comprise up to even half of all processes in a profile. Organic matter decomposition in the profile depends to a large extent on the relief, soil filtration ability and its location.

Discussion
Carbon cycling in ecosystems not only depends on the soil surface stratum but also on deeper strata, with the distribution of the organic matter in the profile playing a very important role. Its distribution can vary widely and depends on many factors like climatic, geologic, and geomorphologic conditions [28] [29] [30]. Since the rise of agriculture, the role of organic matter in the soil has been of key interest. Of particular importance is that soil organic matter is the source of nutrients for both plants and microorganisms. However, over the last years, the organic matter has also been considered to be a substratum with water holding ability. F. J. Veihmeyer was the first scientist to raise this point, at a meeting of the American Society of Agronomy in 1937 [31]. But over the next years, this problem was somehow forgotten. Recently, as water has become increasingly important [32], water holding capacity has become a hot topic all over the world [33] [34] [26] et al. Hudson [31], in a comprehensive paper, extensively debated the role of soil organic matter in water holding capacity, and discussed how widely the results of researchers varied. He assumed that this variation was due to bad research design, and excessive differences in soil texture and insufficient differences in the amounts of organic matter (between studies). In our investigations, texture differentiation (granulometric composition) of the soil horizons was not significantly related to soil organic matter content or water capacity. Instead, the dependence of WHC on SOM was quite significant (R 2 = 0.92) in each tested soil profile. The investigations of Hudson [31] also found that soils with high organic matter content have a much higher water holding ability than soils with a similar texture but lower organic matter content. Hudson suggested that high organic matter content may be due to the hydration of the organic substance by the creation of so-called "water bridges" [25]. In this way, soil organic matter has to be considered not only a nutrient but also an element that widely regulates humidity.
In light of the variation in humidity and organic matter distribution, which are dependent on the type of soil and its location in the relief and ecosystem, it is necessary to take into account the depth of both a soil profile's respiration and fermentation strata during carbon cycle modelling.

Conclusions
Our research on the soil profiles of inland dune peaks led us to conclude the following: • The main factor determining levels of carbon dioxide production in the soil profile was organic matter content, the location of which in the profile influenced carbon cycling. • Oxygen consumption and carbon dioxide excretion decreased exponentially with depth. Oxygen consumption decreased with depth less rapidly than carbon dioxide production.
• Consecutive horizons, from O/A to Es, had statistically different organic matter contents and water holding capacities, whereas horizons B/C and C did not have such differences.
• The effect of increasing humidity on oxygen consumption and carbon dioxide excretion was different for every soil horizon, but the effect was less pronounced in the deeper horizons than in the surface strata.
• In terms of its respiration processes, every profile could be divided into two strata, the surface horizon, where aerobic respiration dominates, and the other horizons, where fermentation processes dominate. In assessing the soil's role in carbon cycling in an ecosystem, it is necessary to consider both the respiratory and fermentation strata because both produce carbon dioxide in considerable quantities.

Conflicts of Interest
The authors declare no conflicts of interest.