AbstractCannabis is an ancient plant that has been used for therapeutic and recreational purposes. Nowadays, industrial hemp, a variety with low concentration of the psychoactive cannabinoid Δ⁹-tetrahydrocannabinol (THC) and high concentration of non-psychoactive cannabinoids, is getting more and more interest in the food, pharmaceutical, and cosmetic industry. However, cannabis not only contains cannabinoids as bioactive components but also other metabolites like terpenes and phenolic compounds, and the content of these interesting secondary metabolites greatly differs with the genetic variety of the plant. Due to the huge complexity of composition of the cannabis matrix, in this work, a comprehensive two-dimensional liquid chromatography (LC × LC) method has been developed as a very power separation technique coupling a pentafluorophenyl (PFP) and a C18 in the first and second dimensions. Two industrial hemp strains (cookie and gelato) were analyzed to determine the difference in their content of cannabinoids and phenolic compounds. To do this, a new demodulation process was applied for the first time to transform 2D raw data into 1D data which allowed carrying out the chemometric analysis needed to determine the statistical differences between the hemp strains. The cookie strain presented a total of 41 cannabinoid markers, while the gelato strain presented more representative phenolic compounds, in total 24 phenolic compounds were detected as potential markers of this sample. These differences in the chemical composition could determine the industrial destiny of the different hemp strains. Graphical abstract