Effect of Forage Plant Mixture and Biostimulants Application on the Yield, Changes of Botanical Composition, and Microbiological Soil Activity

: Recently, an increasing interest in such fertilizers and fertilization methods which not only directly supply nutrients to plants, but also stimulate soil bioactivity is noted. Their effect on both soil microbiota and forage plants has not been fully recognized. The aim of the study was to investigate the combined effect of forage plant mixture type and mineral fertilizers (NPK) with biostimulants based on a marine algae extracts on the botanical composition, yield, the structure of selected taxonomic and trophic groups of soil microorganisms, and the soil enzymatic activity. During the years 2018–2019 a ﬁeld experiment established in split-plot design with two different forage plant mixtures, as a ﬁrst factor, and different fertilization basing on mineral fertilizers amended with biostimulants, as a second factor was conducted. Two types of forage mixtures of sown species were used: grass mixture (GM) and legume-grass mixture (LGM). Every year the following biostimulants were applied: N-14, PinKstart, Physiostart, Physioactive and they were compared with standard NPK fertilisation and no fertilisation as a control. The reaction of forage plant mixtures on applied fertilisation was different. The intensive development of grass species, mainly Lolium perenne , at the expense of Trifolium repens share in LGM was observed. In GM sward dominated Dactylis glomerata . A beneﬁcial effects of biostimulants’ application on the biomass yields of both grass mixtures was observed. The systematic soil acidiﬁcation and a decrease of soil enzymatic activity in result of applied fertilization, except NPK + Physioactive treatment (calcium fertilizer containing 76% calcium carbonate), was noted. Soil reaction to applied fertilisation was dependent on the botanical composition of the sward. The counts of microorganisms in the soil under LGM were almost two times higher than in the soil under GM. The most effective, in reducing the negative effect of nitrogen mineral fertilization on the pH of soil, was fertilization with NPK + Physioactiv.


Introduction
One of the main factors determining the productivity of all agricultural lands, including grasslands, is fertilization, especially mineral fertilization. The main goal of fertilization is to obtain a yield of optimal quality and quantity by supplying various micro-and macronutrients. Common mineral fertilizers such as ammonium nitrate, potassium salt and superphosphate, as well as natural fertilizers such as liquid manure and slurry are applied to the grassland. Fertilizers modify the botanical composition of the sward in The work hypothesis was that biostimulants increase soil microbial activity. These effects can be direct and indirect (soil pH, modification of biomass yield and botanical composition of the sward).
The aim of the study was to investigate how the supplementation of mineral fertilization (NPK) with biostimulants based on a marine algae extract influenced a sward with two types of forage plants mixtures (grasses, and legumes with grasses). The effect of the biostimulants on the botanical composition of the sward, yield, the structure and counts of selected taxonomic and trophic groups of soil microorganisms, and the soil enzymatic activity were investigated.

Study Site and Weather Conditions
The It was located on soil classified into an order of lessive soils, subgroup of typic lessive soils, family of boulder clay, and series of light and heavy loamy sands. According to international World Reference Base [43] classification, the soil was included in Albic Luvisols, and according to Soil Taxonomy in Typic Hapludalfs, regarding granulation in loamy sand underlined by loam [44]. The soil of the experimental field was classified as bonitation class IIIb-IVa, of a very good rye complex. In terms of physic-chemical properties, the soil was characterized by the following parameters: mildly acidic reaction of a soil (pH KCl = 6.5), 1.24% content of humus, 16% share of floated parts. Content of available nutrients in the topsoil, measured each year before fertilization, was high for phosphorus-78 mg P·kg −1 soil, medium for potassium-132 mg K·kg −1 soil (double lactate method), and medium for magnesium-54 mg Mg·kg −1 soil (Schachtschabel method).

Experimental Design
A two-factor experiment in split plot design in triplicate with two different forage plants mixtures, as a first factor, and the different fertilization basing on mineral fertlizers amended with biostimulants, as a second factor was conducted.
The first experimental factor was type of forage plants mixture differing with botanical compositions: grass mixture (GM) and legume-grass mixture (LGM).
Phosphorus-potassium fertilization with granulated triple superphosphate 46% P 2 O 5 and potassium salt containing 60% of K 2 O was applied once a year, in spring, before vegetation in every variant of fertilization in doses of 80 kg·ha −1 P and 80 kg·ha −1 K, apart from an absolute control (unfertilised).

