The Effect of Biogas Slurry Application on Biomass Production and Forage Quality of Lolium Multiflorum

The development of ecological circular agriculture has been highly encouraged to recycle agricultural wastes, reduce mineral fertilizer input, and protect the environment. Biogas slurry (BS), a by-product of biogas production generated from anaerobic digestion of animal waste and crop residues, is often considered a substitute to reduce mineral fertilizer input. Being a cheap source of organic matter and plant nutrients, its application may improve soil fertility and yield quality and quantity. The field experiments were conducted in 2016 and 2017 to study the plant growth responses and forage quality by applying biogas slurry to replace chemical synthetic fertilizer (CSF). Results revealed that biogas slurry combination with chemical synthetic fertilizer significantly (p < 0.05) improved the growth of Italian ryegrass on treatment with T2, and the Italian ryegrass dry matter was increased by more than 9.00%, while the stem-to-leaf ratio was decreased by more than 12% (second cutting), in comparison with only chemical synthetic fertilizer group. In the case of forage quality, the crude protein (CP) and crude fiber (CF) content had a significant difference was observed between the T0 and T2 treatment group. Compare with the chemical synthetic fertilizer group, the CP content improved by 10.35%, and the CF content decreased about 10.00%. Based on these results, it was concluded that the application of 37.5 kg/ha CSF + 100.5 t/ha BS could improve the production of biomass and forage quality in Italian ryegrass.


Introduction
As number of livestock has increased all around the world, the amount of animal manure has increased simultaneously. From 2006 to 2016, the number of total livestock increased from 293.48 to 374.99 million in Bangladesh [1]. In China, the livestock population is estimated to be 9.03 million cattle and buffaloes, 0.61 million horses and donkeys, 44.06 million goats and sheep, and 74.66 million hogs in 2017 [2]. However, the increasing amounts of the livestock wastes, such as feces and urine, were difficulty in disposing of as it induced extensive environmental pollution [3,4]. Besides, economical reason is also important for avoiding over-fertilization. Holm-Nielsen et al. [5] have suggested that increasing animal production areas need suitable manure management and to optimize their recycling. For sustainable development, it is necessary to find out a cost-effective way to dispose of these wastes.
Anaerobic digestion is one of the effective solutions for livestock wastes and satisfies growing concerns on energy supply [6][7][8]. Among the alternative energy generation technologies, anaerobic digestion is an efficient way to transform waste into resources and it benefits by improving fertilizer qualities, reducing odors and pathogens of livestock wastes. The main product of anaerobic digestion is biogas, which is an important and clean energy source. The biogas could be used for human activity like cooking or lighting, thereby Recently, more and more studies have focused on the effects of biogas slurry on crop biomass yield, soil properties, and fertility. Balance of the useage of biogas slurry is not only avoiding environmental risk of over-fertilization, but also for optimum plant growth. However, environmental factors were selected only as their treatment, and the integrated effects that related to crop growth, yield, and quality by applying biogas slurry are rarely investigated. Additionally, previous studies mainly focused on cereal crops or conventional vegetables, research on Italian ryegrass (lolium multiflorum) were reported rarely, which account for a large proportion of the total cool-season forages planting area [34,36]. In brief, there is little information available related to dealing with the use of biogas slurry and synthetic chemical fertilizer as an alternative to only synthetic fertilizer application in Italian ryegrass production. The objective of the study was to evaluate the effects of biogas slurry in combination with synthetic chemical fertilizer on biomass yield production and nutritional status of Italian ryegrass, and find the optimal ratio of synthetic fertilizer and biogas slurry to dissolve baogass slurry and reduce the synthetic fertilizer application. On the other hand, this study can provide basic information about that the balance of Italian ryegrass biomas yield and biogas slurry assimilating capacity.

