NC State
BioResources
Pan, J., Shen, J., Zhou, Z., Xin, Y., Huang, Z., Xiong, J., Liu, Y., Cui, X., and Liu, Y. (2025). "Sustainable management of biogas slurry discharge in biogas engineering: As a chemical fertilizer substitute for garlic cultivation," BioResources 20(1), 790–808.

Abstract

To address the persistent challenge of managing livestock and poultry manure resulting from agricultural intensification and mechanization, the application of biogas engineering has steadily expanded. Biogas slurry, a valuable byproduct abundant in nitrogen and other nutrients, emerges as an attractive alternative to chemical fertilizers. This study investigated the effects of substituting chemical fertilizers with biogas slurry at varying application rates on garlic growth and soil properties. The results indicate that with increasing application of biogas slurry, the promotion of garlic growth and soil nutrients exhibited an initial increase followed by a decrease. Notably, application of biogas slurry with nitrogen content equivalent to chemical fertilizers (G-1 group) gave the most pronounced promotion effect. In the G-1 group, garlic yield and soil organic matter content reached 1.06 kg/m3 and 40.6 g/kg, respectively, representing increases of 11.6% and 17.6%, respectively, compared to the chemical fertilizer group. Furthermore, after the application of biogas slurry, concentrations of both garlic and soil heavy metals remained within standard limits. Biogas slurry can be recommended as an effective substitute for chemical fertilizers, fostering garlic growth, boosting yields, enhancing soil organic matter content, and promoting biological carbon sequestration.


Download PDF

Full Article

Sustainable Management of Biogas Slurry Discharge in Biogas Engineering: As a Chemical Fertilizer Substitute for Garlic Cultivation

Junhui Pan,a,# Jiali Shen,a,# Zixuan Zhou,a Yicong Xin,b Zhenxia Huang,b Jianghua Xiong,b Yuhuan Liu,a Xian Cui,a,* and Yuxin Liu b,*

To address the persistent challenge of managing livestock and poultry manure resulting from agricultural intensification and mechanization, the application of biogas engineering has steadily expanded. Biogas slurry, a valuable byproduct abundant in nitrogen and other nutrients, emerges as an attractive alternative to chemical fertilizers. This study investigated the effects of substituting chemical fertilizers with biogas slurry at varying application rates on garlic growth and soil properties. The results indicate that with increasing application of biogas slurry, the promotion of garlic growth and soil nutrients exhibited an initial increase followed by a decrease. Notably, application of biogas slurry with nitrogen content equivalent to chemical fertilizers (G-1 group) gave the most pronounced promotion effect. In the G-1 group, garlic yield and soil organic matter content reached 1.06 kg/m3 and 40.6 g/kg, respectively, representing increases of 11.6% and 17.6%, respectively, compared to the chemical fertilizer group. Furthermore, after the application of biogas slurry, concentrations of both garlic and soil heavy metals remained within standard limits. Biogas slurry can be recommended as an effective substitute for chemical fertilizers, fostering garlic growth, boosting yields, enhancing soil organic matter content, and promoting biological carbon sequestration.

DOI: 10.15376/biores.20.1.790-808

Keywords: Biogas slurry; Garlic; Alternative fertilizer; Carbon sequestration

Contact information: a: State Key Laboratory of Food Science and Resources, Engineering Research Center for Biomass Conversion, Ministry of Education, Nanchang University, Nanchang, Jiangxi, 330047, China; b: Agricultural Ecology and Resources Protection Station of Jiangxi Province, Nanchang, Jiangxi, 330046, China; *Corresponding author: cuixian@ncu.edu.cn; 15170014031@163.com;

# The first and second author contributed equally to this work.

 

GRAPHICAL ABSTRACT

INTRODUCTION

With the advancement of intensification and mechanization, the scale of livestock and poultry production worldwide is gradually expanding, resulting in an increasingly concentrated production of waste (Jia et al. 2018; Zhang et al. 2021). China stands as one of the largest producers and consumers of livestock and poultry products globally, with an annual production of livestock and poultry manure exceeding 4 billion tons nationwide. Over 40% of this manure remains not fully utilized, leading to serious environmental pollution. Therefore, the effective and comprehensive utilization of livestock and poultry manure is crucial (MARAC 2017; Zhou et al. 2019). Biogas engineering, which is primarily based on the anaerobic digestion of straw and livestock and poultry manure, represents a systematic approach for emissions reduction and carbon sequestration that integrates manure treatment, biogas production, and resource utilization. Biogas engineering not only manages livestock and poultry manure, increases the supply of renewable clean energy, and facilitates carbon peak and carbon neutrality, but it also promotes agricultural recycling through the application of biogas residue to fields, enhancing soil organic matter (Song et al. 2023; Chen et al. 2024). By the end of 2015, the government had subsidized 111,000 large-scale anaerobic digestion plants (NDRC and MOA 2017). Biogas engineering faces challenges in its practical promotion, such as difficulties in biogas slurry utilization and digestion (Gu et al. 2024; Zhou et al. 2024).

