Next Article in Journal
“Purplish Blue” or “Greenish Grey”? Indigo Qualities and Extraction Yields from Six Species
Next Article in Special Issue
Response of Pasture Grasses to Organic Fertilizer Produced from Black Soldier Fly Frass
Previous Article in Journal
Hybridization in the Subtribe Alopecurinae Dumort. (Poaceae) According to Molecular Phylogenetic Analysis: Different Ploidy Level Tells Different Origin of the Groups
Previous Article in Special Issue
Effect of a Fortified Biostimulant Extract on Tomato Plant Productivity, Physiology, and Growing Media Properties
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Effect of Biostimulant, Manure Stabilizer, and Manure on Soil Physical Properties and Vegetation Status

1
Department of Machinery Utilization, Faculty of Engineering, Czech University of Life Sciences Prague, 165 00 Prague, Czech Republic
2
Agrovyzkum Rapotin Ltd., Vyzkumniku 863, 788 13 Rapotin, Czech Republic
*
Author to whom correspondence should be addressed.
Plants 2024, 13(7), 920; https://doi.org/10.3390/plants13070920
Submission received: 27 February 2024 / Revised: 16 March 2024 / Accepted: 18 March 2024 / Published: 22 March 2024
(This article belongs to the Special Issue Biostimulant and Biofertilizer Use in Agricultural Systems)

Abstract

:
Food production sustainability is one of contemporary agriculture’s fundamental challenges. Farmers are currently facing high input prices in crop production and declining organic matter in the soil. For this reason, a field experiment was established to assess the effect of the biostimulant NeOsol (NS), the manure stabilizer Z’fix (ZF), farmyard manure (FM), and their combination in farm practice. In situ measurements provided information on the change in bulk density (BD), unit draft (UD), saturated hydraulic conductivity (SHC), and cone index (CI). Furthermore, the vegetation status was investigated via vegetation indices, and the yield and quality parameters were assessed. Management of the experimental field resulted in an overall decrease in BD over time for the treated variants compared to the control (CL). The decrease with time was also verified in the case of UD and CI at the depth zone of 10–20 cm. Variants FM (by 8.0%), FM_NS (by 7.3%), and FM_ZF_NS (by 3.8%) proved to have lower UD values than CL. An overall increase in SHC and in yield was observed over time. Concerning SHC, only FM (by 58.5%) proved different from CL. The yield of all the treated variants, i.e., NS (by 8.2%), FM (by 10.8%), FM_NS (by 14.1%), FM_ZF (by 17.8%), and FM_ZF_NS (by 20.1%), surpassed CL. Simultaneously, none of the examined treatments proved to have any adverse effect either on soil or on plant-related variables.

1. Introduction

Contemporary agriculture is facing demanding challenges unprecedented in history. Currently, there are eight billion people in the world. The United Nations estimates there will be over 8.5 billion people by 2030 and even more than 11 billion by 2050 [1]. With this population growth, a decrease of 50 million hectares of arable land in developed countries is also expected [2]. Regarding the Green Deal, a decline in livestock production of 10 to 15 percent can be expected [3].
With the anticipated decline in livestock farming, there will also be a decline in farmyard manure production and availability. While farmyard manure has a positive effect on soil fertility and consequently on plant yields, synthetic fertilizers provide large amounts of nutrients to the soil—and hence increase yields—but reduce the carbon content of the soil and do not improve soil health [4]. In terms of sustainability, it is appropriate to integrate plant-based products—such as algae or microorganism-based products—into agricultural practice [5,6]. Some companies may add value to products, e.g., stranded seaweed on beaches, thus facilitating the restoration of ecosystems. This is the case with the French group Olmix [7].
The company produces diverse biostimulants, soil improvers, agents, and stabilizers; this study investigated the following two: NeOsol (formerly known as PRP SOL) and Z’fix (formerly known as PRP FIX). According to the manufacturer, (a) NeOsol stimulates the natural soil microflora, which is involved in the recycling and transforming of crop residues. By supporting humus synthesis, NeOsol enhances the components of soil fertility and soil health. (b) Z’fix reduces ammonia emissions from stable bedding or effluents and moreover improves effluents’ fluidity and homogeneity. This simplifies spreading work and assists in the preservation of nutrients for the crops.
Several studies have been conducted to explore the efficacy of these amendments: (a) Findura et al. [8] observed over a three-year period that biostimulants positively influenced crop yield and soil bulk density and increased oxidizable carbon. Similarly, Šindelková et al. [9] reported a favourable effect on soil physical and chemical properties in a short-term experiment. In addition, the biostimulant improved the C/N ratio [10]. Subsequent short-term research indicated improvements in water retention [11]. The two-year results of the study by Tuba et al. [12] showed a reduction in cone index and also in the energy demand for ploughing. The authors recommended long-term use of the product for cumulative positive effects. A long-term study (10 years) of NeOsol investigated its effect on the chemical composition of soils. Soil pH increased and the content of available forms of magnesium increased; on the other hand, there was a decrease in available forms of potassium [13]. In a two-year study, a biostimulant significantly increased the chlorophyll content of rye leaf blades [14]. NeOsol combined with urea increased soybean yield [15]. The positive effect of the product was also demonstrated in potatoes, where the highest greenness index was observed [16]. The product improved the growth and yield of calendula [17].
(b) The results of the three-year experiment proved that manure treated with Z’fix positively affected the soil conditions and crop status [18]. Slurry stabilizer treatment resulted in higher macronutrient content in the soil [19]. This agent was found to positively influence the fermentation process and the manure quality [20]. Plants usually contained more nitrogen, phosphorus, potassium, magnesium, calcium, and sulphur when treated slurry had been applied [19]. Higher nitrogen content was achieved also in another study [21].
Even with the increasing use of biostimulants in agricultural practice, many scientific community members perceive that biostimulants need peer-reviewed scientific evaluation [22]. The majority of the results described above were based on short-term experiments.
This study aims to verify the effect of the biostimulant, the stabilizer, manure, and their combinations in real conditions of agricultural practice in a long-term experiment carried out on Modal Luvisol of the silt loam soil texture. The hypotheses to be verified are as follows: the use of these improvers would result in (a) a reduction in the unit draft, soil bulk density, and cone index, (b) an increase in saturated hydraulic conductivity, (c) an improvement in crop growth, and (d) an increase in crop yield and improved qualitative parameters.

2. Results

Soil Physical Characteristics