Experimental Materials and Procedures
The experiment was conducted in Sichuan Academy of Agricultural Sciences farm (Hong Ya, Sichuan, China; N 29 • 32 , E 103 • 15 , 442 m) and the soil at the site was identified as viscous yellow soil. The land was prepared well and weeds were carefully removed. The land was divided into 18 plots (each plot was 15 m 2 with 3 m wide and 5 m length) according to the layout of the experiment (Figure 1). The experiment was performed by two growth seasons, the first growth season from October 2016 to April 2017, and the second growth season from October 2017 to April 2018. Italian ryegrass was sowed in September each growth season and the biogas slurry treatment was applied between October and April in each growth season. and has good environmental, ecological, and economic effects [39], which makes it become a good plant with absorbing biogas slurry to reduce environmental pollution.
Recently, more and more studies have focused on the effects of biogas slurry on crop biomass yield, soil properties, and fertility. Balance of the useage of biogas slurry is not only avoiding environmental risk of over-fertilization, but also for optimum plant growth. However, environmental factors were selected only as their treatment, and the integrated effects that related to crop growth, yield, and quality by applying biogas slurry are rarely investigated. Additionally, previous studies mainly focused on cereal crops or conventional vegetables, research on Italian ryegrass (Lolium multiflorum) were reported rarely, which account for a large proportion of the total cool-season forages planting area [34,36]. In brief, there is little information available related to dealing with the use of biogas slurry and synthetic chemical fertilizer as an alternative to only synthetic fertilizer application in Italian ryegrass production. The objective of the study was to evaluate the effects of biogas slurry in combination with synthetic chemical fertilizer on biomass yield production and nutritional status of Italian ryegrass, and find the optimal ratio of synthetic fertilizer and biogas slurry to dissolve baogass slurry and reduce the synthetic fertilizer application. On the other hand, this study can provide basic information about that the balance of Italian ryegrass biomas yield and biogas slurry assimilating capacity.

Experimental Materials and Procedures
The experiment was conducted in Sichuan Academy of Agricultural Sciences farm (Hong Ya, Sichuan, China; N 29°32′, E 103°15′, 442 m) and the soil at the site was identified as viscous yellow soil. The land was prepared well and weeds were carefully removed. The land was divided into 18 plots (each plot was 15 m 2 with 3 m wide and 5 m length) according to the layout of the experiment (Figure 1). The experiment was performed by two growth seasons, the first growth season from October 2016 to April 2017, and the second growth season from October 2017 to April 2018. Italian ryegrass was sowed in September each growth season and the biogas slurry treatment was applied between October and April in each growth season. The experiment was conducted in a randomized block design comprising 6 levels of fertilizer treatment: 112.4 kg/ha chemical synthetic fertilizer (CSF), 74.9 kg/ha CSF + 49.5 t/ha biogas slurry (BS), 37.5 kg/ha CSF + 100.5 t/ha BS, 150 t/ha 2 BS, 300 t/ha BS, 450 t/ha BS, referred to treatments T0, T1, T2, T3, T4, and T5, respectively. The physicochemical properties of the soil, chemical fertilizer, and biogas slurry are given in Table 1. The biogas slurry used for this study was collected from a commercial biogas plant, named Modern animal husbandry Hongya Ranch, which was produced by fermentation of cow dung. The biogas plant has a special pipeline to transport the biogas slurry to the test site. The experiment was conducted in a randomized block design comprising 6 levels of fertilizer treatment: 112.4 kg/ha chemical synthetic fertilizer (CSF), 74.9 kg/ha CSF + 49.5 t/ha biogas slurry (BS), 37.5 kg/ha CSF + 100.5 t/ha BS, 150 t/ha 2 BS, 300 t/ha BS, 450 t/ha BS, referred to treatments T0, T1, T2, T3, T4, and T5, respectively. The physicochemical properties of the soil, chemical fertilizer, and biogas slurry are given in Table 1. The biogas slurry used for this study was collected from a commercial biogas plant, named Modern animal husbandry Hongya Ranch, which was produced by fermentation of cow dung. The biogas plant has a special pipeline to transport the biogas slurry to the test site. For observing the efficiency of biogas slurry, the Italian ryegrass variety 'Splendor' was selected for this experiment, which was normally used for fodder production in China. The Italian ryegrass seeds were sown using a line sowing method, maintaining row-to-row spacing of 30 cm with the long side of plot, and sown 10 line per plot. The seed rate was 30 kg/ha and the seed germination rate was 90%. No irrigation or insecticide was applied during the experimental period. The plant height of Italian ryegrass was measured (cm) before harvesting every time. After second and third cutting of the Italian ryegrass plants each growth season, the stem and leaf weight were measured from 1 kg plants, then the stem-to-leaf ratio was calculated. The Italian ryegrass was harvested when the plants are grown to 50-60 cm. Harvesting was performed manually and samples were kept for further analysis. The biomass yield of the forage from each plot was recorded using a weighing balance, immediately.

