Z-Type Control Methods on a Three-Species Model with an Invasive Prey

: The process by which nature and evolution discriminate between species bound to survive or become extinct has always been an interesting phenomenon. One example of this selection is represented by biological invasions, when alien species spread into new environments causing ecological disruptions, alteration of native population dynamics and ecosystem dysfunctioning. In this paper, we concentrate on the introduction of the Eastern cottontail ( Sylvilagus ﬂoridanus ) into the regions of northern and central Italy. This has inﬂuenced the local predator–prey dynamics of the indigenous European hares ( Lepus europaeus ) and red foxes ( Vulpes vulpes ), consequently producing a “hyperpredation” effect on hares. In this framework, we investigate this scenario with the objective of mitigating the negative effects of the invasion. Speciﬁcally, we devise and apply a Z-control technique in order to drive chosen populations toward a desired state that at least theoretically can ensure the survival of native species. The purpose of such a study is twofold: on one hand, this approach is employed with the aim of reducing the invasive population; secondly, it aims to avoid the extinction of native prey, i


Introduction
Currently, with global communications a phenomenon unknown in centuries before the XXth is occurring, also favored by the ongoing climatic changes.It is represented by the invasions of alien species that find new habitats in regions where they were formerly unknown [1].The negative effects of alien species on indigenous ones can be exerted in various ways [2].For example, demographically, by competition, predation or hybridization, and epidemiologically, via new disease transmission.In particular, during the invasion the alien species often outcompete the native ones, for various reasons, among which we can mention the larger niche occupied by the invaders, their bigger size, absence of specific predators and so on.The alien species may also become a nuisance for agriculture, in addition to the indigenous species [3][4][5].
Modeling pest populations is a fundamental task for biological control, which tries to fight these pernicious agents via their natural enemies, rather than by using chemicals.In the last decades, the latter have indeed often been found to produce counterproductive effects [6].Focusing on insects, spraying also kills useful populations, such as spiders, generalist predators that remove a very sizeable quantity of insects [7][8][9].A second side effect is that upon ingestion by other animals that are not their specific target, these poisons enter into the food web, accumulating in the largest animals.Indeed, the latter have longer lifespans and therefore their intake lasts longer and may ultimately reach dangerous levels in their bodies.Through the passage among various trophic levels, these poisons end up on our tables and finally in our bodies.This advocates the use of alternative methods, such as biological control [10][11][12][13][14].
The output of population model studies is the assessment of the final regime of the dynamical system, e.g., stable equilibria configurations or persistent oscillations.The theoretical analysis allows one to assess under which conditions these configurations occur, in terms of the model parameters.Often, it is the case that these various situations are linked to each other, in the sense that a change of a parameter shifts the ecosystem to a different outcome.The latter can be induced by natural stochastic environmental fluctuations or by temperature or other physical quantity changes brought up by different climatic conditions.It is sometimes possible to obtain a chart in which all the final system configurations appear, with links among "neighboring" equilibria, i.e., those attainable by a parameter change.Such a graph would represent a useful tool for ecosystem managers to drive ecosystems toward a desired state [15].
In this paper, we turn our attention to the introduction of an invasive species into a predator-prey system.This in general can lead to direct or indirect competition between native and exotic species, which in some cases can even lead to the extinction of one or more native species [16].In this regard, let us keep in mind that invasive species are a major cause of animal extinction [17].Several ecological implications of invasions have been analyzed in the literature, for example, in the case of marine invaders in [18] and of terrestrial invaders in [19].
In particular, we consider a specific invasion situation concerning the Eastern cottontail (Sylvilagus floridanus).This lagomorph has shifted the natural equilibrium between red foxes (Vulpes vulpes) and native European hares (Lepus europaeus) in the northern and central regions of Italy [20][21][22].No direct competition seems to occur between the hares and the cottontails.However, indirect competition dynamics exist between the two lagomorphs.In fact, the correlation between fox and hare abundances, which is positive when invasive rodents are few, becomes more and more negative as the latter's population increases [23].A theoretical investigation has allowed an explanation of the phenomenon [24].However, the problem persists.An increase in cottontail abundance would lead to an increase in fox populations, amplifying their predatory impact on hares ("hyperpredation").Alternatively, cottontails attract foxes where they live and where there are also important resting areas for hares, and consequently the increased presence of foxes results in increased predation rates on hares.In order to preserve the native hare, it is necessary to control the invader population.For this task, we apply a relatively novel technique, the Z control [25][26][27][28].
The article is organized as follows.In the next section, the Z-control technique is recalled, the reference model on which to apply this technique is presented and finally, three different indirect control strategies are proposed, respectively, in Sections 2.3-2.5:action on foxes, action on hares and combined action on foxes and hares.Each one of the new models is studied analytically in detail.In particular, the possible equilibrium points are determined with their feasibility and local stability conditions.In the case of action on foxes, a transcritical bifurcation linking the hare-free equilibrium with coexistence has also been identified, analytically demonstrated and graphically represented.In Section 3, several numerical simulations obtained using MATLAB are shown.Further, the biological interpretation of the three proposed controls is discussed.The Conclusion finally debates the best option to manage the ecological control problem.

