Alpine, subalpine and boreal tree species, of low genetic diversity and adapted to low optimal temperatures, are vulnerable to the warming effects of global climate change. The accurate prediction of these species’ distributions in response to climate change is critical for effective planning and management. The goal of this research is to predict climate change effects on the distribution of red spruce (Picea rubens Sarg.) in the Great Smoky Mountains National Park (GSMNP), eastern USA. Climate change is, however, conflated with other environmental factors, making its assessment a complex systems problem in which indirect effects are significant in causality. Predictions were made by linking a tree growth simulation model, red spruce growth model (ARIM.SIM), to a GIS spatial model, red spruce habitat model (ARIM.HAB). ARIM.SIM quantifies direct and indirect interactions between red spruce and its growth factors, revealing the latter to be dominant. ARIM.HAB spatially distributes the ARIM.SIM simulations under the assumption that greater growth reflects higher probabilities of presence. ARIM.HAB predicts the future habitat suitability of red spruce based on growth predictions of ARIM.SIM under climate change and three air pollution scenarios: 10% increase, no change and 10% decrease. Results show that suitable habitats shrink most when air pollution increases. Higher temperatures cause losses of most low-elevation habitats. Increased precipitation and air pollution produce acid rain, which causes loss of both low- and high-elevation habitats. The general prediction is that climate change will cause contraction of red spruce habitats at both lower and higher elevations in GSMNP, and the effects will be exacerbated by increased air pollution. These predictions provide valuable information for understanding potential impacts of global climate change on the spatiotemporal distribution of red spruce habitats in GSMNP.