Stability of nonlinear population systems with individual scale and migration

: In this paper, we study the stability of a nonlinear population system with a weighted total size of scale structure and migration in a polluted environment, where fertility and mortality depend on the density in di ﬀ erent ways. We ﬁrst prove the existence and uniqueness of the equilibrium point via a contraction mapping and give the expression for the equilibrium point. Some conditions for asymptotic stability and instability are presented by means of a characteristic equation. When the e ﬀ ect of density restriction on mortality is not considered, the threshold value of equilibrium stability can be obtained as Λ = 0 . When Λ < 0 , the equilibrium is asymptotically stable, and when Λ > 0 , the equilibrium is unstable. In addition, the upwind di ﬀ erence method is used to discrete the model, and two examples are given to show the evolution of species.


Introduction
In recent years, the status of protecting the environment and population diversity has been improved, and environmental pollution has become the top priority of environmental issues. Our primary concern is the interaction between environmental pollution and biological populations. We need to consider the interdependence between biological populations as well as the interdependence between populations and the environment. In the dynamical model of a biological population, the equilibrium of the model corresponds to the static equilibrium of the biological system, and the stability of equilibrium describes the stability of ecological balance and represents the anti-interference ability of the ecosystem. Therefore, it is paramount to study the stability of the population dynamics model. So far, scholars at home and abroad have established a large number of mathematical models in the process of exploring the evolution of biological populations. Age is an important parameter, affecting the size and behavior of the population. Sharpe and Lotka [1] established a linear model system with an age structure. Since the birth and death rate of individuals are always affected by the number of populations, Curtin and Maccamy [2] established the age-structured population model with density constraint, proved the existence and uniqueness of the model, and obtained that the equilibrium is locally asymptotically stable. Considering the environmental constraints on population growth, a nonlinear age-structured population model was studied [3]. For the population model with age structure, many scholars have studied its optimal control problem [4][5][6][7]. In [8], the nonlinear age-structured population model with nonlocal diffusion and boundary conditions was studied. Lopez et al. [9] studied the discrete nonlinear age-structured population model of environmental changes and obtained conditions for the existence and asymptotic stability of positive equilibrium. Kilic et al. [10] used the forward Euler algorithm to study the discrete-time prey-predator population model with migration, and obtained the asymptotic stability of the unique positive fixed point. Studies by ecologists show that the scale structure is a continuous variable describing the individual characteristics of the population (volume, diameter, length, maturity, etc.), it can better reflect the relationship between groups and individuals in the evolution of biological populations. So it has attracted extensive attention from ecologists and mathematicians [11][12][13]. Dercole et al. [14] explored a nonlinear plant growth model in which the scale variable represented the base diameter of the tree, and obtained the stable properties of the equilibrium. Calsina et al. [15,16] considered environmental factors, they studied the scale structure models when individual growth rate, death rate, and birth rate were related to the environment and the total population. Farkas [17][18][19] successively studied a class of nonlinear and migration models based on scale structure, in which mortality rate and birth rate were related to the total population, and obtained the stability conditions of the equilibrium. Kang et al. [20] studied a nonlinear double physiologically structured population models with two internal variables and obtained the stability of the steady state. Ademosu et al. [21] analyzed the stability and optimal measure for controlling eco-epidemiological dynamics of the prey-predator model. Singh and Emerick [22] studied generalized stability conditions for host-parasitoid population dynamics and identified implications for biological control. Molla et al. [23] analyzed a mathematical model on predator-prey interactions incorporating prey refuge and additive Allee effect on the prey species and obtained the various dynamical behaviors. Wu [24] studied an optimal harvesting problem of a size-structured population model: where p(x, t) represented the number of individuals with scale x at time t, g(x) represented the growth rate, µ, β represented individual natural mortality rate and reproductive rate respectively, p 0 (x) represented the initial distribution of the population, q i (x), i = 1, 2 were different weight functions used to reflect the dependence of the lifetime parameters on density, u(x) represented the harvest intensity of humans to the individual population.
With the development of industry and the intensification of human activities, a large number of toxic substances have been discharged into the ecological environment, destroying the habitat and community structure of species, and seriously threatening the survival of the species exposed to the environment (humans, fish, plankton, etc.). Liu et al. [25] studied the persistence and extinction of a single-species population system with random disturbance and pulse toxicity in a polluted environment. Qi et al. [26] studied the dynamic behavior of a stage-structured population model with transient and non-transient pulse effects in a polluted environment. Wei et al. [27] studied the psychological effect on single-species population model in a polluted environment and obtained the conditions to ensure local extinction and average persistence of the population. Ma et al. [28] studied the stationary distribution and optimal control of a stochastic population model in a polluted environment.
However, there is little literature on the stability of scale-structured populations in a polluted environment. Considering the biological background of the population model, this paper focuses on the external living environment of species and the impact of migration on species survival, so that the bio-mathematics model is more realistic. Therefore, based on reference [24], we study the following population model with scale structure and migration in a polluted environment. Generally speaking, the increase in population has different effects on the birth and death rates of individuals. We consider the density restriction in the process of population evolution with scale structure, and adopt different individual weighting functions to represent the different ways in which life parameters depend on density. The model is as follows: 2) where the fixed constant m > 0, which represents the largest scale of the individual, 0 < T ≤ ∞, T represents the maximum time of the evolution. The variables of the model are described in Table 1.
Note 1: In system (1.1), l represents the individual scale, here denotes as s. In system (1.2), l becomes m, and m represents the largest scale of the individual.
Note 2: We consider the external migration of individuals, which is assumed as f (s) ≥ 0. In a certain environment, due to limited resources, the birth rate will be affected by the proportion of females in the population density γ(s) , that is, the birth rate will decrease with the increase of the population density, and the mortality rate will be affected by the weight function of total population δ(s), that is, the mortality rate increases with the increase of population density. This is reasonable in the actual external environment. individual scale growth rate c 0 (t) the concentrations of the pollutants in organisms c e (t) the concentrations of the pollutants in environment v(t) the pollutants that are input from external to internal environment δ(s) the weight function controlling the total population γ(s) the proportion of females in the population density f (s) the external migration of individual populations p 0 (s) the initial scale distribution k the net absorption rate of toxins in living organisms η the rate of purification of mycomycin by metabolic mechanism in organism k 1 the rate of uptake of mycomycin by the population in the environment g 1 the ratio of pollutants released by organisms into the environment h self-purification rate of mycomycin in the environment P(t) the total weight of the population J(t) the weighted total population at time t R(t) total number of females in the population at time t µ(s, c 0 (t), J(t)) the mortality rate of an individual with scale s and concentration c 0 (t) based on a weighted total J(t) β(s, c 0 (t), R(t)) the birth rate of an individual with scale s and concentration c 0 (t) based on a weighted total R(t) u(s, c 0 (t), J(t)) the harvest rate of an individual with scale s and concentration c 0 (t) based on a weighted total J(t) k, g, η, k 1 , g 1 , h non-negative constant We make the following assumptions: (H 1 ) f (s) is a non-negative continuous function; (H 2 )μ(s, t) := µ(s, c 0 (t), J(t)) ∈ L 1 loc ((0, m) × (0, T )), 0 ≤ µ(s, c 0 (t), J(t)) ≤ µ 0 , ∀x ∈ (0, +∞), t ∈ (0, T ), µ(·, c 0 (t), x) is monotonically increasing with respect to x and locally Lipschitz continuous |µ(s, c 0 (t), x 1 ) − µ(s, c 0 (t), x 2 )| ≤ L(n 1 )|x 1 − x 2 |, L(n 1 ) is the Lipschitz constant; (H 3 )β(s, t) := β(s, c 0 (t), R(t)) ∈ L 1 loc ((0, m) × (0, T )), 0 ≤ β(s, c 0 (t), R(t)) ≤ β 0 , ∀x ∈ (0, +∞), t ∈ (0, T ), β(·, c 0 (t), x) is monotonically decreasing with respect to x and locally Lipschitz continuous |β(s, c 0 (t), This paper mainly studies the stability of equilibrium. Over time, the population either persists or becomes extinct. The extinction of a population means that the total number of individuals in the population tends to zero in a finite or infinite time interval. However, with the continuous development of the population, does the population have one or more equilibria? If so, under what conditions are these equilibria stable or unstable? These problems have guiding significance for the sustainability of ecological balance.
Note 3: In [29], we have established the well-posedness of the system, the system owns a unique solution, which is non-negative, bounded and continuous with respect to the initial distributions. The method of proving the existence and uniqueness of model (1.2) in this paper is similar to that in [29].

The existence and uniqueness of the equilibrium solution
When discussing stability, it is assumed that there is no human interference, namely u ≡ 0.
Since the second and the third equations in system (1.2) are two first order ordinary differential equations, which could be solved by method of variation of constant, we get For system (1.2), when it satisfies the condition 0 < k If there is an equilibrium solution p * (s) in model (1.2), the equilibrium solution satisfies the following equation: From the first equation of (2.1), it can be determined dr , solved by the above formula Substituting the above formula into (2.1), it can be determined The net regeneration number of this population is defined as According to (2.3) , N(c 0 (t), R, J) must be a normal number less than 1.
Next, we discuss the existence and uniqueness of the positive solutions of (2.3).
It can be known by (H 2 ), (H 3 ) and the monotonicity of Then, we obtain by the first formula of (2.3) that From the fourth formula of (2.3), it can be determined that From the fifth formula of (2.3), it can be determined that a 2 , a 3 , a 4 , a 5 ) T ∈ Ω, and define the mapping φ : Firstly, since A ∈ Ω, there is F(s, s, a 2 , a 4 )ds Since 0 ≤ c 0 * ≤ 1, 0 ≤ c e * ≤ 1, therefore, 0 ≤ y 2 ≤ 1, 0 ≤ y 3 ≤ 1.
That is Y = (y 1 , y 2 , y 3 , y 4 , y 5 ) T ∈ Ω. Secondly, it can be proved that φ maps Ω into itself, let It can be determined from the first formula of (2.6) that F(s, s, a 2 2 , a 2 4 )ds| It can be determined from the second formula of (2.6) that It can be determined from the third formula of (2.6) that − F(s, r, a 1 2 , a 2 4 ) + F(s, r, a 1 2 , a 2 4 ) − F(s, r, a 2 2 , a 2 4 ))drds It can be determined from the fourth formula of (2.6) that It can be determined from the fifth formula of (2.6) that It can be known by (2.7)-(2.11) that Since M = max{N 1 , N 2 , N 3 , N 4 , N 5 }, it can be known by the definition of the norm on R 5 that Thus, lemma is proved. Theorem 1. If M = max{N 1 , N 2 , N 3 , N 4 , N 5 }, then system (2.1) has a unique positive solution p * (s).
Proof. From the assumption M < 1, the mapping φ on Ω is a compression map. According to the compression mapping theorem, φ has a unique fixed-point (a * 1 , a * 2 , a * 3 , a * 4 , a * 5 ), which is the unique solution of (2.3). Therefore, (2.1) has a unique solution. It is known that the solution of system (2.3) is positive. The theorem is verified. It can be determined from the first formula of (1.2) that

