In this study, first we expand the cubic population model under Allee effect with a quadratic harvest function that represents harvest effect. Then, using four reaction equations representing the micro-interactions within the population under the influence of demographic noise and harvest, we obtain the mean field equations containing the effect of the harvest from the solution of the master equation. In this way, a relationship is established between the micro and macro parameters of the population. As a result, we calculate Mean Time to Extinction (MTE) by using WKB approximation for single step processes and observe the effect of harvesting.

A detailed examination of the effect of harvesting on a population has been carried out by extending the standard cubic deterministic model by considering a population under Allee effect with a quadratic function representing harvesting. Weak and strong Allee effect transitions, carrying capacity, and Allee threshold change according to harvesting are first discussed in the deterministic model. A Fokker–Planck equation has been obtained starting from a Langevin equation subject to correlated Gaussian white noise with zero mean, and an Approximate Fokker–Planck Equation has been obtained from a Langevin equation subject to correlated Gaussian colored noise with zero mean. This allowed to calculate the stationary probability distributions of populations, and thus to discuss the effects of linear and nonlinear (Holling type-II) harvesting for populations under Allee effect and subject to white and colored noises, respectively.