Highlights on the Z Control
Starting from an autonomous system ẏ = f (y, G) incorporating an unknown control function G = G(t) as input, the central idea of the Z control method is to consider the evolution over time of the tracking error function v 1 (t) = y(t) − y d (t), which measures the discrepancy between the output solution y(t) and the desired outcome y d (t).More in detail, this method is based on forcing the first error v 1 (t) to converge zero exponentially in the following prescribed way: If it is not possible to find an explicit expression for the control function G(t) from Equation (1), the same process is repeated for the subsequent errors v k (t) , recursively defined by until this goal is achieved.This means that, at each step, a new condition is imposed: The control function can involve one or more equations of the original model.We distinguish between indirect and direct controls.In the former case, the equations of the classes to be controlled are not modified and the action is performed on other classes; in the second case, a control function is introduced directly on the classes for which a desired state has been chosen.
In our case, this method will be applied in three different ways, namely, acting on the foxes alone with the aim of reducing the invasive population and possibly to avoid the extinction of native hares; alternatively just acting on the hares to preserve them and possibly eradicate the cottontails; finally, by combining these efforts.

The Reference Model
The starting point of this analysis is a three-species mathematical model (sometimes referred to simply as reference model) which describes the predator-two-prey interaction between the two native species, foxes and hares, and the invasive cottontail species.As already mentioned, the model formulated in [24] has explained the field results obtained in [23].
Formally, the dynamics are described by the following system of first-order nonlinear ordinary differential equations, where the variables V, S and L, respectively, represent the populations of foxes, cottontails and hares: All parameters are nonnegative, with their values being taken from [24], unless otherwise specified.In particular, information from [20][21][22] on the hares' carrying capacity allows the determination of the intraspecific competition rate for this species.These equations essentially model reproduction and intraspecific competition for all species, natural mortality and the hunting of foxes over the two types of prey.These interactions are described via mass action, or Holling type I, terms in view of the low population densities.It is important to remark that direct interspecific competition between prey does not occur, as clearly discussed in [24].
A preliminary analysis of the model leads to eight equilibria E i = (V i , S i , L i ) for i = 0, . . ., 7, whose feasibility conditions and stability properties can be found in [24]; for our discussion, we will just consider (where the names are referred to the ones in [24]) in order to facilitate a comparison with the models in which Z-control is applied.In addition to explaining the "hyperpredation effect" empirically observed, the results of the analysis of [24] show that coexistence of the three species is possible under suitable conditions, but also that the extinction of the hares can occur.In the former case, this means persistence of the invader.Our task here is to try to curb this phenomenon and possibly to enhance the hares' survival.This is performed in three different ways: acting separately on foxes, at first, then on the hare population and finally by a combination of the two.

Invasive Species Removal: Single Indirect Control on Foxes
Here, we apply an indirect control on the cottontail population acting on foxes.This is performed by introducing a control function G(t) in the equation of the foxes.The motivation behind this choice lies in the fact that to curb the invasive cottontails causing the hyperpredation effect on hares, a predator removal is considered, thereby reducing the pressure on hares.The class to be forced is represented by the cottontails S, for which the desired state is here chosen as a nonnegative constant function S d (t) ≡ S d ∀t ≥ 0, with S d ≈ 0.

