Biochar promotes soil aggregate stability and associated organic carbon sequestration and regulates microbial community structures in Mollisols from Northeast China

. Since the 1950s, heavy plowing of Mollisols, combined with a lack of organic matter intake, has resulted in severe soil degradation in Northeast China. The use of biochar in combination with 10 fertilizer is a sustainable method of improving soil quality. In this paper, we conducted field experiments to explore the response of the stability mechanism of the soil aggregate, the dynamic properties of organic carbon, and changes in the microbial community structure to biochar. The biochar input levels were C1, C2, and C3 (9.8, 19.6, and 29.4 Mg C ha -1 , respectively), while the nitrogen (N) fertilizer rates were N1/2 (300 kg N ha -1 ) and N (600 kg N ha -1 ). Results indicated that biochar combined with N fertilizer 15 effectively increases soil carbon storage and aggregates stability ( P < 0.05). And C2N treatment increased the aggregate contents of the > 2 mm and 0.25–2 mm fractions by 56.59% and 23.41%, respectively. The phospholipid fatty acids (PLFA) analysis revealed that microbial community structure was effectively improved with biochar combined with N fertilizer application ( P < 0.05). The F/B ratio increased by 25.22% and the gram-positive (Gm + )/gram-negative (Gm − ) ratio by 4.65% under the 20 C2N1/2 treatment. This study concluded that the response of Mollisols to biochar is primarily determined by the interplay of aggregate, organic carbon

