Global dynamic scenarios in a discrete-time model of renewable resource exploitation: a mathematical study

We consider the two-dimensional map introduced in Bischi et al. (J Differ Equ Appl 21(10):954–973, 2015) formulated as a model for a renewable resource exploitation process in an evolutionary setting. The global dynamic scenarios displayed by the model are not so often encountered in smooth two-dimensional dynamical systems. We explain the occurrence of such scenarios at the light of the theory of noninvertible maps. Moreover, complex structures of basins of attraction of coexisting invariant sets are observed. We analyze such structures by examining stability properties of chaotic sets, in the case in which a non-topological Milnor attractor is present. Stability changes of a chaotic set occur through global bifurcations (such as riddling and blowout) and are detected by means of the study of the spectrum of Lyapunov exponents associated with the set.