The formation of carbonaceous deposits (coke) in zeolite pores during catalysis leads to temporary deactivation of catalyst, necessitating regeneration steps, affecting throughput, and resulting in partial permanent loss of catalytic efficiency. Yet, even to date, the coke molecule distribution is quite challenging to study with high spatial resolution from surface to bulk of the catalyst particles at a single particle level. To address this challenge we investigated the coke molecules in HZSM-5 catalyst after ethanol conversion treatment by a combination of C K-edge X-ray absorption spectroscopy (XAS), 13C Cross polarization-magic angle spinning nuclear magnetic resonance (CP-MAS NMR) spectroscopy, and atom probe tomography (APT). XAS and NMR highlighted the aromatic character of coke molecules. APT permitted the imaging of the spatial distribution of hydrocarbon molecules located within the pores of spent HZSM-5 catalyst from surface to bulk at a single particle level. 27Al NMR results and APT results indicated association of coke molecules with Al enriched regions within the spent HZSM-5 catalyst particles. The experimental results were additionally validated by a level-set-based APT field evaporation model. These results provide a new approach to investigate catalytic deactivation due to hydrocarbon coking or poisoning of zeolites at an unprecedented spatial resolution.