Biochemical Analysis
To determine dry matter (DM) content of Italian ryegrass, the fresh samples were inactivated in oven at 105 • C for 30 min, and then dried at 65 • C to a constant weight. Finally, the dried samples of the first cutting of first growth season were stored in a dry container until used for analysis of crude protein (CP), ether extract (EE), and fiber fractions (crude fiber (CF), neutral detergent fiber (NDF), and acid detergent fiber (ADF)).
The CP content was determined by wet chemistry analysis according to the Kjeldahl method [40] with the KjeltecTM 8400 analyzer unit (FOSS, Hoganas, Sweden). Briefly, the ground samples (0.5 g) were added into a 250 mL TKN digestion tube with 10 mL concentrated sulfuric acid and two digestive tablets (Beijing Jinyuanxingke Technology, Beijing, China). Blanks containing all these reagents were simultaneously processed. All tubes were digested in the preheated digestion block (TecatorTM digestor auto; FOSS, Ho-ganas, Sweden) at 420 • C for 90 min or until the samples were green and clear. The CP content was calculated using the following equation: where: V1 = volume (mL) of standard HCl required for sample; V2 = volume (mL) of standard HCl required for blank; C = molarity of standard HCl; 1.4007 = milliequivalent weight of N × 100; 6.25 = average coefficient of nitrogen conversion into proteins; M = sample weight in grams.
The crude fat (ether extracts, EE) content was determined using the ether extraction method after acid-hydrolysis [40] with a Soxtherm apparatus (Gerhardt, Bonn, Germany). Briefly, a total of 2 g samples were placed in a degreased filter paper package and placed in a flask with carborundum dried to a constant weight (M2) at 70 • C. Then, the extraction was carried out in a Soxtherm apparatus (Gerhardt, Bonn, Germany). After ether extract, the sample package was placed in a flask with carborundum and dried to a constant weight (M1) at 70 • C. The contents of EE were calculated using the following equation: where: M = sample weight; M1 = weight of flask with carborundum; M2 = weight of flask with carborundum and weight of dry residue after the ether extract.
The CF, NDF, and ADF contents were measured by using the methods described by Goering and Van Soest [41] and Van Soest, et al. [42] in an automatic fiber analyzer (ANKOM 2000 Fiber Analyzer, ANKOM Technology, NY, USA). Briefly, for CF analysis, a total of 1 g pulverized samples were placed into the fiber filter bag (F57, ANKOM Technology, NY, USA) and preextracted in a beaker with ether. Then, they were placed into the automatic fiber analyzer for digesting the samples with 0.255 N H 2 SO 4 (FCFA20, ANKOM Technology, NY, USA) and 0.313 N NaOH (FCFB20, ANKOM Technology, NY, USA) solutions. For NDF analysis, a total 0.5 g pulverized samples were placed into the fiber filter bag (F57, ANKOM Technology, NY, USA) and placed into the automatic fiber analyzer for digesting the samples with dilutes solutions of neutral detergent dry concentrate with triethylene glycol (FND20C, ANKOM Technology, NY, USA). The residues obtained after the deter-mination of NDF can be directly digested by the solution of acid detergent concentrate-dry CTAB powder (FAD20C, ANKOM Technology, NY, USA) for ADF determination.The contents of CF, NDF, and ADF were calculated using the following equation: where: M = sample weight; M1 = bag tare weight; M2 = weight of organic matter after extraction by neutral detergent; M3 = weight of organic matter after extraction by acid detergent; C = loss in weight on ignition of bag/sample; C1 = ash-corrected blank bag factor (a running average of the loss of weight after extraction of the blank bag/original blank bag).