Variables describing soil physical characteristics were assessed in three years, i.e., terms (see Materials and Methods), following autumn manure application. Figures present their values relative to the value of the respective control variant (CL) since, in addition to the variants, there may have been other factors affecting the variables’ absolute values, e.g., soil moisture and machinery passes. These factors can be considered though to have had the same effect on all the variants including CL. The tables contain the absolute values with significant differences marked within individual terms. Within individual terms, one-way analyses of variance (ANOVA) provides the same outcome for the absolute values as well as for the relative ones. When using factorial (variant; term) ANOVA, the outcomes of analysis differ depending on whether absolute or relative values are compared. If the latter ones are used, the differences among terms will provide information on the variable development relative to CL.
The bulk density (BD) values related to the CL variant are shown in Figure 1a. Table 1 presents the absolute values and their statistical evaluation analysis for each term. Factorial ANOVA of relative values proved that Term III (95.609%) was significantly different compared to Term II (100.075%) and Term I (102.835%). This suggested a positive effect of the treatments over a longer period of time. Factorial ANOVA of absolute values proved that Term I (1.491 g cm−3) was significantly different compared to Term II (1.328 g cm−3) and Term III (1.353 g cm−3). One-way ANOVA proved that there were significant differences among variants solely in Term I, when standard deviations demonstrated the lowest values. The BD of CL and FM variants attained significantly lower values than that of FM_ZF (see Table 1).
Figure 1b presents the unit draft (UD) values of the variants relative to CL. The factorial ANOVA proved that individual terms significantly differed among themselves, with Term I showing the highest overall value of 99.454%, Term II showing a value of 97.439%, and Term III showing the lowest value of 94.629%. Factorial ANOVA of absolute values proved also significantly different among terms. In that case, the differences should be attributed to other varying factors rather than to the effect of variants, such as soil moisture, working depth, working speed, preceding crop, etc. One-way ANOVA proved that there were significant differences among variants within all three terms. In Term I (see Table 1), variants formed three significantly differing homogenous groups, with FM_ZF forming that with the highest UD and FM_NS forming the group with the lowest UD. In Term II, the average UD of the variants CL, FM_NS, and FM_ZF_NS significantly exceeded the UD of FM and FM_ZF. In Term III, variants formed again three significantly differing homogenous groups with the variants CL, NS, and FM_ZF forming that with the highest UD and FM forming the group with the lowest UD. Based on the analyses described, there does not seem to be any pattern discernible. When performing factorial ANOVA of relative values combined for all three terms though, variants form four homogenous groups with an apparent pattern, where only CL belongs simultaneously to two groups. The group with the highest UD values consists of FM_ZF (101.422%) and CL (100.000%). Then, a group of CL and NS (98.758%) follows, trailed by a group comprising a single variant, FM_ZF_NS (96.231%). Finally, the group with the lowest UD is formed by FM_NS (93.726%) and FM (92.047%).
Saturated hydraulic conductivity (SHC) values relative to the control variant CL are presented in Figure 1c. Factorial ANOVA of relative values proved that Term I (74.658%) was significantly different compared to Term II (127.622%) and Term III (134.559%). This result suggests that the SHC of the treated variants improved over time relative to the control variant. Factorial ANOVA of absolute values proved that there were significant differences among all three terms, with Term I showing the lowest overall value of 5.038 mm h−1, Term III showing a value of 22.982 mm h−1, and Term II showing the highest value of 37.718 mm h−1. One-way ANOVA proved that there were significant differences among variants within all three terms. In Term I (see Table 1), the SHC of CL attained a significantly higher average than those of NS and FM_NS. In Term II, the SHC of FM_ZF_NS attained a significantly higher value than all the other variants except FM. Finally, in Term III, it was the FM variant that demonstrated a significantly higher SHC than all the other variants.
Figure 2 presents the cone index (CI) averages of the variants relative to CL at three soil depth zones of 0 to 10 cm, 10 to 20 cm, and 20 to 30 cm. Their factorial ANOVA proved that there was a significant difference among terms only at the middle zone of 10 to 20 cm, with Term III differing the most with its lowest overall value of 78.742%when compared to the Term I value of 96.234% and the Term II value of 96.825%. With respect to the variants, the factorial ANOVA of relative values found significant differences at the depth zones of 10 to 20 cm and 20 to 30 cm. At the 10 to 20 cm depth, the CI average of FM_NS (82.096%) significantly differed from FM_ZF (97.656%) and from CL (100.000%); the CI average of FM (85.989%) also significantly differed from CL (100.000%). At the 20 to 30 cm depth, the CI average of FM (83.136%) significantly differed from FM_ZF (99.590%) and from CL (100.000%). Factorial ANOVA of absolute values proved that there were significant differences among terms at all three depth zones. At the 0 to 10 cm depth, the CI averages differed significantly among all three terms, with Term III demonstrating the lowest value (0.378 MPa), followed by Term I (0.558 MPa) and, finally, Term II (0.975 MPa). At the 10 to 20 cm depth, the CI averages differed significantly among Term III (again, the lowest value was 0.711 MPa) and both Term I (1.232 MPa) and Term II (1.288 MPa). At the 20 to 30 cm depth, the CI averages differed significantly again among all three terms, with Term III demonstrating the lowest value (1.650 MPa), followed by Term II (1.902 MPa) and, finally, Term I (2.314 MPa). With respect to the variants, the factorial ANOVA of absolute values proved that there were significant differences at the depth zones of 10 to 20 cm and 20 to 30 cm that presented a similar pattern to the factorial ANOVA of relative values. At the 10 to 20 cm depth, the CI average of FM_NS (0.990 MPa) significantly differed from FM_ZF (1.159 MPa) and from CL (1.171 MPa). At the 20 to 30 cm depth, the CI average of FM (1.760 MPa) significantly differed from FM_ZF (2.103 MPa) and from CL (2.105 MPa).
One-way ANOVA proved that there were significant differences among variants only at the soil depth zones of 10 to 20 cm and 20 to 30 cm (see Table 2). At the 10 to 20 cm depth, the CI average of FM_ZF_NS (lowest) significantly differed from FM_ZF (highest) in Term I, FM_NS (lowest) significantly differed from FM_ZF_NS (highest) in Term II, and, finally, FM_NS (lowest) significantly differed from CL (highest) in Term III. At the 20 to 30 cm depth, significant differences among variants were verified solely in Term I. In that case, the CI average of FM_ZF_NS (lowest) significantly differed from FM_ZF (highest), which was the same outcome as for the shallower 10 to 20 cm depth in that term.
Figure 3a presents yields of the variants relative to CL for all the years that the experiment was conducted, i.e., from 2015 (Term I) to 2020 (Term III). Their factorial ANOVA proved that there were significant differences among the years as well as among the variants. Concerning the differences among years, the highest yield relative to CL was reached in 2020 (119.896%) and was significantly higher than in all the other years in question. Moreover, the years 2018 (113.177%) and 2019 (112.578%) demonstrated a significant difference compared to the year 2016 (105.048%). In terms of values relative to CL with all the experimental years combined, the variants FM_ZF_NS (120.130%) and FM_ZF (117.780%) significantly differed from FM (110.836%), NS (108.200%), and CL (100.000%). In fact, the latest mentioned CL proved to be significantly lower than all the other variants. With respect to changing crops, comparing absolute values among individual years was not required. One-way ANOVA of absolute values proved that there were significant differences in average yields among variants only in the years 2018, 2019, and 2020 (see Table 3). In 2018, the average yields of FM_ZF_NS (highest) and FM_ZF significantly differed from that attained by CL (lowest). In 2019, FM_ZF_NS reached again the highest yield that significantly differed from NS and CL (lowest). Also, FM_ZF significantly exceeded CL. In 2020, the yields of all five treated variants significantly surpassed the yield produced by CL.
Figure 3b presents the values of the qualitative parameters of the variants relative to CL for all the years that the experiment was conducted. In all cases but one, the qualitative parameter monitored was crude protein content (see Table 3). The factorial ANOVA of relative values proved that there were significant differences among the years only, not among the variants. The qualitative parameter content relative to CL that was reached in 2016 (113.084%) surpassed significantly the value attained in all the other experimental years.
One-way ANOVA of absolute values proved that there were significant differences in qualitative parameter averages among variants only in the years 2016 and 2019 (see Table 3). In 2016, the average crude protein content in winter wheat with FM_ZF_NS (highest) significantly differed from the one attained by CL (lowest). In 2019, FM_ZF_NS reached again the highest average content that, together with the one attained by FM, significantly exceeded NS (lowest).
Figure 4 summarises the results of the soil- and plant-related variables in question over the whole monitored period and provides outcomes of the factorial analyses of variance (ANOVA). Concerning UD and CI (see Figure 4a), there were improvements in the case of individual variants except for CI at 0–10 cm of depth. FM (by 8.0%), FM_NS (by 7.3%), and FM_ZF_NS (by 3.8%) had a lower UD than CL. Regarding CI, FM_NS (by 17.9%), and FM (by 14%), they demonstrated lower values than CL at 10–20 cm, and FM was lower again (by 16.9%) at 20–30 cm of depth. With regard to BD and quality parameters (Figure 4b), no substantial differences were found. Concerning SHC (Figure 4b), solely FM (by 58.5%) proved different from CL. On the other hand, the yield (Figure 4b) of all the treated variants, i.e., NS (by 8.2%), FM (by 10.8%), FM_NS (by 14.1%), FM_ZF (by 17.8%), and FM_ZF_NS (by 20.1%), was higher compared to CL.
Table 4 presents the correlation matrix of the variables (BD, UD, SHC, CI, yield, qualitative parameters). Their values relative to the value of the respective control variant (CL) were analysed since, in addition to the variants, there may have been other factors affecting the variables’ absolute values, e.g., soil moisture, machinery passes, and different crops grown. These factors can be considered though to have had the same effect on all the variants including CL. A significantly increasing correlation was verified between BD and UD, between BD and CI at 10–20 cm (BD was measured in the topsoil though), between UD and CI at 10–20 cm, between UD and CI at 20–30 cm, and between CI at 10–20 cm and CI at 20–30 cm. On the contrary, a significantly decreasing correlation was confirmed between BD and SHC, between BD and yield, between BD and quality parameters, between CI at 10–20 cm and yield, and between CI at 20–30 cm and yield. In other words, higher bulk density may have adversely affected unit draft, saturated hydraulic conductivity, quality parameters, and yield. The latter was negatively influenced by the cone index at 10–30 cm of depth.
Table S1 (see Supplementary Materials) presents the values of vegetation indices: (a) NDVI, (b) MSI, and (c) GCI. The images originated from Sentinel-2 satellites using the Google Earth Engine platform, where cloud cover did not exceed 10%. Each experimental variant consisted of 92 pixels. These indices provide spatial information on crop conditions in terms of (a) plant health, (b) water content, and (c) chlorophyll content. In addition, the table contains the one-way ANOVA results for each date separately.
From the results of the vegetation indices, the interconnectedness of water stress and plant health can be observed, as the MSI and NDVI results often correspond with each other. The better values for the CL variant at the beginning of the experiment in 2015 may be attributed to later observation dates in terms of the vegetative period. In the 2018 season, the best health status was reached by the FM variant, where the lowest water stress was also observed. In the 2019 season, it was the FM_ZF_NS variant that achieved high NDVI and low MSI values on most dates, while other variants were rated the highest as far as the chlorophyll content according to the GCI index. Regarding the GCI index, the FM_ZF_NS variant reached often the highest values, which corresponded with its yields having been also the highest (see Table 3).

3. Discussion

