Pervasive soil phosphorus losses in terrestrial ecosystems in China.

Future phosphorus (P) shortages could seriously affect terrestrial productivity and food security. We investigated the changes in topsoil available P (AP) and total P (TP) in China's forests, grasslands, paddy fields, and upland croplands during the 1980s-2010s based on substantial repeated soil P measurements (63,220 samples in the 1980s, 2000s, and 2010s) and machine learning techniques. Between the 1980s and 2010s, total soil AP stock increased with a small but significant rate of 0.13 kg P ha-1 year-1 , but total soil TP stock declined substantially (4.5 kg P ha-1 year-1 ) in the four ecosystems. We quantified the P budgets of soil-plant systems by harmonizing P fluxes from various sources for this period. Matching trends of soil contents over the decades with P budgets and fluxes, we found that the P-surplus in cultivated soils (especially in upland croplands) might be overestimated due to the great soil TP pool compared to fertilization and the substantial soil P losses through plant uptake and water erosion that offset the P additions. Our findings of P-deficit in China raise the alarm on the sustainability of future biomass production (especially in forests), highlight the urgency of P recycling in croplands, and emphasize the critical role of country-level basic data in guiding sound policies to tackle the global P crises.

Survey of Heterodera glycines Population Abundance and Virulence Phenotypes During 2021 to 2022 in Heilongjiang Province.

Soybean cyst nematode (SCN), Heterodera glycines, poses a significant threat to global soybean production. Heilongjiang, the largest soybean-producing province in China, contributes more than 40% to the country's total yield. This province has much longer history of SCN infestation. To assess the current situation in Heilongjiang, we conducted a survey to determine the SCN population density and virulence phenotypes during 2021 to 2022 and compared the data with a previous study in 2015. A total of 377 soil samples from 48 counties representing 11 major soybean-planting regions were collected. The prevalence of SCN increased from 55.4% in 2015 to 59% in the current survey. The population densities ranged from 80 to 26,700 eggs and juveniles per 100 cm3 of soil. Virulence phenotypes were evaluated for 60 representative SCN populations using the H. glycines (HG) type test, revealing nine different HG types. The most common virulence phenotypes were HG types 7 and 0, accounting for 56.7 and 20% of all SCN populations, respectively. The prevalence of populations with a female index (FI) greater than 10% on PI 548316 increased from 64.5% in 2015 to 71.7%. However, the FI on the commonly used resistance sources PI 548402 (Peking) and PI 437654 remained low at 3.3%. These findings highlight the increasing prevalence and changing virulence phenotypes of SCN in Heilongjiang. They also emphasize the importance of rotating soybean varieties with different resistance sources and urgently identifying new sources of resistance to combat SCN.

Warming Leads to Changes in Soil Organic Carbon Molecules Due to Decreased Mineral Protection.

Climate change affects the content and composition of soil organic carbon (SOC). However, warming-induced changes in the SOC compounds remain unknown. Using nuclear magnetic resonance spectroscopy, molecular mixing models, and Fourier transform ion cyclotron resonance mass spectrometry, we analyzed the variations and relationships in molecular compounds in Mollisol with 10-56 g C kg-1 soil-1 by translocating soils under six climate regimes. We found that increased temperature and precipitation were negatively correlated with carbohydrate versus lipid and lignin versus protein. The former was consistent across soils with varying SOC contents, but the latter decreased as the SOC content increased. The carbohydrate-lipid correlations were related to dithionite-citrate-extractable Fe, while the lignin-protein correlations were linked to changes in moisture and pyrophosphate-extractable Fe/Al. Our findings indicate that the reduction in the mineral protection of SOC is associated with molecular alterations in SOC under warming conditions.

Pedogenesis of typical zonal soil drives belowground bacterial communities of arable land in the Northeast China Plain.

Belowground bacterial communities play essential roles in maintaining ecosystem multifunction, while our understanding of how and why their distribution patterns and community compositions may change with the distinct pedogenetic conditions of different soil types is still limited. Here, we evaluated the roles of soil physiochemical properties and biotic interactions in driving belowground bacterial community composition across three typical zonal soil types, including black calcium soil (QS), typical black soil (HL) and dark brown soil (BQL), with distinct pedogenesis on the Northeast China Plain. Changes in soil bacterial diversity and community composition in these three zonal soil types were strongly correlated with soil pedogenetic features. SOC concentrations in HL were higher than in QS and BQL, but bacterial diversity was low, and the network structure revealed greater stability and connectivity. The composition of the bacterial community correlated significantly with soil pH in QS but with soil texture in BQL. The bacterial co-occurrence network of HL had higher density and clustering coefficients but lower edges, and different keystone species of networks were also detected. This work provides a basic understanding of the driving mechanisms responsible for belowground bacterial biodiversity and distribution patterns over different pedogenetic conditions in agroecosystems.