Statistical Analysis
All data obtained were analyzed using a one-way analysis of variance (ANOVA) with SPSS 18.0 (IBM, Armonk, NY, USA). A mean comparison was performed using Fisher's least significant difference test (LSD) and the Duncan multiple range test with a significance level of p < 0.05.

The Plant Height of Italian Ryegrass with the Biogas Slurry Applied
Plant height is an important component for biomass production in Italian ryegrass. Table 2 showed that the plant height of Italian ryegrass varied at different harvest times and different levels of biogas slurry treatment. The application of different levels of biogas slurry was non-significantly different at the first cutting in two growth seasons. However, for the second and third cutting, the plant height of Italian ryegrass treated with T2 were 102.52 and 97.00 cm, respectively, which was significantly higher than that treated with Sustainability 2021, 13, 3605 6 of 13 T1, T3, T4, and T5 in the 1st growth season (p < 0.05). In the 2nd growth season, the plant height of Italian ryegrass was non-significantly changed with the increased biogas slurry at second and third cutting, while T2 was obviously higher than T0 which was no biogas slurry applied.  Table 3 showed the fresh weight of biomass yield of Italian ryegrass under CSF with different BS ratios at different harvest times. The significant differences (p < 0.05) were found among all treatment groups. Specially, the Italian ryegrass treated with T2 displayed the highest fresh weight at different harvest times in the 1st growth season, which were 36647.54 kg/ha, 54663.61 kg/ha, and 42956.58 kg/ha, respectively. The recorded improvements were 7.64%, 7.99%, and 6.96% than that of T0 treatment at first cutting, second cutting, and third cutting. In the 2nd growth season, the biomass yield of Italian ryegrass with T2 treatment has higher fresh weight than T0 at different cutting times, and the yields improvements were 4.71%, 23.60%, and 17.70% at first cutting, second cutting, and third cutting, respectively. It was worth mentioning that by increasing the BS ratio to two or three folds in the 2nd growth season, the fresh weight of the second cutting and third cutting increased more significantly than the treatment of T2, which could be the drought in spring at the location.

