Modelling population dynamics in ecological systems reveals properties that are difficult to find by empirical means, such as the probability that a population will go extinct when it is exposed to harvesting. To study these properties, we use an aquatic ecological system containing one fish species and an underlying resource as our models. In particular, we study a class of stage-structured population systems with and without starvation. In these models, we study the resilience, the recovery potential, and the probability of extinction and show how these properties are affected by different harvesting rates, both in a deterministic and stochastic setting. In the stochastic setting, we develop methods for deriving estimates of these properties. We estimate the expected outcome of emergent population properties in our models, as well as measures of dispersion. In particular, two different approaches for estimating the probability of extinction are developed. We also construct a method to determine the recovery potential of a species that is introduced in a virgin environment.