Hybrid zone theory provides a powerful theoretical framework for measuring and testing gene flow and selection. The Senecio aethnensis and Senecio chrysanthemifolius hybrid zone on Mount Etna, Sicily, was investigated to identify phenotypic traits under divergent selection and to assess the contributions of intrinsic and extrinsic selection against hybrids to hybrid zone maintenance. Senecio samples from 14 sites across Mount Etna were analyzed for 24 quantitative traits classified into four groups (QTGs), six allozymes and seven simple sequence repeat (SSR) loci to describe patterns of variation throughout the hybrid zone. Narrower cline widths or shifts in cline centre position were observed for three QTGs relative to the molecular clines, indicating that these traits are likely to be under extrinsic environmental selection. Altitude was key to describing species distributions, but dispersal and intrinsic selection against hybrids explained patterns at smaller spatial scales. The hybrid zone was characterized by strong selection against hybrids, high dispersal rates, recent species contact and few loci differentiating QTGs based on indirect measures. These results support the hypothesis that extrinsic and intrinsic selection against hybrids maintains the hybrid zone and species distinctiveness despite gene flow between the two Senecio species on Mount Etna.