The Stem-to-Leaf Ratio of Italian Ryegrass with the Biogas Slurry Applied
Stem-to-leaf ratio is a component of canopy architecture that determines diet selection and forage intake in forages. Increased stem-to-leaf ratio is a major cause of the decline in forage quality, and also the loss in quality under adverse hay curing conditions [31]. Table 4 illustrates the ratio of the stem and leaf in Italian ryegrass under CSF and different BS ratios at different cutting times. By increasing the biogas slurry fertilizer ratios to 2/3 (T2 treatment), the ratio of stem and leaf reduced to the minimum value at second cutting and third cutting in the 1st and 2nd growth season, respectively. Mean values with different letters in the same row differ significantly, p < 0.05. Figure 2 shows the dry matter (DM) accumulation in above-ground parts of Italian ryegrass plants in two growth seasons. The DM accumulation of Italian ryegrass with T2 treatment was significantly higher than other treatments in the 1st growth season (p < 0.05). However, in the 2nd growth season, the Italian ryegrass DM accumulation treated with T2 occupied third place in all treatments, and the treatment of T4 and T5 were significantly higher than others (p < 0.05). Analyzing the total DM yield, contributing characters according to the DM accumulation of Italian ryegrass at every cutting time, it could be found that the Italian ryegrass hay yield of the third cutting contributed a lot to the DM accumulation in the 1st growth season, and the same results were also found in the 2nd growth season. tion and forage intake in forages. Increased stem-to-leaf ratio is a major cause of the decline in forage quality, and also the loss in quality under adverse hay curing conditions [31]. Table 4 illustrates the ratio of the stem and leaf in Italian ryegrass under CSF and different BS ratios at different cutting times. By increasing the biogas slurry fertilizer ratios to 2/3 (T2 treatment), the ratio of stem and leaf reduced to the minimum value at second cutting and third cutting in the 1st and 2nd growth season, respectively. Mean values with different letters in the same row differ significantly, p < 0.05. Figure 2 shows the dry matter (DM) accumulation in above-ground parts of Italian ryegrass plants in two growth seasons. The DM accumulation of Italian ryegrass with T2 treatment was significantly higher than other treatments in the 1st growth season (p < 0.05). However, in the 2nd growth season, the Italian ryegrass DM accumulation treated with T2 occupied third place in all treatments, and the treatment of T4 and T5 were significantly higher than others (p < 0.05). Analyzing the total DM yield, contributing characters according to the DM accumulation of Italian ryegrass at every cutting time, it could be found that the Italian ryegrass hay yield of the third cutting contributed a lot to the DM accumulation in the 1st growth season, and the same results were also found in the 2nd growth season.

Nutritional Status of Italian Ryegrass with the Biogas Slurry Applied
The nutritional status of Italian ryegrass in the different biogas slurry treatment groups was presented in Table 5. Crude protein (CP) content of T2 (25.06% DM) was significantly higher (p < 0.05) than T0 (22.71% DM), T1 (23.85% DM), T3 (23.11% DM), T4, and T5 (both with 23.68% DM), and CP content gradually increased in T1 and T2, but decreased in T3, T4, and T5. Ether extract (EE) content was higher for T2 (6.40% DM) than for the other treatment groups, but there was no significant difference except T5. In the case of CF content, significant differences (p < 0.05) were between the treatment groups. The lowest CF content was in T2 (15.50% DM), and CF gradually decreased in T1 and T2, but increased in T3, T4, and T5. For the neutral detergent fiber (NDF) and acid detergent fiber (ADF) content, there were no statistically significant differences between treatment groups, indicating that biogas slurry had no effect on the NDF and ADF.