Characteristics of Fertilizers with Biostimulants
N-14-Timac Agro Poland (Roullier Group France) nitrogen fertilizer contained 14% of nitrogen (7% in an ammonia form and 7%-in an amide form), 22% of CaO, 2% of MgO, 28% of SO 3 and Pheoflore complex. The complex is based on the extract of marine algae, rich in carbohydrates and polypeptides and, according to the producer, it enhances growth of bacterial biota, which may allow for the increase in nitrogen and phosphorus available in soil.
PinKstart ® -Timac Agro Poland (Roullier Group France) is the starter fertilizer, containing 48% of CaCO 3 , 4.5% of SO 3 , 28% of P 2 O 5 , 5% of K 2 O and Physio+ complex based on the extract of marine algae, rich in carbohydrates and polypeptides. In this additive, cooperation of aminopurine and Mezolac (highly reactive calcium carbonate) was used. PinKstart stimulates intensive growth of plants' root systems (especially of root hair).
Physioactiv (the current name is Physiomax 975 physio+)-Timac Agro Poland (Roullier Group France) is a calcium fertilizer which contains 76% of CaCO 3 in the form of Mezocalc, 3% of MgO and aminopurine-a natural extract from marine algae which provides bio-stimulation of plants' root systems.
According to the producer, starter fertilizers application results in more effective intake of calcium and intensive development of a root system from the beginning of vegetation. It also prevents plants from the effect of "a lazy root" as plants build a root system despite favourable conditions connected with an optimal availability of nutrients [45,46].

Botanical Composition and Biomass Yield
The botanical composition of the sward was assessed three times during the growing season (before the harvest of each regrowth) with the botanical-weight method developed by Stebler and Schröter and modified by Filipek [47] and Novak [48]. With this method, the percentage share of each species in the sward was assessed.

Soil Sampling and Analysis
In the third (2018) and fourth year (2019) of sward utilisation soil samples were collected for microbiological analysis and pH measurement. The samples were collected two times, on the dates of the harvest of the first (May 2018, June 2019) and third (October 2018 and 2019) of regrowth of the sward. The samples were collected with a soil sampler probe (20 punctures per plot), from the top layer of the soil profile (0-15 cm). In total, 42 soil samples were collected each sampling.
Microbiological analyses were conducted to determine the counts of selected groups of soil microorganisms (the total counts of heterotrophic bacteria, oligotrophic bacteria, copiotrophic bacteria, and fungi) with the pour plate method on appropriate agar mediums. Five replications of each analysis were made. The mean count of colonies was expressed as colony forming units (CFU)·g −1 DM of soil.
The total count of heterotrophic bacteria was determined on nutrient agar (3.0 g yeast extract; 5.0 g peptone from casein; 5.0 g sodium chloride; 12.0 g agar, 1.0 dm 3 H 2 O) [49], after 5-6 days of incubation at the temperature of 28 • C. Oligotrophic bacteria were counted on diluted nutrient agar (0.1 g peptone, 0.1 g beef extract, 0.05 g sodium chloride, 20.0 g agar, 1 dm 3 H 2 O) after 21 incubation days at 28 • C [50]. Copiotrophic bacteria were determined on nutrient broth medium (10.0 g peptone, 10.0 g beef extract, 5.0 g sodium chloride, 20.0 g agar, 1.0 dm 3 H 2 O) after 7 days of incubation at 28 • C [50]. Number of fungi (yeasts and moulds) was determined using Martin agar (1.0 g KH 2 PO 4 , 0.5 g MgSO 4 , 5.0 g peptone, 10.0 g glucose, 3.3 mL Rose bengal, 0.1 g chlortetracycline, 25.0 g agar, 1.0 dm 3 H 2 O), counted after 5 days of incubation at 28 • C [51]. The mean number of colonies was converted into soil dry matter on the basis of used dilution of soil solution and moisture of the soil sample.
The activity of acid (ACP) and alkaline phosphatases (ALP) was determined with pNPP (p-nitrophenyl phosphate) solutions (pH 6.5 for ACP and pH 11.0 for ALP) as buffered substrates. After 1-h incubation at 37 • C p-NP (p-nitrophenol) was produced. It was extracted and stained with sodium hydroxide. Its amount was measured spectrophotometrically at a wavelength of 400 nm. The enzyme activity was expressed as µmol pNP · g −1 DM of soil · h −1 [53].
The urease activity (URE) was also determined spectrophotometrically, using urea as a substrate. After 1-h incubation at 37 • C the amount of non-hydrolysed urea was measured at a wavelength of 410 nm. The enzyme activity was expressed as µg N-NH 4 + ·g −1 DM of soil·18 h −1 [54].
Catalase activity (CAT) in the soil was determined by means of titration (permanganometry). The soil with 0.3% H 2 O 2 solution was incubated for 20 min (temp. 20 • C) and then 1.5 M H 2 SO 4 was added. The resulting solution was titrated with 0.02 M KMnO 4 . The catalase activity was expressed as µmol H 2 O 2 ·g −1 DM of soil·min −1 [55]. The biological index of fertility (BIF) was also determined on the basis of the dehydrogenase and catalase activities [56].

Statistical Analyses
The results were tested by using multifactorial analyses of variance (ANOVA). Mean separations were made for significant effects with Tukey tests at the probability of p ≤ 0.05. Principal Component Analysis (PCA) [57], and Pearson correlation coefficient, was used to illustrate the relationship between the variables. All statistical analyses were carried out with the Statistica 13.1 software.