Z-Controlled Model Design on Foxes
The first modified model is obtained acting on the first equation of the reference model ( 2) by introducing a forcing function G pred (t), while the other two equations remain unchanged: From now on, we will refer to it as G(t) for simplicity.In general, the input function could assume both positive and negative values, but what is considered relevant in this context is G(t) > 0. This function could possibly represent a culling or removal of predators followed by their transfer elsewhere.
Introducing now a design parameter λ > 0 , we define the first error function We apply the Z-control method by imposing v 1 (t) to approach 0 exponentially fast for t → +∞.This corresponds to setting In this case, we will also need a further condition involving the second error function, defined as so that an explicit expression of the control function G(t) can be found in the following way.Imposing once more and substituting in it (5), we are led to and finally to It is now sufficient to rewrite (6) in terms of the original variables, omitting for convenience the dependence on t, to obtain and rearranging By letting and observing that we thus obtain We now solve in terms of G Equation (8) to get the control function Note that in so doing we must assume S, V = 0.This is the reason why we cannot impose that S d = 0. From (7), we finally observe that for S d (t) Note that the equation describing the dynamics of foxes in the model controlled in this way explicitly depends not only on V but also on S.
Theorem 1.For any positive initial state (V 0 , S 0 , L 0 ) and for any continuously differentiable and bounded desired state S d (t), the tracking error function v 1 (t) converges to zero exponentially.
Proof.We observe that setting u 1 = v 1 and u 2 = v1 , ( 6) is equivalent to the following system of first-order equations This system has the general solution where, using the initial conditions Applying Lemma 1 of [29], it is possible to control with an exponential function the term c 2 te −λt appearing in the expression of v 1 (t).Choosing C, Λ > 0 appropriately, we can then bound v 1 (t) from above From this, the exponential convergence of the error to zero with rate λ is proved.

Equilibria
We now look for the equilibrium points of the controlled system and study their feasibility conditions, excluding the cases in which either V or S vanish, as stated above.The only population that is theoretically allowed to disappear is thus L. The equilibrium equations lead only to the two possible equilibria Observe that S = S d when G d (t) ≡ λ 2 S d , which is consistent with the control requesting S(t) to achieve a stationary state.Furthermore, the value of V at both equilibria is the same.It depends only on the choice of S d .Hence, the only difference between these points is determined by the survival or extinction of the hare population.As far as these two equilibria are concerned, the first one does not reflect an optimal ecological situation, as the hares disappear while foxes, being generalist predators, would be supported not only by cottontails but also by the alternative feeding resources modeled via the carrying capacity in System (2).Instead, the second one is preferable because "near extinction", if S d ≈ 0, only affects the invasive species.There is the possibility then that the pristine invader-free ecosystem could be restored.
We rewrite here for convenience the actual populations at equilibrium with the specific choice of S d (t) Finally, the equilibria have biological sense if the populations are nonnegative.Thus, the following feasibility conditions for E L must hold: while those for E c are If not otherwise specified, from now on we will assume they are verified.As already remarked earlier, an important point is that S d = 0.

Local Stability Analysis
For a local stability analysis, we consider the matrix of the linearized system in a neighborhood of each equilibrium point, i.e. the Jacobian matrix evaluated at each equilibrium.Some derivatives are indeed immediate, since F 1 , f 2 and f 3 do not depend on L and S, respectively.
The Jacobian matrix of the Z-controlled model is so that its characteristic polynomial factorizes to immediately provide an explicit real eigenvalue µ 1 = d f 3 dL , while the remaining ones come from the roots of the quadratic In addition to the negativity of µ 1 , the local stability for the equilibria is obtained by applying the Routh-Hurwitz (RH) criterion.Proposition 1.For the Z-controlled Models (2) and (3):

