Critical loads (CLs) are frequently used to quantify terrestrial ecosystem impacts from nitrogen (N) deposition using ecological responses such as the growth and mortality of tree species. Typically, CLs are reported as a single value, with uncertainty, for an indicator across a species' entire range. Mediating factors such as climate and soil conditions can influence species' sensitivity to N, but the magnitudes of these effects are rarely calculated explicitly. Here, we quantify the spatial variability and estimation error in N CLs for the growth and survival of 10 different tree species while accounting for key environmental factors that mediate species sensitivity to N (e.g., soil characteristics). We used a bootstrapped machine learning approach to determine the level of N deposition at which a 1% decrease occurs in growth rate or survival probability at forest plot locations across the United States. We found minimal differences (<5 kg N ha−1 year−1) when comparing a single species' CLs across climatic regimes but found considerable variability in species' local N CLs (>8.5 kg N ha−1 year−1) within these regimes. We also evaluated the most important factors for predicting tree growth rates and mortality and found that climate, competition, and air pollution generally have the greatest influence on growth rates and survival probability. Lastly, we developed a new probability of exceedance metric for each species and found high likelihoods of exceedance across large portions (46%) of some species' ranges. Our analysis demonstrates that machine learning approaches provide a unique capability to: (1) quantify mediating factor influences on N sensitivity of trees, (2) estimate the error in local N CL estimates, and (3) generate localized N CLs with probabilities of exceedance for tree species.