The Effect of Long-Term Farmyard Manure and Mineral Fertilizer Application on the Increase in Soil Organic Matter Quality of Cambisols

: Soil organic matter (SOM) quantity and quality are important factors that signiﬁcantly inﬂuence soil fertility. SOM quality indicators change throughout time. In this study, long—term ﬁeld experiments (22–50 years in duration) on a Cambisol at four sites in the Czech Republic were selected. Seven crops were successively rotated in the sequence: clover, winter wheat, early potato, winter wheat, spring barley, potato, and spring barley with interseeded clover. Five treatments were investigated, including an unfertilized treatment, farmyard manure, and various combinations of farmyard manure and mineral fertilization. A total of 40 t ha − 1 of farmyard manure was applied to the early potato and potato crops. Combining organic and mineral fertilizers increased soil organic matter quality and quantity over unfertilized or organic only treatment. The highest intensity of mineral fertilizers in our trials elevated the mean of carbon sequestration efﬁciency to 45.6% in comparison to pure manure treatment which reached only 22.9% efﬁciency. A strong correlation was established between the total glomalin content and soil organic matter carbon, fulvic acid, humic acid, carbon hot water extraction, potential wettability index (PWI), and aromaticity index. The PWI was also strongly correlated with these indicators. The E4/E6 ratio indicator was shown to be a much less sensitive method for reﬂecting the change in soil organic matter quality.


Introduction
The quantity and quality of soil organic matter (SOM) constitute fundamental factors for soil fertility.A great deal of attention has been paid to soil organic matter and its stabilization, particularly concerning its role in mitigating climate change.Soil carbon sequestration can be a part of the solution in tempering climate change processes.
The pivotal role of soil organic matter in various soil characteristics, such as sorption properties and soil structural stability, is widely acknowledged.This, in turn, influences hydrological and filtering characteristics, as well as resistance to water erosion.SOM quantity is a very important factor, but even more important are SOM quality indicatorsthese significantly influence soil fertility.SOM quality indicators change throughout time as a result of soil organic matter complexities and the evolution of analytical procedures.
Prior notions of humic substances (C HS ), fulvic acids (C FA ), and humic acids (C HA ) characterized as stable, robust macromolecules with substantial molecular weights [1] have been supplanted by novel concepts.Humic substances were better described as supramolecular associations rather than macromolecular polymers [2].Supramolecular associations can be defined as assemblies of relatively small, chemically diverse organic molecules within the humic structure, linked together by hydrophobic interactions and hydrogen bonds.The accumulation of humus and the creation of microaggregates, characterized by increased hydrophobic traits, result from the concentration of hydrophobic humic molecules in the soil.This, in turn, enhances the stability of soil aggregates.This concept became the base for describing natural humic molecules through sequential fractionation without disrupting the C-C bonds [3].Infrared spectroscopy in the mid-infrared spectrum can identify SOM composition without physical or chemical extraction or separation from soil.Although it has its limits due to the complex nature of the soil, this method can detect some functional groups within soil organic matter and offer insights into relative differences in SOM composition across distinct soil horizons or soil types [4].
Additionally, the potential wettability index (PWI) defines SOM by evaluating the ratio of aliphatic (C-H) to carboxylic (C-O) bonds.Elevated PWI values point towards reduced aggregate wettability [5].A correlation between the PWI and C SOM content has also been noted [6].A robust correlation was observed between the the PWI and the total glomalin content (TG) in a Luvisol [7] but was not confirmed in Chernozem soil type [8].
To further investigate the soil characteristics, the aromaticity index (iAR) can be determined based on the reflectance values of aliphatic and aromatic bands [9].Increased values of the iAR can be the result of the increased formation of aliphatic compounds within aggregates due to soil organic matter mineralization [10].
Conventional methodsfor determining SOM quality indicators are based on extraction principles.The soil organic matter can be categorized into labile and stable fractions.The stable fraction is represented by humic acids (C HA ), fulvic acids (C FA ), and humins [1,11].The labile forms are represented by potentially mineralizable carbon (C HWC ) [12] or dissolved organic carbon (C DOC ) [7,8].
The amount glomalin (GRSP) in the soil can be regarded as a reliable indicator of soil fertility [13,14] due to the positive correlation between glomalin and C SOM content [15].In our previous study [7], we identified a positive correlation between GRSP content and humic acid content (C HA ), as well as a positive correlation between GRSP content and the ratio of the humic acid to fulvic acid (C HA /C FA ) in the context of long-term maize monoculture production on a Luvisol In addition, a positive correlation of glomalin content with the potential wettability index (PWI) was found.Glomalin content in the soil can be increased by through the long-term application of manure [16,17].
The E4/E6 ratio exhibits an inverse correlation with the degree of condensation of the aromatic network in humic substances.A diminished E4/E6 ratio suggests high degree of condensation among aromatic constituents, whereas a heightened ratio signifies smaller degree of aromatic condensation and the prevalence of comparatively large proportions of aliphatic structures [18].Carboxylic groups contribute to a higher amount of acidity than phenolic groups in fulvic and humic acids using the E4/E6 ratio measurements of these fractions [19].The E4/E6 ratio has been recorded as higher in fulvic acids than in humic acids.No substantial relationship between (C SOM ) and the E4/E6 ratio was discerned in long-term experiments involving maize monoculture [8].
This study aims to: (i) assess the alterations in both the quality and quantity of theSOM in Cambisol under the influence of long-term exposure to organic and mineral fertilizers.;(ii) evaluate the viability of the glomalin soil extraction method as an indicator of SOM quality; (iii) evaluate the suitability of the PWI method as indicator of the SOM quality; iv) establish, whether the E4/E6 ratio is sensitive enough to reflect the changes in theSOM quality and quantity.We hypothesize that: (a) long-term fertilization will have a significant effect on the SOM quality and quantity; (b) the glomalin soil extraction method will be a suitable method for determining the SOM quality; (c) the PWI determination will be a suitable method to evaluate the SOM quality; and (d) theE4/E6 ratio will be sensitive to changes in the SOM quality and quantity.
For the purpose of this study, four long-term (22 to 50 years) field experiment sites on Cambisol soils in the Czech Republic were selected.We attempted to verify the influence of long-term farmyard manure application as a practice often mentioned in relation to carbon sequestration in agricultural soils.In temperate climate zones, Cambisols characterized by high base saturation are recognized as some of the most fertile and productive soils on Earth.This soil type has an area of 15 million km 2 worldwide [20].