Weather Condition
The weather data for characteristic of weather conditions in years of study came from the measurement station located in Brody. The local climate is classified as temperate. It is characterised by mutual influence of maritime and continental climates. It has a seasonal, transitional, and changeable nature, especially in summer [58]. The average air temperature for the period 1960-2017 was 8.4 • C, while the average rainfall was 603.1 mm·m 2 . It was noted that the weather conditions were diverse in the study period, especially regarding the amount and distribution of rainfall [ Figure 1]. In the first year (2018) mean air temperature was 9.8 • C and was higher by 1.4 • C than the long-term average temperature. Annual sum of precipitation in 2018 was 444.1 mm·m 2 and was lower by 161.7 mm·m 2 than sum from the long-term period. The amount and distribution of rainfall in the growing season 2019 was more favourable for plant growth compared to the previous year.
Catalase activity (CAT) in the soil was determined by means of titration (permanganometry). The soil with 0.3% H2O2 solution was incubated for 20 min (temp. 20 °C) and then 1.5 M H2SO4 was added. The resulting solution was titrated with 0.02 M KMnO4. The catalase activity was expressed as µmol H2O2·g −1 DM of soil·min −1 [55].
The biological index of fertility (BIF) was also determined on the basis of the dehydrogenase and catalase activities [56].

Statistical Analyses
The results were tested by using multifactorial analyses of variance (ANOVA). Mean separations were made for significant effects with Tukey tests at the probability of p ≤ 0.05. Principal Component Analysis (PCA) [57], and Pearson correlation coefficient, was used to illustrate the relationship between the variables. All statistical analyses were carried out with the Statistica 13.1 software.

Weather Condition
The weather data for characteristic of weather conditions in years of study came from the measurement station located in Brody. The local climate is classified as temperate. It is characterised by mutual influence of maritime and continental climates. It has a seasonal, transitional, and changeable nature, especially in summer [58]. The average air temperature for the period 1960-2017 was 8.4 °C, while the average rainfall was 603.1 mm·m 2 . It was noted that the weather conditions were diverse in the study period, especially regarding the amount and distribution of rainfall [ Figure 1]. In the first year (2018) mean air temperature was 9.8 °C and was higher by 1.4 °C than the long-term average temperature. Annual sum of precipitation in 2018 was 444.1 mm·m 2 and was lower by 161.7 mm·m 2 than sum from the long-term period. The amount and distribution of rainfall in the growing season 2019 was more favourable for plant growth compared to the previous year.

Botanical Composition
Lolium perenne developed intensively in all the plots with LGM. The highest share of this grass species was observed in the biomass sward fertilized with NPK mineral fertilizers (41.0% in 2018 and 42% in 2019) and NPK + N-14 fertilizer applied at dose of 300 kg·ha −1 (40.6 and 42.0%) and NPK + N-14 fertilizer applied at dose of kg·ha −1 900 kg·ha −1 (34.4 and 41.0%).
The share of this species at the other treatments ranged from 24.3-27.3% (control) to about 38% (NPK + Physiostart). Among the other grasses sown in this mixture, Festuca arundinacea had a significant share-in 2019 it was 24.6% at the plots fertilised with NPK + PinKstart. Festuca pratensis was less abundant-its share ranged from 1.3% (control) Agronomy 2021, 11, 1786 7 of 22 to 7.6% (NPK + PinKstart). The share of the other grass species sown in the mixture (Lolium westerwoldicum and Phleum pratense) amounted to only about a few per cent. Among the unsown, spontaneously occurring species in the sward, there was a relatively high share of Dactylis glomerata, which appeared in the sward fertilised with: NPK (7.0-10.3%), NPK + Physioactiv (3.3-7.3%), and NPK + PinKstart (3.3-8.6%). There were also other unsown grass species appearing spontaneously in the sward of most of the plots, Bromus inermis being the dominant species. Although grasses were dominant on the majority of treatments there were some treatments where Trifolium repens had a considerable share in the yielded biomass.
In 2019 the share of this species in the sward of the control, where no fertilizers had been applied, was as high as 52%. Among the fertilised plots, the highest share of Trifolium repens (nearly 40%) was observed in the sward treated with the calcium fertilizer combined with the biostimulant, i.e., NPK + Physioactiv. The share of Trifolium repens was much lower in the other treatments ( Figure 2). The lowest share of this species, i.e., only about a dozen or so percent, was found in the sward where standard NPK fertilizers had been applied and in the combinations fertilised with NPK + Physiostart and NPK + PinKstart. Meadow herbs also appeared in the LGM sward in most of the fertilizer combinations. The dominant herb species was Taraxacum officinale. Its share was the highest in the sward fertilised with NPK + Physiostart, but it was absent from the swards fertilised with the N-14 at a dose of 300 kg·ha −1 and with the calcium fertilizer NPK + Physioactiv.
The share of this species at the other treatments ranged from 24.3-27.3% (control) to about 38% (NPK + Physiostart). Among the other grasses sown in this mixture, Festuca arundinacea had a significant share-in 2019 it was 24.6% at the plots fertilised with NPK + PinKstart. Festuca pratensis was less abundant-its share ranged from 1.3% (control) to 7.6% (NPK + PinKstart). The share of the other grass species sown in the mixture (Lolium westerwoldicum and Phleum pratense) amounted to only about a few per cent. Among the unsown, spontaneously occurring species in the sward, there was a relatively high share of Dactylis glomerata, which appeared in the sward fertilised with: NPK (7.0-10.3%), NPK + Physioactiv (3.3-7.3%), and NPK + PinKstart (3.3-8.6%). There were also other unsown grass species appearing spontaneously in the sward of most of the plots, Bromus inermis being the dominant species. Although grasses were dominant on the majority of treatments there were some treatments where Trifolium repens had a considerable share in the yielded biomass.
In 2019 the share of this species in the sward of the control, where no fertilizers had been applied, was as high as 52%. Among the fertilised plots, the highest share of Trifolium repens (nearly 40%) was observed in the sward treated with the calcium fertilizer combined with the biostimulant, i.e., NPK + Physioactiv. The share of Trifolium repens was much lower in the other treatments ( Figure 2). The lowest share of this species, i.e., only about a dozen or so percent, was found in the sward where standard NPK fertilizers had been applied and in the combinations fertilised with NPK + Physiostart and NPK + PinKstart. Meadow herbs also appeared in the LGM sward in most of the fertilizer combinations. The dominant herb species was Taraxacum officinale. Its share was the highest in the sward fertilised with NPK + Physiostart, but it was absent from the swards fertilised with the N-14 at a dose of 300 kg·ha −1 and with the calcium fertilizer NPK + Physioactiv. Year 2018 Year 2019

