This paper proposes a system of difference equations to model mosquito population. In this study, we develop and analyze the stage-structured models which consist of four distinct mosquito metamorphic stages: eggs, larvae, pupae, and adults. First, a model with constant birth rate is studied, and the inherent net reproduction number 0 of the model is derived. If 0 < 1, the extinction equilibrium is globally asymptotically stable. If 0 > 1, there exists a unique positive equilibrium which is uniformly persistent. When breeding is seasonal for a special case, it indicates that a unique globally asymptotically stable periodic solution is admitted when the net reproductive number is larger than one. When this value is less than one, the mosquito population goes to extinction. Finally, numerical simulations to demonstrate our findings and brief discussion are also provided.At last, we proceed to verifying the last inequality (v). For convenience, we denote m = -bs 3 J 21 J 32 , n = -J 44 . Then (v) is equal to 1b 2 s 2 3 J 2 21 J 2 32 2b 2 s 2 3 J 2 21 J 2 32 J 2 44 --bs 3 J 21 J 32 J 2 44