•
If feasible, the hare-free equilibrium E L is locally asymptotically stable if and only if • The coexistence equilibrium E c , when feasible, is instead locally asymptotically stable for every choice of parameters, including λ.
Proof.The asymptotic local stability for the equilibria is therefore equivalent to requiring the RH conditions on the principal minor of (12), tr( J) < 0 and det( J) > 0, which are the same for both equilibria and unconditionally satisfied, as for all λ > 0 Hence, the stability is determined by the sign of the first eigenvalue, for which we distinguish the two cases.For E L µ (L) and for E c which holds unconditionally.
In conclusion, stability for E c is unconditional when it is feasible, namely, (11) holds, while for E L in addition to feasibility (10) it is also necessary to require (13).
From the two previous properties, the bifurcation scenarios may be discussed, taking as bifurcation parameter, for instance, b.Choosing for b a value much larger than the one of a could represent a hyperpredation situation.Such a choice may produce stability for E L , as the hare population is subject to a higher pressure.
Assume that all ecological parameters are fixed so that the feasibility conditions of E L are satisfied.Studying the sign change of a(u Note that some values of b do not guarantee the feasibility of E c .Specifically, (11) are not satisfied if b ≥ b * .

Existence of the Transcritical Bifurcation
We now actually prove the existence of the transcritical bifurcation for which E L emanates from E c as soon as b increases past b * , using the classical tool of Sotomayor's theorem [30].
Evaluating in E L the Jacobian matrix J of the Z-controlled model with S d (t) ≡ S d = const., we obtain Focusing attention on the parameter b, we can observe that the third element on the diagonal of J(E L ) vanishes when b = b * , and in view of the structure of J(E L ), it is in fact an eigenvalue.Right and left eigenvectors of Furthermore, denoting by F the system's right-hand side, we find Consequently, we have Thus, Sotomayor's theorem allows us to conclude that if s = n, there is a transcritical bifurcation from This transcritical bifurcation is graphically represented in Figure 1.The initial conditions are while the ecological parameter values are those of [24], except b, which are

Action on Hares
In view of the hyperpredation effect on hares exerted by foxes in the presence of the invading cottontails [23,24], it is imperative to find ways of curbing this negative effect.
For this second scenario, the goal remains the same, but in this case instead of culling the predators, the importation of hares from outside is assumed.This sustains the hare population, although it also provides more resources for foxes.However, one should keep in mind that the introduction of alien hares may have a negative effect on the native population, with the possible emergence of a new hybrid species and the pushing of the native one toward smaller and more inhospitable areas, as discussed in [31].
We now consider the second approach for indirect control of cottontail population.We introduce an unknown input function G prey (t), to modify the reference Model (2) so that now the last equation becomes while the other two equations are the same as in (2).From now on we will refer to G pred (t) as G(t) for simplicity.The idea behind this choice is that we consider an additional "immigration" term for hares in order to preserve them from extinction.We apply the Zcontrol method to drive S(t) to the desired state S d .Since G(t) is part of the third equation, we need to impose error conditions mathematically until the derivative of L with respect to time is involved.In particular, we prove that by repeating the Z-control process up to the third error, an explicit expression for G(t) can be obtained.Recalling ( 4) and ( 5), we define the third error function as follows and impose once again dv 3 dt (t) = −λv 3 (t).
Back substitution now leads to Rewriting it in terms of S, we obtain from which an explicit expression G prey (t) is obtained for SVL = 0 where We can therefore rewrite the right-hand side of the third equation more conveniently as:

Equilibria
For this model as well, we search for possible stationary configurations.Solving the equilibrium equations leads to a unique coexistence equilibrium E = ( V, S, L) since all the populations cannot be zero.Assuming G d (t) to be constant, the coexistence population values are Once again, we observe that by forcing the solution S(t), at equilibrium we obtain a value of S = S d , which cannot be zero.The feasibility conditions for V and S are the same as (11), while for the third one L > 0, we now obtain

Local Stability Analysis
For the local stability analysis, proceeding similarly as in the previous section, we consider the Jacobian matrix: , where letting H(V, S) = (c SS S) 2 − 3λc SS S + 3λ 2 − ea 2 VS , the derivatives of the third row at equilibrium have the following expressions: Proposition 2. The coexistence equilibrium E of the Z-controlled Model (2) together with (16) is locally asymptotically stable if and only if where these quantities are defined in the proof, (20) and (21).
Proof.The eigenvalues of J( E) are solutions of the characteristic equation with coefficients Considering the minors involved in the RH criterion applied to the monic polynomial in the brackets of ( 19) E is asymptotically stable if and only if (18) holds.

Combined Control Over Two Species 2.5.1. Direct and Indirect Control Combination Providing Cottontails' Extinction and Hares' Survival
The results of the previous sections indicate that a single indirect control, by acting on a species other than the invasive one, does not produce optimal results, because the hare-free equilibrium may arise, entailing the hares' extinction.
We thus try a new approach, where we study the effect produced by a simultaneous control of classes S and L.More precisely, the idea is to choose two desired states, S d (t) and L d (t), in order to drive both prey populations to a predefined configuration.To this end, we will consistently choose S d (t) ≡ S d ≈ 0 and L d (t) ≡ L d > 0 ∀t > 0. Introducing two different control functions for foxes and hares, respectively, G V (t) and G L (t), we modify the first and third equations of the reference Model (2), obtaining We separately apply a Z-type control method to the two populations starting from two error functions We are thus producing an indirect control for the cottontail population and a direct one for the hares, because G L (t) directly affects their dynamics.
The approach we propose has similarities with the Z-control method.The difference between this approach and the conventional one lies in defining the two error components (23).However, we cannot impose joint design conditions on the components of the vector-valued error function v 1 (t) = (v 1,S , v 1,L ).This is because we need to act separately on each one of them, as explained in the following.
To find an explicit expression for G V (t), we have to impose different conditions for the two parts of the error, v 1,S (t) and v 1,L (t).Indeed, note that for the latter only one condition is required, since the hares' dynamics are directly controlled, namely, Instead, because the only way to find G V is to manipulate an equation containing a derivative of V with respect to time, for the former we need to use both v 1,S (t) and v 2,S (t) as is done in Section 2.3.1.Now, the conditions are respectively equivalent to As was previously done, we will restrict the study to the case in which G d,S (t) ≡ λ 2 S d ≈ 0 and G d,L (t) ≡ λL d = 0.
First, assuming S, V = 0, from (25), we obtain a general explicit expression for the control G V (t) From (26) taking L = 0, we obtain instead In Theorem 1, it has already been shown that S(t) → S d as t → +∞.Therefore, we only need to prove the convergence of L(t) to the chosen state.Theorem 2. For any positive initial state (V 0 , S 0 , L 0 ) and for any continuously differentiable and bounded desired state L d (t), the tracking error function v 1,L (t) converges to zero exponentially.

How to Choose the Hares' Desired State L d
In the context of the optimal choice for L d , it is important to make two remarks: 1.
L d should be chosen in such a way that L d ≤ K (L) .Here, K (L) denotes the carrying capacity related to this class, as well as the population level attained by hares in the absence of the other two species, i.e., the value L 4 in the two-vanishing populations equilibrium E 4 = (0, 0, L 4 ) [24].Consistently with the previous sections, in the simulations we will use a value for the carrying capacity within the range of approximately 26-40 hares/km 2 [32], and so K (L) = 30 hares/km 2 .
2. It would be advisable to drive L(t) to the hares value attained at the native species coexistence equilibrium.This value corresponds to the one at the cottontail-free equilibrium E 5 of the reference model ( 2) [24].

Equilibria
In this scenario with both foxes and hares subject to forcing, the hare-free equilibrium does not exist, in contrast to the case of single action on foxes.The Z-controlled model (22) constructed in this way exhibits a unique coexistence equilibrium.Taking again V, S, L = 0 for an arbitrary choice of constant desired state functions, it is Since V, S, L = 0, we cannot take S d = 0. Furthermore, E c is biologically consistent if the general feasibility conditions hold

Local Stability Analysis
The Jacobian matrix J of the system is almost the same as the one in Section 2.3.1, with only a change in the third row; indeed, here derivatives with respect to V and L of G L (t) also need to be taken.Proposition 3.For any parameter choice, including the design parameter λ > 0, the equilibrium E c , when feasible, is unconditionally locally asymptotically stable.
Proof.The Jacobian evaluated at E c is It thus factorizes to give the negative eigenvalue −λ < 0 which has no influence on the equilibrium stability.The remaining ones are the roots of the quadratic equation where A = −2λ < 0 and B = λ 2 > 0. The latter represent the RH conditions, which are thus unconditionally satisfied and the claim follows.

Results
In this section, we show several numerical simulations to support the theoretical results previously obtained.In all simulations, we use the ecological parameters of [24] and the following initial conditions for all models

Numerical Simulations for the Action on Foxes
Suppose we want to reduce the invasive population to a lower value, so that S d is set to, say, 0.5.We distinguish two cases for the choice of the parameter b.A low value, specifically when b = 0.3, simulates a situation where foxes hunt the two prey almost indistinguishably, at the same rates b ∼ a.On the contrary, the larger value b = 3 represents a scenario of hyperpredation on hares.We observe that b * ≈ 0.4076 and therefore these two different situations could indeed occur.Figure 2 shows the model with no control, respectively, with b = 0.3 and b = 3.In the first case, the cottontail population naturally decreases and with this choice of parameters the equilibrium E 5 of the uncontrolled reference model (2) is feasible and locally stable.The equilibrium values for the native species are, respectively, V 5 ≈ 3.5529 and L 5 ≈ 7.6017 (27) and these levels are reached relatively soon, after about 3 years.This scenario would not require an excessive control, as the invasive species naturally tends to disappear.However, we discuss how the application of Z-control can affect this process and which scenario may be more advantageous.Figures 3 and 4 show the population behavior in the controlled model with b = 0.3, as well as the control function G pred (t) with λ = 0.01, λ = 0.025, λ = 0.05, λ = 0.1 and λ = 0.7, respectively.We observe that, since b is low, a control with a small λ, for example λ = 0.01 or λ = 0.025, should already produce a good result.In this case, however, a weak control initially produces an increase in the cottontail population and thus slows down the natural process for cottontails' extinction, so its employment does not add any advantages.Additional forcing means a higher convergence speed to the desired state, as expected from a theoretical point of view, and a gradual lowering of the cottontails' curve.Furthermore, we observe that when a faster convergence is required, as in the case of λ = 0.1, the control function G pred (t) indicates that a considerably high effort is necessary in order to guide the populations towards the configuration E c .The stable equilibrium of the controlled model is E c = (3.5028,0.5, 7.9177); comparing these values with those of E 5 , compare (27), we observe that the control does not result in a significantly better situation than the uncontrolled one.Focusing now on the second scenario, represented in Figures 5 and 6, we observe that in the uncontrolled model (b = 3), the cottontails' invasion leads to the extinction of the native prey L and the local feasibility and stability conditions for equilibrium E 3 = (3.1779,9.7277, 0), where the predator and invasive species coexist, are satisfied.In this situation, the Z-control action plays a crucial role to eradicate the invasive species but it is not sufficient to ensure the survival of hares since, with b > b * , the only admissible stationary configuration is represented by E L = (3.5028,0.5, 0).In particular, E L is locally asymptotically stable.As in the previous case, as the value of λ increases, the convergence of S to S d strengthens, and the effort employed by G pred (t) increases as well.

Numerical Simulations for the Action on Hares
We now discuss the second approach (16) employed for reducing the invasive species, where the desired state is assumed to be S d = 0.5.For this model, the equilibrium E = (3.5028,S d , L) is always feasible with the choice made on parameters, but its local stability depends on the λ values.We also observe that the hares value L at the equilibrium depends on b: according to (17), it is inversely proportional to the hunting parameter.In order to investigate the stability conditions (18) for this situation, we can graphically observe in Figure 7 that D < 0, i.e., −λ 3 + (c VV V + c SS S d − 3λ)(c VV + ea)(aS d ) V < 0, if λ > 0.6668.Moreover, the numerical real solution resulting from the intersection of the cubic polynomials 3λ B(λ) and D(λ) involved in (18), i.e., the real solution of 7).Therefore, local stability for E is guaranteed if λ > 0.6668 and does not depend on the choice of b.In line with the previous section, we perform simulations to explore the system's behavior which may lead to instability scenarios, across a range of values for λ together with the two previously used values of b.In these situations, the convergence speed to the equilibrium is quite high and, once again, it increases for higher values of λ.
To observe a different situation, let us consider, for example, the case where b = 3 and λ = 0.025, illustrated in Figure 10: the control function in the bottom panel is unbounded and, after approximately 5 years, the introduction of hares into the territory grows indefinitely, therefore influencing the native population's behavior.This fact could be interpreted as a restoration of the coexistence of native species, as the growth is attributed to a continuous introduction of prey into the territory.Indeed, the cottontail population decrease with time, leading to extinction.