Numerous studies have demonstrated a correlation between soil organic matter content [23,24] or organic amendments [25,26] and bulk density [24]. According to the USDA, the ideal bulk density for crop growth of silt loam soils is <1.4 g cm−3, the root growth is adversely affected at the level of 1.55 g cm−3, and root growth is restricted at the level >1.65 g cm−3 [27]. In this respect, root growth was not restricted throughout the experiment. In Terms II and III, all variants reached the ideal bulk density, except for the CL variant in Term III (see Table 1). The lowest BD values were achieved in Terms II and III for the FM variant, which was in line with the results of Yu et al., who found a reduction in soil bulk density after cow manure application [28]. The FM_ZF variant reached higher values than FM alone; this finding was supported by the results of Čermáková et al. [29]. The significant reduction in BD between Term I and the following terms indicated a positive effect of long-term application of organic matter to soil, as confirmed by Fu et al. [30].
Tillage is one of the most energy-intensive operations in agriculture [31]; high amounts of energy are devoted to cutting, breaking, and possibly overturning soil layers, reducing the size of clods and rearranging soil aggregates [32]. The results of this study for the NS variant, namely the 1.2% reduction in UD (see Table 1), though not significant, support the findings of the study by Urbanovičová et al. [33], who reported that biostimulants had a positive effect on UD reduction. In that study, UD was reported to decrease by 5.71% when the biostimulant was applied. The decrease in UD with biostimulant treatment was also confirmed by the results of a study conducted by Tuba et al. [12], who reported a 9% reduction in UD in the first year after application. The FM_ZF variant had overall higher UD values than the FM variant; this finding is not consistent with the authors’ previous study [18]. The combination of stabilized manure and biostimulant led to a 3.8% reduction in UD compared to CL. This finding is consistent with a previous study by Šařec et al. [34], where a reduction in UD was also observed after a four-year experiment. Overall, the FM variant presented the highest reduction in UD, which was consistent with a study [35] where plough draft was reduced by 27–38% after eight years of manure application, although the reduction was not as high in our study. The effect of organic matter addition on the reduction in the plough draft was also confirmed in another study [36]. In our study, UD was reduced by almost 8% for the FM variant compared to CL. This reduction could lead to fuel savings of approximately 0.8 L ha−1 (assuming an average power delivery efficiency of around 50% and a fuel demand for soil tillage at the level of 20 L ha−1). In terms of sustainability, it is necessary to consider the high energy requirements of tillage not only from the agronomic perspective but also from the economic and ecological ones [37]. The carbon dioxide (CO2) emitted by the combustion of one litre of diesel oil is approximately equivalent to 2.6–2.8 kg [38,39]; an additional 0.5 kg of CO2 is emitted to produce and distribute this one litre of diesel oil. In total, an emission rate of around 3.2 kg of CO2 per litre of combusted diesel oil can be assumed [40]. For the FM variant, the CO2 emissions coming from diesel combustion are therefore 5.12 kg ha−1 lower. Although the reduction in CO2 emissions can be ensured by reduced tillage [41], it is necessary to address this issue also in the case of conventional technology. CO2 emissions in agriculture, however, are a very complex issue that needs to be considered in a broad ecosystem cycle.
The cone index (CI), also often known as penetration resistance, is considered one of the crucial indicators of soil compaction [42]. Cone index values between 2 and 3 MPa, measured at field capacity, are considered upper limits for root growth [43,44,45]. Exceeding these values can lead to a restriction in root growth and a reduction in plant biomass and yield [46]. CI values up to the depth of 30 cm did not exceed 2.5 MPa (see Table 2). From this point of view, there was no significant restriction of root growth. Manure treated with a stabilizer reached lower values in the upper soil layer than untreated manure, this finding is in accordance with the results of previous studies [18,29]. Celik et al. [47] reported a positive influence of organic amendments on CI; this was confirmed by our results in depth zones of 10–20 cm and 20–30 cm, where FM yielded 14% lower values compared to CL. Tuba et al. [12] reported a 17–23% reduction in penetration resistance when the biostimulant was used; our results did not confirm this finding.
Saturated hydraulic conductivity (SHC) is an important soil hydraulic parameter that affects water flow and dissolved solute transport [48]. SHC is strongly correlated with the soil’s physical properties [49] and can be used as an indicator of water erosion [50]. This issue is currently very significant, as soil loss due to water erosion is expected to increase by 13–22.5% in the EU and UK by 2050 [51]. According to estimates by the Czech Ministry of Agriculture, water erosion causes a total damage of CZK 4.3 billion annually in the Czech Republic. This figure represents solely soil loss and does not include any other damages, e.g., to properties [52]. The SHC values can present a relatively large variability for the Simplified Falling Head Method, especially when a smaller ring diameter is used [53]. However, all outliers in our study were excluded based on Dixon’s Q test. The significant overall increase in SHC over the years suggests a positive effect of field management on soil infiltration (Table 1). Despite the large variability in results, very similar SHC values were measured for the FM and FM_ZF variants in Terms II and III. The positive effect of manure on SHC was reported by Miller et al. [54], who found an increase in SHC from 76 to 128%. This is in accordance with the findings of our study, where there was an increase of 77% in Term II and 308% in Term III. However, treatment of manure with a stabilizer did not lead to an increase in SHC compared to untreated manure, which is inconsistent with the authors’ previous findings [18]. The use of the biostimulant did not lead to an improvement in SHC, which is in accordance with the results of the previous study [55].
Currently, it is necessary to increase crop yields in order to ensure food security related to the demands and needs of a growing population [56]. While organic fertilizer application improves soil health [24,57], the inappropriate use of industrial fertilizers, especially nitrogen and phosphate fertilizers, has a harmful effect on soil health and soil-related ecosystem services [58]. The application of the biostimulant resulted in a higher yield each year compared to the CL variant, though significantly in 2018 and 2019 only (Table 3). This result was consistent with the study by Sulewska et al. [59], who found that phosphorus–potassium fertilization could be replaced with a biostimulator without yield loss. Manure treatment with stabilizer resulted in increased yield each year compared to untreated manure, though again significantly in 2018 and 2019 only. This result was consistent with the study [18], where an increase in sugar beet yield was reported, and with another study [60], where an increase in silage maize yield was observed. The increasing yields of the treated variants compared to CL over the years could be attributed to an extent to an improvement in the soil environment.
Fertilization using manure, whether treated or not, led to an increase in quality parameters, though significantly in 2016 and 2019 only (Table 3). This finding is consistent with studies that show that fertilizer management affects the crude protein content of maize [61,62,63], winter wheat [64,65], and barley [66]. Nitrogen application also influenced crude fat content in canola [67]. Sulewska et al. [68] noted that biostimulant application led to a similar fat content in dry weight of rape seeds as fertilization with conventional P and K. This was in line with the results of our study, where no significant difference was proven between NS and CL variants. Zielewicz et al. [69] observed a 7% increase in crude protein in ryegrass sward when a biostimulant was used, but no evidence for this increase was found in our results, although the crops were different. In addition, Możdżer et al. [70] found that treatment of slurry with a stabilizer and PK fertilization increased the macro-elements’ concentrations in winter rape and spring wheat grain; however, the significant qualitative parameter improvement was not confirmed in our results when the stabilizer was used.
Remote sensing data provide a cheap and relatively accurate source of information that can be used directly in agricultural practice. However, augmenting these data with other in situ data alleviates correct interpretation. In this respect, for instance, an interesting relation could be observed between the SHC values and MSI index in Term I (Table S1), where the lowest water stress was achieved by the CL variant, which presented the highest SHC. In Term II, the variant with the highest SHC had the highest GCI, suggesting the highest chlorophyll content. The results of other researchers indicated a positive effect of the biostimulant on chlorophyll content [14,16,71,72]. These findings partially agree with our results, where the NS variant achieved higher GCI values on most dates compared to the CL variant.

4. Materials and Methods

4.1. Biostimulant NeOsol

NeOsol, which the Olmix Group produces in Brégan, France, is a granular biostimulator designed to enhance essential soil functions. Utilizing the patented Mineral Inducer Process (MIP) technology, this soil amendment uses the bioactive properties of minerals and specific trace elements to activate biological reactions in the soil. Crucial ingredients such as iron, manganese, copper, and boron stimulate enzymatic reactions essential for the decomposition of raw organic matter, particularly aiding humification processes (e.g., α-glucosidase, β-glucosidase, etc.). The NeOsol complements the MIP technology with SEAweed DRY algae extracts, which provide rich nutrients that stimulate the soil biota.
The composition of this biostimulant is reported to include 28.0% w/w of CaO, 17.0% w/w of MgO, and 98.9% w/w of dry matter, with combustible substances constituting 7% w/DMw. With a pH ranging from 8 to 10, it is highly alkaline.
Depending on the specific crop and local soil conditions, the recommended dosage varies between 110 and 220 kg ha−1. NeOsol is typically applied post-harvest and is spread on the soil surface.

4.2. Stabilizer Z’fix

Z’fix, which the Olmix Group produces in Brégan, France, is used as an activator of biological transformation in stables, aiming to improve the quality of bedding by controlling the fermentation process of organic matter. The main emphasis lies in improving animal welfare. In addition, the manufacturer mentions a secondary benefit: its ingredients control fermentation in organic matter to preserve the maximum amount of fertilising elements. Z’fix is produced in the form of granules that contain calcium and magnesium carbonates, together with a mixture of micro- and macro-elements such as potassium, sodium, sulphur, iron, and manganese. These ingredients are designed to regulate the fermentation processes in manure, slurry, and compost. The Z’fix also uses MIP technology. Table 5 shows the selected chemical parameters of treated and non-treated FM. The FM used originated from cattle stabled on deep straw bedding. Z’fix was applied at the rate of 1 kg head−1 week−1 directly to deep bedding.
The composition of Z’fix is organic matter—5% w/w, CaO—37.5% w/w, MgO—4.4% w/w, Na2O—3.9% w/w, SO3—0.7% w/w, K2O—0.5% w/w, P2O5—0.1% w/w, Fe—2000 ppm, Mn—150 ppm, and Zn—30 ppm.

4.3. Site and Crop Management

The six-year field experiment was conducted near the town of Lázně Bohdaneč in the Královéhradecký region of the Czech Republic (50°27′15.2″ N 15°34′12.5″ E, 365 masl.). According to the Czech Hydrometeorological Institute, this site is located in an area with an annual average temperature of 8.2 °C and an annual rainfall of 732 mm (normal between 1991 and 2020). Table 6 represents the site’s soil conditions at the beginning of the experiment. The experiment was performed on high-quality agricultural soil, namely Modal Luvisol, which, according to the United States Department of Agriculture (USDA) texture triangle, belongs to the silt loam category.
Experimental field of the agricultural company ZEPO Bělohrad a.s. contained six plots, i.e., variants differing by fertilizer and biostimulant applied: farmyard manure stabilized by Z’fix and combined with biostimulant NeOsol (FM_ZF_NS); farmyard manure stabilized by Z’fix (FM_ZF); conventional farmyard manure combined with NeOsol (FM_NS); farmyard manure (FM); biostimulant NeOsol (NL); and a control variant (CL), where only NPK fertilizer was applied at the rate according to the farm standards, respecting crop requirements. In addition to CL, the NPK fertilizer was applied also within the other variants at a rate securing the same overall nitrogen dose., e.g., in the case of NL, the NPK rate did not differ from CL; when complementing farmyard manure, the NPK rate calculated acknowledged manure’s multi-year decomposition. Each variant plot was a rectangle of 30 m × 300 m (0.9 ha), and the positioning on the experimental site was chosen with respect to the headlands. Table 7 provides information on crop rotation and application rates of farmyard manure, both stabilized by Z’fix and conventional manure, and of the biostimulant NeOsol for the respective variants. For the purpose of the article, soil physical characteristics were assessed each season (term) following the autumn manure application, i.e., in the years 2015, 2017, and 2020 (see Table 7).

4.4. Data Acquisition and Processing

