A multimedia fate model with spatially resolved air and soil phases was developed and evaluated. The model was used for calculation of phenanthrene concentrations in air, water, soil, and sediment in Tianjin area and transport fluxes between the adjacent bulk phases under steady-state assumption. Both air and soil phases were divided into 3113 individual compartments of 4 km2 each to assess the spatial variation of phenanthrene concentrations and fluxes. Independently measured phenanthrene concentrations in air, water, and soil were used for model validation. The spatial variation in soil was validated using a set of measured phenanthrene concentrations of 188 surface soil samples collected from the area. Most data used either for model calculation or for model validation were collected during the last 5 years. As the results of the model validation, the calculated mean values for phenanthrene concentrations in various bulk phases are in fair agreement with those independently observed and are very close to those calculated using the model without spatial variation. The absolute difference between the calculated and the measured mean concentrations are 0.14, 0.48, and 0.13 log-units (mol/m3) for air, water, and soil, respectively. The spatial distribution patterns of phenanthrene in both air and soil were well modeled. Spatially, however, the model overestimated the soil phenanthrene level at low concentration range and underestimated it at high concentration range. The calculated distribution of phenanthrene in the air matches well with the emission from fossil fuel combustion, while the calculated distribution pattern in the soil is similar to that observed.