Numerical Simulations for the Combined Action on Foxes and Hares
We finally consider the model (22) with action on both foxes and hares choosing as desired states the values S d = 0.5 and L d = L 5 , as first trial.We will then discuss the case L d < K (L) and compare which choice may be better.For this controlled model, the equilibrium is E * c = (3.5028,S d , L d ).We observe that L 5 , whose expression is presented again in the next equation, depends on the choice of b and in particular its value decreases as b increases (assuming the other parameters are fixed) Therefore, in order to have a feasible equilibrium, b has to be chosen so that and, furthermore, if b is larger than a certain threshold, the use of this method could drive the class L towards a desired state L 5 that could represent extinction.Figure 11 shows the simulations in the case in which ( 28) is satisfied, particularly with b = 0.3 (for different values of λ) and b = 1.In these two cases, L 5 = 7.6017 and L 5 = 0.3664, respectively.From Figure 11, we observe that an increase in the value of λ leads to a faster convergence.In this case, we also note that, as this parameter increases, the need for adding hares becomes lower and lower and consequently this is reflected in the decreasing of the function G L (t).This fact can be explained by observing that higher values of λ better reflect the natural undisturbed system behavior.Instead, lower ones impose a convergence speed slower than the one which is naturally exhibited; counterintuitively, it produces a general slowdown of the convergence process.This point is also apparent in the sequence of panels of Figure 3. On the other hand, it appears that above a certain threshold, the system is forced to converge at a faster rate than its natural behavior, resulting in greater external effort.Figure 12 shows the changes in the cottontails' curve for different values of λ in the case of the controlled model with action on foxes (left panel) compared to the double action on V and L (right panel).In particular, smaller values of λ show a peak in the cottontail population, which decreases with the increase in this parameter and thus better approaches the natural trajectory.However, for early times, the best situation appears to be given by an intermediate value, such as λ = 0.25.As observed in previous cases, in the uncontrolled system with b = 0.3, the invasive species naturally goes extinct, and the dual action through the Z-control works from a mathematical perspective but is not necessary in practice.On the other hand, by increasing b, e.g., b = 1, the class L is indeed driven towards a lower value (for instance L 5 = 0.3664).If b does not satisfy (28), hares are driven to extinction because L 5 turns out to be negative.Figure 13 shows the behavior of the class L and the control functions for b varying between 0.3 and 1.8.To conclude this section, the case L d < K (L) is considered.For example, set L d approximately equal to the value achieved at the equilibrium configuration E 5 when b = 0.3 as calculated in (27), so L d = 7.6.An instance of this situation is presented in Figure 14.