Biogas slurry (BS) is the liquid product of the anaerobic digestion of organic materials such as straw and livestock and poultry manure in biogas engineering (Mohanty et al. 2024). There are two main methods for treating BS: standard discharge and return to fields. Due to the high cost of BS treatment, its application to fields is the primary method of treatment. Nitrogen fertilizer is one of the indispensable important fertilizers in agricultural production, and it plays a vital role in the growth and development of crops (Tei et al. 2020; Zhang at el. 2022). Among them, nitrate (NO3) and ammonium nitrogen (NH4+) are the main inorganic nitrogen sources of N (Ye et al. 2022). Nitrated nitrogen can be leached due to its high mobility, resulting in nitrogen loss. In addition, inorganic nitrogen is lost to the environment by participating in the volatilization of NH2 and nitrification into forms that are not usable by plants, such as NO (Lyu et al. 2024). Chemical fertilizers are currently the most commonly used inorganic fertilizers due to their ease of use. However, the long-term application of inorganic nitrogen fertilizers accelerates soil compaction and acidification, affecting crop growth. Moreover, excessive use of chemical fertilizers leads to significant greenhouse gas emissions and water pollution issues, thus hindering the achievement of carbon peak and neutrality (Nishita et al. 2017; Iqbal et al. 2020; Ma et al. 2024). To develop sustainable low-carbon agriculture, using of nitrogen-efficient fertilizers is essential (Ma et al. 2024).

Compared to other organic fertilizers, BS contains abundant nitrogen, phosphorus, and other nutrients, providing sufficient nutrition for crop growth. BS also contains a rich array of nutrients, such as various trace elements (iron, zinc, manganese, etc.), amino acids, B vitamins, humic acid, plant growth hormones, and dozens of active substances that prevent crop pests and diseases, making it a high-quality organic fertilizer (Liang et al. 2023a; Nascimento et al. 2023). The application of BS has a good regulatory effect on the soil environment. Nutrient elements can improve soil fertility, while organic matter and biologically active substances can regulate soil physicochemical properties and physical structure, enhancing soil water and nutrient retention, thereby increasing crop yields (Zhang et al. 2018; Ibrahim et al. 2020). Soil total nitrogen includes nitrate nitrogen, ammonium nitrogen, and alkali-hydrolyzed nitrogen. The content of soil available nitrogen, namely nitrate nitrogen and ammonium nitrogen, reflects the nitrogen supply capacity of soil and nitrogen absorption by plants. In addition, alkali-hydrolyzed nitrogen is not easy to leach, and it is sensitive to soil nitrogen application, which is an important index to evaluate soil nutrients (Chen et al. 2016).

Garlic is an essential vegetable, seasoning, and medicinal material, with China ranking first globally in terms of garlic planting area, total production, and export volume. Allicin, a compound found in garlic, has been demonstrated to have negative effects on the growth of tumor cells (Brugnoli et al. 2021; Melguizo-Rodríguez et al. 2022). Liang et al. (2023b) showed that the combined application of biogas slurry and biochar improves soil health and increased soil carbon content. Li et al. (2023) found that a 3% addition of biogas slurry was most beneficial for the growth and disease resistance of cucumbers, without the risk of soil salinization. However, at present, the research on the effect of biogas slurry on garlic is still blank.

Biogas slurry is an organic alternative to chemical fertilizers, and its application could efficiently utilize resources and maintain soil health in future biogas engineering endeavors. Therefore, the objectives of this study are: (i) to elucidate the influence of biogas slurry application rates on garlic growth parameters and nutritional quality, (ii) to investigate the effects of different fertilization treatments on soil nutrients and heavy metal content, and (iii) to analyze the feasibility of substituting chemical fertilizers with biogas slurry for garlic cultivation to enhance crop yields while improving soil quality.

EXPERIMENTAL

Biogas Slurry

