Properties of humic acids depending on the land use in different parts of Slovakia

Many studies report organic carbon stabilization by clay minerals, but the effects of land use and soil type on the properties of humic acids (HAs) are missing. The aim of the paper is to determine the effects of land use and soil types on the characteristics of HAs, which have a considerable influence on organic matter quality. It was hypothesised that the effect of the land use on HAs properties depends on the particular size distribution. The research was performed in three ecosystems: agricultural, forest, and meadow, located in Slovakia. From each of them, the samples of 4 soil types were taken: Chernozem, Luvisol, Planosol, and Cambisol. The soil samples were assayed for the content of total organic carbon (TOC) and the particle size distribution. HAs were extracted with the Schnitzer method and analysed for the elemental composition, spectrometric parameters in the UV-VIS range, and hydrophilic and hydrophobic properties, and the infrared spectra were produced. The research results have shown that the properties of HAs can be modified by the land use and the scope and that the direction of changes depends on the soil type. The HAs of Chernozem and Luvisol in the agri-ecosystem were identified with a higher “degree of maturity”, as reflected by atomic ratios (H/C, O/C, O/H), absorbance coefficients, and the FT-IR spectra, as compared with the HAs of the meadow and forest ecosystem. However, as for the HAs of Cambisol, a higher “degree of maturity” was demonstrated for the meadow ecosystem, as compared with the HAs of the agri- and forest ecosystem. The present research has clearly identified that the content of clay is the factor determining the HAs properties. Soils with a higher content of the clay fraction contain HAs with a higher “degree of maturity”. Supplementary Information The online version contains supplementary material available at 10.1007/s11356-021-14616-9.