Discussion
With the growing energy demand and inadequate natural resource availability, biogas technology has been considered a new renewable energy alternative, obtained by using cattle dung, poultry litter, and other wastes [43]. However, with the synthesis of biogas through anaerobic digestion, solids and other wastes are released as slurry, which is commonly considered a waste product and potential caused environmental pollution. An optimized level of BS fertilizer is a good source of plant nutrients and can significantly improve soil properties and crop productivities [44], providing an opportunity to treat and reutilize various organic wastes and minimize greenhouse gas emissions [14,45].
Analysis of variance (ANOVA) is the most efficient method available for the analysis of experimental data, and it is a technique for analyzing the way in which the mean of a variable is affected by different types and combinations of factors. One-way analysis of variance is the simplest form. It is an extension of the independent samples t-test [46] and can be used to compare any number of groups or treatments. This method is widely used in agricultural research. Previous study showed that ANOVA was used to compare changes in the phenological properties, yield, and phytochemical composition of guar (Cyamopsis tetragonoloba L.) landaraces under various irrigation regimes and planting dates [47]. In the present study, we designed a randomized block experiment to compare the effects of CSF + BS on the biomass yield and forage quality of Italian ryegrass with ANOVA, which could have statistical meaning for agriculture application.
It has been proved that the application of biogas slurry as fertilizer could effectively stimulate the growth of plants and is the ideal substitute of chemical synthetic fertilizer. As a kind of biofertilizer with abundant nutrients, bioslurry has been reported to have significantly increased yields in S. tuberosum L. [48], Spinacia oleracea L. [49], Capsicum annuum L. [49], Abelmoschus esculentus (L.) Moench [50], and Lactuca sativa L. [51] compared with only the application of synthetic fertilizer. However, in order to obtain optimized growth and be environment friendly, the suitable application rate of biogas slurry should be study. Based on this study, approximately 37.5 kg/ha CSF + 100.5 t/ha BS is the optimum level for Italian ryegrass growth, and the yield greater than the yield of traditional farmer practice (only used the chemical synthetic fertilizer), which was consistent with previous study in tomato by Ferdous et al. [52]. Marino et al. [53] reported that the nitrogen fertilizer significantly increased the forage dry matter and N accumulation. Although Wentzel and Joergensen [54] reported the positive effect of biogas slurry fertilizer on Italian ryegrass plant biomass, the biomass yield was non-significantly increased with the excess use of biogas slurry. Similarly, in order to obtain the maximum growth of maize fodder, an optimal level of biogas slurry (70 kg of slurry N ha −1 ) was suggested to be applied, while excessive use beyond this optimal level might induce adverse effects [27]. Other factors like genetic makeup of plant, soil fertility, and environment conditions were all important to plant height of forage.
It has been demonstrated that BS might increase the availability of macro-and micronutrients, particularly nitrogen, which in turn promotes meristematic growth, leading to higher bioactive compounds like anthocyanin and vitamin. BS is the rich source of nitrogen which holds prime importance since it is the key element in all proteins including enzymes which regulate all the metabolic processes in plants [27]. Nitrogen also has a direct correlation with the important metabolic pathways such as growth, photosynthesis, and consequently affects the light energy efficiency and yield because it is the main element of leaf chlorophyll content [55,56]. The plant photosynthetic rate may be linear in correlation to leaf nitrogen content but the continuous lack of nitrogen reduces photosynthesis rate, soluble protein, enzymes, and chlorophyll [57]. The observed increase in forage yield in response to slurry nitrogen might be owing to its positive effect on cell elongation, cell division, and increasing photosynthetic surface, resulting in more production and accumulation of photosynthetic compounds.
Usually, crops can directly use soil-available NH4 + -N, P, and K in biogas slurry [8,58]. Further, it has been reported that the use of biogas slurry in combination with chemical fertilizer could effectively increase soil organic matter, nutrient concentrations, and N use efficiency [59]. The increase in nitrogen use efficiency could promote root and shoot growth, and increase the crop biomass yield [52]. Katuwal and Bohara [60] reported an increase in vegetable crop yields among 68% of farming households, along with an increase in revenue among 42% households after use of biogas slurry. Muhmood et al. [49] reported that synthetic fertilizer in combination with biogas slurry applied could increase economic returns from spinach (Spinacia oleracea) and chilli (Capsicum frutescens) production. In this study, the results showed that the optimum level of biogas slurry for Italian ryegrass could benefit plant growth, and the yield, which was similar to those reported by Daudén and QuíLez [61], who conducted a maize yield experiment using different levels of slurry. However, these studies mainly explored the effect of BS on yield, rather than both the yield and quality of productions.
Forage DM content is also an indicator of its quality according to variety, stage of forage maturity, season, temperature, and environment conditions. In the current study, due to soil fertilizer residues in the 1st growth season and the improvement of soil structure by biogas slurry, dry matter production in the 2nd growth season was higher than that of the previous year [62]. Interestingly, the DM content of the second and third cutting were higher than other treatments results where the DM production of T4 and T5 in 2017 were higher than others. When the Italian ryegrass suffered drought stress during the second and third growth, the treatment of T4 and T5 can provide enough water for Italian ryegrass to grow. Besides, the CP and CF content of forage is the most important indicator contributing to its forage quality. Generally speaking, it has a direct effect on milk production in dairy animals and body growth in cattle. Sarker [63] reported that the increased use of nitrogen fertilizer could effectively increase the CP content in Zamboo grass (Hybrid jowar). The previous study of mixed sown pasture in North Wales indicated that application of digestate leads to a similar response in pasture yield as when compound inorganic fertilizer is applied, and better than straight N organic fertilizer [64]. In this study, application of 37.5 kg/ha chemical synthetic fertilizer and 100.5 t/ha biogas slurry significantly increased the CP content. In addition, the CF content of Italian ryegrass decreased when the biogas slurry was used in combination with chemical fertilizer. The reason may be that biogas slurry had the potential to increase Italian ryegrass yield and improve the forage quality. Additionally, use of chemical fertilizer with biogas slurry could be a viable input for production of high yield and high quality forage. Similar results were also reported in switchgrass by the University of Life Sciences in Lublin where the biogas digestate can increase the biomass yield, the content of protein, P, and Mg of switchgrass [65]. The present study provides guidance for the use of BS to grow Italian ryegrass, especially in the dairy farm and the developed area of animal husbandry. Italian ryegrass can provide feed for cow, and the cow dung is used to produce biogas to provide energy, and BS as a by-product can provide nutrients for Itanlian ryegrass. Due to an 800 m 3 volume, the biogas system needs to discharge 15 tons of biogas slurry daily; this circulation model can be applied on a large scale to reduce environmental pollution. The biogas slurry produced by a scale farm can be supplied to Italian ryegrass planting for animal feed. However, the biogas slurry transportation seems troublesome, and it is best to have pipeline laying as a cost-effective way, which sometimes requires government financial supply.
Biogas slurry used in crop cultivations is an agricultural eco-circulatory system for the reuse of livestock and poultry waste and reduction of chemical syntheic fertilizer input. Although this study explored the appropriate amount of biogas slurry to improve the yield of Italian ryegrass, there is much more to be learned. Uncertain chemical characteristics of biogas slurry including liquidity, high pH, high ammonium (NH4 + ) content, and potentially high heavy metal content [20,29,31], which are very different with traditional mineral fertilizers, could induce biology security risks. Thus, necessary evaluating studies must be performed to elucidate the effects of biogas slurry on crop development and environment before its extensive application. Especially, the accumulation of heavy metals in plants and soil to prevent them from affecting human health through the food chain.