Materials and Methods
Between 1972 and 2000, the Central Institute for Supervising and Testing in Agriculture established long-term on-farm trials on Cambisols soil.Table 1 outlines the characteristics of these experimental sites.A rotation of seven crops was implemented in the following sequence: clover, winter wheat, early potato, winter wheat, spring barley, potato, and spring barley with interseeded clover.C SOM -soil organic matter carbon content; 1 long-term (30 years) annual average; 2 according to NRCS USDA [21]; * values in this row list the following treatment order Con, F, F+M1, F+M2, F+M3; Con-control, F-farmyard manure, F+M-farmyard manure with increasing amounts of mineral fertilizer.
The experiment follows a randomized complete block design across all sites, with three replications of blocks.Specifics regarding plot sizes on each site are outlined in Table 1.The study investigates five treatments: (1) an unfertilized control (Con); (2) farmyard manure only treatment (F); (3) farmyard and mineral fertilizers (F+M1, 57, 13, 33 kg N, P, K ha −1 , respectively); (4) F+M2 (85, 26, 66 kg N, P, K ha −1 , respectively); and (5) F+M3 (114, 53, 133 kg N, P, K ha −1 , respectively).Plant residues are left unincorporated into the soil.A mass of 40 t ha −1 of farmyard manure is applied to early potato and potato, occurring twice during a crop rotation.Carbon input over the experiment's duration is calculated from a dry matter content of 23.0% and a content of carbon 27.9% in dry matter.Values derived from long-term monitoring on the sites.Mature farmyard cattle manure (at least six months of storage before application) was used.Average N content was 0.48% in fresh matter, C:N ratio was 13.5:1.