Biological Interpretation for Indirect Control on Foxes
For the first model, external intervention is made through predator removal or reintroduction into the territory.Numerical results enable us to primarily compare two distinct situations.For the first one, where b = 0.3, one may expect that the application of Z-control speeds up the process of extinction of the cottontails.This is possible from a theoretical point of view, since the cottontails' class converges to the desired state, but it may be resource-demanding from a practical perspective.For instance, the first row of Figure 4 illustrates that around half a month after the invasion, the number of individuals to be removed or relocated exceeds 100.Moreover, the values attained at the equilibrium reflect the scenario that would naturally occur without any intervention.As a further example, looking at the second row of Figure 4, we observe in the graph of G pred that predator insertion and removal should alternate, presenting a more practical approach.However, non-intervention might be preferable.
On the other hand, the second scenario, where b = 3, represents an invasive phenomenon where the invasive species settles in the territory, negatively affecting the survival of the native prey.Unlike the previous case, here both prey populations go extinct.Hence, the natural pre-invasion dynamics between native species are disrupted and predators may be forced to seek new resources elsewhere.

Biological Interpretation for Indirect Control on Hares
Among the two methods indirectly applied to control the cottontails' curve, the second one may yield better results.In this situation, the behavior of the function G prey represents the introduction or removal of native hares, depending on its sign.From an analysis of Figures 8-10, we can observe that, in addition to the removal of the invasive species, the survival of hares can also be assured without being imposed.Furthermore, even for those values of λ that do not make the unique equilibrium locally stable, situations of coexistence among the native species can still be achieved.Hence, a restoration of the dynamic interaction between foxes and hares becomes feasible for these simulations.