∂[g(s)p(s, t)]
∂s Let d(s, t) = p(s, t) − p * (s), that is p(s, t) = d(s, t) + p * (s). Then, by substituting that into above formula, it can be determined that Substituting the first formula of (2.1) and set Performing a Taylor expansion of β(s, c 0 (t), R(t)) with respect to the third variable R(t), there is substituting the above formula into the fourth formula of (2.1) That is , Substituting D(s) into D 1 , it can be determined that where Substituting D(s) into D 2 , it can be determined that It can be determined from (3.7)-(3.9) that The second equation of (3.10) is multiplied by A 33 (λ) and used as input into the third equation,   There exists a nonzero solution to equation (3.12) if and only if the determinant of its coefficient is 0. From the above discussion, the following theorem can be obtained. Theorem 2. Linear systems (3.1) and (3.2) have unique formal solutions d(s, t) = D(s)e λt , if and only if λ is the root of the characteristic equation

Equilibrium stability
The following theorems can be obtained from the linear stability theory.
For the reader convenience we recall that any solution of the first-order linear differential equation system dX dt = AX with constant coefficients, it can be expressed as In order to obtain a specific discriminant condition, the death rate µ is considered to be independent of weighted size J(t) . Corollary 1. Let µ be independent of J(t), g(0) = 1, then (1) when Λ < 0, the equilibrium of (2.1) p * (s) is asymptotically stable; (2) when Λ > 0, the equilibrium of (2.1) p * (s) is unstable, where It can be determined that lim λ→+∞ K(λ) = 0 and K (λ) < 0.

