This study presents a description of the emission, transport, dispersion, and deposition of heavy metals contained in atmospheric aerosols emitted from a large industrial complex in southern Spain using the HYSPLIT model coupled with high- (MM5) and low-resolution (GDAS) meteorological simulations. The dispersion model was configured to simulate eight size fractions (<0.33, 0.66, 1.3, 2.5, 5, 14, 17, and >17 μm) of metals based on direct measurements taken at the industrial emission stacks. Twelve stacks in four plants were studied and the stacks showed considerable differences for both emission fluxes and size ranges of metals. We model the dispersion of six major metals; Cr, Co, Ni, La, Zn, and Mo, which represent 77% of the total mass of the 43 measured elements. The prediction shows that the modeled industrial emissions produce an enrichment of heavy metals by a factor of 2–5 for local receptor sites when compared to urban and rural background areas in Spain. The HYSPLIT predictions based on the meteorological fields from MM5 show reasonable consistence with the temporal evolution of concentrations of Cr, Co, and Ni observed at three sites downwind of the industrial area. The magnitude of concentrations of metals at two receptors was underestimated for both MM5 (by a factor of 2–3) and GDAS (by a factor of 4–5) meteorological runs. The model prediction shows that heavy metal pollution from industrial emissions in this area is dominated by the ultra-fine (<0.66 μm) and fine (<2.5 μm) size fractions. ► We modeled the dispersion of heavy metals originating from industrial complex. ► The model was refined by configuration of eight size fractions of metals. ► High- and low-resolution meteorological models were compared. ► Model driven by MM5 had reasonable prediction of metal pollution. ► Industrial emissions enhanced metal pollution levels by a factor of 2–5.