ABSTRACTIn this paper, the bifurcations of a discrete stage-structured population model with self-limitation between the two subgroups are investigated. We explore all possible codimension-one bifurcations associated with transcritical, flip (period doubling) and Neimark-Sacker bifurcations and discuss the stabilities of the fixed points in these non-hyperbolic cases. Meanwhile, we give the explicit approximate expression of the closed invariant curve which is caused by the Neimark-Sacker bifurcation. After that, through the theory of approximation by a flow, we explore the codimension two bifurcations associated with 1:3 strong resonance. We convert the nondegenerate condition of 1:3 resonance into a parametric polynomial, and determine its sign by the theory of complete discrimination system. We introduce new parameters and utilize some variable substitutions to obtain the bifurcation curves around 1:3 resonance, which are returned to the original variables and parameters to express for easy verification. By using a series of complicated approximate identity transformations and polar coordinate transformation, we explore 1:6 weak resonance. Moreover, we calculate the two boundaries of Arnold tongue which are caused by 1:6 weak resonance and defined as the resonance region. Numerical simulations and numerical bifurcation analyzes are made to demonstrate the effective of the theoretical analyzes and to present the relations between these bifurcations. Furthermore, our theoretical analyzes and numerical simulations are explained from the biological point of view.