Conclusions
Biogas slurry, a by-product of biogas production generated from anaerobic digestion of animal waste and crop residues, is often considered a substitute to reduce mineral fertilizer input. The results of the present study showed that biogas slurry in combination with chemical synthetic fertilizer significantly (p < 0.05) improved the growth of Italian ryegrass on treatment with T2, and the Italian ryegrass dry matter increased by more than 9.00%, while the stem-to-leaf ratio decreased by more than 12% (second cutting), in comparison with only chemical synthetic fertilizer group (T0). The forage quality index of CP and CF also had a significant difference and was observed between the T0 and T2 treatment group. Compare with the chemical synthetic fertilizer group, the CP content improved by 10.35%, and the CF content decreased by about 10.00%. Based on these results, it was concluded that the application of 37.5 kg/ha CSF + 100.5 t/ha BS could improve the production of biomass and forage quality in Italian ryegrass, but also replace part of the chemical fertilizer applied in agriculture.
Author Contributions: W.X. and C.L., Conceived the idea, analyzed, and interpreted the data, drafted and revised the manuscript; X.W., H.W. and L.Y., Collected and analyzed the data; Y.Z. and L.J., revised the manuscript. All authors have read and agreed to the published version of the manuscript. Data Availability Statement: All data has been shared in this manuscript.

Conflicts of Interest:
The authors declare no conflict of interest.
Ethical Approval: This article does not contain any studies with human participants or animals performed by any of the authors.

Abbreviations
The abbreviation were used in this study.