Human activity is changing air quality and climate in the US Pacific Northwest. In a first application of non-metric multidimensional scaling to a large-scale, framework dataset, we modeled lichen community response to air quality and climate gradients at 1416 forested 0.4 ha plots. Model development balanced polluted plots across elevation, forest type and precipitation ranges to isolate pollution response. Air and climate scores were fitted for remaining plots, classed by lichen bioeffects, and mapped. Projected 2040 temperatures would create climate zones with no current analogue. Worst air scores occurred in urban-industrial and agricultural valleys and represented 24% of the landscape. They were correlated with: absence of sensitive lichens, enhancement of nitrophilous lichens, mean wet deposition of ammonium >0.06 mg l(-1), lichen nitrogen and sulfur concentrations >0.6% and 0.07%, and SO(2) levels harmful to sensitive lichens. The model can detect changes in air quality and climate by scoring re-measurements.