Global dynamics in a stage-structured discrete-time population model with harvesting

► We study the stability of equilibria in the system using analytical methods. ► We apply an automated rigorous numerical method to classify global dynamics. ► We obtain a comprehensive overview of the dynamics of the system. ► We analyze the hydra effect-when harvesting may increase population abundance. ► We show a bubble effect-when population is destabilized and stabilized again.