Algorithms and examples
Given the degree of complexity of the characteristic equation G(λ) = 0, it is difficult to analyze the distribution of roots. We use numerical analysis to study the stability of the equilibrium The model is a class of first-order hyperbolic equations with nonlinear variable coefficients whose solutions have local dependence on the initial values. Therefore, a windward differential format, with good approximation accuracy, is chosen. The upwind differential format for system (1.2) can be summarized as follows: where a is the number of evolution time steps and τ is the time step. The scale interval is taken as [0, 1] and equally divided into N segments. The scale step size h = 1 N , in order to make the difference format stable, according to the freezing coefficient method and the Fourier analysis, the following two steps can be determined τ h max j |g j | ≤ 1.

Example 1. Population persistence
The initial distribution function of the population is assumed as In this example, the scale step size is h = 0.01, the time step is τ = 0.02, and the time iteration step is a = 120. Figure 1 is obtained from the calculation results. It can be seen that as time goes by, the density function of the population gradually moves away from the level surface, the scale distribution interval of the individual population becomes longer and longer, and the peak value increases. That is, the sustainable viability of the population is enhanced. In this example, the scale step size is h = 0.025, the time step is τ = 0.01, and the time iteration step is a = 50. Figure 2 is obtained from the calculation results. It can be seen that over time, the density of the population gradually approaches zero, that is, the population tends to become extinct.