Conversion of steppe to cropland increases spatial heterogeneity of soil functional genes.

The microbiome function responses to land use change are important for the long-term prediction and management of soil ecological functions under human influence. However, it has remains uncertain how the biogeographic patterns of soil functional composition change when transitioning from natural steppe soils (NS) to agricultural soils (AS). We collected soil samples from adjacent pairs of AS and NS across 900 km of Mollisol areas in northeast China, and the soil functional composition was characterized using shotgun sequencing. AS had higher functional alpha-diversity indices with respect to KO trait richness and a higher Shannon index than NS. The distance-decay slopes of functional gene composition were steeper in AS than in NS along both spatial and environmental gradients. Land-use conversion from steppe to farmland diversified functional gene profiles both locally and spatially; it increased the abundances of functional genes related to labile carbon, but decreased those related to recalcitrant substrate mobilization (e.g., lignin), P cycling, and S cycling. The composition of gene functional traits was strongly driven by stochastic processes, while the degree of stochasticity was higher in NS than in AS, as revealed by the neutral community model and normalized stochasticity ratio analysis. Alpha-diversity of core functional genes was strongly related to multi-nutrient cycling in AS, suggesting a key relationship to soil fertility. The results of this study challenge the paradigm that the conversion of natural to agricultural habitat will homogenize soil properties and biology while reducing local and regional gene functional diversity.

Structural convergence of maize and wheat straw during two-year decomposition under different climate conditions.

Straw decomposition plays an important role in soil carbon sequestration. Litter quality and climate condition are considered to be key factors that regulate straw decomposition. This study investigated the decomposition characteristics of wheat and maize straw under cold temperate, warm temperate, and midsubtropic climate conditions, and examined whether the chemical structures of straw residues became similar during decomposition under different climate conditions. Straws were put in 0.074-mm-mesh size litter bags to exclude soil fauna and buried in black soil plots at three experimental stations located in the aforementioned climate regions to rule out the impact of soil type. The decomposition rate constants of wheat straw and maize straw increased linearly with temperature, and the former was more sensitive to temperature. Climate conditions and straw quality had marked effects on the residual material structure in the first half year of decomposition, but then decreased. Wheat and maize straw showed common decomposition characteristics with a decrease of O/N-alkyl carbons and di-O-alkyls, and a simultaneous increase of alkyl carbons, aromatic carbons, aromatic C-O groups, and COO/N-C ═ O groups. Overall, the results indicated that the chemical compositions of the two types of straw became similar after 2-year decomposition under different climate conditions.

Soil Organic Carbon Mineralization and Its Temperature Sensitivity under Different Substrate Levels in the Mollisols of Northeast China.

Soil organic carbon (SOC) mineralization plays an important role in global climate change. Temperature affects SOC mineralization, and its effect can be limited by the substrate available. However, knowledge of the effects of temperature and substrate quality on SOC mineralization in the Mollisols of Northeast China is still lacking. In this study, based on a spatial transplant experiment, we conducted a 73-day incubation to examine the effects of temperature on SOC mineralization and its temperature sensitivity under different carbon levels. We found that the SOC content, incubation temperature and their interaction had significant effects on SOC mineralization. A higher SOC content and higher incubation temperature resulted in higher SOC mineralization. The temperature sensitivity of SOC mineralization was affected by the substrate quality. The temperature sensitivity of SOC mineralization, showed a downward trend during the incubation period, and the range of variation in the Q10 declined with the increment in the SOC content. The study suggested that there was a higher SOC mineralization in high levels of substrate carbon when the temperature increased. Further, SOC mineralization under higher SOC contents was more sensitive to temperature changes. Our study provides vital information for SOC turnover and the CO2 sequestration capacity under global warming in the Mollisols of Northeast China and other black soil regions of the world.

Vulnerability and driving factors of soil inorganic carbon stocks in Chinese croplands.