The original anaerobic digestion slurry used in the experiment was sourced from the discharge station of a biogas project in Ganzhou City, Jiangxi Province. Biogas slurry was used for this study after anaerobic fermentation (15 days) – photosynthetic bacteria purification (15 days) – struvite precipitation (15 days) – nitrifying bacteria treatment (15 days). The physicochemical properties of the treated biogas slurry are detailed in Table 1.

Table 1. Properties of Biogas Slurry

Experimental Design

The garlic field experiment took place on a farm located in Ganzhou City, Jiangxi Province (24°90′N, 115°04′E), situated within a subtropical monsoon climate zone characterized by mild weather conditions. The area experiences an annual average temperature of 19.8°C and an average annual precipitation of 1320 millimeters.

The garlic field experiment involved five treatment groups: (a) G-0: Organic fertilizer was applied as a one-time base fertilizer during the embryonic stage, and no fertilizer was applied during the growing stage. (b) G-F: Organic fertilizer was applied as a one-time base fertilizer during the embryonic phase, and urea (about 46.6% nitrogen) was applied exclusively was the growth phase. (c) G-0.5: Organic fertilizer was applied as a disposable base fertilizer during the embryonic period, and anaerobic slurry (about 23.3% nitrogen) was applied alone during the growth period. (d) G-1: Organic fertilizer was applied as a disposable base fertilizer during the embryonic phase and anaerobic slurry (about 46.6% nitrogen) was applied alone during the growth phase. (e) G-2: Organic fertilizer was applied as a disposable base fertilizer during the embryonic phase and anaerobic slurry (nitrogen content of about 93.2%) was applied alone during the growth phase. Garlic received a total of three topdressings throughout the entire growth period. The first topdressing was applied during the embryonic stage (4 to 7 days after colonization). The second and third topdressings were applied during the bolting period, with a 15-day interval between each fertilization. The specific application rates for each treatment group are detailed in Table 2. Each treatment group comprised three replicates, totaling 15 plots. Garlic seeds were planted with a spacing of 8 cm between plants. The plots, measuring 4 m in length and 4 m in width, with an area of 16 m², were arranged in a north-south direction. A 0.5 m wide buffer zone was established between plots. Plot allocation followed a random principle, and on-site planting involved randomizing the layout through a drawing process. The fertilizer used for topdressing is sprayed as foliar fertilizer.

Table 2. Experimental Group Settings of Garlic (L/plot (16 m2))

Garlic Analysis

After the initial fertilization and at harvest, garlic samples were randomly collected. Sampling was conducted within each experimental plot by selecting three random sample plots, each measuring 1 m × 1 m. Within each selected sample plot, the overall growth of garlic was observed, including plant height and leaf counts. Upon collection, the samples were transported to the laboratory, where the root soil was cleaned. Subsequently, the whole fresh weight and edible fresh weight of each garlic plant were measured. The relative chlorophyll content in garlic leaves was determined using a chlorophyll content meter (SPAD 502, Konica Minolta, Japan). Subsequently, a portion of the harvested garlic samples was ground. The vitamin C content was determined using high-performance liquid chromatography, following the method described by Wu et al. (2023), while the allicin content was determined using the mercury nitrate method. The remaining samples were dried at 60 °C, ground, and sieved for protein, soluble sugar, dietary fiber, and heavy metal content analysis. The protein content was analyzed using an automated Kjeldahl nitrogen analyzer (K9860, Hanon, China). The soluble sugar content was determined using the anthrone-sulfuric acid colorimetric method, and the dietary fiber content was determined as described by Li and Cardozo (1994).

The heavy metal content in garlic was determined using the microwave digestion method (Yan et al. 2024). A 0.1 g sifted garlic sample was placed into a polytetrafluoroethylene digestion tubes, to which 5 mL of nitric acid and 1 mL of sulfuric acid were added. The digestion process was conducted in a microwave digestion system (MDS-6G, Sineo, China), with three stages of temperature increase: the first stage at 150 °C for 10 min, the second stage at 180 °C for 5 min, and the third stage at 200 °C for 20 min. Subsequently, the solution was cooled to 60 °C and diluted with deionized water. After filtration through a 0.22 μm membrane, the heavy metal content was analyzed using an inductively coupled plasma mass spectrometer (OPTIMA 8000, PerkinElmer, USA). The plant heavy metal enrichment factor was calculated using Eq. 1,

(1)

where Cgalic is the concentration of a specific heavy metal in plant leaves (mg/kg), and Csoil is the concentration of the same heavy metal in the soil (mg/kg).

Soil Analysis