Conclusions
Many achievements have been made in an age-structured model. In the growth process of populations, individuals of different ages can exist in the same scale stage, and scale is a better evaluation of population growth than age. The population model with scale structure proposed in this paper still has many problems to be solved, for example, the continuous variation of individual scales directly leads to mathematical models in the form of partial differential-integral equations, usually with global feedback boundary conditions, which are not a standard mathematical object and cannot be directly applied to the existing research results of mathematical theory. When the scale structure is considered, the growth rate of the population is no longer synchronized with time, which greatly increases the difficulty of analyzing the model. Through careful analysis, we obtained the existence and uniqueness condition of the equilibrium, derived the characteristic equation, gave the stability rule for judging the equilibrium, and obtained the gradual change of the equilibrium solution of the system without considering the effect of density restriction on mortality. Considering the nonlinear characteristics of the state system, it cannot be solved accurately. Therefore, an approximation algorithm is given in this paper, and the evolutionary behavior of the population is shown through two examples.
From the above stability conditions, it can be clearly seen that when the population has continuous immigration from the outside world, the stability conditions of its equilibrium are different from the situation in [30] that does not consider migration. This has some reasonableness in the actual situation. For example, for endangered species, due to the continuous immigration, the population can maintain a stable distribution under condition Λ < 0, and condition Λ < 0 also provides a reference for protecting endangered species. When µ is related to J(t), its stability analysis is difficult and needs further study. In this paper, stability conditions are used to describe the process of population development, which plays an important role in bio-diversity conservation and species management.