The long-term stability of soil inorganic carbon (SIC) and its minimum contribution towards global C cycle has been challenged, as recent studies have showed rapid decreases in SIC stocks in intensive agricultural systems. However, the extent of SIC losses and its driving factors remains unclear. Here, we compared changes in SIC density (SICD) in Chinese croplands between the 1980s and 2010s. The SIC contents in 1980s were obtained from second national soil survey (n = 949) and published studies (n = 47). The SIC contents in 2010s were based on resampling of soil profiles from the same locations during 2019 and 2020 (n = 30), as well as data from published studies and national soil survey (n = 903). We found that Chinese croplands have lost 27-38% of SICD from the 0-40 cm soil layer and that the soil pH has decreased by 0.53 units over the past 30 years. These SIC losses increased with the ratio of precipitation (P) to potential evapotranspiration (PET) and most notably with nitrogen (N) fertilization. The SICD decreased greatly in humid and semiarid regions, and these losses were enhanced by high N fertilization rates; however, the SICD increased in very arid regions. This analysis demonstrates that the water balance and N fertilization are major drivers leading to dramatic losses of SICD in croplands and, consequently, to decreases in soil fertility and functions.

Significant loss of soil inorganic carbon at the continental scale.

Widespread soil acidification due to atmospheric acid deposition and agricultural fertilization may greatly accelerate soil carbonate dissolution and CO2 release. However, to date, few studies have addressed these processes. Here, we use meta-analysis and nationwide-survey datasets to investigate changes in soil inorganic carbon (SIC) stocks in China. We observe an overall decrease in SIC stocks in topsoil (0-30 cm) (11.33 g C m-2 yr-1) from the 1980s to the 2010s. Total SIC stocks have decreased by ∼8.99 ± 2.24% (1.37 ± 0.37 Pg C). The average SIC losses across China (0.046 Pg C yr-1) and in cropland (0.016 Pg C yr-1) account for ∼17.6%-24.0% of the terrestrial C sink and 57.1% of the soil organic carbon sink in cropland, respectively. Nitrogen deposition and climate change have profound influences on SIC cycling. We estimate that ∼19.12%-19.47% of SIC stocks will be further lost by 2100. The consumption of SIC may offset a large portion of global efforts aimed at ecosystem carbon sequestration, which emphasizes the importance of achieving a better understanding of the indirect coupling mechanisms of nitrogen and carbon cycling and of effective countermeasures to minimize SIC loss.

High stability and metabolic capacity of bacterial community promote the rapid reduction of easily decomposing carbon in soil.

Irreversible climate change alters the decomposition and sequestration of soil carbon (C). However, the stability of C components in soils with different initial organic matter contents and its relationship with the response of major decomposers to climate warming are still unclear. In this study, we translocated Mollisols with a gradient of organic matter (OM) contents (2%-9%) from in situ cold region to five warmer climatic regions to simulate climate change. Soil C in C-rich soils (OM >5%) was more vulnerable to translocation warming than that in C-poor soils (OM ≤ 5%), with a major loss of functional groups like O-alkyl, O-aryl C and carboxyl C. Variations of microbial ß diversity with latitude, temperature and precipitation indicated that C-rich soils contained more resistant bacterial communities and more sensitive fungal communities than C-poor soils, which led to strong C metabolism and high utilization ability of the community in C-rich soils in response to translocation warming. Our results suggest that the higher sensitivity of soils with high organic matter content to climate change is related to the stability and metabolic capacity of major bacterial decomposers, which is important for predicting soil-climate feedback.

Microbial metabolism and necromass mediated fertilization effect on soil organic carbon after long-term community incubation in different climates.

Understanding the effects of changing climate and long-term human activities on soil organic carbon (SOC) and the mediating roles of microorganisms is critical to maintain soil C stability in agricultural ecosystem. Here, we took samples from a long-term soil transplantation experiment, in which large transects of Mollisol soil in a cold temperate region were translocated to warm temperate and mid-subtropical regions to simulate different climate conditions, with a fertilization treatment on top. This study aimed to understand fertilization effect on SOC and the role of soil microorganisms featured after long-term community incubation in warm climates. After 12 years of soil transplantation, fertilization led to less reduction of SOC, in which aromatic C increased and the consumption of O-alkyl C and carbonyl C decreased. Soil live microbes were analyzed using propidium monoazide to remove DNAs from dead cells, and their network modulization explained 60.4% of variations in soil labile C. Single-cell Raman spectroscopy combined with D2O isotope labeling indicated a higher metabolic activity of live microbes to use easily degradable C after soil transplantation. Compared with non-fertilization, there was a significant decrease in soil α- and ß-glucosidase and delay on microbial growth with fertilization in warmer climate. Moreover, fertilization significantly increased microbial necromass as indicated by amino sugar content, and its contribution to soil resistant C reached 22.3%. This study evidentially highlights the substantial contribution of soil microbial metabolism and necromass to refractory C of SOC with addition of nutrients in the long-term.

