We construct a topological invariant that classifies density matrices of symmetry-protected topological orders in two-dimensional fermionic systems. As it is constructed out of the previously introduced Uhlmann phase, we refer to it as the topological Uhlmann number nU. With it, we study thermal topological phases in several two-dimensional models of topological insulators and superconductors, computing phase diagrams where the temperature T is on an equal footing with the coupling constants in the Hamiltonian. Moreover, we find novel thermal-topological transitions between two nontrivial phases in a model with high Chern numbers. At small temperatures we recover the standard topological phases as the Uhlmann number approaches to the Chern number. Although these studies provide a successful picture for quantum systems in pure states, typically the ground state, very little is known about the fate of those topological phases of matter when the system is in a mixed quantum state represented by a density matrix. In fact, the correct understanding of this situation becomes particularly relevant in order to address unavoidable thermal effects on topological phases and nonequilibrium dynamics under dissipation [19].Over the last few years, there have been some results establishing the absence of stable topological phases subject to thermal effects, like for TOs with spins [20][21][22], for SP-TOs with spins [23], or SP-TOs with fermions under certain conditions [24]. However, in a recent work [25] we have shown that it is possible to characterize a thermal topological phase for topological insulators and superconductors in one-dimensional systems. This thermal topological phase, classified by Uhlmann holonomies [26,27], is separated from a trivial phase by a critical finite temperature, at which the topological phase abruptly disappears. This result paves the way towards the characterization of SP-TOs with fermions in thermal states, or more general density matrices.In this paper, we have achieved this goal for twodimensional fermion systems, either insulating or superconducting. This extension from one to two-dimensional systems is nontrivial in the sense that a direct generalization of the topological invariants for pure states (Chern numbers) to density matrices via the Uhlmann approach leads to trivial results. However, we have succeeded in circumventing this problem and introducing a suitable notion of topological Uhlmann numbers. These are gauge invariant and observable quantities which allows for a classification of topological phases of density matrices in two-dimensional (2D) quantum systems. Specifically, we