Introduction
The role of organic matter (OM), including humic substances (HS), in the soil formation and development is unquestionable. OM is also an important indicator of soil fertility, and it plays an essential role in the efficiency of agriculture, the quality of the environment and global climate changes (Zech et al. 1997;Lal 2004;Schmidt et al. 2011;Tan 2014).
The content and properties of organic matter are conditioned by the physicochemical properties of the parent rock and climate conditions. However, the land use change is the most dynamic factor of soil organic matter (SOC) changes (Guo and Gifford 2002;Poeplau et al. 2011;Jonczak 2013;Viscarra-Rossel et al. 2014;Kukuls et al. 2019;Kunlanit et al. 2019). In the temperate climate zone, the capacity of soil for SOC storage increases in the following order: cropland < forest < grassland (Meersmans et al. 2008;Martin et al. 2011). Leifeld et al. (2005) and Lettens et al. (2005) claim that the forest SOC stocks can be higher than grassland stocks, and a change in use from forest to farmland results in significant declines in SOC stocks (Guo and Gifford 2002;Murty et al. 2002;Wei et al. 2014;Kunlanit et al. 2019). Wei et al. (2014) report on the largest decrease in the SOC stock in topsoil (0-30 cm) in temperate regions -52% followed by tropical regions -41% and boreal regions -31% (the values are the average decrease in the SOC stock for 3 periods -≤10, 11-50, > 50 years which lapsed since the woodland was transformed into arable land).
One of the key SOC-stabilizing mechanisms is the soil interaction with mineral particles (Sollins et al. 1996;Von Lutzow et al. 2006). Due to a strong correlation between the SOC stocks and the content of clay presented in numerous research reports, the soil texture can be a promising factor applied as the SOC storage indicator (Hassink et al. 1997;Kaiser and Guggenberger 2000;Arrouays et al. 2006;Zinn et al. 2007).
As already mentioned, one of the most important organic carbon (OC) reservoirs is the soils of grasslands and forest soils. For that reason, the meadow and forest ecosystems are a land use preventing OC losses in soils. The OC losses triggered by soil processes are gas emissions of CO 2 and CH 4 as well as soluble organic carbon leaching to the groundwaters. It has been demonstrated that meadow ecosystems help limiting the OC losses from soil, and they are essential for the sustainable management of that element in the environment (Minami et al. 1993;Conant et al. 2001;Mannetje 2002;Lal 2011;Kampf et al. 2016).
Many reports point to progressive decreases in HS concentrations in soils converted from forest to arable land (Spaccini et al. 2006;Gonzalez-Perez et al. 2007;Guimaraes et al. 2013; Barancikova et al. 2016;Kunlanit et al. 2019), and as claimed by Feller and Beare (1997), Watanabe et al. (2001), and Jonczak (2013), the changes in land use can alter the chemical properties of soil HS. HS (fulvic acids, HAs, and humins) are the main fractions of SOM making its specific contribution to soil fertility. HS constitute a considerable part of the resources of OC and nitrogen (N) (Lal 1994;Milori et al. 2002;Zhang et al. 2019). Those compounds take part in all the processes which occur in soil and which affect its physical, chemical, and biological properties. HS improve the soil buffering capacity, supply the plants with available micronutrients, and immobilize organic contaminants and metals (McCarthy 2001;Yamashita et al. 2008;Canellas et al. 2010;Lanyi 2010). HS also determine the soil production potential, and by participating in the global carbon cycling, they play environmental functions (Hayes and Clapp 2001;Piccolo 2001;Lal 2006). For that reason, it is important to understand the nature, composition, and the dynamics of HS.
The molecules of HAs, with a humification progress, show a growing stability, e.g., in terms of elemental composition and spectroscopic properties (Dergacheva et al. 2012;Zhang et al. 2019). For that reason, one of the basic indicators applied to evaluate the properties of HAs is the elemental composition and, determined with it, values of the H/C, O/C, and N/C atomic ratios. The numerical values of atomic ratios facilitate an approximation of the structure of the molecules of humic acids by evaluating the degree of condensation of aromatic rings (H/C ratio) and the degree of maturity (O/C, O/H, N/C ratios) (Rice and MacCarthy 1991;Fuentes et al. 2007;Canellas et al. 2010;Trubetskaya et al. 2013;Boguta et al. 2016).
An important criterion characterizing the HAs in terms of the molecular composition and their origin is the values of absorbance of their solutions in the UV-VIS range: A 280 , A 465 , A 665 , and the coefficients of absorbance A 2/4 , A 2/6 , A 4/6 , and ΔlogK. Coefficients A 2/4 , A 2/6 , A 4/6 , and ΔlogK are important indices of the degree of advancement of the humification of organic materials and the characteristics of the HS produced, as well as changes in the properties of the HAs which occur due to various anthropogenic factors (Kumada 1987;Chin et al. 1994;Tan 1998;Gonet and Debska 1999;Chen et al. 2002;Weishaar et al. 2003;Moran Vieyra et al. 2009;Polak et al. 2011;Rodriguez et al. 2016). For example, the research of the spectrometric parameters of the HAs of soils collected from primary forests, secondary forests, coffee plantations, and cultivated lands performed by Watanabe et al. (2001) suggests that the degree of HAs humification varied from site to site.
Infrared spectroscopy (FTIR) has also been a technique used for the structural characterization of HS, especially for the identification of functional groups within the humus macromolecule (Tan 1998;Cocozza and Miano 2002;Pajaczkowska et al. 2003;Polak et al. 2011;Kukuls et al. 2019).
Some essential information on the transformation of HAs is also provided by the results of high-performance liquid chromatography. Woelki et al. (1997), Preuße et al. (2000), Banach-Szott and Debska (2006), Sierra et al. (2006), and Debska et al. (2010), using the chromatographic analysis, separated the hydrophilic (HIL) and hydrophobic fractions (HOB-1 and HOB-2). The ratios of those fractions affect the solubility of HAs and, as a result, their migration deep down the soil profile. According to  and Debska and Gonet (2007), with an increase in the degree of humification, the share of the HIL fraction in the molecules of HAs increases and the share of HOB fractions decreases. As a result, the HAs molecules with a higher "degree of maturity" showed a higher value of the HIL/ΣHOB ratio.
The aim of the paper has been to present the land use and soil types effect on the characteristics of HAs determining the OM quality in the soil of various types. There are many results on the processes of OC stabilization by binding clay; however, some results on the changes in elemental composition and the other properties of HAs depending on the land use and soil types are missing. It was hypothesised that the effect of the land use on the properties of HAs will depend on the particle size distribution (soil types). The research has been performed in three ecosystems: agricultural, forest, and meadow, located in various parts of Slovakia. The properties of HAs have been determined from the assay of the elemental composition, spectroscopy in the UV-VIS, and IR ranges and the HIL-HOB properties applying the high-performance liquid chromatography.
The areas are located in different parts of Slovakia (Fig. 1). Chernozem (Ch) comes from the location of Voderady (48°16′N, 17°34′E), found on the northern border of the Danube Basin. The geological structure shows the Neogene strata, which consist mainly of claystones, sandstones, and andesites, covered with younger quaternary rocks represented by different fluvial and Aeolian sediments (Sajgalik et al. 1986). It is located in a warm climate region with an average annual temperature of 9.6°C and the total annual precipitation of 560 mm (Korec et al. 1997).
Planosol (Pl) comes from the location of Vavrecka (49°38′ N, 19°47′E), found in Horna Orava, which is part of the Flysch Belt. The geological structure shows the alternation of clay slates, sandstones, and conglomerates in the layers of various depth. It is located in a cold climate region with an average annual temperature of 4.6°C and the total annual precipitation of 1010 mm (Sajgalik et al. 1986).
Cambisol (Ca) comes from the location of Prietrz (48°40′ N, 17°26′E), on the Myjava Hills. The geological structure includes the Carpathian flysch, which consists mainly of marine deposits of claystones, shales, and sandstones (Sajgalik et al. 1986). It is located in a temperate climate region with an average annual temperature of 8.0°C and the total annual precipitation of 681 mm (Korec et al. 1997).
Luvisol (Lu) comes from the location of Plave Vozokany (48°06′N, 18°46′E), in the northeast of the Danubian Upland. The geological structure involves clay, gravels, and sands covered with quaternary sediments (Sajgalik et al. 1986). It is located in a warm climate region with an average annual temperature of 10.4°C and the total annual precipitation of 589 mm (Korec et al. 1997).
The agri-ecosystem included four crop rotations for each soil and land use ( Table 1). The fields in the agri-ecosystems were located on different farms under real production conditions. The forest ecosystems were natural forests with human control, in all the cases they were hundreds of years of age. The dominant species in the stands were on Chernozem -Quercus, on Luvisol -Alnus, on Planosol -Picea, and on Cambisol -Fagus, Quercus, and Carpinus. The meadow ecosystems were created by man at least 30 years ago.