Organism body size structures the soil microbial and nematode community assembly at a continental and global scale.

Body size is a key life-history trait that influences community assembly by affecting how ecological processes operate at the organism level. However, the extent to which the relative influences of ecological processes mediate the assembly of differentially sized soil organisms is still unknown. Here, we investigate the community assembly of differentially sized soil microorganisms and microfauna using a continental-scale sampling effort combined with a global-scale meta-analysis. Our results reveal a general relationship between organism body size and the stochastic-deterministic balance operating on community assembly. The smallest microorganisms (bacteria) are relatively more influenced by dispersal-based stochastic processes, while larger ones (fungi, protists and nematodes) are more structured by selection-based deterministic processes. This study elucidates a significant and consistent relationship between an organism life-history trait and how distinct ecological processes operate in mediating their respective community assemblages, thus providing a better understanding of the mechanisms supporting soil biodiversity.

[Interactions of Fungi Community and Relationship with the Carbon Structure in Arable Mollisols with Gradient Organic Matter Content].

Fungi play an important role in the accumulation and transformation of soil organic matter (SOM) and nutrient cycling. To investigate the relationship between the fungal community and soil organic carbon functional groups under gradient SOM contents in arable mollisols, arable mollisols with 2%-9% SOM content were collected in Northeast China. Solid-state 13C-NMR technology was used to explore the differences in the functional group structure of SOM, and ITS high-throughput sequencing was used to investigate the fungal community structure. The potential interactions between different taxonomic groups of soil fungal community and their associations with organic carbon molecular structures were compared by constructing molecular ecological networks under low SOM (2%-5%) and high SOM (7%-9%) conditions. The 13C-NMR results indicated an increase in the relative abundance of Alkyl C (25.8% to 35.9%). The decrease in Alkyl C/O-Alkyl C indicated a smaller degree of decomposition in high SOM soils. Sordariomycetes and Mortierellomycotina dominated the fungal community and their relative abundance increased with the SOM gradient (P<0.05) from 14.33% to 28.17% and from 7.32% to 23.14%, respectively. The network analysis showed simpler ecological topological properties of the fungal community in low SOM soils, with lower numbers of nodes, edges, and average clustering coefficients than those in high SOM soils. A closer relationship between fungi and organic carbon functional groups, especially LOC, was observed in low SOM soils. The random forest model showed that LOC had the largest amount for fungal interactions in low SOM soils (10%), followed by recalcitrant organic carbon (ROC). In comparison, LOC contributed less to the variations in fungal interactions in high SOM soils (7.4%). With globally increasing soil carbon loss, the limition of the carbon resources, especially the reduction of LOC, may reduce the stability and ecological functions of soil fungal communities.

The Diversity and Geographic Distribution of Cultivable Bacillus-Like Bacteria Across Black Soils of Northeast China.

Bacillus-like species are gram-positive bacteria that are ubiquitous in soils. Many of Bacillus-like bacteria are demonstrated as beneficial microbes widely used in industry and agriculture. However, the knowledge related to their diversity and distribution patterns in soils is still rudimentary. In this study, we developed a combined research method of using culture-dependent and high-throughput sequencing to investigate the composition and diversity of cultivable Bacillus-like bacterial communities across 26 soil samples obtained from the black soil zone in northeast China. Nearly all bacterial 16S rDNA sequences were classified into the order Bacillales. Fifteen genera were detected, with Bacillus, Paenibacillus, and Brevibacillus being the three most abundant genera. Although more than 2,000 OTUs were obtained across all samples, 33 OTUs were confirmed as the abundant species with a relative abundance over 5% in at least one sample. Pairwise analysis showed that the diversity of Bacillus-like bacterial communities were significantly and positively correlated with soil total carbon contents and soil sampling latitudes, which suggests that a latitudinal gradient diversity of Bacillus-like bacterial communities exists in the black soil zone. The principal coordinates analysis revealed that the Bacillus-like bacterial communities were remarkably affected by soil sampling latitudes and soil total carbon content. In general, this study demonstrated that a distinct biogeographic distribution pattern of cultivable Bacillus-like bacterial communities existed in the black soil zone, which emphasizes that the strategy of local isolation and application of beneficial Bacillus-like strains is rather important in black soil agriculture development.

Biogeographic Distribution Patterns of the Archaeal Communities Across the Black Soil Zone of Northeast China.

Although archaea are ubiquitous in various environments, the knowledge gaps still exist regarding the biogeographical distribution of archaeal communities at regional scales in agricultural soils compared with bacteria and fungi. To provide a broader biogeographical context of archaeal diversity, this study quantified the abundance and community composition of archaea across the black soil zone in northeast China using real-time PCR and high-throughput sequencing (HTS) methods. Archaeal abundances across all soil samples ranged from 4.04 × 107 to 26.18 × 107 16S rRNA gene copies per gram of dry soil. Several soil factors were positively correlated with the abundances including soil pH, concentrations of total C, N, and P, and available K in soil, and soil water content. Approximately 94.2, 5.7, and 0.3% of archaeal sequences, and 31, 151, and 3 OTUs aligned within the phyla Thaumarchaeota, Euryarchaeota, and Crenarchaeota, respectively. Within the phylum of Thaumarchaeota, group 1.1b was a dominating genus accounting for an average of 87% archaeal sequences and phylogenetically classified as Nitrososphaera, a genus of ammonia oxidizing archaea. The response of dominating OTUs to environmental factors differed greatly, suggesting the physiological characteristics of different archaeal members is diversified in the black soils. Although the number of OTUs was not related with any particular soil parameters, the number of OTUs within Thaumarchaeota and Euryarchaeota was marginally related with soil pH. Archaeal community compositions differed between samples, and a Canonical correspondence analysis (CCA) analysis indicated that soil pH and the latitude of sampling locations were two dominating factors in shifting community structures. A variance partitioning analysis (VPA) analysis showed that the selected soil parameters (32%) were the largest drivers of community variation, in particular soil pH (21%), followed by geographic distances (19%). These findings suggest that archaeal communities have distinct biogeographic distribution pattern in the black soil zone and soil pH was the key edaphic factor in structuring the community compositions.

Ammonia-Oxidizing Archaea Show More Distinct Biogeographic Distribution Patterns than Ammonia-Oxidizing Bacteria across the Black Soil Zone of Northeast China.

Black soils (Mollisols) of northeast China are highly productive and agriculturally important for food production. Ammonia-oxidizing microbes play an important role in N cycling in the black soils. However, the information related to the composition and distribution of ammonia-oxidizing microbes in the black soils has not yet been addressed. In this study, we used the amoA gene to quantify the abundance and community composition of ammonia-oxidizing archaea (AOA) and ammonia-oxidizing bacteria (AOB) across the black soil zone. The amoA abundance of AOA was remarkably larger than that of AOB, with ratios of AOA/AOB in the range from 3.1 to 91.0 across all soil samples. The abundance of AOA amoA was positively correlated with total soil C content (p < 0.001) but not with soil pH (p > 0.05). In contrast, the abundance of AOB amoA positively correlated with soil pH (p = 0.009) but not with total soil C. Alpha diversity of AOA did not correlate with any soil parameter, however, alpha diversity of AOB was affected by multiple soil factors, such as soil pH, total P, N, and C, available K content, and soil water content. Canonical correspondence analysis indicated that the AOA community was mainly affected by the sampling latitude, followed by soil pH, total P and C; while the AOB community was mainly determined by soil pH, as well as total P, C and N, water content, and sampling latitude, which highlighted that the AOA community was more geographically distributed in the black soil zone of northeast China than AOB community. In addition, the pairwise analyses showed that the potential nitrification rate (PNR) was not correlated with alpha diversity but weakly positively with the abundance of the AOA community (p = 0.048), whereas PNR significantly correlated positively with the richness (p = 0.003), diversity (p = 0.001) and abundance (p < 0.001) of the AOB community, which suggested that AOB community might make a greater contribution to nitrification than AOA community in the black soils when ammonium is readily available.

Responses of Labile Organic Nitrogen Fractions and Enzyme Activities in eroded Mollisols After 8-year Manure Amendment.

Soil erosion will cause a degradation in soil nitrogen supplying capacity (SNSC) and manure amendment is an effective way to restored eroded soils. Both labile fractions of soil organic N (SON) and N transformation enzymes are indicators for SNSC, but the effect of manure amendments on labile SON fractions and the relationship between labile SON fractions and enzyme activities remains unclear. In this study, five degrees of erosion were simulated in Mollisols (removal of 0, 5, 10, 20 and 30 cm of topsoil) to analyse the changes in labile SON fractions and nitrogen transformation enzyme activities after 8-year manure amendment. We found that soil total N (TN), labile SON fractions and enzyme activities all increased after manure amendments. The largest labile SON fraction was particle organic nitrogen (POM-N) and the second was light fraction organic nitrogen (LFOM-N), which accounted >60% for TN in total. Correlation analysis showed that both urease and protease activities were significantly correlated with POM-N, LFOM-N, microbial biomass N and dissolvable organic N, indicating that both urease and protease activities can be used to predict labile SON pools and enzyme activities worked similarly in indicating SNSC with labile SON fractions. Altogether, 8-year manure amendment could recover SNSC of lightly eroded Mollisols to natural levels, i.e. erosion depths at 5 cm and 10 cm; however, it is not able to recover SNSC in Mollisols suffering severe erosion.

Soil type recognition as improved by genetic algorithm-based variable selection using near infrared spectroscopy and partial least squares discriminant analysis.

Soil types have traditionally been determined by soil physical and chemical properties, diagnostic horizons and pedogenic processes based on a given classification system. This is a laborious and time consuming process. Near infrared (NIR) spectroscopy can comprehensively characterize soil properties, and may provide a viable alternative method for soil type recognition. Here, we presented a partial least squares discriminant analysis (PLSDA) method based on the NIR spectra for the accurate recognition of the types of 230 soil samples collected from farmland topsoils (0-10 cm), representing 5 different soil classes (Albic Luvisols, Haplic Luvisols, Chernozems, Eutric Cambisols and Phaeozems) in northeast China. We found that the PLSDA had an internal validation accuracy of 89% and external validation accuracy of 83% on average, while variable selection with the genetic algorithm (GA and GA-PLSDA) improved this to 92% and 93%. Our results indicate that the GA variable selection technique can significantly improve the accuracy rate of soil type recognition using NIR spectroscopy, suggesting that the proposed methodology is a promising alternative for recognizing soil types using NIR spectroscopy.

[Effect of temperature, rainfall and soil properties on farmland soil nitrification].

Climate conditions, soil properties and management practices control soil nitrification process which affects nitrogen cycling and balance in agro-ecosystems. The interaction of temperature, rainfall, soil type and fertilization on the soil nitrification process was studied by a soil transplantation experiment installed in 3 experiment stations of Chinese Ecological Research Network, i.e., Hailun, Fenqiu and Yingtan Agroecological Experiment Station, which represents middle temperature, warm temperature and middle subtropical zone, respectively. Three types of cropland soils were selected, i.e., neutral black soil (Phaeozem), alkaline Chao soil (Cambisol) and acidic red soil (Acrisol). Then one-meter depth soil profiles for each soil were transplanted in 3 stations to build the field experiment. The two-year experimental results (2006-2007) showed soil nitrification intensity (SNI) changed with the temperature and rainfall during the maize tasseling stage. From Hailun to Yingtan, with an increase of monthly average temperature from 22.3 degrees C to 26.8 degrees C and the monthly rainfall from 100.8 mm to 199.6 mm, SNI decreased by 64.2%-67.2% for black soil, 52.1%-52.5% for Chao soil, and 41.7%-75.2% for red soil, respectively. There were significant negative correlations between SNI and temperature and rainfall, with a correlation coefficient of r = -0.354 (p < 0.01) and r = -0.290 (p < 0.01), respectively. The total number of soil nitrobacteria and the intensity of soil nitrification was affected by soil types, which increased in a sequence of Chao soil > black soil > red soil. Among soil properties, pH affected SNI significantly, with a correlation coefficient of r = 0.551 (p < 0.01). In generally, climate condition (temperature and rainfall), soil type and fertilization present an integrated impact on soil nitrification process, and there were significant interactions of climate x soil type, climate x fertilization, soil type x fertilization, and climate x soil type x fertilization.