Aromaticity is an established and widely used concept for the prediction of the reactivity of organic molecules. However, its role remains largely unexplored in on-surface chemistry, where the interaction with the substrate can alter the electronic and geometric structure of the adsorbates. Here we investigate how aromaticity affects the reactivity of alkyne-substituted porphyrin molecules in cyclization and coupling reactions on a Au(111) surface. We examine and quantify the regioselectivity in the reactions by scanning tunnelling microscopy and bond-resolved atomic force microscopy at the single-molecule level. Our experiments show a substantially lower reactivity of carbon atoms that are stabilized by the aromatic diaza[18]annulene pathway of free-base porphyrins. The results are corroborated by density functional theory calculations, which show a direct correlation between aromaticity and thermodynamic stability of the reaction products. These insights are helpful to understand, and in turn design, reactions with aromatic species in on-surface chemistry and heterogeneous catalysis. While aromaticity is a useful concept for assessing the reactivity of organic compounds, the connection between aromaticity and on-surface chemistry remains largely unexplored. Now, scanning probe experiments on cyclization reactions of porphyrins on Au(111) show that the peripheral carbon atoms outside of the aromatic 18-& pi; electron pathway exhibit a higher reactivity.