After garlic harvesting, soil samples were collected from each plot area, with three random samples taken per plot. Upon mixing, approximately 1 kg of each sample was retained using the quadrat method. Stones, roots, and other debris were removed, and the samples were transported to the laboratory for air-drying. Subsequently, the samples were ground through a 50-mesh sieve and sealed in bags for nutrient and heavy metal analysis. Soil heavy metal content was determined as described previously. Total nitrogen content was determined using the Kjeldahl method, while total phosphorus content was measured using the molybdenum antimony anti-colorimetric method (Kowalenko and Babuin 2007). Total potassium content was determined using flame photometry. Available nitrogen was determined using the diffusion absorption method, and available phosphorus was determined using the sodium bicarbonate extraction-molybdenum antimony anti-spectrophotometric method. Available potassium content was determined by ammonium acetate-flame photometry. Organic matter was analyzed using a total organic carbon/total nitrogen analyzer (Multi N/C 3100+HT 1300, Analytik Jena, Germany).

Data Analysis

All data were statistically analyzed and calculated using Excel 2016 (Microsoft, Redmond, WA, USA). Figures and tables were generated using Origin 2021 software (OriginLab, Northampton, MA, USA). Data variability analysis and principal component analysis (PCA) were performed using IBM SPSS Statistics 22.0 software (SPSS, Inc, Chicago, IL, USA).

RESULTS AND DISCUSSION

Garlic Growth

The growth parameters of garlic following different fertilization treatments are shown in Fig. 1, both after the initial fertilization and at harvest. After the initial fertilization, compared to the chemical fertilizer (G-F group), the growth indexes (plant height, leaf number, SPAD value, fresh weight and edible fresh weight) of all treatment groups were higher than those of blank group (G-0), except the fresh weight and edible fresh weight of the G-F and G-2 group. Compared with the G-F group, G-0.5 and G-1 group improved the growth indexes of garlic to different degrees.

Fig. 1. Effects of different fertilization treatments on (a) plant height, (b) number of leaves, (c) SPAD value in leaves, (d) fresh weight, (e) fresh weight of edible parts, (f) yield and (g) growth photograph of garlic after the first fertilization and after harvest

Among these, the growth rate of G-1 surpassed that of the G-0.5 group (excluding plant height), whereas the G-2 group exhibited only greater plant height compared to the G-F group, with other indicators being lower than the G-F group (Fig. 1). This suggests that, during garlic growth, the application of biogas slurry was more advantageous for promoting garlic growth compared to chemical fertilizer. The initial fertilization coincides with the vigorous growth stage of garlic, a critical period for its development. Biogas slurry, being richer in nutrients compared to chemical fertilizer, facilitates increased absorption of nitrogen, micronutrients, and other nutrients by garlic, thereby enhancing its growth (Cui et al. 2021).

At harvest, the plant heights of the G-0, G-0.5, and G-2 group were 60.57, 65.79, and 61.95 cm, and the average leaf numbers were 5.06, 5.52 and 5.98, respectively. The plant height and leaf number of garlic were significantly increased by nitrogen fertilizer application during growth period. As the slurry application increased, the garlic plant height decreased, while the leaf count gradually increased (Fig. 1a and Fig.1b). This study employed foliar spraying of biogas slurry, allowing for rapid nutrient absorption by the leaves to support garlic growth. Thus, foliar spraying can swiftly supplement the nutrients lacking in garlic. The nutrient content in slurry exceeded that in chemical fertilizers, and slurry application promoted garlic’s nutrients absorption. The G-2 group exhibited the highest leaf numbers, as sufficient nitrogen serves as the basis for leaf formation and improved photosynthetic efficiency.

Varying concentrations of biogas slurry had no significant effect on increasing chlorophyll in mature garlic (Fig. 1c). After harvesting, the order of single plant fresh weight for each group was as follows: G-1 > G-2 > G-0 > G-0.5 > G-F (Fig. 1d), and the order of edible fresh weight for each group was as follows: G-2 > G-1 > G-0 > G-0.5 > G-F (Fig. 1e). The garlic yield of the G-1 group was 1.06 kg/m3, representing increases of 29.3% and 11.6% compared to the G-0 and G-F groups, respectively. The highest garlic yield was observed in the G-2 group, reaching 1.12 kg/m3, representing increases of 36.6% and 17.9% compared to the G-0 and G-F group, respectively (Fig. 1f).