Lolium perenne Dactylis glomerata Festuca pratensis Festuca arundinacea
Other grasses Trifolium repens Herbs In the third and fourth years after sowing of GM Dactylis glomerata was the dominant grass species in many of the experimental plots, except for the control plot ( Figure 3). In 2018 the share of this species in the sward of the fertilized plots ranged from 53% (NPK + Physioactiv) to 67% (NPK + Physiostart). In 2019 it was even greater, i.e., from 58% in the sward fertilised with the NPK + Physioactiv to 74% in the plots treated with the NPK + Physiostart fertilizers. On the control plots, where no fertilizers had been applied, the share of Dactylis glomerata in the sward was about a dozen or so per cent. The second most numerous species in the sward of this mixture was Lolium perenne. The highest share of this species, i.e., over 60%, was found in the sward of the control plot. In the other plots, which had been treated with the fertilizers, its share was lower-it ranged from 18% in the plots treated with NPK + Physiostart to 32% in the ones fertilised with NPK + Physioactiv. The share of the other two species sown in the mixture, i.e., Festuca pratensis and Phleum pratense, was low or they were not found in the sward. Apart from the sown species, the sward of all the plots had some spontaneously occurring species of dicotyledons such as Trifolium repens and other species commonly classified as weeds (Capsella bursa-pastoris and Taraxacum officinale).
sward fertilised with the NPK + Physioactiv to 74% in the plots treated with the NPK + Physiostart fertilizers. On the control plots, where no fertilizers had been applied, the share of Dactylis glomerata in the sward was about a dozen or so per cent. The second most numerous species in the sward of this mixture was Lolium perenne. The highest share of this species, i.e., over 60%, was found in the sward of the control plot. In the other plots, which had been treated with the fertilizers, its share was lower-it ranged from 18% in the plots treated with NPK + Physiostart to 32% in the ones fertilised with NPK + Physioactiv. The share of the other two species sown in the mixture, i.e., Festuca pratensis and Phleum pratense, was low or they were not found in the sward. Apart from the sown species, the sward of all the plots had some spontaneously occurring species of dicotyledons such as Trifolium repens and other species commonly classified as weeds (Capsella bursapastoris and Taraxacum officinale). LGM swards), other grasses (representing mid-and tallgrasses, understood as all grasses without Lolium perenne) and Trifolium repens (the only representative of the legumes). The significance of the influence of the investigated factors on the botanical composition of the sward was assessed regarding the share of groups shown above. The analysis of variance (Table 1) showed a significant influence of the applied experimental factors (sown mixture of seeds and type of fertilization) and the effect of interactions between them in relation to all analyzed variables characterizing the botanical composition of the sward. Year 2018 Year 2019

Lolium perenne Dactylis glomerata Festuca pratensis Phleum pratense
Other grasses Trifolium repens Herbs The starting point for further analysis of the research results was the distinction within the botanical composition dominant taxonomic-functional groups of species, characterized by specific biological and utility properties and fulfilling specific ecological functions in the biocenosis: Lolium perenne (the only representative of shortgrasses, distinguished by a large share in the sward of all experimental combinations of both-GM and LGM swards), other grasses (representing mid-and tallgrasses, understood as all grasses without Lolium perenne) and Trifolium repens (the only representative of the legumes). The significance of the influence of the investigated factors on the botanical composition of the sward was assessed regarding the share of groups shown above. The analysis of variance (Table 1) showed a significant influence of the applied experimental factors (sown mixture of seeds and type of fertilization) and the effect of interactions between them in relation to all analyzed variables characterizing the botanical composition of the sward.

