Difference equations are used in order to model the dynamics of populations with non-overlapping generations, since the growth of such populations occurs only at discrete points in time. In the simplest case such equations have form $N_{t+1}= F(N_t)$, where $N_t >0$ is the population size at a moment of time $t$, and $F$ is a smooth function. Among such equations the discrete Skellam model are most often used in practice. In the given paper the Skellam model of the form $N_{t+1}=a N_t/ (1+b N_t)$, $N_{t+1}=a N_t^2/ (b^2+ N_t^2)$, $N_{t+1}=a N_t/ (1+b N_t^2)$ is considered, where the parameters $a,b>0$ with taking an effect of harvesting. Positive equilibrium points and conditions for their stability for these equations were found. It is shown in analytical form that these equations do not have periodic solutions with period $T=2$, which means, according by the Sharkovskii theorem, periodic solutions of any periods. In the model with harvesting, only regimes with monotonic stabilization of the population size are observed. Therefore, in all models of Skellam with harvesting, the existence of a critical conception is show, beyond which the population will be completely destroyed. For practice it is important to know the permissible limits of harvesting intensity, which are found in this paper.