The formulation of a kinetic model for a complex reaction network typically yields reaction rates which vary over orders of magnitude. This results in time scale separation that makes the model inherently stiff. In this work, a graph-theoretic framework is developed for time scale decomposition of complex reaction networks to separate the slow and fast time scales, and to identify pseudo-species that evolve only in the slow time scale. The reaction network is represented using a directed bi-partite graph and cycles that correspond to closed walks are used to identify interactions between species participating in fast/equilibrated reactions. Subsequently, an algorithm which connects the cycles to form the pseudo species is utilized to eliminate the fast rate terms. These pseudo-species are used to formulate reduced, non-stiff kinetic models of the reaction system. Two reaction systems are considered to show the efficacy of this framework in the context of thermochemical and biochemical processing. (C) 2016 Elsevier Ltd. All rights reserved.