In situ measurements were always considered for the years following manure applications, i.e., Term I (2015), Term II (2017), and Term III (2020). With regard to the more homogenous moisture distribution in the soil profile, the following variables were measured in spring: cone index, bulk density, and saturated hydraulic conductivity.
Cone index (CI) is used as a crucial indicator of pedocompaction. Its advantage is the relatively rapid and uncomplicated measurement over the entire soil profile—depending on the moisture content. CI is, in fact, a measurement of the soil’s resistance against a cone with a well-defined geometric characteristic (angle, area). The PN70 penetrometer used in this study was developed at the Czech University of Life Sciences Prague (Prague, the Czech Republic). This tailor-made instrument fulfils all the requirements of the ASAE S313.3—the American Society of Agricultural and Biological Engineers (ASABE) agricultural standard. CI was measured in ten repetitions per variant.
Furthermore, bulk density (BD) was determined as another indicator of soil compaction since it can influence all the key soil processes (infiltration, rooting depth, available water capacity, soil porosity, availability of nutrients, activity of soil micro-organisms, etc.) and the productivity of the soil. A soil sample ring kit (Eijkelkamp, Giesbeek, Zevenaar, The Netherlands), where the volume of each ring was 100 cm3, was used to collect undisturbed soil samples. BD was sampled in six repetitions per variant, in the topsoil (0–5 cm), without the soil crust—if present. Subsequent analysis of the BD soil samples was carried out in the laboratories of the Faculty of Engineering, Czech University of Life Sciences Prague, according to the national standard CSN EN ISO 17892-2 [73].
The unit draft (UD) of a tillage implement was measured in autumn—after the crop was harvested. The method of one tractor towing another tractor using a bar equipped with a dynamometer, particularly with strain gauges S-38/200 kN/(LUKAS, Prague, Czech Republic), was used. The towed tractor was mounted with the implement. The towing bar required a horizontal positioning between the two tractors, and it was essential that the working speed and working depth were as constant as possible. In the case of this method, several passes of the above-described gang of machines per one variant were necessary. The passes differed by the implement either working (at least four passes per variant) or not working (at least two passes per variant). In the first case, the dynamometer measured the overall draft of the towed tractor and the working implement, i.e., the overall rolling resistance of the towed tractor and the implement and the implement’s functional draft. In the latter case of the implement not working, the dynamometer measured only the overall rolling resistance of the towed tractor and the implement. Their difference thus provided the value of a functional draft of the implementation. Also, the direction of passes, i.e., potential slope-induced downhill and uphill movement, had to be allowed for. Otherwise, the functional draft might have been biased by the tangential component of the machinery’s gravity. After each pass, the working depth of the implement was measured at both sides of the working width in order to enable unit draft calculation. The tillage depth was set to the maximum depth that the traction conditions and the engine power allowed. The tillage depth was 12.2 ± 0.7 cm in Term I, 22.4 ± 1.2 cm in Term II, and 13.6 ± 1.0 cm in Term III. In all three terms, the Farmet Hurikán tine cultivator was used—in Terms I and III, the towing tractor was a New Holland T8030 (see Figure 5), and in Term II, it was a John Deere 8320R. An NI CompactRIO controller (National Instruments Corporation, Austin, TX, USA) with relevant modules including GPS was employed for data acquisition with the frequency set to 0.1 s. Trimble Business Center 5.70 (Trimble, Westminster, CO, USA) appointed attained data to specific experimental variants.
Saturated hydraulic conductivity (SHC) was determined by a Simplified Falling-Head Technique described by Bagarello et al. [45]. This method is rapid, cost-effective, and does not demand as much water as others. The SHC was measured in fifteen repetitions per variant. Circular infiltrometers (in this study, 0.15 m diameter) and a known quantity of water (0.5 L) were used. The water was poured on the soil surface in the area of the infiltrometer, and the time of infiltration was measured. Before and after the water application, volumetric soil moisture was measured using a Theta Probe (Delta-T Devices, Ltd., Cambridge, UK).
The above-mentioned variables were used for the subsequent calculation of the SHC, resp. the Kfs value (Equation (1)):
K f s = ( Δ θ ) ( 1 Δ θ ) t α D Δ θ D + 1 α ( 1 Δ θ ) ln 1 + 1 Δ θ D Δ θ ( D + 1 α )   ,
where Δθ = difference between initial soil moisture content and saturated soil moisture content; D = the ratio of V (volume of water) and A (area of a cylinder), which is the water level corresponding to the water volume; tα = infiltration time; and α* = constant according to Elrick et al. [74]; in this study, α* = 12 m−1, which is appropriate for the majority of structured soils.

4.5. Crop Status

The European Space Agency’s Sentinel-2 satellites were used to acquire the remote sensing data. Their use is recommended for agricultural practice due to the free high-resolution and multi-spectral data (10 m per px) [75]. The data were obtained using Google Earth Engine (Google LLC, San Francisco, CA, USA) and subsequently processed using qGIS (Open Source Geospatial Foundation, Beaverton, OR, USA) and STATISTICA (TIBCO, Palo Alto, CA, USA).
Vegetation status was evaluated using the following indices: the Normalized Difference Vegetation Index (NDVI) [76], one of the most common vegetation indices used for local-scale management purposes as a direct indicator of plant health and growth [77,78,79]; the Moisture Stress Index (MSI) [80], which is highly sensitive to moisture stress in plants and correlates with soil moisture [81]; and the Green Chlorophyll Index (GCI) [82], which is used to estimate the chlorophyll content of various plant species. The chlorophyll content reflects the physiological condition of vegetation; thus, the GCI can be used as an indicator of vegetation health [83]. Comprehensive results of remote sensing observations can be found in Supplementary Materials Table S1.
Grain yield was based on three combine harvester passes per variant. The grain weight of each pass was determined using a trailer positioned on a DINI ARGEO WWSB 16t portable static axle scale (DINI ARGEO S.r.l., Modena, Italy). Mix samples for qualitative analysis were taken from each crossing (three repetitions per variant) and subsequently analysed in laboratories.
Data were processed using MS Excel (Microsoft Corp., Redmond, WA, USA), MS Access (Microsoft Corp., Redmond, WA, USA), and STATISTICA 14 (TIBCO, Palo Alto, CA, USA). Significant differences in the above-mentioned variables among investigated variants were determined through factorial (factors: variant; term) analyses of variance (ANOVA) and Tukey’s HSD post hoc tests with a 95% confidence interval. Correlation analysis of the examined variables (BD, UD, SHC, CI, yield, qualitative parameters) was performed at a probability level of 0.05.

5. Conclusions

This study aimed to verify the effect of the biostimulant NeOsol (NS), the manure stabilizer Z’fix (ZF), farmyard manure (FM), and their combinations in real-life conditions of agricultural practice in a six-year experiment. First of all, these treated variants were supposed to reduce soil bulk density (BD), unit draft (UD), and cone index (CI). On the other hand, they should have increased saturated hydraulic conductivity (SHC), yield, and quality parameters. When all three treatments were joined, i.e., farmyard manure stabilised by Z’fix and combined with NeOsol (FM_ZF_NS), only UD decreased significantly and yield increased. The variant of farmyard manure stabilised by Z’fix (FM_ZF) proved solely to increase yield. Conventional farmyard manure combined with NeOsol (FM_NS) significantly decreased UD and CI at 10–20 cm and increased yield. Farmyard manure alone (FM) significantly reduced UD and CI at 10–20 cm and at 20–30 cm. On the contrary, SHC and yield augmented significantly with conventional farmyard manure application. The application of NeOsol significantly raised yield alone. None of the examined treatments proved any adverse effect either on soil or on plant-related variables. Additionally, they are likely to have improved other variables that were not monitored, e.g., ammonia emissions in animal housing, plant canopy cover, etc. A more environmentally sustainable agricultural system is unlikely to be harmful to the environment; it rather generates benefits and protects the land against possible risks.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/plants13070920/s1, Table S1: Means and standard deviations of vegetation indices: Normalized Difference Vegetation Index (NDVI), Moisture Stress Index (MSI), and Green Chlorophyll Index (GCI) for different variants of farmyard manure (FM), farmyard manure stabilized by Z’fix (FM_ZF), biostimulant NeOsol (NS), and control variant (CL) application in individual dates; lowercase letters denote significant differences at probability level of 0.05 within each date separately.

Author Contributions

Conceptualization, V.N., P.Š., and O.L.; methodology, V.N. and P.Š.; software, V.N. and P.Š.; validation, P.Š.; formal analysis, V.N.; investigation, V.N. and P.Š.; resources, V.N. and P.Š.; data curation, V.N. and P.Š.; writing—original draft preparation, V.N. and P.Š.; visualization, P.Š.; supervision, P.Š.; project administration, O.L.; funding acquisition, O.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Technology Agency of the Czech Republic, grant number TH02030169, and by the Czech University of Life Sciences, Faculty of Engineering, in the frame of the internal project IGA 2022: 31180/1312/3106.

Data Availability Statement

The data presented in this study are available from the corresponding author upon reasonable request.

Acknowledgments

The authors would like to thank all colleagues that took part in the research for their support and cooperation.

Conflicts of Interest

