The MM5-CMAx-PSAT modeling approach was presented to identify the variation of emission contribution from each modeling grid to regional and urban air quality per unit emission rate change. The method was applied to a case study in Tangshan Municipality, a typical industrial region in northern China. The variation of emission contribution to the monthly atmospheric SO2concentrations in Tangshan from each modeling grid of 9 Ć 9ākm per 1000āt/yr of emission rate change was simulated for four representative months in 2006. It was found that the northwestern part of Tangshan region had the maximum contribution variation ratio (i.e., greater than 0.36%) to regional air quality, while the lowest contribution variation ratio (i.e., less than 0.3%) occurred in the coastal areas. Principal component analysis (PCA), canonical correlation analysis (CCA), and Pearson correlation analysis indicated that there was an obvious negative correlation between the grid-based variation of emission contribution to regional air quality and planetary boundary layer height (PBLH) as well as wind speed, while terrain data presented insignificant impacts on emission contribution variation. The proposed method was also applied to analyze the variation of emission contribution to the urban air quality of Tangshan (i.e., a smaller scale).