Thermodynamic properties of chiral spin liquids are investigated for a variant of the Kitaev model defined on a decorated honeycomb lattice. Using the quantum Monte Carlo simulation, we find that the model exhibits a finite-temperature phase transition associated with the time reversal symmetry breaking, in both topologically trivial and nontrivial regions. While changing the exchange constants, the phase transition changes from continuous to discontinuous one, apparently correlated with the change in the excitations from Abelian to non-Abelian anyons. We show this coincidence by computing the topological quantities: the Chern number and the thermal Hall conductivity. In addition, we find, as a diagnostic of the chiral spin liquids, successive crossovers with multi-stage entropy release above the critical temperature, which indicates that the hierarchical fractionalization of a quantum spin occurs differently between the two regions.PACS numbers: 75.10. Kt,75.10.Jm,75.30.Et Understanding of quantum spin liquids (QSLs) in magnets, where strong quantum fluctuations suppress magnetic ordering even at the lowest temperature (T ), has been one of the most challenging subjects in strongly correlated electron systems [1,2]. Among many possible realizations of QSLs, the chiral spin liquid (CSL), in which the time reversal symmetry is broken, has attracted considerable attention in not only condensed matter physics but also quantum information. This is because it may have the excitations obeying the nonAbelian anyon statistics, which are utilized as the operators in topological quantum computing [3]. To explore this exotic state, quantum spin systems on geometrically frustrated lattices have been intensively studied thus far [4]. For instance, the possibility of CSLs has been studied theoretically in the Heisenberg model on a kagome lattice [5][6][7]. Experimentally, a possible CSL was discussed for a metallic pyrochlore compound Pr 2 Ir 2 O 7 [8].Besides the analyses of geometrically-frustrated quantum magnets, a class of the models that have the exact CSL ground states has opened a new avenue in the study of CSLs [3,9,10]. One of them was originally suggested by A. Kitaev [3] and studied in detail by H. Yao and S. Kivelson [9]. This model is a variant of the honeycomb Kitaev model, which is defined on a decorated honeycomb lattice, composed by extending each honeycomb lattice site to a triangle. The exact solution of this model shows that the ground state accommodates two different types of the CSLs accompanied with Abelian and nonAbelian anyons as the elementary excitations. Interestingly, the non-Abelian CSL has topologically nontrivial Majorana fermion bands with a chiral edge mode.Although the exact solutions for the ground states serve as good references in the exploration of CSLs, it is a crucial issue how the CSLs behave against thermal fluctuations at finite T . As the discrete chiral symmetry is broken in the ground state, one expects a phase transition to the CSL at a finite T , even in two dimensions. An ...