Seabirds breeding on islands are vulnerable to introduced predators, such as rats and cats, and the removal of such predators is generally viewed as a priority for seabird conservation and restoration. However, multiple invasive mammal species interacting may generate unexpected outcomes following the removal (eradication) of one species. Generally these indirect interactions are not well understood or demonstrated. We propose and study a prey (seabird)- mesopredator (rat)-superpredator (cat) model, taking into account the juvenile stages in the prey population, in order to direct conservation management for seabird conservation. We give a more biologically realistic differential system than those studied before (Courchamp et al. [1999]; Fan et al. [2005]), in particular for long-lived seabird species. We present a theoretical study and show existence and uniqueness of a positive solution as well as a qualitative study of the equilibria that may appear. Because standard numerical methods, usually implemented in scientific softwares, can fail to give the right biological approximations (Anguelov et al. [2009]), we propose a reliable algorithm that preserves most of the qualitative properties of the continuous system, using the theory of nonstandard finite difference methods. Finally, we use biologically realistic parameters available for the representative Barau's petrel (Pinet et al. [2008]), an endemic species from Réunion island, to present numerical simulations that support the theoretical study. Cats play the major role in seabird prey population dynamics. Seasonality in seabird breeding delays but does not prevent extinction. In all scenarios, cat control (or preferably eradication) is imperative to prevent extinction of vulnerable long-lived seabirds, like the Barau's petrel.