Biological Interpretation for Combined Indirect Control on Foxes and Hares
By combining the previous approaches, it becomes possible to impose two conditions to ensure both the removal of the invasive prey and the survival of the native one.In principle, this technique is expected to be more resource-intensive since it involves two parallel efforts.However, it may also be more effective, as the controlled classes converge to the desired states.Specifically, forcing the hares' population towards a level L d , it is possible to achieve favorable outcomes with control functions that do not entail excessive effort, as shown in Figure 11.Therefore, this method can provide, for the situation examined in this work, a balanced strategy to address the invasion issue of the pre-existing ecosystem.

Conclusions
In this paper, the Z-type control method was applied in an ecological context involving a native predator, the fox, and two preys, one native and one invasive, hares and cottontails, respectively.Our goal was to explore some situations in which the invasive species could be kept at low values while preserving the critically endangered native prey, in order to possibly give indications to ecosystem managers.Specifically, we considered three possible alternatives to act on the ecosystem.At first, we considered predator removal in order to curb the invasive cottontails causing the hyperpredation effect on hares.We applied an indirect control on the cottontails acting on foxes to reduce the pressure on hares.Then, we considered the possibility of importing new individuals of the indigenous prey population from outside, employing a direct control on this class.This choice sustains the hare population although it also provides more resources for foxes.However, this may entail other problems due to the emergence of hybrid species, a problem addressed in [31].Finally, we combined the two strategies initially proposed: direct control on hares and indirect control on the cottontails acting on foxes.In the first case, the controlled model can evolve toward two equilibrium points, coexistence and hare-free equilibrium.Both are found to be conditionally admissible and locally asymptotically stable.Moreover, a transition from one to the other is possible through a transcritical bifurcation as the value of the hunting rate of foxes on hares varies.In the last two cases, instead, the controlled system can only converge to a conditionally admissible and stable coexistence equilibrium.Consequently, with a view to preserving the hare population, we can conclude that the best solutions among the three proposed to manage the cottontail invasion problem would seem to be the last two.Indeed, by applying indirect control on cottontails acting on foxes, the hare-free equilibrium may arise, entailing the hares' extinction.This cannot happen in the remaining two cases because the corresponding controlled models can only evolve toward the coexistence of the three populations involved, where the cottontail population can be kept as small as desired.In this regard, however, we stress that choosing either of the last two proposed strategies could give rise to other secondary problems.Considering the numerical results, comparisons between the three models and the original unmodified one have been made.Specifically, the attention was focused on different scenarios relying on the choice of the parameter b.The situations in which it was assumed that hares are hunted at rate b = 3, approximately 10 times higher than the value chosen for a, may be more interesting and significant from a biological point of view.Instead, the case in which b was kept lower, e.g., b = 0.3, could be more interesting from a theoretical perspective, since the system naturally evolves to the desired state.For the specific situation where b = 3, our conclusions indicate that it is not necessary to apply the Z-control method, mainly for two reasons.On one hand, the cottontails' curve undergoes a significant slowdown in the convergence to zero, since the choice of λ greatly affects the dynamics of the system; on the other one, in some cases, requiring higher convergence speed to the desired solution involves a great effort that could become not practicable.In conclusion, our goal was primarily to show that the method is at least theoretically viable, a task that has been achieved, and to explore different ways of implementing the Z-control approach in a real-world scenario.