Biomass Production
An impact of an applied fertilization on the biomass yields in both years of study was significant ( Table 2). The effect of grass mixture on sward yield was proven only in 2019. First of all, very visible differences between years were observed ( Figure 4). In the first year of study significantly the lowest yields (8.1 t·ha −1 of FM) were stated on control plots. The average biomass yields for individual fertilization treatments were higher, ranged from 12.1 t·ha −1 of FM to 13.7 t·ha −1 of FM, but did not differ significantly. No significant effect of grass mixture was proved. Such extremely low biomass yielding in 2018 was due to weather conditions i.e., low amount of precipitation in the vegetation season. In the second year noticed yields on average exceeded the level of 20 t·ha −1 of FM. Both type of plants mixtures and fertilization treatment significantly influenced the sward yielding. Biomass yielding of GM was significantly higher than LGM. Regardless of the mixture sward yields ranged from 15.5 t·ha −1 of FM (control plots) to 24.0 t·ha −1 of FM (NPK + Physioactiv).

Soil pH and Counts of Microorganisms
Both types of grass mixtures and fertilization treatment significantly influenced the evaluated soil parameters (Table 3). Only the soil pH (F = 1.97) was not significantly influenced by the plant mixture. While the counts of heterotrophic and oligotrophic bacteria, were influenced by type of grass mixture and fertilization. Table 3. Two-way ANOVA analyses of influence of fertilization and plant mixture on soil pH and counts of microorganisms in 2018-2019 (n = 168).

Mixture
Fertilization

Soil pH and Counts of Microorganisms
Both types of grass mixtures and fertilization treatment significantly influenced the evaluated soil parameters (Table 3). Only the soil pH (F = 1.97) was not significantly influenced by the plant mixture. While the counts of heterotrophic and oligotrophic bacteria, were influenced by type of grass mixture and fertilization. The average soil pH values for individual fertilization treatments ranged from 3.5 to 6.5 and were similar in both grass mixtures ( Figure 5). The fertilization had a highly significant effect on soil pH, which was not proved for the type of sward (Table 3, Figure 4). The highest pH value (about 6.5) was noted in the treatment where the NPK + Physioactiv was used. The second highest pH value was stated in the control plot (pH ≈ 6.0). All the other fertilization treatments did not differ significantly in pH values, which ranged from 5.3 to 5.5. The effect of the sward type and fertilization on the counts of studied soil microorganisms was significant (Table 3). It is noteworthy that the counts of microorganisms in the soil under LGM were almost twice higher compared to the values in the soil under the GM. This dependence was particularly noticeable for all groups of bacteria, and much less marked for fungi ( Figure 6).
The analysis clearly showed that the applied fertilizers differently affected the counts The effect of the sward type and fertilization on the counts of studied soil microorganisms was significant (Table 3). It is noteworthy that the counts of microorganisms in the soil under LGM were almost twice higher compared to the values in the soil under the GM. This dependence was particularly noticeable for all groups of bacteria, and much less marked for fungi ( Figure 6).
The analysis clearly showed that the applied fertilizers differently affected the counts of soil bacteria and fungi ( Figure 6). Similar patterns were obtained for bacteria groups, especially heterotrophic and oligotrophic bacteria. The control plot and the plots fertilised with the N-14 applied at a dose of 900 kg·ha −1 were distinguished by the count of these microorganisms. Copiotrophic bacteria were the most numerous in the treatment with the N-14 applied at a dose of 900 kg·ha −1 . This dependence was particularly noticeable in GM. The smallest counts of all bacterial groups were found in both mixtures in soil from plots fertilised with the NPK + PinKstart and NPK + Physiostart fertilizers. In most cases the sward type and applied fertilizers did not interact with each other to take effect on the count of soil microorganisms. However, such interaction was observed in the case of heterotrophic and oligotrophic bacteria (Table 3). In comparison with the other fertilization treatments, the NPK + PinKstart treatment affected the counts of both groups of soil microorganisms and reduced differences between GM and LGM.

Soil Enzymes Activity
The influence of the experimental factors on the activity of soil enzymes is shown in Figure 7. A significant difference between the experimental treatments and significant interactions between the experimental factors were stated ( Table 4). The exception was the urease activity, which was not significantly influenced by both experimental factors. The wide confidence intervals point to a large discrepancy in the results within the individual experimental variants. In most cases the sward type and applied fertilizers did not interact with each other to take effect on the count of soil microorganisms. However, such interaction was observed in the case of heterotrophic and oligotrophic bacteria (Table 3). In comparison with the other fertilization treatments, the NPK + PinKstart treatment affected the counts of both groups of soil microorganisms and reduced differences between GM and LGM.

