1. Introduction
The stability of an ecosystem has always been the focus of biologists and it involves many fields, such as the protection of the ecosystem, economic development, political management and so on. Traditionally, scholars are interested in the asymptotic stability of ecological models and have conducted several qualitative analyses and quantitative calculations. In previous studies, it is known that asymptotic stability can be characterized by resilience. In 1977, Pimm and Lawton [
1] proposed the definition of resilience in a biological sense for the first time, which means the asymptotic decay rate of perturbations. Given a linearized (or linear) system
The resilience is calculated as , is the eigenvalue of A with the largest real part. The greater the resilience, the faster the decay rate of the disturbance in the long run and the more stable the system. Although the long-term asymptotic dynamic behavior is significant, the short-term transient dynamic behavior also cannot be ignored. Some systems are ultimately stable; however, the trajectory of these systems may rapidly move away from the equilibrium in the early stage after being disturbed. For example, the pests that will eventually become extinct can surge in the beginning. The infectious diseases, which will be solved, can cause a large number of deaths at first. Therefore, if we only consider the final state and ignore the temporary behavior of the system, we will suffer serious losses.
In the past few years, scholars have made remarkable achievements in transient dynamics for the ecosystem. In 1997, Neubert and Caswell [
2] introduced reactivity as a supplement to resilience, which means the maximum initial rate that a small perturbation can grow. It characterizes the instantaneous responses (
). The reactivity is denoted by
,
is the Hermitian part of the Jacobian matrix
A. The equilibria is called reactive if
. The larger the reactivity, the more sensitive the system and the faster a perturbation can grow. When
and
, the system is eventually stable but reactive. This situation is possible and not accidental.
Resilience and reactivity describe two limits of solution behaviors
and
, but neither describes the transient behaviors between zero and infinity. Therefore, Neubert and Caswell [
2] defined the amplification envelope
. It is the maximum possible amplification that any initial-value perturbation could induce in a solution at time
t, given by
, where
is the Euclidean matrix norm of
. In addition, the maximum value of
is denoted by
. The time at which
occurs is denoted as
.
Based on these conclusions, scholars further studied the resilience, reactivity and envelope amplification of ecosystems. Regarding resilience, Folke, Carpenter and Walker [
3] had a heated discussion around resilience and put forward several significant concepts related to ecological resilience, such as adaptability, transformability, etc. In the same year, Hughes and Graham et al. [
4] found that resilience played an important role in the management of marine ecosystems (such as coral reef systems), and proposed that we should learn how to avoid unnecessary phase transformations of ecosystems. In 2019, Wright [
5] proposed a model of social resilience and sustainable development in the 21st century from the lessons of historical ecology. For reactivity, Snyder [
6] analyzed what makes the ecosystem reactive. Then, Vesipa and Ridolfi [
7] also focused on the short-term dynamics of long-term stable and reactive ecosystems in 2017. Taking a seasonal forced prey–predator model as an example, it is proved that seasonality can greatly affect the transient dynamics of the system.
In 2014, Evelyn Buckwar and Conall Kelly [
8] extend the above measure indicators to stochastic cases for a perturbed predator–prey model, which are called root-mean-square resilience, root-mean-square reactivity and root-mean-square amplification envelope. Since ecosystems exist in a real world and are affected by the external environment continuously, it is more reasonable to consider a random biological model. Compared with the traditional predator–prey model, the ratio-dependent predator–prey model proposed by Arditi and Ginzburg [
9] in 1989 avoids the limitation that the trophic function of the traditional predator–prey model ignores the influence of predators. Therefore, it is more meaningful to consider that the stochastic ratio-dependent predator–prey model, whose characteristic is that a functional response depends on the ratio of prey to predator density.
To the best of our knowledge, Berezovskaya, Karev and Arditi [
10] made a complete parameter analysis of the stability properties and dynamics of the deterministic ratio-dependent predator–prey model. Recently, some scholars have made careful studies on the dynamic behaviors of the deterministic predator–predator models, such as Chen and Zhang [
11] in 2021 and Molla, Sahabuddin and Sajid [
12], in 2022. Furthermore, Bandyopadhyay and Chattopadhyay [
13] discussed the deterministic and stochastic behaviors and stochastic stability near the internal equilibrium of the ratio-dependent predator–predator model. They also gave the conditions for the deterministic model to enter the Hopf bifurcation. Wu, Huang and Wang [
14] researched the dynamical behavior of a stochastic ratio-dependent predator–prey system. They obtained the conclusions of stochastically ultimate boundedness and stochastic permanence, among which persistence is one of the most important characteristics of the biological model. In addition, in 2016, Fan, Teng and Muhammadhaji [
15] analyzed the global dynamics of a stochastic ratio-dependent predator–prey system. They gave results for the existence, uniqueness and global attractivity of the solution for the stochastic model. Ji, Jiang and Li [
16] described the connection and difference between a deterministic and random ratio-dependent predator–predator model. As mentioned in the article, the previous literature had shown that the solution of the deterministic model is uniformly bounded. In addition, the author had obtained that the solution of the random model is uniformly bounded in mean, which is similar to the deterministic system. They also studied the persistence and extinction of model populations. The conditions under which the stochastic model species is not permanent but the deterministic model species may be permanent are given.
In this paper, we will study the root-mean-square resilience, root-mean-square reactivity and root-mean-square amplification envelope of stochastic ratio-dependent predator–prey model. The layout of the paper is as follows. In
Section 2, we briefly introduce the studied model and some basic concepts. Then, in
Section 3, we compute the mean-square stability matrix
S of stochastic ratio-dependent prey–predator model by the Kronecker product. Meanwhile, we use the numerical method to simulate parameter graphs of the three measure indicators: root-mean-square resilience, root-mean-square reactivity and root-mean-square amplification envelope. Finally, the analysis conclusions are obtained by comparing the numerical figures in
Section 4.
4. Numerical Simulation and Analysis
First, choose the appropriate parameter values
for the drawing of root-mean-square resilience, root-mean-square reactivity and root-mean-square amplification envelope numerical figures. (The selection of parameters refers to aticle [
13]). When
and
, the model system has a stable Hopf-bifurcating periodic solution and all other trajectories around the limit cycle eventually approach it.
4.1. Root-Mean-Square Resilience and Root-Mean-Square Reactivity
Figure 1 shows the variation of the root-mean-square resilience and root-mean-square reactivity of the studied system within the parameter range when there is no random disturbance. It can be used as a comparison figure for analyzing the influence of random environmental fluctuations on the system. In
Figure 1, the parameter ranges of the mean-square stable stage for this system are
, respectively. Resilience shows a trend of monotonous increase first and then decrease. If
b and
c are changing parameters, the interior equilibria of the system are both reactive (reactivity
).
The first column of
Figure 2 shows the perturbation only applied to the prey, the disturbance intensity is 0.2 and 0.4, correspondingly. The variable parameter is a. It can be seen that the mean-square stable area of the system is reduced from
to
and
. The maximum value of resilience (representing the system’s self-recovery ability; the greater the resilience, the more stable the system) is reduced from 0.357 to 0.313 and 0.275. The second column is the disturbance to predator and, similarly, the mean-square stable area is reduced from
to
and
. The maximum value of resilience is reduced from 0.357 to 0.298 and 0.210. Regarding the reactivity, there is no obvious change when the disturbance is applied to the prey, but when it acts on the predator, the initial value of the reactivity becomes larger as the intensity of the disturbance increases. Summarizing the above conclusions, the system is more sensitive to predator disturbance and reacts more strongly.
In
Figure 3, the range of the change parameter b is
. When the perturbation is only on the prey, the mean-square stable regions of the system are
and
, respectively. The maximum values of resilience are 0.198 and 0.163. When
and
, the model is mean-square stable if the predator is disturbed. Comparing the four sub-figures in
Figure 3, it is found that increasing the intensity of the disturbance increases the reactivity, the resilience decreases and the mean-square stable area shrinks. Similar to
Figure 2, it is still the predator disturbance that has a stronger impact on the system.
The trends of resilience and reactivity are opposite, which also shows that the more stable the system, the stronger the resilience; the more sensitive the system, the stronger the reactivity.
Increasing the random perturbation to any species (predator or prey) will reduce the range of variable parameters of the system in the mean-square asymptotical stable state (where the values of root-mean-square resilience are positive).
In
Figure 2,
Figure 3 and
Figure 4, we can clearly see that when perturbation is added to any population, the maximum root-mean-square resilience is visibly reduced. In contrast, when the random perturbation is increased for any species, the value of root-mean-square reactivity increases but the increase is not obvious.
The same random disturbance intensity acts on predator and prey, respectively, the system will be more strongly affected when the predator is disturbed.
4.2. Root-Mean-Square Amplification Envelope (Hopf-Bifurcating Periodic Solution)
Figure 5 shows the root-mean-square amplification envelope of the positive equilibrium of (
11) with
. The positive equilibrium is a stable Hopf-bifurcating periodic solution.
Figure 5 illustrates the type of this solution because the change of disturbance is also periodic.
It can be found that adding perturbation to any population at the Hopf-bifurcating equilibrium (even if the perturbation is small) will make the system unstable. In comparison, perturbation to the predator makes the system more sensitive.
Conclusions can be drawn from all the above analyses, and some appropriate suggestions can be made to managers for this ecosystem. For example, when regulating the ecosystem, pay attention to controlling the parameters in the mean-square asymptotically stable state of system (the value of the root-mean-square resilience is positive). Second, according to the numerical graph of root-mean-square amplification envelope, managers can clearly know the maximum value of disturbance magnification: and the time when the maximum disturbance appears: . If the manager sets a small disturbance to the ecosystem, it needs to consider that the perturbation amplification should be controlled within the range that the system can bear, that is, the system will be mean-square asymptotically stable. gives critical information, and decision makers should give orders and actively implement them before , so as to avoid crucial losses and save money and manpower.
In order to make the reader better understand the significance of root-mean-square amplification envelope, the following numerical example is listed.
4.3. Root-Mean-Square Amplification Envelope (Locally Stable Positive Equilibrium)
In this example, select the parameter values
. According to literature [
13], the interior equilibrium is locally stable in this case. From
Figure 6, we find that when the predator is perturbed, and the perturbation intensity reaches 0.4, the system can still restore a stable state through its own regulation. However, the perturbation can no longer disappear and even shows a constant increasing trend when the prey is perturbed with the perturbation intensity 0.4. At the same time, we can see that no matter which species is disturbed, the greater the perturbation intensity, the larger the amplification envelope.