Methods
The soil samples were collected with the Egner stick in three replications down to the depth of 0.30 m in the forest and meadow ecosystems. In the agri-ecosystems, collections were made in different fields with different crop rotations and different applications of farmyard manure (in 12 replications, the depth of 0.30 m). The distance between the replications was 20 m in equilateral triangle. The distance between fields differed, depending on the locations; however, it was not more than 2 km. The distance between the locations ( Fig. 1) was not more than 200 km. The sample replications for chemical analyses were combined and carefully mixed, dried in room temperature, and sieved (2 mm).

Basic parameters of soils
The particle size distribution was determined after dissolving CaCO 3 with 2 M HCl and the oxidation of OM-with 30% H 2 O 2 . Silt, sand, and clay fractions were assayed applying the pipette method (Van Reeuwijk 2002). In the soil samples, the TOC was assayed applying wet combustion (Orlov and

Extraction of humic acids
HAs were extracted and purified according to the following procedure (Debska et al. 2010): -Decalcification (24 h) with 0.05 M HCl (1:10 w/v). After centrifugation the residue was rinsed with distilled water till neutral. -The extraction (24 h) of the remaining solid with 0.5 M NaOH (1:10 w/v), with occasional mixing, followed by centrifugation. -The precipitation (24 h) of HAs from the resulting alkaline extract with 2 M HCl to pH=2 and centrifugation. -The purification of the resulting HAs; the HAs residue was treated with a mixture of HCl/HF (950 mL H 2 O, 5 mL HCl, 5 mL HF) over a 24-h period, followed by centrifugation. This procedure was repeated three times. The HAs residue was treated with distilled water until a zero reaction to chloride was achieved.
The preparations were lyophilised and powdered in the agate mortar. The ash content in the HAs preparations was lower than 2%.

Characteristics of humic acids
The HAs separated were analysed for: -Elemental composition (Perkin Elmer, series II CHN analyser, Shelton, USA). The H/C, O/C, O/H, and N/C atomic ratios and ω (internal oxidation degree) were calculated using the following formula: where O, N, H, C-content in atomic %.

Statistical analysis
The significance of differences of the parameters between land uses (agri-ecosystem (AE), forest (FE), and meadow ecosystems (ME)) within a soil type was evaluated with Duncan's post hoc test at p < 0.05. The effect of the soil type and the land use on the properties of HAs was defined with cluster analysis. The method involves dividing the data set into groups to produce clusters in which the elements are similar to one another and, at the same time, different from the elements from the other groups. The groups of similar treatments are presented in a form of dendrogram. In a given group, the smaller the Euclidean distance, the more similar the objects. Data clustering was performed with the Ward method. The analysis was made after data standardization. The cluster analysis was performed based on the elemental composition (C, H, N, O, H/C, O/H, O/C, N/C, ω), spectrometric parameters (A 280 , A 465 , A 665 , A 2/4 , A 2/6 , A 4/6 , and ΔlogK), as well as the HIL-HOB properties. The method involves dividing the data set into groups to produce clusters where the elements are similar to one another and, at the same time, different from the elements of the other groups. The relationships between the share of clay, silt, and sand fraction and the basic parameters of HAs were defined using the Pearson's correlation coefficients (P ≤ 0.05). The above relationships were determined using statistics software STATISTICA MS 12.

Basic parameters of soils
One of the key factors controlling the rate of the SOM circulation is field aggregates produced by combining clay particles with organic particles (mineral-organic particles). The soils differed significantly in terms of the particle size distribution ( Table 2). The highest content of silt and clay fractions and the lowest of the sand fraction were recorded for Luvisol. Chernozem showed a similar content of sand and silt fractions (34.9 and 36.6%, respectively). Cambisol demonstrated a similar content of sand and clay fraction (23.0, 24.5%, respectively) and Planosol a similar content of silt and clay. The lowest content of TOC was reported for Cambisol and the highest for Planosol. The analysis of correlation did not show a significant dependence between the content of clay and the content of TOC; however, as reported in literature, one cannot exclude the effect of clay on the stability of SOM (Wiesmeier et al. 2019). As evident from the literature reports (Six et al. 2002;Wiesmeier et al. 2015), it is not just the quantity of the fine fraction but also its quality which drives SOC retention in soils. Liang et al. (2009) andWiesmeier et al. (2015) show that the correlation between SOC and clay strongly depends on climate conditions, land use, and clay type. However, as reported in literature (Wiesmeier et al. 2019), despite a lack of significant correlations, one cannot exclude the effect of clay on SOM stability.

Properties of humic acids
According to De Moraes et al. (2011), the changes in land use practices can alter the chemical properties of soil HS. In this paper, the changes in land use practices are clearly reflected in the elemental composition (in atomic %) of HAs (Table 3). The content of C in the molecules of HAs ranged from 33.38 to 39.55%, H from 35.57 to 43.81%, N from 2.55 to 4.21%, and O accounted for 18.10-22.80%. The highest share of C was recorded for the molecules of the HAs of Luvisol sampled from the forest ecosystem. The lowest share of that element as well as N was found for the HAs of Planosol (FE). The lowest content of H and the highest content of O were recorded in the HAs of Chernozem from the agri-ecosystem (AE). Table 3 also presents the values of the C, H, O, and N atomic ratios. With the values of H/C atomic ratio, one can determine the degree of condensation of aromatic rings and, as a result, the "degree of maturity" of the molecules of HAs. The values of the H/C ratio for all the variants ranged from 0.91 to 1.31 (Table 3), which shows that the HAs analysed contained the aromatic systems coupled with the aliphatic chain containing up to 10 atoms of C (Van Krevelen 1950). The parameter used to describe the advancement of the process of humification is also the degree of internal oxidation (ω) of the HAs molecules and the O/C, O/H, and N/C ratios. Higher ω, O/C, and O/H values and lower H/C values correspond to the HAs with a higher "degree of maturity" (Sanchez-Monedero et al. 2002;Gonet et al. 2007). Dergacheva et al. (2012) stress that the value of the H/C ratio depends on the conditions the HAs originated from. As seen from the data presented (Table 3)   Of the HAs of Luvisol the highest content of C and the lowest of H were noted for the HAs of the forest soil, whereas the HAs of the arable soil showed the highest value of ω and the O/C ratio. The HAs quality parameters in meadow soil, in general, did not demonstrate any significant differences, as compared with the molecules of the HAs of forest ecosystems.
Of the HAs of Planosol, the highest content of C and the lowest content of H and, as a consequence, the lowest value of the H/C ratio and the highest values of O/H and ω were recorded for the HAs of meadow soil. The HAs of the agriecosystem showed the lowest value of the O/C ratio.
As for the HAs of Cambisol, the lowest content of C and O and the highest content of H and the highest value of the H/C ratio and the lowest values of the O/H ratio and ω were noted for the HAs of the arable soil. Interestingly, the significant differences were reported in the value of the O/H ratio and parameter ω between the HAs of the soils of all the ecosystems investigated.
The impact of the land use on the properties of HAs is confirmed by the diagram of the relationships between the values of the H/C atomic ratio and parameter ω (Fig. 2). The molecules of the HAs of meadow soils showed similar values irrespective of the soil type. As for the HAs of the other ecosystems, that soil type was an additional factor determining their properties (e.g. HAs of Luvisol). Generally, the HAs of Luvisol showed the highest contents of C and lower contents of H, as compared with the HAs of the other soils (Table 3). The highest C content in the HAs of Luvisol of the forest ecosystem resulted in a low content of the H/C ratio and a high value of ω, as compared with the HAs of the forest ecosystem, of the other soil types.
The existing research results have indicated some dependencies between the values of absorbance at the wavelength of 280, 465, and 665 nm and the structure of HAs. The lowest values of absorbance coefficients were found for the HAs of Luvisol irrespective of the land use and for the molecules of the HAs of Chernozem sampled from the agri-ecosystem (Table 4). As reported by, e.g., Kumada (1987), Howard et al. (1998), Tinoco et al. (2015, and Filcheva et al. (2018), lower values of absorbance and higher values of the following coefficients, A 2/4 , A 2/6 , A 4/6 , and ΔlogK, point to a chemical "young age" of HAs. Young humic acids show a lower degree of condensation of aromatic structures and a lower molecular weight, as compared with the HAs with a high degree of humification. As results from the data presented in Table 4, the HAs differed in the degree of maturity. Of all the HAs of Luvisol, higher values of the coefficients of absorbance were recorded for the HAs of meadow soil, as compared with the HAs of forest soil and under agricultural use. Similarly, for the HAs of meadow soil of Planosol, there were found higher values of A 2/6 and ΔlogK, as compared with the HAs of the agri-ecosystems. However, the HAs of Cambisol under meadow ecosystem use showed, in general, lower values of the coefficients of absorbance as compared with the HAs of the other variants. The values of the coefficients of absorbance recorded for the HAs of Cambisols confirm the dependencies reported by Jonczak (2013), who for the HAs of soils under different land uses (forest, meadow, arable field, fallow, postarable afforestation) noted a higher maturity of humus in stands with grass vegetation (meadow, fallow, afforestation with birch with dense grass cover in forest floor) in relation to other stands. Watanabe et al. (2001) have demonstrated that the relationship between the land use and the degree of HAs humification differed depending on the soil sampling site. On the one hand, the degree of humification of HAs in the topsoil was greatest under secondary forest and least under coffee plantation, whereas for the soil sampled from another site, the degree of humification of HAs was greater under agriculture than under forests. According to the authors, the increase in the degree of humification with decreasing amount of HAs due to changing land use suggests selective decomposition of the molecules or moieties of humic acids with low degrees of humification.
HAs are built from structures with both hydrophobic and hydrophilic properties. By applying the HPLC method, one can divide the molecules of HAs into fractions which are HIL and HOB in nature (Fig. 3). The fractions of the retention time of 4.0-7.0 min show greater HIL properties, whereas the fractions of the retention time between 14.0 and 25.0 min become more and more HOB (Woelki et al. 1997;Preuße et al. 2000;Debska et al. , 2010Debska et al. , 2012. The share of HIL fractions in the HAs molecules of the soil analysed was lower than the share of HOB fractions, and it ranged from 15.90 (Lu FE) to 29.46% (Ch FE) (Table 5). With the patterns of chromatograms for the HOB fractions, for most HAs, there were separated 3 smaller fractions: HOB-1, HOB-2, and HOB-3 (Fig.  3). Irrespective of the factors discussed (land use, soil type), of those hydrophobic fractions, in general, the HAs recorded the highest share of the HOB-2 fraction. It should be noted that the biggest differences between the shares of fractions HOB- Fig. 2 Relationship between the H/C atomic ratio values and parameter ω, where H, Ccontent in atomic %, ωinternal oxidation degree, AEagri-ecosystem, FEforest ecosystem, MEmeadow ecosystem 1, HOB-2, and HOB-3 (Table 5) were identified in the molecules of the HAs of the meadow ecosystem (6.22-8.94 pp). The highest share of the ΣHOB fraction was recorded for the HAs of the forest and agri-ecosystem of Luvisol and the lowest for the HAs of the forest and meadow ecosystem of Chernozem. Of the Chernozem HAs, the highest share of ΣHOB was recorded for the HAs of the agri-ecosystem and of the Planosol HAs and the Cambisol HAs in the meadow ecosystem. The results point to a strong effect of the soil type on the ΣHOB properties of HAs. The changes in the share of the respective fractions resulted in changes in the value of the HIL/ΣHOB ratio. The HAs isolated from the Luvisol, except for ME, showed the lowest and the HAs extracted from the Chernozem, except for AE, the highest values of the HIL/ ΣHOB ratio. The values of the HIL/ΣHOB ratio were similar for the HAs of Cambisol of the agri-and forest ecosystem and for the HAs of Chernozem of the forest and meadow ecosystem. For the HAs of Chernozem, no significant differences were found between the values of HIL/ΣHOB in the forest and meadow ecosystem and as for the HAs of Planosol between the agri-and meadow ecosystem. The share of HIL and HOB fractions and parameter HIL/ΣHOB is connected with the degree of OM humification. The values of HIL/ΣHOB of the parameter increase with an increase in the degree of maturity of HAs molecules (Woelki et al. 1997;Preuße et al. 2000;Debska et al. , 2010Debska et al. , 2012. It should also be emphasized that the proportions of both fractions determine the solubility of HAs and, as a result, their migration deep down the soil profile (Woelki et al. 1997;. Figure 4 provides sample FTIR spectra of the HAs. All the spectra were identified with a presence of the same absorption bands the ranges of which are presented in Table 6. As a result of a detailed spectral analysis, it was found that the intensity of the absorption bands in the molecules of HAs depended on the land use. The HAs of Chernozem and Luvisol under agricultural use showed a higher band intensity in the range of 1730-1710 cm −1 , as compared with the molecules of the HAs of the forest and meadow ecosystem. The intensity of the  To acquire complete information on the differences (similarities) in the chemical composition of HAs, depending on the soil type and the land use, the cluster analysis was applied based on the elemental composition, spectrometric parameters, as well as the HIL-HOB properties dividing the HAs into two groups (Fig. 5). In the first one, two subgroups were identified. The first subgroup included mostly the HAs of the meadow ecosystem of the following soil types: nd non-detectable (below detection limit) Fig. 4 FT-IR spectra of humic acids of Planosol (Pl) showing the dependence of transmittance (T) on the wavenumber; a in the range from 3600 to 2500 cm −1 , b in the range from 1800 to 1000 cm −1 (AE -agri-ecosystem, FEforest ecosystem, MEmeadow ecosystem). Assignment of peaks (bands), see Table 5 Planosol and Cambisol as well as the HAs of the forest ecosystem of Cambisol and HAs Planosol under agricultural use.
In the second subgroup, the most similar properties were found for the molecules of the HAs of Chernozem in the forest and meadow ecosystem. In the second group, similar properties were noted for the HAs isolated from Luvisol, irrespective of the ecosystem, and the HAs extracted from Chernozem sampled from the agri-ecosystem. As seen from the diagram, both the soil type and the land use are the factor which determines the properties of HAs. With some approximation, one can assume that Luvisol is the type of the soil which is modified even least considerably by the land use.  Table 7 presents the correlation coefficients (significant for p ≤ 0.05) between the particle size distribution and the properties of HAs. The HAs with a higher "degree of maturity" (a negative significant correlation between the clay content and content of H in HAs, the values of H/C, A 2/4 , A 2/6 , A 4/6 and positive between the content of C and parameter ω and the values of absorbance) are characteristic for the soils with a higher clay fraction content and a lower sand fraction content. Even earlier Jindaluang et al. (2013) noted the effect of the particle size composition, especially the content of the clay fraction, on the content and the properties of SOM. As seen from the research, the particle size distribution, especially the content of the clay fraction, is the factor affecting also the properties of HAs, mainly their elemental composition and spectrometric properties in the UV-VIS range.

Conclusions
From H/C, O/C, O/H, and ω parameters and absorbance coefficients and the FT-IR spectra, it has been found that the HAs of Chernozem and Luvisol in the agri-ecosystem show a higher "degree of maturity", as compared with the HAs of the meadow and forest ecosystem. However, for the HAs of Cambisol, a higher "degree of maturity" was demonstrated for the meadow ecosystem, as compared with the HAs of the agriand forest ecosystem.
The research results have demonstrated unambiguously that the properties of HAs can be modified by the land use and the scope and that the direction of changes depends on the soil type. The research has identified that the content of clay affects the properties of HAs. Soils with a higher content of clay fraction include HAs with a higher "degree of maturity". It is seen from the positive (significant) values of  the coefficient of correlation between the clay fraction content and the content of C in the molecules of HAs, the values of the degree of internal oxidation and absorbance values, and the negative (significant) correlations between the content of the clay fraction and the content of H and the values of the H/C, A 2/4 , A 2/6 , and A 4/6 ratios.
Acknowledgements We thank the Ministry of Science and Higher Education for their support.
Availability of data and materials Data sharing is not applicable to this article. Funding The research has been made as part of BS-3/2017 research project, financed by the Ministry of Science and Higher Education.

Declarations
Ethics approval and consent to participate Not applicable.

Consent for publication
The authors confirm that the final version of the manuscript has been reviewed, approved, and agreed for publication by all authors.

Competing interests The authors declare no competing interests.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.