Soil Enzymes Activity
The influence of the experimental factors on the activity of soil enzymes is shown in Figure 7. A significant difference between the experimental treatments and significant interactions between the experimental factors were stated ( Table 4). The exception was the urease activity, which was not significantly influenced by both experimental factors. The wide confidence intervals point to a large discrepancy in the results within the individual experimental variants. The dehydrogenase activity (Figure 7) ranged widely from about 4 to 14 µmol TP·kg −1 DM of soil·24 h −1 . It was the highest in the control variant and in the plots treated with the calcium fertilizer-Physioactiv. The differences between the mixtures were statistically insignificant, except for the control plots and the plots with the standard NPK fertilization.
The influence of the experimental treatments on the acid phosphatase activity was inconclusive. There was a noticeable interaction between both factors-type of fertilization and the type of the sward (Table 4, Figure 7). The lowest activity of acid phosphatase was below 0.13 µmol PNP·g −1 DM of soil·h −1 . However, this value was noted only in the soil under GM treated with traditional NPK and NPK + PinKstart fertilizers. The acid phosphatase activity in the soil under LGM was much higher (0.17-0.18 µmol PNP·g −1 DM of soil·h −1 ) both in these two fertilization treatments and in the other ones, except for the control and the NPK + Physioactiv treatment. In these fertilization treatments the acid phosphatase activity in the soil under both plant mixtures was almost identical, but this enzyme exhibited the lowest activity in LGM and the highest in GM (apart from the fertilization with the N-14 applied at a dose of 900 kg·ha −1 ).
The alkaline phosphatase activity was much more unequivocal. In both mixtures it was the highest in the control and the NPK + Physioactiv treatment, where its level was very similar, i.e., about 0.12 µmol PNP·g −1 DM of soil·h −1 . The differences between the plant mixtures were particularly noticeable in the other fertilization treatments. The alkaline phosphatase activity in these treatments was close to the highest values in LGM, whereas in the GM it was much lower-below 0.07 µmol PNP·g −1 DM of soil·h −1 .
In case of catalase activity, the highest values (close to 8.0) were noted in the control variant and the treatment with the NPK + Physioactiv fertilizer. However, there was a noticeable interaction between both factors-the effect of the fertilization depended on the plant mixture (Table 4). However, the application of the N-14 fertilizer at a dose of 900 kg·ha −1 in both mixtures resulted in high catalase activity.
The diagram showing the mean values of urease activity is to some extent similar to the diagrams showing the activities of the other enzymes. The highest urease activity was noted in the control, NPK + N-14 applied at a dose of 900 kg·ha −1 , and the NPK + calcium fertilizer-Physioaktiv. However, the differences between the experimental treatments were statistically insignificant, mainly due to a large spread in the results between the first and second term of the analyses. The dehydrogenase activity (Figure 7) ranged widely from about 4 to 14 µmol TP·kg −1 DM of soil·24 h −1 . It was the highest in the control variant and in the plots treated with the calcium fertilizer-Physioactiv. The differences between the mixtures were statistically insignificant, except for the control plots and the plots with the standard NPK fertilization.  The value of BIF, which to some extent results from the activity of dehydrogenases and catalase, exposed the interactions between the analyzed experimental factors. It turned out to be significantly the highest in the control, but only in LGM and in the treatments with Physioactiv-only in the GM.