Oldřich Látal was employed by the company “Agrovyzkum Rapotin Ltd.”. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

  1. United Nations. World Population Prospects: The 2017 Revision, Key Findings and Advance Tables; United Nations: New York, NY, USA, 2017. Available online: https://population.un.org/wpp/Publications/Files/WPP2017_KeyFindings.pdf (accessed on 1 January 2024).
  2. FAO. How to Feed the World in 2050; FAO: Rome, Italy, 2009; Available online: https://www.fao.org/fileadmin/templates/wsfs/Summit/WSFS_Issues_papers/WSFS_feeding_E.pdf (accessed on 1 January 2024).
  3. Jongeneel, R.; Silvis, H.; Martinez, A.G.; Jager, J. The Green Deal: An Assessment of Impacts of the Farm to Fork and Biodiversity Strategies on the EU Livestock Sector; Wageningen Economic Research: Wageningen, The Netherlands, 2021. [Google Scholar]
  4. Cai, A.; Xu, M.; Wang, B.; Zhang, W.; Liang, G.; Hou, E.; Luo, Y. Manure Acts as a Better Fertilizer for Increasing Crop Yields than Synthetic Fertilizer Does by Improving Soil Fertility. Soil Tillage Res. 2019, 189, 168–175. [Google Scholar] [CrossRef]
  5. Weisser, M.; Mattner, S.W.; Southam-Rogers, L.; Hepworth, G.; Arioli, T. Effect of a Fortified Biostimulant Extract on Tomato Plant Productivity, Physiology, and Growing Media Properties. Plants 2023, 13, 4. [Google Scholar] [CrossRef] [PubMed]
  6. Soares, C.R.F.S.; Hernández, A.G.; da Silva, E.P.; de Souza, J.E.A.; Bonfim, D.F.; Zabot, G.L.; Ferreira, P.A.A.; Brunetto, G. Applications and Market of Micro-Organism-Based and Plant-Based Inputs in Brazilian Agriculture. Plants 2023, 12, 3844. [Google Scholar] [CrossRef]
  7. Bourgougnon, N.; Burlot, A.S.; Jacquin, A.G. Algae for Global Sustainability? 1st ed.; Elsevier Ltd.: Amsterdam, The Netherlands, 2021; Volume 100, ISBN 9780128216927. [Google Scholar]
  8. Findura, P.; Šindelková, I.; Rusinek, R.; Karami, H.; Gancarz, M.; Bartoš, P. Determination of the Influence of Biostimulants on Soil Properties and Field Crop Yields. Int. Agrophysics 2022, 36, 351–359. [Google Scholar] [CrossRef]
  9. Šindelková, I.; Badalíková, B.; Kubíková, Z. The Soil Biostimulant Usage Effect on Soil Properties in Dry Area. In Proceedings of the International Multidisciplinary Scientific GeoConference Surveying Geology and Mining Ecology Management; SGEM, 2019; Volume 19, pp. 561–568. Available online: https://www.sgem.org/index.php/elibrary-research-areas?view=publication&task=show&id=5766 (accessed on 1 January 2024).
  10. Šindelková, I.; Badalíková, B. Use of the Biostimulants in a Dry Area within South Moravia. Zahrad. Vědecká Příloha 2019, XVIII, 103–113. [Google Scholar]
  11. Rumeur, B. Le Trace Elements and Seaweeds Extracts, Tools for a Sustainable Management of Soil Functions and Plat Production. In Proceedings of the Agriculture in Nature and Environment Protection, Osijek, Republic of Croatia, 7–9 September 2020; Available online: https://www.hdpot.hr/images/files/Vukovar%20zbornici/Zbornik%20Agro%202020.pdf (accessed on 1 January 2024).
  12. Tuba, G.; Kovács, G.; Sinka, L.; Nagy, P.; Rivera-Garcia, A.; Bajusová, Z.; Findura, P.; Zsembeli, J. Effect of Soil Conditioning on Soil Penetration Resistance and Traction Power Demand of Ploughing. Agriculture 2021, 67, 113–123. [Google Scholar] [CrossRef]
  13. Spychalski, W.; Sulewska, H.; Ratajczak, K.; Kaczmarek, T. Composition of Soil Solution after 10 Years of PRP SOL Fertilization Based on the Selected Chemical Properties. J. Res. Appl. Agric. Eng. 2017, 62, 139–143. [Google Scholar]
  14. Swędrzyńska, D.; Zielewicz, W.; Swędrzyński, A. Comparison of Soil Bioconditioners and Standard Fertilization in Terms of the Impact on Yield and Vitality of Lolium Perenne and Soil Biological Properties. Open Life Sci. 2019, 14, 666–680. [Google Scholar] [CrossRef] [PubMed]
  15. Jukić, G.; Sunjić, K.; Varnica, I.; Gašo, D.; Labudović, B. Effect of different kinds of fertilizer on soybean yield. In Proceedings of the 8th International Scientific/Professional Conference, Agriculture in Nature and Environment Protection, Vukovar, Croatia, 1–3 June 2015; pp. 110–114. [Google Scholar]
  16. Sulewska, H.; Niewiadomska, A.; Majchrzak, L.; Panasiewicz, K. Potatoes Reaction on PRP SOL Fertilisation. J. Res. Appl. Agric. Eng. 2012, 57, 116–121. [Google Scholar]
  17. Ratajczak, K.; Sulewska, H.; Szymańska, G.; Wolna-Maruwka, A.; Faligowska, A. The Effect of Soil Type and Soil Additives on the Selected Growth Parameters and Yield of Flowerheads of Calendula officinalis L. Herba Pol. 2016, 62, 17–30. [Google Scholar] [CrossRef]
  18. Novák, V.; Šařec, P.; Křížová, K.; Novák, P.; Látal, O. Soil Physical Properties and Crop Status under Cattle Manure and Z’Fix in Haplic Chernozem. Plant Soil Environ. 2021, 67, 390–398. [Google Scholar] [CrossRef]
  19. Mozdzer, E.; Chudecka, J. Impact of Natural Fertilization Using PRP Fix on Some Soil Fertility Indicators. J. Ecol. Eng. 2017, 18, 137–144. [Google Scholar] [CrossRef]
  20. Látal, O.; Sedláčková, I.; Pozdíšek, J.; Farmyard Manure Quality Optimisation by Using the Z`fix Activator Application in Cattle Farms. Obilnářské Listy 2022. Available online: https://www.vukrom.cz/userfiles/files/obilnarske_listy/2022/obil_listy_3-4_2022.pdf (accessed on 1 January 2024).
  21. Sedláčková, I.; Látal, O.; Hambálková, M.; Polišenská, I.; Spring Barley Quality in Kroměříž Biological Transformed Organic Matter Trial. Obilnářské Listy 2021, 77–81. Available online: https://www.vukrom.cz/userfiles/files/obilnarske_listy/2021/Obil_listy_3_2021.pdf.pdf (accessed on 1 January 2024).
  22. Calvo, P.; Nelson, L.; Kloepper, J.W. Agricultural Uses of Plant Biostimulants. Plant Soil 2014, 383, 3–41. [Google Scholar] [CrossRef]
  23. Dhaliwal, S.S.; Naresh, R.K.; Mandal, A.; Walia, M.K.; Gupta, R.K.; Singh, R.; Dhaliwal, M.K. Effect of Manures and Fertilizers on Soil Physical Properties, Build-up of Macro and Micronutrients and Uptake in Soil under Different Cropping Systems: A Review. J. Plant Nutr. 2019, 42, 2873–2900. [Google Scholar] [CrossRef]
  24. Rayne, N.; Aula, L. Livestock Manure and the Impacts on Soil Health: A Review. Soil Syst. 2020, 4, 64. [Google Scholar] [CrossRef]
  25. Chen, X.; Xi, K.; Yang, Z.; Lu, J.; Zhang, Q.; Wang, B.; Wang, K.; Shi, J. Long-Term Increases in Continuous Cotton Yield and Soil Fertility Following the Application of Cotton Straw and Organic Manure. Agronomy 2023, 13, 2133. [Google Scholar] [CrossRef]
  26. Dudka, O.; Tanchyk, S.; Pavlov, O. Bulk Density of Soil Depending on the Farming Systems and Soil Tillage for Growing Spring Wheat in the Right Bank Forest-Steppe of Ukraine. Tech. Rom. J. Appl. Sci. Technol. 2022, 4, 31–37. [Google Scholar] [CrossRef]
  27. United States Department of Agriculture. Soil Bulk Density/Moisture/Aeration; United States Department of Agriculture: Washington, DC, USA, 2019. Available online: https://www.nrcs.usda.gov/sites/default/files/2022-10/Soil%20Bulk%20Density%20Moisture%20Aeration.pdf (accessed on 1 January 2024).
  28. Yu, F.; Zhao, S.; Zhao, Y.; Wang, Y.; Zhai, C.; Zhong, R.; Zhang, J.; Meng, Q. Long-Term Cattle Manure Application to Saline-Sodic Soil Increases Maize Yield by Decreasing Key Obstacle Factors in the Black Soil Region of Northeastern China. Int. J. Agric. Biol. Eng. 2023, 16, 176–183. [Google Scholar] [CrossRef]
  29. Čermáková, N.; Šařec, P.; Látal, O. Impact of Manure and Selected Conditionners on Physical Properties of Clay Soil. In Proceedings of the TAE 2019—Proceeding of 7th International Conference on Trends in Agricultural Engineering 2019, Prague, Czech Republic, 17–20 September 2019; Czech University of Life Sciences Prague: Prague, Czech Republic, 2019; pp. 93–98. [Google Scholar]
  30. Fu, Y.; De Jonge, L.W.; Moldrup, P.; Paradelo, M.; Arthur, E. Geoderma Improvements in Soil Physical Properties after Long-Term Manure Addition Depend on Soil and Crop Type. Geoderma 2022, 425, 116062. [Google Scholar] [CrossRef]
  31. Forgó, Z.; Tolvaly-Rosca, F.; Pásztor, J.; Kovari, A. Energy Consumption Evaluation of Active Tillage Machines Using Dynamic Modelling. Appl. Sci. 2021, 11, 6240. [Google Scholar] [CrossRef]
  32. Choudhary, S.; Upadhyay, G.; Patel, B.; Naresh; Jain, M. Energy Requirements and Tillage Performance Under Different Active Tillage Treatments in Sandy Loam Soil. J. Biosyst. Eng. 2021, 46, 353–364. [Google Scholar] [CrossRef]
  33. Urbanovičová, O.; Krištof, K.; Findura, P.; Mráz, M.; Jobbágy, J.; Križan, M. The Effect of Soil Conditioner on the Spatial Variability of Soil Environment. Agron. Res. 2018, 16, 2197–2210. [Google Scholar] [CrossRef]
  34. Šařec, P.; Látal, O.; Novák, P.; Holátko, J.; Novák, V.; Dokulilová, T.; Brtnický, M. Changes in Soil Properties and Possibilities of Reducing Environmental Risks Due to the Application of Biological Activators in Conditions of Very Heavy Soils. Agron. Res. 2020, 18, 2581–2591. [Google Scholar] [CrossRef]
  35. Mclaughlin, N.B.; Gregorich, E.G.; Dwyer, L.M.; Ma, B.L. Effect of Organic and Inorganic Soil Nitrogen Amendments on Mouldboard Plow Draft. Soil Tillage Res. 2002, 64, 211–219. [Google Scholar] [CrossRef]
  36. Liang, A.; Mclaughlin, N.B.; Ma, B.L.; Gregorich, E.G.; Morrison, M.J.; Burtt, S.D.; Patterson, B.S.; Evenson, L.I. Changes in Mouldboard Plough Draught and Tractor Fuel Consumption on Continuous Corn after 18 Years of Organic and Inorganic N Amendments. Energy 2013, 52, 89–98. [Google Scholar] [CrossRef]
  37. Rusu, T. Energy Efficiency and Soil Conservation in Conventional, Minimum Tillage and No-Tillage. Int. Soil Water Conserv. Res. 2014, 2, 42–49. [Google Scholar] [CrossRef]
  38. Moradi, E.; Miranda-Moreno, L. Vehicular Fuel Consumption Estimation Using Real-World Measures through Cascaded Machine Learning Modeling. Transp. Res. Part D Transp. Environ. 2020, 88, 102576. [Google Scholar] [CrossRef]
  39. Wilson, D.L.; Talancon, D.R.; Winslow, R.L.; Linares, X.; Gadgil, A.J. Avoided Emissions of a Fuel-Efficient Biomass Cookstove Dwarf Embodied Emissions. Dev. Eng. 2016, 1, 45–52. [Google Scholar] [CrossRef]
  40. Martelli, A.; Cachiba, S. Mitigating Measures Carried out by the Environment Secretariat of the Municipality of Itapira-SP in Relation to Greenhouse Gases. South Fla. J. Dev. 2020, 1, 11–20. [Google Scholar] [CrossRef]
  41. Mühlbachová, G.; Kusá, H.; Růžek, P.; Vavera, R. CO2 Emissions in a Soil under Different Tillage Practices. Plant Soil Environ. 2022, 68, 253–261. [Google Scholar] [CrossRef]
  42. Jiang, Q.; Cao, M.; Wang, Y.; Wang, J. Estimating Soil Penetration Resistance of Paddy Soils in the Plastic State Using Physical Properties. Agronomy 2020, 10, 1914. [Google Scholar] [CrossRef]
  43. Pardo, A.; Amato, M.; Chiarandà, F.Q. Relationships between Soil Structure, Root Distribution and Water Uptake of Chickpea (Cicer arietinum L.). Plant Growth and Water Distribution. Eur. J. Agron. 2000, 13, 39–45. [Google Scholar] [CrossRef]
  44. Otto, R.; Silva, A.P.; Franco, H.C.J.; Oliveira, E.C.A.; Trivelin, P.C.O. High Soil Penetration Resistance Reduces Sugarcane Root System Development. Soil Tillage Res. 2011, 117, 201–210. [Google Scholar] [CrossRef]
  45. Bengough, A.G.; McKenzie, B.M.; Hallett, P.D.; Valentine, T.A. Root Elongation, Water Stress, and Mechanical Impedance: A Review of Limiting Stresses and Beneficial Root Tip Traits. J. Exp. Bot. 2011, 62, 59–68. [Google Scholar] [CrossRef]
  46. Kuhwald, M.; Hamer, W.B.; Brunotte, J.; Duttmann, R. Soil Penetration Resistance after One-Time Inversion Tillage: A Spatio-Temporal Analysis at the Field Scale. Land 2020, 9, 482. [Google Scholar] [CrossRef]
  47. Celik, I.; Gunal, H.; Budak, M.; Akpinar, C. Effects of Long-Term Organic and Mineral Fertilizers on Bulk Density and Penetration Resistance in Semi-Arid Mediterranean Soil Conditions. Geoderma 2010, 160, 236–243. [Google Scholar] [CrossRef]
  48. Zeleke, T.B.; Si, B.C. Scaling Relationships between Saturated Hydraulic Conductivity and Soil Physical Properties. Soil Sci. Soc. Am. J. 2005, 69, 1691–1702. [Google Scholar] [CrossRef]
  49. Chyba, J.; Kroulik, M.; Kristof, K.; Misiewicz, P.A. The Influence of Agricultural Traffic on Soil Infiltration Rates. Agron. Res. 2017, 15, 664–673. [Google Scholar]
  50. Jadczyszyn, J.; Niedźwiecki, J. Relation of Saturated Hydraulic Conductivity to Soil Losses. Pol. J. Environ. Stud. 2005, 14, 431–435. [Google Scholar]
  51. Panagos, P.; Ballabio, C.; Himics, M.; Scarpa, S.; Matthews, F.; Bogonos, M.; Poesen, J.; Borrelli, P. Projections of Soil Loss by Water Erosion in Europe by 2050. Environ. Sci. Policy 2021, 124, 380–392. [Google Scholar] [CrossRef]
  52. The Ministry of Agriculture of the Czech Republic. Water Erosion. Available online: https://eagri.cz/public/portal/en/mze/soil (accessed on 1 January 2024).
  53. Bagarello, V.; Baiamonte, G.; Caia, C. Variability of Near-Surface Saturated Hydraulic Conductivity for the Clay Soils of a Small Sicilian Basin. Geoderma 2019, 340, 133–145. [Google Scholar] [CrossRef]
  54. Miller, J.J.; Sweetland, N.J.; Chang, C. Hydrological Properties of a Clay Loam Soil after Long-Term Cattle Manure Application. J. Environ. Qual. 2002, 31, 989–996. [Google Scholar] [CrossRef] [PubMed]
  55. Novák, V.; Šařec, P.; Křížová, K.; Novák, P.; Látal, O. Potential Impact of Biostimulator NeOsol and Three Different Manure Types on Physical Soil Properties and Crop Status in Heavy Soils Conditions. Sustainability 2021, 14, 438. [Google Scholar] [CrossRef]
  56. Pandey, P.C.; Pandey, M. Highlighting the Role of Agriculture and Geospatial Technology in Food Security and Sustainable Development Goals. Sustain. Dev. 2023, 31, 3175–3195. [Google Scholar] [CrossRef]
  57. Hernandez, T.; Berlanga, J.G.; Tormos, I.; Garcia, C. Organic versus Inorganic Fertilizers: Response of Soil Properties and Crop Yield. AIMS Geosci. 2021, 7, 415–439. [Google Scholar] [CrossRef]
  58. Krasilnikov, P.; Taboada, M.A. Amanullah Fertilizer Use, Soil Health and Agricultural Sustainability. Agriculture 2022, 12, 462. [Google Scholar] [CrossRef]
  59. Sulewska, H.; Szymańska, G.; Śmiatacz, K.; Koziara, W.; Niewiadomska, A. Effect of Prp Sol Application in Maize Grown for Grain. J. Res. Appl. Agric. Eng. 2013, 58, 161–166. [Google Scholar]
  60. Šařec, P.; Novák, V.; Křížová, K. Effect of organic fertilizers, biochar and other conditioners on modal luvisol. In Proceedings of the 7th International Conference on Trends in Agricultural Engineering, Prague, Czech Republic, 17–20 September 2019; pp. 494–499. [Google Scholar]
  61. Shynkaruk, L.; Lykhochvor, V. Influence of Fertilization and Foliar Feeding on Maize Grain Qualitative Indicators. Ukr. J. Ecol. 2021, 11, 113–116. [Google Scholar]
  62. Da Silva, P.R.F.; Strieder, M.L.; Coser, R.P.S.; Rambo, L.; Sangoi, L.; Argenta, G.; Forsthofer, E.L.; Silva, A.A. Grain Yield and Kernel Crude Protein Content Increases of Maize Hybrids with Late Nitrogen Side-Dressing. Sci. Agric. 2005, 62, 487–492. [Google Scholar] [CrossRef]
  63. Wang, H.; Zhang, X.; Zhang, G.; Fang, Y.; Hou, H.; Lei, K.; Ma, Y. Regulation of Density and Fertilization on Crude Protein Synthesis in Forage Maize in a Semiarid Rain-Fed Area. Agriculture 2023, 13, 715. [Google Scholar] [CrossRef]
  64. Litke, L.; Gaile, Z.; Ruža, A. Effect of Nitrogen Fertilization on Winter Wheat Yield and Yield Quality. Agron. Res. 2018, 16, 500–509. [Google Scholar] [CrossRef]
  65. Tsvey, Y.; Ivanina, R.; Ivanina, V.; Senchuk, S. Yield and Quality of Winter Wheat (Triticum aestivum L.) Grain in Relation to Nitrogen Fertilization. Rev. Fac. Nac. De Agron. Medellín 2021, 74, 9413–9422. [Google Scholar] [CrossRef]
  66. Klem, K.; Křen, J.; Šimor, J.; Kováč, D.; Holub, P.; Míša, P.; Svobodová, I.; Lukas, V.; Lukeš, P.; Findurová, H.; et al. Improving Nitrogen Status Estimation in Malting Barley Based on Hyperspectral Reflectance and Artificial Neural Networks. Agronomy 2021, 11, 2592. [Google Scholar] [CrossRef]
  67. Groth, D.A.; Sokólski, M.; Jankowski, K.J. A Multi-Criteria Evaluation of the Effectiveness of Nitrogen and Sulfur Fertilization in Different Cultivars of Winter Rapeseed—Productivity, Economic and Energy Balance. Energies 2020, 13, 4654. [Google Scholar] [CrossRef]
  68. Sulewska, H.; Ratajczak, K.; Niewiadomska, A.; Panasiewicz, K. Impact of Four-Year Period Application of Prp Sol and Prp Ebv on the Yield and Quality of Winter Oil Seed Rape. J. Res. Appl. Agric. Eng. 2016, 61, 176–181. [Google Scholar]
  69. Zielewicz, W.; Wróbel, B.; Stopa, W.; Paszkiewicz-Jasińska, A.; Szulc, P.; Klarzyńska, A. Effect of Soil Biostimulant Application on Protein and Carbohydrate Content in Per-Ennial Ryegrass (Lolium perenne L.). J. Res. Appl. Agric. Eng. 2023, 68, 36–48. [Google Scholar] [CrossRef]
  70. Możdżer, E.; Meller, E.; Sammel, A. The effect of organic fertilisation of liquid manure and the prp fix preparation on the yield and chemical composition of winter rape seeds and spring wheat grain. J. Ecol. Eng. 2017, 18, 139–146. [Google Scholar] [CrossRef]
  71. Borowiak, K.; Niewiadomska, A.; Sulewska, H.; Szymanska, G.; Gluchowska, K.; Wolna-Maruwka, A. Effect of PRP SOL and PRP EBV Nutrition on Yield, Photosynthesis Activity and Soil Microbial Activity of Three Cereal Species. Fresenius Environ. Bull. 2016, 25, 2026–2035. [Google Scholar]
  72. Šařec, P.; Novák, P.; Kumhálová, J. Impact of Activators of Organic Matter on Soil and Crop Stand Properties in Conditions of Very Heavy Soils. In Proceedings of the Engineering for Rural Development, Jelgava, Latvia, 24 May 2017; Volume 16, pp. 486–491. [Google Scholar]
  73. CSN EN ISO 17892-2; Geotechnický Průzkum a Zkoušení—Laboratorní Zkoušky Zemin—Část 2: Stanovení Objemové Hmotnosti. Available online: https://www.technicke-normy-csn.cz/csn-en-iso-17892-2-721007-218314.html (accessed on 1 January 2024).
  74. Elrick, D.E.; Reynolds, W.D.; Tan, K.A. Hydraulic Conductivity Measurements in the Unsaturated Zone Using Improved Well Analyses. Groundw. Monit. Remediat. 1989, 9, 184–193. [Google Scholar] [CrossRef]
  75. Darra, N.; Espejo-Garcia, B.; Kasimati, A.; Kriezi, O.; Psomiadis, E.; Fountas, S. Can Satellites Predict Yield? Ensemble Machine Learning and Statistical Analysis of Sentinel-2 Imagery for Processing Tomato Yield Prediction. Sensors 2023, 23, 2586. [Google Scholar] [CrossRef]
  76. Rouse, R.W.H.; Haas, J.A.W.; Deering, D.W. Monitoring Vegetation Systems in the Great Plains with ERTS. In Third Earth Resources Technology Satellite-1 Symposium- Volume I: Technical Presentations; NASA, 1974; pp. 309–317. Available online: https://ntrs.nasa.gov/citations/19740022614 (accessed on 1 January 2024).
  77. Huang, S.; Tang, L.; Hupy, J.P.; Wang, Y.; Shao, G. A Commentary Review on the Use of Normalized Difference Vegetation Index (NDVI) in the Era of Popular Remote Sensing. J. For. Res. 2021, 32, 1–6. [Google Scholar] [CrossRef]
  78. McVeagh, P.; Yule, I.; Grafton, M. Pasture Yield Mapping from Your Groundspread Truck. In Proceedings of the Advanced Nutrient Management: Gains from the Past—Goals for the Future; 2012; pp. 1–5. Available online: http://tur-www1.massey.ac.nz/~flrc/workshops/12/Manuscripts/McVeagh_2012.pdf (accessed on 1 January 2024).
  79. Tunca, E.; Köksal, E.S.; Taner, S.Ç. Silage Maize Yield Estimation by Using Planetscope, Sentinel-2A and Landsat 8 OLI Satellite Images. Smart Agric. Technol. 2023, 4, 100165. [Google Scholar] [CrossRef]
  80. HUNTJR, E.; ROCK, B. Detection of Changes in Leaf Water Content Using Near- and Middle-Infrared Reflectances. Remote Sens. Environ. 1989, 30, 43–54. [Google Scholar] [CrossRef]
  81. Welikhe, P.; Quansah, J.E.; Fall, S.; McElhenney, W. Estimation of Soil Moisture Percentage Using LANDSAT-Based Moisture Stress Index. J. Remote Sens. GIS 2017, 6, 1–5. [Google Scholar] [CrossRef]
  82. Gitelson, A.A.; Vina, A.; Arkebauer, T.J.; Rundquist, D.C.; Keydan, G.; Leavitt, B. Remote Estimation of Leaf Area Index and Green Leaf Biomass in Maize Canopies. Geophys. Res. Lett. 2003, 30, 4–7. [Google Scholar] [CrossRef]
  83. Kumar, V.; Sharma, A.; Bhardwaj, R.; Thukral, A.K. Comparison of Different Reflectance Indices for Vegetation Analysis Using Landsat-TM Data. Remote Sens. Appl. Soc. Environ. 2018, 12, 70–77. [Google Scholar] [CrossRef]
Figure 1. Values of soil physical characteristics: (a) bulk density (BD); (b) unit draft (UD); and (c) saturated hydraulic conductivity (SHC) relative to the respective control variant (CL) for different variants of farmyard manure (FM), farmyard manure stabilized by Z’fix (FM_ZF), and biostimulant NeOsol (NS) application in individual terms; vertical bars denote 0.95 confidence intervals.
Figure 1. Values of soil physical characteristics: (a) bulk density (BD); (b) unit draft (UD); and (c) saturated hydraulic conductivity (SHC) relative to the respective control variant (CL) for different variants of farmyard manure (FM), farmyard manure stabilized by Z’fix (FM_ZF), and biostimulant NeOsol (NS) application in individual terms; vertical bars denote 0.95 confidence intervals.
Plants 13 00920 g001
Figure 2. Values of cone index (CI) for different variants of farmyard manure (FM), farmyard manure stabilized by Z’fix (FM_ZF), and biostimulant NeOsol (NS) application in individual terms; vertical bars denote 0.95 confidence intervals.
Figure 2. Values of cone index (CI) for different variants of farmyard manure (FM), farmyard manure stabilized by Z’fix (FM_ZF), and biostimulant NeOsol (NS) application in individual terms; vertical bars denote 0.95 confidence intervals.
Plants 13 00920 g002
Figure 3. Average values of (a) yield and (b) given qualitative parameter relative to the respective control variant (CL) for different variants of farmyard manure (FM), farmyard manure stabilized by Z’fix (FM_ZF), and biostimulant NeOsol (NS) application in individual years; vertical bars denote 0.95 confidence intervals.
Figure 3. Average values of (a) yield and (b) given qualitative parameter relative to the respective control variant (CL) for different variants of farmyard manure (FM), farmyard manure stabilized by Z’fix (FM_ZF), and biostimulant NeOsol (NS) application in individual years; vertical bars denote 0.95 confidence intervals.
Plants 13 00920 g003
Figure 4. Change in (a) soil-related variables (unit draft, cone index) and (b) soil- and plant-related variables (bulk density, saturated hydraulic conductivity, yield, qualitative parameters) relative to the respective control variant (CL) for different variants of farmyard manure (FM), farmyard manure stabilized by Z’fix (FM_ZF), and biostimulant NeOsol (NS) application over the monitored period; lowercase letters denote significant differences at probability level of 0.05.
Figure 4. Change in (a) soil-related variables (unit draft, cone index) and (b) soil- and plant-related variables (bulk density, saturated hydraulic conductivity, yield, qualitative parameters) relative to the respective control variant (CL) for different variants of farmyard manure (FM), farmyard manure stabilized by Z’fix (FM_ZF), and biostimulant NeOsol (NS) application over the monitored period; lowercase letters denote significant differences at probability level of 0.05.
Plants 13 00920 g004
Figure 5. Unit draft (UD) measurement using the towing tractor New Holland T8030 and the implement Farmet Hurikán tine cultivator.
Figure 5. Unit draft (UD) measurement using the towing tractor New Holland T8030 and the implement Farmet Hurikán tine cultivator.
Plants 13 00920 g005
Table 1. Means and standard deviations of soil bulk density (BD), unit draft (UD), and saturated hydraulic conductivity (SHC) for different variants of farmyard manure (FM), farmyard manure stabilized by Z’fix (FM_ZF), and biostimulant NeOsol (NS) application in individual terms; lowercase letters denote significant differences at probability level of 0.05 within each term separately; numbers in bold indicate the most favourable value.
Table 1. Means and standard deviations of soil bulk density (BD), unit draft (UD), and saturated hydraulic conductivity (SHC) for different variants of farmyard manure (FM), farmyard manure stabilized by Z’fix (FM_ZF), and biostimulant NeOsol (NS) application in individual terms; lowercase letters denote significant differences at probability level of 0.05 within each term separately; numbers in bold indicate the most favourable value.
VariableVariantTerm ITerm IITerm III
Bulk density
(g·cm−3)
FM_ZF_NS1.517 ab ± 0.0251.285 a ± 0.0111.298 a ± 0.126
FM_ZF1.553 b ± 0.0351.374 a ± 0.0111.382 a ± 0.053
FM_NS1.483 ab ± 0.0421.327 a ± 0.0621.323 a ± 0.091
FM1.457 a ± 0.0211.304 a ± 0.0251.318 a ± 0.019
NS1.487 ab ± 0.0121.349 a ± 0.0691.381 a ± 0.020
CL1.450 a ± 0.0261.327 a ± 0.0241.415 a ± 0.026
Unit Draft
(kN·m−2)
FM_ZF_NS95.1950 b ± 10.50838.662 b ± 4.81587.699 b ± 8.954
FM_ZF105.784 c ± 7.37736.538 a ± 3.33793.372 c ± 10.549
FM_NS88.856 a ± 6.66238.730 b ± 4.64286.271 b ± 6.012
FM94.670 b ± 8.70636.089 a ± 4.76481.004 a ± 6.764
NS96.625 b ± 13.55437.882 ab ± 2.80994.208 c ± 7.613
CL97.755 b ± 7.85238.773 b ± 3.60794.812 c ± 6.999
Saturated
hydraulic
conductivity
(mm·h−1)
FM_ZF_NS5.324 ab ± 3.48172.178 b ± 35.34418.816 a ± 21.883
FM_ZF5.451 ab ± 2.82925.238 a ± 17.72525.033 a ± 20.985
FM_NS3.355 a ± 1.88931.189 a ± 21.1986.546 a ± 3.372
FM5.589 ab ± 2.84852.283 ab ± 37.83152.630 b ± 26.385
NS3.859 a ± 1.70315.206 a ± 11.37415.736 a ± 8.318
CL6.748 b ± 3.00829.555 a ± 42.19817.080 a ± 15.432
Table 2. Means and standard deviations of cone index (CI) for different variants of farmyard manure (FM), farmyard manure stabilized by Z’fix (FM_ZF), biostimulant NeOsol (NS), and control variant (CL) application in individual terms; lowercase letters denote significant differences at probability level of 0.05 within each term separately, numbers in bold indicate the most favourable value.
Table 2. Means and standard deviations of cone index (CI) for different variants of farmyard manure (FM), farmyard manure stabilized by Z’fix (FM_ZF), biostimulant NeOsol (NS), and control variant (CL) application in individual terms; lowercase letters denote significant differences at probability level of 0.05 within each term separately, numbers in bold indicate the most favourable value.
Depth ZoneVariantCone Index (MPa)
Term ITerm IITerm III
0–10 cmFM_ZF_NS0.515 a ± 0.1181.040 a ± 0.3380.380 a ± 0.083
FM_ZF0.590 a ± 0.1680.950 a ± 0.2280.345 a ± 0.110
FM_NS0.525 a ± 0.1770.885 a ± 0.2430.363 a ± 0.076
FM0.605 a ± 0.1610.970 a ± 0.2360.400 a ± 0.103
NS0.550 a ± 0.1610.965 a ± 0.2810.400 a ± 0.086
CL0.550 a ± 0.2541.040 a ± 0.2410.380 a ± 0.062
10–20 cmFM_ZF_NS1.043 a ± 0.3581.403 b ± 0.2890.727 ab ± 0.328
FM_ZF1.437 b ± 0.4661.270 ab ± 0.1840.770 ab ± 0.518
FM_NS1.140 ab ± 0.3581.197 a ± 0.2480.593 a ± 0.319
FM1.237 ab ± 0.3811.223 ab ± 0.2900.627 ab ± 0.236
NS1.253 ab ± 0.4361.303 ab ± 0.2470.637 ab ± 0.346
CL1.280 ab ± 0.4451.330 ab ± 0.2600.903 b ± 0.495
20–30 cmFM_ZF_NS1.985 a ± 0.6312.055 a ± 0.5511.430 a ± 0.571
FM_ZF2.575 b ± 0.5881.875 a ± 0.4301.860 a ± 0.737
FM_NS2.210 ab ± 0.4121.930 a ± 0.5421.611 a ± 0.768
FM2.185 ab ± 0.4851.605 a ± 0.4541.490 a ± 0.620
NS2.485 ab ± 0.7121.915 a ± 0.3821.665 a ± 0.595
CL 2.445 ab ± 0.7172.030 a ± 0.5911.840 a ± 0.426
Table 3. Means and standard deviations of yield from standard dry matter (DM) content and of a given qualitative parameter for different variants of farmyard manure (FM), farmyard manure stabilized by Z’fix (FM_ZF), biostimulant NeOsol (NS), and control variant (CL) application in individual years; lowercase letters denote significant differences at probability level of 0.05 within each year separately, numbers in bold indicate the most favourable value..
Table 3. Means and standard deviations of yield from standard dry matter (DM) content and of a given qualitative parameter for different variants of farmyard manure (FM), farmyard manure stabilized by Z’fix (FM_ZF), biostimulant NeOsol (NS), and control variant (CL) application in individual years; lowercase letters denote significant differences at probability level of 0.05 within each year separately, numbers in bold indicate the most favourable value..
Year:2015 Term I20162017 Term II201820192020 Term III
Crop:Silage maizeWinter wheatWinter barleyOilseed rapeWinter wheatWinter wheat
Standard moisture (% DM):38.014.514.59.014.514.5
VariantYield at standard moisture (t ha−1)
FM_ZF_NS35.200 a ± 3.6817.870 a ± 0.2147.780 a ± 0.7234.040 b ± 0.4127.860 c ± 0.6027.380 b ± 0.478
FM_ZF33.500 a ± 3.3417.850 a ± 0.1887.660 a ± 0.4403.990 b ± 0.3147.620 bc ± 0.4147.320 b ± 0.672
FM_NS33.800 a ± 2.3567.680 a ± 0.4157.440 a ± 0.3813.660 ab ± 0.2597.390 abc ± 0.5177.130 b ± 0.406
FM32.000 a ± 2.3147.630 a ± 0.5407.410 a ± 0.4583.400 ab ± 0.1307.300 abc ± 0.5096.950 b ± 0.319
NS31.800 a ± 0.9357.310 a ± 0.6667.200 a ± 0.7763.440 ab ± 0.3266.660 ab ± 0.4107.020 b ± 0.533
CL30.100 a ± 1.2897.230 a ± 0.4566.580 a ± 0.4033.200 a ± 0.2646.400 a ± 0.3635.780 a ± 0.174
Qualitative parameter and its value (% DM)
Crude proteinCrude proteinCrude proteinCrude fatCrude proteinCrude protein
FM_ZF_NS8.472 a ± 0.54812.900 b ± 0.98710.200 a ± 0.82943.000 a ± 3.17813.600 b ± 0.8769.700 a ± 0.769
FM_ZF8.493 a ± 0.51012.200 ab ± 0.91210.100 a ± 0.89644.830 a ± 2.48012.800 ab ± 1.0879.800 a ± 0.685
FM_NS8.361 a ± 0.76312.300 ab ± 1.08210.900 a ± 0.27343.120 a ± 4.14612.800 ab ± 0.6759.700 a ± 0.628
FM8.976 a ± 0.43811.900 ab ± 0.69710.200 a ± 0.80943.920 a ± 1.06713.100 b ± 1.0059.700 a ± 0.823
NS8.509 a ± 0.54212.600 ab ± 0.7379.500 a ± 0.67143.770 a ± 2.77111.100 a ± 0.7529.500 a ± 0.747
CL8.575 a ± 0.54510.700 a ± 0.78210.100 a ± 0.51243.990 a ± 2.10312.500 ab ± 0.7989.400 a ± 0.494
Table 4. Correlation matrix of Pearson correlation coefficients (r) and probabilities (p) for soil bulk density (BD), unit draft (UD), saturated hydraulic conductivity (SHC), cone index (CI) at three depth zones, yield, and quality parameters; * denotes significant differences at probability level of 0.05.
Table 4. Correlation matrix of Pearson correlation coefficients (r) and probabilities (p) for soil bulk density (BD), unit draft (UD), saturated hydraulic conductivity (SHC), cone index (CI) at three depth zones, yield, and quality parameters; * denotes significant differences at probability level of 0.05.
BDUDSHCCI atYieldQuality
0–10 cm10–20 cm0–30 cm Parameter
BD r = 0.6008 *r = −0.4802 *r = −0.0700r = 0.6366 *r = 0.4597r = −0.5221 *r = −0.4951 *
p = 0.0084 p = 0.0437 p = 0.783p = 0.0045p = 0.0549p = 0.0263p = 0.0367
UDr = 0.6008 * r = −0.3177r = 0.1233r = 0.6863 *r = 0.7554 *r = −0.4133r = −0.1499
p = 0.0084 p = 0.1989p = 0.6260p = 0.0017p = 0.0003p = 0.0882p = 0.5527
SHCr = −0.4802 *r = −0.3177 r = 0.1612r = −0.1181r = −0.1903r = 0.2672r = 0.3161
p = 0.0437 p = 0.1989 p = 0.5229p = 0.6407p = 0.4495p = 0.2837p = 0.2013
CI atr = −0.0700r = 0.1233r = 0.1612 r = 0.0899r = 0.0637r = −0.1855r = −0.0579
0–10 cmp = 0.783p = 0.6260p = 0.5229 p = 0.7227p = 0.8018p = 0.4612p = 0.8196
10–20 cmr = 0.6366 *r = 0.6863 *r = −0.1181r = 0.0899 r = 0.6730 *r = −0.6607 *r = −0.3527
p = 0.0045 p = 0.0017 p = 0.6407p = 0.7227 p = 0.0022p = 0.0028p = 0.1512
20–30 cmr = 0.4597r = 0.7554 *r = −0.1903r = 0.0637r = 0.6730 * r = −0.4737 *r = −0.1719
p = 0.0549p = 0.0003 p = 0.4495p = 0.8018p = 0.0022 p = 0.0470p = 0.4952
Yieldr = −0.5221 *r = −0.4133r = 0.2672r = −0.1855r = −0.6607 *r = −0.4737 * r = 0.3516
p = 0.0263 p = 0.0882p = 0.2837p = 0.4612p = 0.0028p = 0.0470 p = 0.1525
Qualityr = −0.4951 *r = −0.1499r = 0.3161r = −0.0579r = −0.3527r = −0.1719r = 0.3516
parameterp = 0.0367 p = 0.5527p = 0.2013p = 0.8196p = 0.1512p = 0.4952p = 0.1525
Table 5. Chemical analysis of farmyard manure treated by Z’fix (FM_ZF) and of conventional farmyard manure (FM).
Table 5. Chemical analysis of farmyard manure treated by Z’fix (FM_ZF) and of conventional farmyard manure (FM).
VariantDry Matter %N kg t−1C:NP2O5K2OCaOMgOpH
FM_ZF237.115.1:14.27.96.441.999.5
FM235.922.3:13.56.34.901.608.9
Table 6. Chemical and physical soil analysis undertaken on 13 August 2014.
Table 6. Chemical and physical soil analysis undertaken on 13 August 2014.
Soil Profile Depth (m)Unit
0.00–0.300.30–0.65
Clay (<0.002 mm)1521% w/w
Silt (0.002–0.05 mm)6667% w/w
Very fine sand (0.05–0.10 mm)31% w/w
Fine sand (0.10–0.25 mm)1611% w/w
Bulk density1.561.52g·cm−3
Porosity41.9742.64% w/w
Hummus content1.810.58% w/w
Cation exchange capacity110120mmol·kg−1
Volumetric moisture31.523.7% v/v
pH (KCl)4.995.09
Table 7. Crop rotation and field management.
Table 7. Crop rotation and field management.
SeasonCropVarietySowing DateHarvest DateFarmyard Manure (t ha−1)NeOsol (kg ha−1)
2014Winter wheat---50200
2015 Term IMaizeSY Kairo27 April 20159 September 2015-200
2016Winter wheatTurandot7 October 201517 August 201630200
2017 Term IIBarleyTitus20 October 201619 July 2017-150
2018Rape seedInspiration15 August 201716 July 2018-150
2019Winter wheatReform22 September 20189 August 201930150
2020 Term IIIWinter wheatReform24 September 20199 August 2020-150
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Novák, V.; Šařec, P.; Látal, O. Effect of Biostimulant, Manure Stabilizer, and Manure on Soil Physical Properties and Vegetation Status. Plants 2024, 13, 920. https://doi.org/10.3390/plants13070920

AMA Style

Novák V, Šařec P, Látal O. Effect of Biostimulant, Manure Stabilizer, and Manure on Soil Physical Properties and Vegetation Status. Plants. 2024; 13(7):920. https://doi.org/10.3390/plants13070920

Chicago/Turabian Style

Novák, Václav, Petr Šařec, and Oldřich Látal. 2024. "Effect of Biostimulant, Manure Stabilizer, and Manure on Soil Physical Properties and Vegetation Status" Plants 13, no. 7: 920. https://doi.org/10.3390/plants13070920

APA Style

Novák, V., Šařec, P., & Látal, O. (2024). Effect of Biostimulant, Manure Stabilizer, and Manure on Soil Physical Properties and Vegetation Status. Plants, 13(7), 920. https://doi.org/10.3390/plants13070920

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop