Abstract
Global patterns in soil microbiomes are driven by non‐linear environmental thresholds. Fertilization is known to shape the soil microbiome of terrestrial ecosystems worldwide. Yet, whether fertilization influences global thresholds in soil microbiomes remains virtually unknown. Here, utilizing optimized machine learning models with Shapley additive explanations on a dataset of 10,907 soil samples from 24 countries, we discovered that the microbial community response to fertilization is highly dependent on environmental contexts. Furthermore, the interactions among nitrogen (N) addition, pH, and mean annual temperature contribute to non‐linear patterns in soil bacterial diversity. Specifically, we observed positive responses within a soil pH range of 5.2–6.6, with the influence of higher temperature (>15°C) on bacterial diversity being positive within this pH range but reversed in more acidic or alkaline soils. Additionally, we revealed the threshold effect of soil organic carbon and total nitrogen, demonstrating how temperature and N addition amount interacted with microbial communities within specific edaphic concentration ranges. Our findings underscore how complex environmental interactions control soil bacterial diversity under fertilization.