Discussion
The fertilization treatments used in the experiment aimed to show whether the formula of biostimulants resulted in an increased soil biological activity compared to standard mineral NPK fertilization. As the botanical composition of sward has great influence on the function of the soil microbiome [59][60][61][62], the influence of the biostimulants on the populations of selected taxonomic and trophic groups of soil microorganisms and the soil enzymatic activity was analysed in two forage plants mixtures (GM and LGM) differing with share of functional groups of grasses and legumes.
The microbiological analyses of the soil were conducted in the third and fourth year of mowing utilisation. Therefore, the share of individual plant species in the sward compared with the share of the species in the sown seed mixtures changed significantly (Figures 2 and 3). These changes were significantly influenced by the fertilization treat-ments used in the experiment (Table 1). This effect was caused by several mechanisms. Firstly, fertilization, especially with nitrogen, promoted development of nitrophilous, highly competitive grass species. This phenomenon has been described by other authors [63,64]. The most abundant grass species in LGM were Lolium perenne and Festuca arundinacea (Figure 2). In GM it was Dactylis glomerata-a tall grass species, which turned out to be more competitive in these conditions than the photophilic, but low grass-Lolium perenne [65]. Lolium perenne was the dominant species in GM, but only in the control plot ( Figure 3). Luxuriantly developing grasses displaced legume species-Trifolium repens from the sward. It was easily noticeable both in LGM, where the highest share of Trifolium repens was found in the non-fertilised control plots (Figure 2), and in GM, where Trifolium repens appeared spontaneously and most abundantly in the non-fertilised plots (Figure 3).
It is well known that fertilization with nitrogen stimulates grasses development, increasing their competitiveness against legumes in legume-grass mixtures [66]. The growth of legume plants can be enhanced with P fertilizers in pastures [67]. On the other hand, some studies revealed that the combined applications of N and P fertilizer resulted in the improvement of grass growth in pasture swards but not in legume species [68].
The negative effect of nitrogen fertilization on the competitiveness of Trifolium repens against grasses is also caused by the limitation of the nitrogen fixation capacity of rhizobia, which live in symbiosis with legumes [69]. A large amount of nitrogen ions in the soil environment deactivates the enzymatic complex responsible for nitrogen fixationnitrogenase [70][71][72]. At the same time, the opposite situation, i.e., the lack of nitrogen fertilization, significantly increases the ability of legumes to compete with grasses. In our study this effect was clearly noticeable in the control plots in both plant mixtures (Figures 2  and 3).
It is noteworthy that in our study the competitiveness of Trifolium repens, an extremely photophilic species typically found in pastures [73], against grasses was reduced by the relatively low sward mowing frequency (three times a season), which promoted tall species [74]. Similarly, the competitiveness of the short Lolium perenne against the tall Dactylis glomerata was also reduced [75].
The influence of the fertilization treatments on the botanical composition of the sward also depended on the way they modified the soil pH. Mineral fertilizers used in all treatments (except the control plot) caused systematic soil acidification. This effect is particularly noticeable for fertilizers based on nitrogen and potassium [76,77]. At the same time, some of the applied biostimulants contained calcium carbonate, which has deacidifying properties. However, only one of the fertilizers (Physioactiv) contained as much as 76% CaCO 3 , so it was mostly a calcium fertilizer. The plots treated with the NPK + Physioactiv fertilizer had the highest pH (higher than in the control variant) ( Figure 5). A higher soil pH (6-6.5), noted in the control and the plots treated with the Physioactiv fertilizer, can stimulated the rhizobia growth thus promoted the development of legumes [78]. This dependence was noticeable in our study, because in both plant mixtures these treatments (the control plots and NPK + Physioactiv) were characterised by the highest share of Trifolium repens in the sward. Also another calcium fertilizer, containing CaCO 3 from ocean deposits, also had a similar, positive effect on the share of white clover in the sward and its durability [79].
As a result of these observations, the Physioactiv fertilizer significantly reduced the negative effect of nitrogen fertilization on the capacity of Trifolium repens to compete with grasses. The most likely mechanism of this interaction was pH optimisation. A comparative study on other calcium fertilizers might provide an answer to the question whether the specific origin of calcium carbonate or the effect of other components of this fertilizer had significant influence on the competing capacity of Trifolium repens.
Principal component analysis (PCA) (Figure 8) and Pearson's correlation coefficient (Table 5)  As can be seen in the diagram in Figure 8, there are four clearly d groups of related variables. Most of the soil biological parameters, especi of all trophic groups of bacteria and the catalase activity, were strongly rela of Trifolium repens in the sward. Although this is not evidence of a cause pendence, most reports indicate it. Legume forage species are known to in mass of the soil microbial community, because they enrich soil with organ root exudates, dying tissues, including root nodules, etc. [80]. Helios [81] [82] observed most of the fragments of leaves and other aerial organs w harvested by mowing come from Trifolium repens rather than tall grasses. T ments form a specific layer of felt (mulch), which increases the count of so isms. Residue quality, especially C:N ratio, is one of the factors that regul tion and, therefore, persistence time of plant biomass in the soil [83]. The l decomposition of grasses, in comparison with legumes, may be explained ratio [84]. Other factors that affect persistence of residues in the soil are t (micro) organisms, soil chemical and physical characteristics, and enviro tions [85]. Legumes include several pioneer species, able to colonize margi improve the nutritional status and organic matter content (and thus waterity) of the soil [86][87][88]. Moreover, it was shown that apart from rhizobia, n inhabited by many other non-rhizobial bacterial taxa, e.g., Bacillus sp., P Stenotrophomonas maltophilia, Micrococcus luteus, Erwinia persicina, and C lathyri [89,90].
It should be recalled at this point that 70-90% of microorganisms inha ments are non-cultivable [91,92]. Therefore, the CFU method we used to counting from soil samples for the quantification of bacteria and fungi offe As can be seen in the diagram in Figure 8, there are four clearly distinguishable groups of related variables. Most of the soil biological parameters, especially the counts of all trophic groups of bacteria and the catalase activity, were strongly related to the share of Trifolium repens in the sward. Although this is not evidence of a cause-and-effect dependence, most reports indicate it. Legume forage species are known to increase the biomass of the soil microbial community, because they enrich soil with organic matter from root exudates, dying tissues, including root nodules, etc. [80]. Helios [81] and Wolinska [82] observed most of the fragments of leaves and other aerial organs which cannot be harvested by mowing come from Trifolium repens rather than tall grasses. These dead fragments form a specific layer of felt (mulch), which increases the count of soil microorganisms. Residue quality, especially C:N ratio, is one of the factors that regulate decomposition and, therefore, persistence time of plant biomass in the soil [83]. The lowest biomass decomposition of grasses, in comparison with legumes, may be explained by higher C:N ratio [84]. Other factors that affect persistence of residues in the soil are the presence of (micro) organisms, soil chemical and physical characteristics, and environmental conditions [85]. Legumes include several pioneer species, able to colonize marginal soils and to improve the nutritional status and organic matter content (and thus water-holding capacity) of the soil [86][87][88]. Moreover, it was shown that apart from rhizobia, nodules are also inhabited by many other non-rhizobial bacterial taxa, e.g., Bacillus sp., Pseudomonas sp., Stenotrophomonas maltophilia, Micrococcus luteus, Erwinia persicina, and Chryseobacterium lathyri [89,90]. It should be recalled at this point that 70-90% of microorganisms inhabit all environments are non-cultivable [91,92]. Therefore, the CFU method we used to determine the counting from soil samples for the quantification of bacteria and fungi offers only limited amount of information. Despite its limitations, this method is still used in studies of soil microbial ecology, as it provides a reliable picture of their response to modifications in the soil environment [93].
Most of the studied enzymes: dehydrogenase, alkaline phosphatase, and urease were strongly related to the soil pH ( Table 5). The enzymatic activity of soil usually increases with increasing soil pH [94,95]. In our study the highest pH value was noted in the control plot (pH 5.8-6.0) and the NPK + Physioactiv combination (pH ≈ 6.5). It was optimal for the dehydrogenase activity [96,97] and alkaline phosphatase, which is responsible for the transformation of organic phosphorus into phosphates [98][99][100]. The share of Lolium perenne in the sward was slightly less related to this group ( Figure 8). However, it is noteworthy that the share of this species in the sward differed considerably from the share of the other grasses, which formed a separate group of variables in the diagram as they exhibited antagonism to almost all other parameters. This dependence confirmed the special bioceonotic role of Lolium perenne. This short grass species perfectly harmonises with the high share of Trifolium repens in the sward [101][102][103][104]. Moreover, it creates completely different conditions for the soil bioactivity compared to the tall Dactylis glomerata and Festuca arundinacea. Lolium perenne, has a very strong, shallowly spaced root system, which creates a much more compact turf and denser sward [105] that better protects the soil from overheating and drying [106]. Moreover, this grass has high sugar content and grows and develops very quickly, which translates into the amount and chemical characteristics of root exudates [6]. All these factors significantly influence the microbiome and soil enzyme activity [107,108].
The biomass yield shows a weak correlation with other variables and in the graph (Figure 8) it occupies a separate position, between Lolium perenne, which determines the yield of LGM, and other grasses, which determine the yield of the GM.
There may be some doubts about the relatively weak dependence between the count of soil bacteria and the activity of dehydrogenases, which are commonly regarded as enzymes of microbial origin [109][110][111]. However, it is noteworthy that dehydrogenases are usually the most positively correlated with the count of soil actinobacteria [112], which were not analysed in our study.
Fungi were a separate group showing very weak positive relations with the other parameters. They were characterised by a noticeably clear negative correlation with the soil pH. Soil fungi generally exhibit wider pH ranges for optimal growth, than bacteria [108], and tolerates lower pH better [1,100,108].
Apart from the urease activity, the other microbiological parameters almost always had higher values in the soil under LGM. This mixture was characterised by much higher shares of Trifolium repens and Lolium perenne in the sward. The influence of the fertilization treatments on the soil microbiological properties and its metabolism largely coincided with their influence on the soil pH and the sward botanical composition. It clearly resulted from the PCA and the high similarity of the diagrams shapes showing the influence of the fertilizer combinations on pH and on the activity of dehydrogenases and alkaline phosphatase. They had the highest values in the control plots and the treatments with the Physioactiv fertilizer. The shape of the diagrams showing the count of fungi ( Figure 6) was exactly opposite and these two combinations were characterised by the lowest values. Many authors stress the fact that even small differences in pH translate into the size and structure of the soil microbiome [113].
Apart from dehydrogenases and to some extent urease (the effect of the latter proved to be statistically insignificant in our study), catalase is the enzyme whose activity is the most strongly related to the soil microbiome, the intensity of changes occurring in it, and plants' viability [114]. The diagrams illustrating the counts of bacteria and the catalase activity also showed its highest values in the control plot and in the NPK + Physioactiv treatment. It was also high in the treatment with the N-14 fertilizer applied at a dose of 900 kg·ha −1 (Figures 6 and 7). Why was the effect of the N-14 biostimulant so clearly pronounced despite the low pH? It is most likely that this effect was caused by high share of Trifolium repens in the sward.

1.
Traditional NPK fertilization without biostimulants decreased the counts of most of the groups of soil microorganisms under analysis and the activity of soil enzymes, as compared with the non-fertilised absolute control variant. This effect was changed only in two cases, where the mineral NPK fertilization was accompanied by biostimulants, i.e., the N-14 at a dose of 900 kg·ha −1 and Physioactiv. 2.
The applied fertilizers had multidimensional and largely indirect influence on the soil microbiome and the activity of soil enzymes. It was mainly caused by the modification of the share of Trifolium repens in the sward and the pH of the soil environment. The effect of additional substances contained in the biostimulants seemed to be significant only at very high doses of these fertilizers-the N-14 applied at a dose of 900 kg·ha −1 was effective but it was ineffective at a dose of 300 kg·ha −1 .

3.
The optimisation of the soil pH with CaCO 3 applied at a dose increasing its value from 5.5 to 6.5 may reduce the negative effect of intensive nitrogen fertilization on the competitiveness of Trifolium repens against grasses. In our experiment this effect was observed after the application of the Physioactiv biostimulant. It is necessary to check whether the same effect can be observed after the application of other fertilizers containing calcium in the carbonate form.

4.
The effect of biostimulants on forage plants and soil microflora is not well understood yet. The investigation of new solutions and combinations of mineral nutrients with new biostimulants in fertilizers, which will affect plants and soil not only by optimizing soil pH are still needed.