We develop a systematic theory of symmetry fractionalization for fermionic topological phases of matter in (2+1)D with a general fermionic symmetry group G f . In general G f is a central extension of the bosonic symmetry group G b by fermion parity, p´1q F , characterized by a nontrivial cohomology class rω2s P H 2 pG b , Z2q. We show how the presence of local fermions places a number of constraints on the algebraic data that defines the action of the symmetry on the supermodular tensor category that characterizes the anyon content. We find two separate obstructions to defining symmetry fractionalization, which we refer to as the bosonic and fermionic symmetry localization obstructions. The former is valued in H 3 pG b , KpCqq, while the latter is valued in either H 3 pG b , A{t1, ψuq or Z 2 pG b , Z2q depending on additional details of the theory. KpCq is the Abelian group of functions from anyons to Up1q phases obeying the fusion rules, A is the Abelian group defined by fusion of Abelian anyons, and ψ is the fermion. When these obstructions vanish, we show that distinct symmetry fractionalization patterns form a torsor over H 2 pG b , A{t1, ψuq. We study a number of examples in detail; in particular we provide a characterization of fermionic Kramers 1. Z 2 pG b , Z 2 q fermionic obstruction, Υ ψ locality-violating, and SOp3q 3 26 2. H 3 pG b , A{t1, ψuq fermionic obstruction, Υ ψ locality-preserving, and Spp2q 2 b t1, ψu 27 VII. Discussion 28 VIII. Acknowledgements 29 A. Characters of the fusion algebra of super-modular tensor categories 29 B. Comments on gauge transformations of U 30 C. Proof of when Υ ψ respects locality 32 D. Proof of Eq. 132 33 References 34 F .