Accurate discrimination of metal species is a significant analytical challenge. Herein, we propose a novel methodology based on liquid-phase cyclic chemiluminescence (CCL) for the identification of cobalt speciation. The CCL multistage signals ( I n ) of the luminol-H 2 O 2 reaction catalyzed by different cobalt species have different decay coefficients k . Thereby, we can facilely identify various cobalt species according to the distinguishable k values, including the complicated and structurally similar cobalt complexes, such as analogues of [Co(NH 3 ) 5 X] n + (X = Cl - , H 2 O, and NH 3 ), Co(II) porphyrins, and bis(2,4-pentanedione) cobalt(II) derivatives. Especially, the number of substituent atoms also influences the k value greatly, which allows excellent discrimination between complexes that only have a subtle difference in the substituent group. In addition, linear discriminant analysis based on I n provides a complementary solution to improve the differentiating ability. We performed density functional theory calculations to investigate the interaction mode of H 2 O 2 over cobalt species. A close negative correlation between the adsorption energy and the k value is observed. Moreover, the calculation of energy evolutions of H 2 O 2 decomposition into a double hydroxide radical shows that a high level of consistency exists between the activation energy barrier and the k value. The results further demonstrate that the decay coefficient of the CCL multistage signal is associated with the catalytic reactivity of the cobalt species. Our work not only broadens the application of chemiluminescence but also provides a complementary technology for speciation analysis.