Soil Analysis
Soil samples were collected in 2022 after the harvest of spring barley with interseeded clover.The experimental sites were instituted from 1972 to 2000, employing a crop rotation cycle involving seven crops (as detailed earlier).Soil sampling took place post the barley harvest in 2022, marking the conclusion of the crop rotation, utilizing a soil probe at a depth of 30 cm.From every experimental plot fifteen soil samples were collected and combined.Combined samples were subjected to air-drying at 25 • C, followed by homogenization and sieving through a 2 mm sieve.Furthermore, a 0.4 mm soil size fraction was also prepared for the C SOM determination.Soil was analyzed using following procedures: The Soil organic carbon (C SOM ) content was assessed through oxidation using the CNS Analyzer Elementar Vario Macro (Elementar Analysensysteme, Hanau near-Frankfurt am Main, Germany).
The humic substance fractionation (C HS ) followed the procedure outlined by Kononova [11] to acquire the pyrophosphate-extractable fraction, encompassing the cumulative carbon in humic acids (C HA ) and fulvic acids (C FA ).
Humus quality, as indicated by the E4/E6 ratio, was assessed using the spectrophotometric method.Soil samples underwent extraction with sodium pyrophosphate (0.05 M Na4P2O7) and were then measured for an absorbance ratio at 400 and 600 nm [22] (Lambda 25 UV/Vis (Perkin Elmer, Waltham, MA, USA).
Extractable organic carbon was determined using CaCl 2 and hot water extraction.
For CaCl 2 extraction (C DOC ), the extraction was performed according to Houba et al., [23].The C DOC content was determined in soil samples using segmental flow analysis with infrared detection using a Skalarplus System (Skalar, Breda, The Netherlands).
Hot water extraction (C HWC ) was used to assess extractable soil organic carbon [12].The C HWC was determined using a segmental flow analysis with infrared detection using a Skalarplus System (Skalar, Breda, The Netherlands).

PWI = A/B
The aromaticity index was calculated based on the reflectance values of aliphatic bands ranging from 3000 to 2800 cm −1 (AL) and the aromatic band at 1520 cm −1 (AR) [24].
Humification indices were calculated according to Raiesi [25] and Iqbal et al., [26]: where C FA is the fulvic acid carbon, C HA is the humic acid carbon, and C SOM is the total organic carbon in soil.
The carbon sequestration efficiency (CSE) was calculated as follows: C SOMtreatment represents the quantity of carbon (C) present in the soil under the fertilized treatment.C SOMunfert signifies the amount of C in the soil under the unfertilized Con treatment.TCI denotes the total carbon input (t ha −1 ) applied through organic fertilizers throughout the course of individual experiments [27].

Statistical analysis
The obtained results underwent assessment through ANOVA analysis and Pearson's correlation coefficient using the Statistica program ver.12.3 (TIBCO, Palo Alto, CA, USA).One-way ANOVA statistical analysis, complemented by Tukey's test, was employed to evaluate treatment and site effects (p < 0.05).Pearson's correlation coefficients were utilized to scrutinize relationships among the studied variables.A significance level of p < 0.05 or lower was deemed statistically significant.

Results
Statistical evaluation encompassed all the measured values of soil organic matter quality indicators from each treatment, as detailed in Table 2.The impact of site was a major factor in these results and caused statistically significant differences among all observed the SOM quality indicators.To assess the impact of the treatments, the existing variable values were substituted with relative ones.The relative values were computed as follows: V treatment /V site -average, where V treatment is represented the value of each treatment, and Vsite-average denoted the average value of a specific site across all the treatments (Table 3).This allowed us to dispose of site influence and focus purely on the influence of individual treatments (Table 3).Along with the relative values, the average values from all four sites are also shown.
The average content of C SOM was 1.51%.The combination of farmyard manure (F) with mineral fertilizers (F+M) led to a significant increase in soil organic carbon content (C SOM ) compared to the unfertilized treatment (Con).

Carbon Sequestration Efficiency
The carbon sequestration efficiency (CSE) was determined by assessing the difference (C SOM ) between the fertilized and unfertilized treatments relative to the total carbon applied in the manure.The application of balanced mineral fertilizer in conjunction with farmyard manure resulted in an increased soil carbon sequestration-specifically, by 22.9% in the farmyard manure (F) treatment and by 45.6% in the F+M3 treatment.

Carbon Fractions and E4/E6 Ratio
In order to estimate and evaluate the content of C HS , C FA , and C HA , the fractionation of humic substances was performed.The acquired results facilitated the computation of the humification degree (HA), humification rate (HR), and humification index (HI).Notably, there were no significant differences observed among the treatments for the HA, HR, or HI.
Furthermore, the quality C SOM quality was also assessed through hot water extraction (C HWC ) and 0.01 M L −1 of CaCl 2 extraction (C DOC ).The application of intensive mineral fertilizer treatments (F+M2 and, F+M3) resulted in higher C HWC values compared to the unfertilized treatment (Con) and other treatments.To obtain better insight into these results, the C HWC /C SOM and C DOC /C SOM ratios were also calculated.The averages for the C HWC /C SOM and C DOC /C SOM ratios were 0.030 and 0.003, respectively.In light of these findings, CHWC and CDOC can be characterized as potentially mineralizable carbon and easily extractable carbon, respectively.
The optical density ratio of the humic acids (E4/E6) serves as an additional qualitative indicator of humic substances, providing insights into the humification process in the soil.The average E4/E6 ratio value was 6.46.A significant increase in the values of this indicator were observed with the F+M2 treatment.To establish and analyze a relationship between the indicators, Pearson's correlation coefficients were used (Table 4).The C SOM was significantly correlated to the C HS , C FA , C HA , and C HWC .There was no correlation between C SOM and E4/E6 ratio.No relationship (Tables 4 and 5) of the E4/E6 ratio with the C SOM was present.

The Potential Wettability Index
Using the DRIFT spectra, the potential wettability index (PWI) and index of aromaticity (iAR) were determined.Significant differences were recorded between the F + M3 and the Con, the F treatments and their PWIs, and the iAR parameters.A significant relationship between the PWI and the C SOM (Table 4) was established and further confirmed with regression analysis (Table 5).No relationship between the E4/E6 ratio and the PWI or iAR was established.The correlation results indicate a strong relationship of thepotential wettability index and aromaticity index with the C HS , C FA , and C HA, as well as the C WHC .

The Glomalin-Related Soil Protein (GRSP) Content
The assessment of organic matter quality involved the utilization of indicators such as easily extractable glomalin (EEG) content and total glomalin (TG).content There were significant differences in the EEG and TG among the observed treatments.The average EEG/TG ratio was 0.265.The average values for the EEG/C SOM and TG/C SOM were 0.050 and 0.186, respectively.There were no significant differences among the treatments with respect to these indicators.The EEG and TG exhibited strong correlations with the C SOM content, as well as the C FA and C HA content.Both C HA and GRSP content emerged as crucial contributors to the stability of soil organic matter.Notably, a more robust relationship was identified in the correlation with the G content as opposed to the EEG content.Additionally, a significant and strong relationship was established between the PWI, iAR, and GRSP content, encompassing both EEG and TG contents.Regression analyses further investigated these relationships.A significant influence of the C SOM , C HS , and C HA over the values of the EEG and TG contents was discovered (Table 5).

Discussion
The focus of this study was the Cambisol soil type of the Czech Republic.The combination of the location of the Czech Republic incentral Europe and the fact, that Cambisol soil covers a vast area of the world [20] may provide useful insight for similar conditions.Older sources in the literature describe Cambisols as a soil type with a low proportion of humin acids with aromatic compounds and high fulvic acid content; in addition, these soils have been described as containing more aliphatic compounds [1,28].
Cambisols typically have a lower organic matter quality index and decomposition index in comparison to Chernozems and Luvisols [4].However, the stability of soil organic matter is influenced not only by the composition of organic matter, its aromaticity, and nitrogen content but also by other factors.The presence of iron oxide is an important factor, with polysaccharides and proteins selectively preserved in organo-metallic complexes [29,30].Other studies showed that microbiologically originated biomolecules that are rich in nitrogen, such as amino acids, are also selectively bound to mineral surfaces [31].In Cambisols, the stabilization of soil organic matter (SOM) occurs not solely through the formation of complexes with the mineral fraction but also through the preservation of stable aggregates, aided by the positive influence of various forms of iron on aggregate stability [32].

Carbon Sequestration Efficiency
As per Wang et al., [27] carbon sequestration efficiency is predominantly associated with soil fertility.On average, 22.9% and 45.6% of the carbon input in the farmyard manure in our experiment was transformed to the soil organic carbon content in the F and F+M3 treatments, respectively.These results are close to the ones reported by Sedlář et al., [33], where an increase in the carbon sequestration in the soil was recorded following mineral fertilizer treatments.The F+M2 treatment in our experiment showed the highest content of C HA and the highest C HA /C FA ratio.In addition, the F+M2 treatment produced the highest humification index (HI).An elevated humification index fosters more favorable conditions for carbon sequestration [34,35].

Carbon Fractions and the E4/E6 Ratio
Based on the results of the E4/E6 ratio (Table 2), the Cambisols of the highest quality were at the U. Ostroh and Horažd'ovice sites.The relationship between the E4/E6 ratio and the C HA /C FA has been found to be site-specific [36,37].There were no statistically significant differences in the C HA /C FA ratios between the treatments.On the other hand, significant differences were established between the F+M2 and unfertilized Con treatments for the E4/E6 ratio.Similar results of increasing E4/E6 ratios in farmyard manure treatments were also published by Song et al., [38] and Galantitni and Rossel [39], although they reported increased contents of fulvic acids resulting fron farmyard manure treatments.Higher aliphatic and phenolic -OH group contents are generally found post organic fertilizer application [39].The impact of fertilization on carbon fractions (C FA , C HS , HI, HR, HA) was smaller compared to impact on the E4/E6 ratio.Similarly, Gerzabek et al., [40] and Oktaba et al., [41] reported significant changes in theE4/E6 ratio caused by fertilization, while there was no effect on the C HA /C FA ratio.The mineral fertilizer treatments exhibited a higher E4/E6 ratio compared to the Con or F treatments.Notably, the F+M3 treatment displayed a significantly distinct E4/E6 ratio in comparison to the Con (refer to Table 3).A parallel finding was observed in another study, wherein treatments involving mineral + organic fertilizer yielded a higher E4/E6 ratio than the control or pure organic fertilizers [42].The elevated E4/E6 ratio observed in the F+M2 treatment in our experiment is indicative of a low aromaticity of humic substances and the prevalence of substantial quantities of aliphatic structures.No significant relationship and no significant regression were present between the C SOM and E4/E6 ratio (Tables 4 and 5).

The Potential Wettability Index
The potential wettability index (PWI) can be used to describe the quality of the soil organic matter.Lower aggregate wettability is indicated by high PWI values [5].Significant differences in the PWI values between F+M2 and other treatments were observed.In our previous work, a significant influence of farmyard manure fertilization was reported [7].In line with this, Demyan et al., [6] reported a similar increase in the PWI following farmyard manure application.The application of manure led to an augmentation of the hydrophobic particles and played a role in the development of larger soil aggregates in these treatments.Secondary metabolites generated during the organic matter decay process can possess hydrophobic characteristics [43].Additionally, mineral fertilization has been shown to contributed to heightened root biomass production, an increase in root exudate production, and the formation of stable aggregates [5].
A significant correlation between the PWI and the C SOM , C HS , C FA , C HA , C HWC and C DOC was observed (Table 4).The relationship between PWI and C SOM was further investigatedusing regression analysis and confirmed (Table 5).Literature also confirms this robust correlation [44].In our experiment, The TG and PWI exhibited a strong correlation, as depicted in Table 4.This correlation was also identified in our prior study involving a maize monoculture on a Luvisol [7].Given that Glomalin-Related Soil Protein (GRSP) is a temperature-stable, adhesive, hydrophobic glycoprotein [13], the rise in glomalin content resulted in a higher ratio of hydrophobic particles, leading to an observable increase in the PWI.The calculation of the aromaticity index (iAR) was based on the reflectance of aliphatic and aromatic bands [9].Notably, significant differences in the iAR were observed among the treatments, as detailed in Table 3.Furthermore, a robust correlation was identified between the iAR and the C SOM , along with its quality indicators (C HS , C FA , C HA , and C HWC ).

Relationship of GRSP and Organic Matter Quality Indicators
Glomalin estimation is a very specific extraction method.Glomalin-related soil proteins (GRSPs) are a mixture of humic substances and heat-stable glycoproteins [45].A strong correlation was established between the glomalin content (EEG and TG contents) and the C SOM , C HA, C FA , and C HWC contents (Table 4).Regression analysis further confirmed these relationships (Table 5).A positive correlation was also observed between the GRSP content and the C SOM , C HA , and C HA /C FA ratios in our previous work on long-term experiments with a maize monoculture on a Luvisol [7].Similarly, other studies in the literature [46] established a correlation between humic and fulvic acid contents, C OX , glomalin content, and the C/N ratio.Correlation coefficient values were greater for theTG content than for the EEG content.The repeated extraction used for the TG estimation could be the cause for a higher correlation for the TG content.In addition to GRSP, there is an additional release of humic substances.It's worth noting that the Bradford assay can exhibit a cross-reaction with humic acids, polyphenolic compounds, sugars, and lipids [47], and this can interfere with glomalin determination.Likewise, Li et al., [15] highlighted that the total glomalin content in soil exhibits strong linkages with the overall content of the soil organic matter.In our study, a significant increase in the glomalin content after treatment with farmyard manure was observed.Similarly, Bertagnoli et al., [48] and Zhang et al., [16] reported a heightened the TG content following the application of farmyard manure.The synergistic use of organic and mineral fertilizers led to an enhanced TG content, particularly evident in the F+M3 treatment.The contents of TG and EEG in the F+M3 treatment were 25% and 24% higher than in the Con treatment, respectively.The F+M3 treatment yielded a high biomass, closely linked to elevated post-harvest residues (stubble) and root biomass production.This likely resulted in an increased exudate production, manifesting as a rise in soil organic carbon (C SOM ) content.Therefore, higher arbuscular mycorrhiza activity does not have to be the primary contributor to the increase in the TG content.More clarity into this issue can be provided by a different study with an 80-year fallow [49].The authors conclude that recycling of the SOM continually produces soil protein.28.3% of TG content was founded by the EEG content.Interestingly, a very similar proportion (28.8%) was estimated on a Luvisol at the Červený Újezd experimental site [7].The ratio of TG in C SOM was found to be 0.186, on average.This is a significant proportion, and it is obvious that glomalin has a role in soil carbon sequestration.Glomalin plays an important role in promoting the stability of soil organic carbon, in terms of the slow decay of glomalin (and a low soil turnover rate) [50].

Conclusions
Drawing from the outcomes of long-term field experiments conducted on Cambisols at four sites in the Czech Republic, it can be asserted that: (i) The amalgamation of organic and mineral fertilizers enhances both the quantity and quality of soil organic matter when compared with unfertilized conditions or treatments involving solely organic fertilizer.The highest intensity of mineral fertilizers (F+M3) increased the average value of carbon sequestration efficiency (CSE) to 45.6% in comparison with pure manure treatment (F) which reached only 22.9%.(ii) We established a robust correlation between total glomalin (TG) content and indicators of both the quality as well as quantity of soil organic matter.(C SOM , C FA , C HA , C HWC , and C SOM /N t ).On the other hand, the correlation between easily extractable glomalin (EEG) and the same indicators of soil organic matter quality and quantity mentioned above was weaker.(iii) We demonstrated a strong correlation between the PWI, and soil organic matter quality as well as quantity indicators (C SOM , C FA , C HA , C HWC , C DOC, and C SOM /N t ) TG content and EEG content.The data regarding the (PWI) can be employed for the investigation of organic matter quality.The E4/E6 ratio indicatorwas found to be a considerably less sensitive method for reflecting alterations in both the quality and quantity of soil organic matter.

Table 1 .
Experimental site characteristics and their founding years.

Table 2 .
Carbon content and qualitative parameters of soil organic matter content at the individual experimental sites.

Table 3 .
The impact of fertilization on both soil carbon content and the qualitative parameters of soil organic matter.
Rows represent relative values of monitored variables.Different letters are employed to denote statistically significant differences among treatments.Tukey HSD test, p < 0.05.n = 12.

Table 5 .
Linear regression equation of C SOM and selected indicators of soil organic matter quality.