Figures 8 and 9
Figures 8 and 9 show the controlled model with λ = 0.7 and λ = 1.4,both for b = 0.3 and b = 3, where the values for L are L (0.3) ≈ 7.1191 and L (3) ≈ 0.7119, respectively.In these situations, the convergence speed to the equilibrium is quite high and, once again, it increases for higher values of λ.To observe a different situation, let us consider, for example, the case where b = 3 and λ = 0.025, illustrated in Figure10: the control function in the bottom panel is unbounded and, after approximately 5 years, the introduction of hares into the territory grows indefinitely, therefore influencing the native population's behavior.This fact could be interpreted as a restoration of the coexistence of native species, as the growth is attributed to a continuous introduction of prey into the territory.Indeed, the cottontail population decrease with time, leading to extinction.

Figure 6 .
Figure 6.(Left): Population behavior and control function G pred (t) for parameters b = 3 and λ = 0.1.(Right): zoom of the control function near the origin.

Figure 10 .
Figure 10.Population behavior (left) and control function G prey (t) (right) as function of time for parameters b = 3 and λ = 0.025.

Figure 12 .
Figure 12.Comparison between the cottontails' curve of the uncontrolled Model (2) for b = 0.3 (continuous line identified with nc in the legend) and the controlled models with single action on foxes (Model (3), (left) panel) and action on foxes and hares (Model (22), (right) panel) for different values of λ.

Figure 14 .
Figure 14.Population behavior ((left) panel) and control functions G V (t) and G L (t) ((right) panel, respectively, in magenta and green) as a function of time for parameters b = 1 and λ = 0.05.