technology for recycling straw while also potentially improving agricultural soils (Islam et al., 2021).
Biochar can enhance SOC storage, soil granular structure, and cation exchange capacity (Hu et al., 2020, Mete et al., 2015. For example, Wang et al. (2019) discovered that biochar improved the structural stability of Latosols in southern China. The aggregate mean weight diameter (MWD) and geometric mean diameter (GMD) was improved by 36.3 and 28.3%, respectively. Furthermore, Xiu et al. (2019) 60 investigated the effect of corn stalk biochar application dose on Albic soils in northern. They discovered that a high biochar application level reduced the bulk density of Albic soils by 9.93% while increasing the pH value. Biochar was also found to significantly improve soil granular structure and organic carbon aggregation . Thus, biochar had a favorable influence on soil quality and aggregation in these acidic soils, which could be attributed to the liming activity of biochar treatments on those acidic 65 soils and the neutralization of the soil pH, which consequently had a significant effect on soil aggregation (Islam et al., 2021). However, there are still some researchers who believe that there is no significant effect (Zhang et al., 2015). Due to the low quantity of biochar minerals and inorganic nitrogen, several studies have indicated that only combination application with other fertilizers can improve soil fertility (Song et al., 2020).  proposed that an 8-year manure amendment could recover soil 70 nitrogen supplying capacity of lightly eroded Mollisols to natural levels. Fungo et al. (2017) conducted a two-year field trial in the impoverished Ultisol of western Kenya and found that biochar combined with urea increased MWD by 13%, whereas biochar alone was less effective. Biochar combined with an organic/inorganic fertilizer has the potential to improve soil fertility (Li et al., 2020), promote plant growth (Aneseyee et al., 2021;Mete et al., 2015), and carbon storage potential (Wang et al., 2019). 75 Despite these benefits, a quantitative understanding is scarce of how combined effects between biochar and nitrogen fertilizer contribute to soil fertility by modifying microbe-soil interactions in agroecosystems.
Principal ecological activities in soil, including organic matter formation and decomposition, nutrient cycling, and soil aggregate size redistribution are all controlled by soil microbial populations (Chen et 80 al., 2022;Trivedi et al., 2017). Phospholipid fatty acids (PLFAs) are the main components of living cell membranes, which play an important role in maintaining cellular fluids, nutrient transportation, elimination of metabolites, etc. Changes in their components can more accurately express the response of soil microbial biomass and community structure to environmental disturbances (Zhang et al., 2013).
The structure of the microbial community is closely related to the change in soil function (Ng et al., 85 2014). The higher the ratio of soil fungal to bacterial fatty acids, the more sustainable and stable the soil ecosystem . High Gm + /Gm − bacterial ratios facilitate soil organic carbon accumulation. Soil total nitrogen (TN) content is the main driver of variations in the community composition . Wang et al. (2021) discovered that after using biochar in rice fields, the abundance of bacteria (B) and fungi (F) increased by 102 and 178%, respectively, which was likely 90 related to an increase in soil total organic carbon (TOC), TN, and rice biomass. According to the study of , the improvement of microbial community structure by biochar was determined by the Gm + / Gm − and F/B in the paddy soil of central-southern China. In addition, Tian et al. (2016) investigated the mechanism of combined effects between biochar and fertilizer addition on microbial community and soil organic matter cycling in heavy loam soil. It was found that the addition of biochar 95 alone did not significantly improve microbial community structure and that its effect on microbial community structure was dependent on fertilization. The ability of biochar and nitrogen fertilizer to stimulate microbial activity is regulated by the soil conditions and application rates (Palansooriya et al., 2019).
Soil organic carbon sequestration and microbial activity are critical for soil health and quality regulation. 100 However, the beneficial effects of biochar on soil aggregate, associated SOC, and microbial activity have been observed primarily in nutrient-poor acidic soils (e.g. Ultisol and Albic soils), and relevant studies on Mollisols in Northeast China have been limited. Furthermore, studies on the combined application of biochar and nitrogen fertilizer are insufficient, limiting the scope of production practice and theory. Therefore, this study using the northeast Mollisols as a pilot, the objectives are to (1) explore the effects 105 of three biochar gradients combined with N fertilizer on the size, proportion, stability, and carbon content of Mollisols aggregates; (2) explore the influence mechanism of biochar on microbial population structure and identify the major determinants for microbial community composition changes; (3) develop scientific and effective field management measures for Mollisols by improving the structure of soil aggregates and microbial communities. 110 2 Materials and methods

Site description
The field experimental site was located at the test base of the Northeast Institute of Geography and Agroecology, Jilin Province (43° 59' 51" N, 125° 24' 5" E). The annual average temperature is 4.6°C , the precipitation is 600-700 mm, and the frost-free period during the whole year is 140-150 d. For many 115 years, continuous maize cropping has been carried out in conventional tillage patterns. The soil of the field was classified as Mollisols (Mei et al., 2021). The experimental soil pH was approximately 6.06, TN was 1.26 g kg -1 , available phosphorus was 26.78 mg kg -1 , available potassium was 133.54 mg kg -1 , and organic matter was 26.72 g kg -1 . The biochar was created by pyrolyzing corn straw at 400-500°C for 4 h under anaerobic conditions. The biochar had a mean particle diameter of 0.003-3.5 mm, a surface 120 area per volume of 0.7 m 2 g -1 , and an ash concentration of 45% (Biochar particles need to pass through a 2 mm sieve before application). Also, the biochar with a pH of 9.16, the total carbon content was 62.64%, and the C/N was 39.08. The fertilizer was high-quality urea that was produced by Erdos Yi Ding Ecological Agriculture Development Co. Ltd., the TN was ≥ 46%, and the particle size range was 1.18-3.35 mm. 125

Field experimental design
A split zone design was adopted for the field experiment and three biochar input levels were set: 9.8 Mg ha -1 (C1), 19.6 Mg ha -1 (C2), and 29.4 Mg ha -1 (C3). Nitrogen was applied as a basal fertilizer at rates of 300· kg N ha -1 (N1/2) and 600 kg N ha -1 (N). The CK treatment was used as a control. In total, ten treatments were studied: CK, C1, C2, C3, C1N1/2, C2N1/2, C3N1/2, C1N, C2N, and C3N. Each 130 treatment was performed on a plot with the dimensions 3.9 × 6.5 m, and each treatment plot had a 1 m buffering zone. A randomized block design was used to conduct the three replicate plots. Biochar with N fertilizer was applied to the soil in April 2013 and 2021, and corn was sown in May 2013 and 2021.

Soil bulk density and water content
On October 29, 2021, after the corn harvest was complete, soil samples were obtained from each plot 135 using the five-point sampling method, which involved taking 1 kg of soil samples from each plot.
Undisturbed soil columns were collected using a soil drill and were placed into ziplocked bags after the removal of plant and animal residues. Some of the soil was promptly refrigerated at 4°C for PLFA measurement. A 5 mm mesh screen was used to remove the water-stable soil aggregates from the rest of the sample, which was then allowed to dry naturally. For the determination of the bulk TOC, subsamples 140 of 2 mm soil particles were passed through a 0.15 mm filter after being air-dried. The TOC in the aggregate fractions was determined by K2Cr2O7 titration . Next, the surface (0-10 cm) and bottom (10-20 cm, 20-40 cm) soils were sampled with a cutting ring (V = 100 cm 3 ) and dried at 105 °C for 24 h to measure the soil bulk density and water content using the following formulae: where X is the field water holding capacity (%), ρb is the soil bulk density (g cm -3 ), m is the dry soil weight (g), v is the cutting ring volume (cm 3 ), m2 is the total weight of the cutting ring and soil after 2 h on dry sand, and m1 is the total weight of the cutting ring and soil after drying.

Soil water-stable aggregate analysis and calculation 150
In this experiment, the soil aggregates were fractionated utilizing a modified version of the wet sieving method which was given by Zhang et al. (2018). The dry soil sample (100 g) was uniformly coated on automatic vibrating sleeve screens of 2, 0.25, and 0.053 mm in diameter.
The formula for calculating the mass fraction of the water-stable aggregates is as follows: (3) 155 where Wt is the percentage of the component weight of the ith-sized aggregate.
The MWD and GMD represent the size distribution of the soil aggregates. The larger the values, the more stable the aggregates. The formulas are as follows: where j is the aggregate size, Xj is the average diameter of the particle size, Wj is the ratio of the aggregate sample weight of each particle size on the screen, Xi is the average diameter of a size i aggregate, Mi is the weight of a size i aggregate, and Mt is the total weight of all the aggregates. The aggregate content was determined as follows: where R0.25 is the aggregate content (%) with an aggregate size of > 0.25 mm, Mr > 0.25 is the weight of the soil aggregate that is> 0.25 mm, and MT is the total weight of all the aggregate fractions.
The formula for the soil carbon contribution rate of each aggregate grain size is as follows: where Cc represents the contribution rate of each particle size aggregate to the carbon level in the soil 170 sample, wi is the weight percent (%) of the i-sized aggregate component, Ci is the organic carbon content of the soil aggregate at size i, and Cs represents the soil TOC content.
The concentration of the target PLFAs in the sample was calculated as follows: 180 where CPLFA is the concentration of the target PLFA (nmol g -1 ), FPLFA is the peak area for the PLFAs, FIS is the area of the internal standard peak, CIS is the internal standard concentration (25 ng μl -1 ), MPLFA is the molecular weight of the target PLFA, V is the sample dissolution volume (120 μl), and m is the soil weight (4 g). 185

Statistical analyses
IBM Statistics SPSS 22.0 software was used to test the data normality, homogeneity and principal component analysis (PCA). Significant differences among treatments were calculated by Duncan's Multiple Range test with P < 0.05. If the data did not meet the criteria, a nonparametric Kruskal-Wallis test was performed to determine the statistical significance. Two-way ANOVA was conducted on a subset 190 of treatments to examine the effects of biochar and organic fertilizer on bulk density, soil water content, soil aggregates, C storage and microbial community abundance. Canoco 5 (Windows Release 5.02 trial version) software was used for redundancy analysis (RDA).

Soil physical properties 195
The biochar had a substantial impact on the soil (0-10 cm) bulk density (P < 0.05; Fig. 1), but its coupling effect with N fertilizer was not significant. Also, soil bulk density showed distinct regularities in all profiles and increased with soil depth. The C2N1/2 treatment had the greatest improvement effect of all treatments, and the soil bulk densities of the 0-10, 10-20, and 20-40 cm layers decreased by 13.00%, 8.62%, and 2.94%, respectively. The surface soil (0-10 cm) had the highest moisture content under the 200 CK treatment, while the 10-20 cm soil had the lowest water content. Additionally, there was a substantial positive relationship between biochar application amount and the soil water content in the profile (P < 0.01; Fig. 1), with the C3 treatment improving the most when compared to the CK. Furthermore, the soil moisture content increased by 15.12%-35.39%. The two-factor ANOVA (Table S1) showed that biochar significantly improved soil water content (P < 0.01) and that biochar contributed significantly to soil 205 bulk density and water content.
First, the number of large aggregates fractions was lower in the 20-40 cm soil layer than in the 0-10 cm soil layer. Second, the biochar considerably increased the percentage of large aggregates (11.59%-50.40%) while decreasing the percentage of < 0.053 mm aggregates (5.12%-38.66%). Third, biochar 215 combined with N fertilizer had a synergistic effect, and the proportion of large aggregates continued to increase (38.98%-56.59%) before stabilizing.
According to the interactive analyses, N fertilizer had a greater effect on the fraction of large aggregates in the profile (Table S2). The C2N treatment increased the > 2 and 0.25-2 mm fractions of soil aggregates by 56.59% and 23.41%, respectively. Furthermore, the proportions of aggregates 0.053-0.025 and < 220 0.053 mm decreased by 4.09% and 43.64%, respectively. The C2N treatment had the highest growth rate of large aggregates within the 0-10 cm layer, which was 3.66 and 20.16% higher than that of the C2N1/2 and C2 treatments, respectively. Furthermore, as soil depth increased, the water-stable aggregates were gradually replaced with 0.053-0.25 mm sized aggregates (35.95%-46.42%).

Total organic carbon distribution in the bulk soil and aggregate fractions
The average TOC content of the surface layer was 20.26% higher than that of the 20-40 cm soil layer (Fig 4). The TOC content was significantly correlated with the application rates of the biochar and nitrogen fertilizer (P < 0.01). Among all the treatments, the C3N treatment in comparison to the CK 240 resulted in the greatest increase in organic carbon content, and the TOC increased by 35.59%, 30.62%, and 29.53% in the soil profile from surface to bottom. The TOC was significantly associated with aggregate fractions of > 2 mm and 0.25-2 mm but inversely 245 associated with fractions of 0.25-0.053 mm and 0.053 mm aggregates (Fig. 5). We also compared the TOC of the particle size components of the various aggregates under different biochar treatments (Fig. 6 a, b, and c) and found that large aggregates had higher carbon content than microaggregates. The C3N1/2 treatment increased the TOC content in the > 2 mm, 2-0.25 mm, 0.25-0.053 mm, and < 0.053 mm fractions by 36.89%, 20.39%, 15.41%, and 16.14% respectively (P < 0.05). Furthermore, the 0.25-2 mm 250 aggregate fractions contributed the most to TOC, followed by the > 2 mm fractions (Fig. 6 d, e, and f).
There was no significant change in the contribution of biochar with full N fertilizer (N) to TOC compared to biochar with reduced N fertilizer (N1/2) treatment. The bars indicate the standard error.

Microbial community structure
The PLFAs of microorganisms (i.e., bacteria, fungi, actinomycetes, Gm + bacteria, and Gm − bacteria) in the soil were identified (Fig. 7). The biochar treatment resulted in the highest increases in F/B and 265 Gm + /Gm − proportions of 28.17 and 7.91%, respectively (Fig. 7 g and h). Also, the two-factor ANOVA The RDA was performed to determine the relationship between soil environmental change and the PLFA response variables (Fig. 8). The two RDA axes were significant (P < 0.05), accounted for 94.12% of the overall variation in the soil microbial characteristics. Soil bulk density was the most significant variable, accounting for 62.61% of the microbial community characteristics, followed by MWD, soil moisture, TOC, R0.25, and GMD, all of which were significantly correlated with the microbial community 275 composition and explained 15.90%, 13.42%, 4.01%, 2.83%, and 1.28% of the various rates of microbial PLFAs, respectively. The PCA was used to evaluate the effects of various treatments on the soil traits in Northeast China (Table   1, Table S4). The results showed that the first 3 principal components (F1-F3) explained 90.13% of the 290 total variance. The higher the F value, the better the improvement effect, and the C2N1/2 treatment was optimal.

The effects of the biochar and nitrogen fertilizer treatments on soil physical properties 305
Soil compaction often makes soil structure destroying and soil degrading which is not conducive to the transport of water, air and solutes in the soil and the propagation of soil microorganisms (Gaia et al., 2020). Our results demonstrated that soil bulk density of Mollisols soil had a negative correlation with biochar application rate. The C3 treatment reduced soil bulk density by up to 12.69% (Fig. 1). The bottom soil bulk density was on average 18.88% higher than that of the surface soil by the biochar amendment, 310 though the improvement in the bottom soil bulk density was not significant (P > 0.05), which are consistent with Xiu et al. (2019). This trend might be due to the unique properties of biochar, such as complex microporous structure, large specific surface area, and light texture, et al. (Zhang et al., 2015).
The biochar had a slow and gradual effect on the soil improvement. According to Chaganti et al. (2015), the biochar in the soil will gradually migrate to the lower soil over time due to natural factors and human 315 activities. Also, Luo et al. (2020) concluded that biochar was often applied to the surface layer, resulting in a greater decline in the bulk density of the surface soil than the underlying soil. This suggests that biochar has a great benefit in ameliorating soil compaction problems in modern agriculture. Our study also found a considerably strong correlation between the soil water content of the Mollisols and the amount of biochar applied, particularly in the surface soil. The soil water content gradually increased 320 with increasing rate of biochar application rate. This improvement was the largest during the single application of biochar, with an average increase of 18.07%. The two-factor ANOVA showed that the increase in soil water content was mainly attributed to biochar, though there was also a synergistic effect of biochar and nitrogen fertilizer on the increase in soil moisture content (Table S1). An et al. (2022) explained this phenomenon by CT scanning techniques, where the addition of biochar reduced soil 325 porosity, reduced pore size and increased water retention, meaning that water was stored in the smaller pores of the soil and drainage was delayed. The porosity, hydrophilic domains, and huge specific surface area of biochar may be favorable the improvement of water retention (Leonard et al., 2014). However, some studies contradicted this study, and found either reduced water retention capacity (Madari et al., 2017) or no effect (Baiamonte et al., 2015) after biochar application. The variation in the 330 actions may be attributed to biochar properties, soil texture type, climate change, and experimental design and duration. The input of biochar combined with N fertilizer actively participates in the formation of soil large aggregates, which enhances the soil water-holding capacity and alleviates soil erosion by improving the soil aggregate structure (Gaia et al., 2020;Islam et al., 2021).

The effects of biochar and nitrogen fertilizer on soil aggregate distribution and stability 335
Soil aggregate is essential for the performance of soil functions and is primarily responsible for the formation of the soil structure . In this study, biochar increased the formation of large aggregates (> 0.25 mm), especially small aggregates (0.25-2 mm), but decreased the number of microaggregates in Mollisols. Grunwald et al. (2016)

also confirmed this point by treating Haplic
Phaeozem and Gleyic Luvisol with biochar in field experiments. Our findings showed that when biochar 340 was combined with N fertilizer, the fraction of large aggregates steadily increased while the content of the microaggregates and clay particles decreased (Fig. 2). This suggests that surface hydrophobichydrophilic interactions between clay minerals and biochar particles, as well as the biochar ability to integrate with the soil biota, and labile carbon, may all contribute to soil aggregation (Joseph et al., 2010).
Furthermore, surface area, microporous structure, and O/C ratio are key biochar features for binding to 345 organo-mineral complexes, an initial stage in aggregate formation and stability (Du et al., 2017).
According to the findings of this study, the soil aggregate stability increased by 10.9-23.49%, which is consistent with the findings of a meta-analysis (Peng et al., 2015). The initial low TOC level (26.72 g kg -1 ) and protracted field experiments (8 years) could explain this. In a laboratory incubation experiment, biochar application increased soil aggregate stability by 2.33%-6.55% in Albic soil soils of northeast 350 China (Xiu et al., 2019). Our data showed that biochar improves the agglomeration of Mollisols better than Albic soil. According to the MWD (Fig. 3), increased TOC and microbial biomass (Fig. 7) were responsible for the significant increase in aggregation caused by biochar addition. This was also found to be the case in other studies, which found that biochar served as a cementing material, assisting more microaggregates, silt, and clay components to cement together into large aggregates (Xu et al., 2019). 355 The N input enhancing aggregation possible is interpreted in three aspects. Firstly, N fertilizer increased macro-aggregates and MWD stemmed from N-stimulation of root growth (Bai et al., 2021). Roots are a major driver of soil aggregation as they not only constitute the primary source of SOC, but also function as physical binding agents for aggregation (Sokol and Bradford, 2019); Second, the N-increase of plant photosynthesis generally allows plants to allocate more photosynthates to the roots and their associated 360 mycorrhizal fungi, which may physically facilitate the formation of soil macroaggregate through fine root and hyphal enmeshment of microaggregates (Miller and Jastrow, 1990;Bai et al., 2021); Third, the active functional groups in the N fertilizer adsorb soil Ca 2+ and then tend to combine with clay minerals to form a clay-humus complex (Gao et al., 2019). Thus, biochar and N fertilizer application have a longerterm positive influence on aggregate stability. At the same time, it provides a theoretical basis for 365 preventing the thinning of the humus layer and reducing surface runoff and soil erosion (Grunwald et al., 2016). Our findings were in contrast with those of Zhou et al. (2019), who discovered neutral or even antagonistic effects on soil aggregate formation and stabilization due to fewer binding agents produced during the decomposition of recalcitrant biochar. Therefore, the response of soil aggregates to biochar is influenced by the initial TOC, clay content, biochar attributes, application rate, etc. (Peng et al., 2015). 370 As a result, the evaluation results should be thoroughly examined, considering these factors as well as the effect of time in the field.
According to a two-factor ANOVA (Table S2), there was a synergistic effect of biochar blend with nitrogen fertilizer on increasing the proportion of soil large aggregates, which suggests that biochar with organic fertilizer is beneficial for the stability of soil aggregates. Previous research proposed that biochar 375 combined with N fertilizer promotes crop root growth, improves crop root fungi reproductive capacity.
Fungal hyphae and root secretions promote soil aggregate by binding and entangling mineral particles to form large aggregates structures (Islam et al., 2021). Consistently, our data suggest that the increased stability of aggregate results from increased root activity and the significant role of exogenous carbon as a binder for soil particles (Wang et al., 2019). 380

The effects of biochar combined with nitrogen fertilizer on the total organic carbon
In this investigation, the TOC level of the Mollisols increased significantly following biochar application, which is consistent with the results of Dong et al. (2016). More recently, Shi et al. (2020) proposed that the combined application of biochar and nitrogen fertilizer was conducive to soil carbon sequestration, with the cumulative mineralization rate of TOC decreasing by 0.6-1.1% when compared to the CK 385 treatment. These findings can be interpreted in three ways. First, the use of biochar increased soil microbial activity (Fig. 7) and crop yields (Lin et al., 2020), thereby promoting further degradation and transformation of the plant residues, increasing TOC. Second, when added to the soil, biochar with a high organic carbon concentration (34.9%) directly improved the soil's organic matter content. Xiu et al. (2019) found similar results in Albic soil. Third, the enrichment phenomenon of the organic carbon occluded 390 within the large aggregates and small aggregates (Fig.5, 6) were higher than that in the microaggregates.
Increasing the proportion of large aggregates promoted soil aggregate carbon fixation .
The fourth explanation is that biochar has a high inert carbon content, which increased the Gm + /Gm − (Fig. 7) in the decomposition of persistent and complex substrates, indicating that carbon accumulation was greater than carbon decomposition (Dong et al., 2020). Thus, biochar effectively prevented the bulk 395 TOC in the Mollisols from decreasing.
In this study, the TOC concentration was positively correlated with the proportion of large aggregate size ( Fig. 5), which is consistent with the aggregate hierarchy model proposed by Tisdall (1982). Figure 6 shows that the > 0.053 mm fractions had a much higher carbon content than silt and clay, especially in the 0.25-2 mm fraction. Our findings confirmed those of Du et al. (2017) and Dong et al. (2016). 400 These results showed that the combined application of biochar and half amount of N fertilizer was more economically efficient. Under the treatment, the carbon of the < 0.053 mm aggregates in the 0-20 and 20-40 cm soil layers decreased significantly, which could be explained by the finding of Ying (2018) that N fertilization promoted the mineralization rates of primary organic carbon by affecting the soil microbial community. Overall, there was no advantage in increasing the organic carbon content of the 405 soil aggregates from the increased N fertilizer ration. This could be due to the high N content, which caused an imbalance in the soil C/N ratio, affecting the breakdown and turnover of soil organic matter (Kimetu et al., 2010). Therefore, biochar in combination with nitrogen fertilizer as amendments effectively improve the soil aggregates and carbon sequestration.

The effects of biochar combined with nitrogen fertilizer on microbial community biomass and 410
structure Biochar can alleviate the negative effects of soil structure and function degradation on soil microbial activities, particularly when applied in conjunction with nitrogen fertilizer (Oksana et al., 2022).
According published research biochar addition alone did not change the microbial community structure in spring maize fields or rice paddy fields, but when combined with fertilizer, the structure was changed 415 (Luo et al., 2017;Tian et al., 2016). These findings are consistent with our experimental results. Soil F/B and total PLFA contents were significantly increased following biochar and N fertilizer treatments, which may be accompanied by increased C and N cycling and mineralization rates (Khadem et al., 2021). The higher the ratio of PLFA of soil fungi to bacteria, the more stable the soil ecosystem (Thiet et al., 2006).
Compared to Gm − bacteria, Gm + bacteria generally possess a greater proportion of peptidoglycan, which 420 is a relatively decay-resistant soil organic matter (Zhang et al., 2013). The high Gm + /Gm − bacteria ratio means that SOC accumulation is higher than mineralization . Therefore, the effect of biochar and organic fertilizer application on microbial community structure may be more inclined to the retention of easily decomposed organic carbon in northeast Mollisols (Jiang et al., 2016). In this study, C2N1/2 treatment had the best improvement effect on microbial community structure. The N source 425 always an important source of microbial nutrients. However, when N fertilizer application exceeds crop absorption capacity and soil retention capacity, the excess N may be leached to deep soil and pollute groundwater. Meanwhile, the unbalanced C/N ratio also became the main factor for the significant decline of soil biodiversity (Yuan et al., 2017).
The RDA showed that the biomass of fungi, bacteria, actinomycetes, Gm + , and Gm − was positively 430 related to the fraction of large aggregates and negatively linked to the soil bulk density. The increased fungal abundance has been proposed as an important biological factor in soil aggregate formation (Yuan et al., 2015;Zheng et al., 2020). Previous research has shown that aggregate stability and carbon storage are important prerequisites for enhancing microbial communities . Our results showed that the mutual effects of biochar and reduced N fertilizer could effectively affect the abundance of microorganisms, which is attributed to the increased soil C/N content as a result of the applied N fertilizer providing more N sources for microbial decomposition and organic matter utilization (Jia et al., 2020). These findings were consistent with those of Zhang et al. (2021), who discovered that combining biochar with fertilizer significantly increased microbial abundance in the soil sample, implying that the addition of inorganic fertilizer reduced crop N limitation and microbial N immobilization. Furthermore, 440 the TOC and C/N affected the fungal community composition, most likely because fungi were the primary decomposers of TOC (Chen et al., 2013). This conclusion is further confirmed by Sekaran et al. (2019), who found that the amount of soil microbial PLFAs and the ratio of soil carbon to nitrogen were strongly and positively correlated, but biochar and a full dose of N fertilizer had little effect. Based on the sequestration of SOC and the sustainability and stability of the ecosystem, we selected C2N1/2 as the 445 most reasonable biochar ratio.

Conclusion
The field experiments showed that the porous structure of biochar and its carbon source can effectively improve soil structure and carbon storage. The proportion of soil large aggregates and the stability of aggregates were significantly increased. Biochar combined with N fertilizer provided an abundance of 450 living space and nutrients for soil microorganisms, but microbial activity and abundance were limited by carbon input and soil nitrogen availability. The effect of excessive N application was unsatisfactory, which affects the further improvement of soil microbial abundance. This study highlighted that biochar combined with N fertilizer application could be a potential option for mitigating soil degradation, reasonable applicating N, and enhancing soil carbon storage, which would support sustainable use of 455 Mollisols. In the future, we will further investigate the long-term effects of biochar application on soil C and N cycles in agroecosystem.

Credit authorship contribution statement
Jing Sun: Investigation, Experimentation, Data collection and analysis, manuscript writing. Xinrui Lu and Guoshuang Chen: Revise the manuscript. Xiujun Li: Concept and design, Project administration, 460 Funding acquisition. Nana Luo and Qilin Zhang: Investigation, Material preparation, and Experimentation. The final manuscript was read and approved by all authors.