Compared with the G-0 group, the yield of all treatment groups was increased to different degrees, which indirectly indicated the importance of nitrogen for garlic growth. Compared with chemical fertilizer, the growth advantage of the G-1 and G-2 groups of garlic was more obvious, and the yield increase was larger. Additionally, garlic yield increased with increasing slurry concentration. Apart from plant height, the other indicators of the G-1 group were the highest among the treatment groups. This indicates that the application of biogas slurry with the equivalent nitrogen content as chemical fertilizers was more suitable for garlic growth, with the highest nitrogen utilization in the slurry at this stage. The G-F group exhibited the lowest yield (0.95 kg/m3), possibly due to the prolonged single application of chemical fertilizer in the soil, leading to soil degradation, severe compaction, inhibition of plant nutrient absorption, and antagonistic effects on soil micronutrients, resulting in incomplete expression of its immediate effect (Iqbal et al. 2020). Following biogas slurry application, the nutrients supply becomes more adequate and balanced, facilitating the adjustment of soil nutrient supply. Moreover, the antibiotics and bioactive substances in biogas slurry can inhibit the growth of pathogens and pests on leaf surfaces, aiding in disease and pest prevention and control, thereby further increasing garlic yield. However, it is worth noting that nitrogen is a major influencing factor in plant growth, and excessive field irrigation with biogas slurry, resulting in nitrate nitrogen exceeding the plant’s absorption capacity, may adversely affect crop growth (Chen et al. 2013).

Garlic Nutrients

Figure 2 illustrates the impact of different fertilization methods on the nutritional components of garlic. Compared with the G-0 group, the nutritional composition of garlic in other treatment groups did not increase significantly (p > 0.05), on the contrary, the vitamin C content and dietary fiber content in the G-0.5 and G-2 group were significantly decreased (p < 0.05), which may be due to the influence of heavy metal content in biogas slurry. As the application rate of biogas slurry increased, the contents of garlic protein, vitamin C, dietary fiber, reducing sugars, and allicin exhibited a trend of initially increasing and then decreasing. This observation suggests that the accumulation of nutritional components within garlic plants is dependent on the concentration of biogas slurry. Compared to the G-0.5 and G-2 groups, the G-1 group exhibited the most favorable conditions for the accumulation of nutritional components within garlic plants, albeit less effective than the G-F group. The relatively low protein content in the G-0.5 group, measuring at only 5.81 ug/g, may be attributed to the reduced application of nitrogen, insufficient to support extensive protein synthesis in garlic.

Fig. 2. Effects of different fertilization treatments on the content of nutrients in garlic: (a) protein and vitamin C; (b) soluble sugar, dietary fiber and allicin

Similarly, the protein content in the G-2 group, also registering at 5.84 ug/g, could be a consequence of excessive nitrogen application leading to the accumulation of nitrate ions in plant tissues, inhibiting the protein conversion process (Htwe and Ruangrak 2021). Additionally, copper (Cu) in biogas slurry is crucial for protein synthesis, and higher Cu content may exert toxic effects on garlic (Shin et al. 2012). The application of different concentrations of biogas slurry has no significant difference in the accumulation of soluble sugars in garlic plants (p < 0.05). Compared to the fertilizer group, the content of soluble sugars in garlic decreases after applying biogas slurry, and the reasons for this decrease require further investigation. While the application of an appropriate amount of biogas slurry to garlic may reduce dietary fiber content while enhancing its palatability, the resultant decrease is insufficient to impact its nutritional quality significantly. Compared with the G-0 group, application of nitrogen fertilizer significantly increased allicin content in garlic. Compared to the G-F group with the equivalent nitrogen content, the content of allicin in the G-1 group garlic increased by 26.0% (Fig. 2). Allicin, a sulfur-containing compound, is generated through the decomposition of alliin by alliinase in garlic. The synthesis of alliin necessitates amino acids and sulfur elements. Given the complex composition of biogas slurry, containing higher levels of sulfur elements and amino acids compared to chemical fertilizers, biogas slurry can promote the synthesis of alliin in garlic, with the application of biogas slurry at the equivalent nitrogen as chemical fertilizers being most conducive to the synthesis and accumulation of alliin in garlic.

The heavy metal content in garlic leaves serves as an indicator reflecting both plant physiology and soil conditions, playing a crucial role in assessing food safety and environmental risks. Table 3 shows the content of various metals in garlic leaves following different fertilization treatments. Compared with the G-0 group, the application of nitrogen fertilizer promoted the accumulation of Mg and Zn in garlic leaves, while the application of biogas slurry only promoted the accumulation of Zn. Compared with G-F group, the application of biogas slurry promoted the absorption and accumulation of Ca and K in garlic leaves but it was detrimental to the accumulation of Mg and Zn elements. With an increase in the application level of biogas slurry, the Mg, Ca, and Zn contents showed a trend of initially increasing and then decreasing, peaking at 546, 2450, and 18.8 mg/kg DM in the G-1 group, respectively, while the K content increased with the increasing application of biogas slurry. Biogas slurry had a significant impact on the absorption and accumulation of minerals in garlic leaves, and as the application of biogas slurry increased, the influence gradually diminished. This phenomenon may stem from the dilution effect of biogas slurry on minerals, leading to decreased absorption and accumulation by plants (Gupta et al. 2016). For heavy metals, compared to the G-F group, the application of biogas slurry reduced the accumulation of Cr and Pb in garlic, while significantly increasing the content of Cu, As, Se, and Cd in garlic leaves. Furthermore, with an increase in the application of biogas slurry, the content of Cu, Cr, and As in garlic leaves increases, reaching maximum values in the G-2 group of 4.95, 0.42, and 0.4 mg/kg DM, respectively (Table 3). The higher the dosage of biogas slurry, the higher the heavy metal content of garlic, which may be due to the high heavy metal content of biogas slurry itself. After applying biogas slurry, the content of Cr, Zn, As, Se, Cd, Hg, and Pb in garlic leaves does not exceed the limits specified in the national standard GB 2762-2022, indicating that the application of biogas slurry is safe for garlic cultivation. However, long-term application warrants further investigation into the accumulation of heavy metals.

Table 3. Effects of Different Fertilization Treatments on Metal Content in Garlic Leaves

The average heavy metal enrichment coefficient of the G-1 group was 4.66, while the enrichment coefficient of G-F with the equivalent nitrogen content was 2.16 (Fig. 3).

Fig. 3. Enrichment coefficient of heavy metals in garlic leaves under different fertilization treatment

The higher enrichment coefficient associated with biogas slurry application compared to the fertilizer group suggests that biogas slurry facilitated the absorption of heavy metals from the soil by garlic, thereby aiding the transfer of heavy metals from the soil to the edible part of garlic (i.e., the leaves), particularly for Cd and As metals. The enrichment coefficients of Pd and Cr in different treatment groups were all below 1%, significantly lower than the other heavy metals, while the soil contained higher concentrations of Pb and Cr. This suggests that elevated levels of heavy metals in the soil can inhibit their migration from the soil to vegetables. Moreover, as the application rate of biogas slurry increased, the average enrichment coefficient of heavy metals decreased, possibly due to the high water content of biogas slurry, which dilutes the concentration of heavy metals. The enrichment coefficient of the G-F group was lower than that of the G-0 group, which may be due to a small amount of heavy metals contained in the fertilizer itself.

Heavy Metals and Nutrients in Soil

The experimental results from earlier sections indicate that compared to the application of chemical fertilizers, the application of biogas slurry exerted a more substantial promoting effect on the growth and development of garlic, culminating in a more pronounced increase in yield. However, there was a potential risk of heavy metal accumulation in the soil from biogas slurry application. To ensure the safety of biogas slurry application on soil, this study investigated the heavy metal content (Pb, Cd, Cr, Cu, Hg, and As) in the soil after different fertilization treatments, as shown in Table 4. Compared with the G-0 group, there was no significant change in soil heavy metal content in the G-F group (p < 0.05), yet there remained a risk of heavy metal accumulation with long-term application. Compared to the G-F group, the application of biogas slurry at the G-1 group significantly increased the As and Cu content by 210% and 27.7%, respectively. Se and Cd increased from “undetected” to 0.27 and 0.09 mg/kg DM, respectively, while the Cr content decreased by 27.4%. The differences in the content of other heavy metals were not significant (p > 0.05), and all heavy metal content remained within the limits stipulated by the national standard GB 15618-2018. Consequently, biogas slurry can effectively replace chemical fertilizers and applying appropriate amounts of anaerobic slurry can increase the yield of garlic and other vegetable crops.

Table 4. Content of Heavy Metals in Soil After Different Fertilization Treatments

The nutrient content of soil after different fertilization treatments is shown in Fig. 4. Compared with the soil background, the content of soil organic matter in all treatments increased significantly after application of nitrogen fertilizer. The soil organic matter content in the G-1 group was the highest, reaching 40.6 g/kg, marking a 34.7% increase over the G-0 group and a 17.6% increase over the G-F group (Fig. 4a).