Dynamics of a new stage-structured population model with transient and nontransient impulsive effects in a polluted environment

*Correspondence: jiaojianjun05@126.com 1School of Mathematics and Statistics, Guizhou University of Finance and Economics, Guiyang, 550025, P.R. China Abstract In this paper, we consider a new stage-structured population model with transient and nontransient impulsive effects in a polluted environment. By using the theories of impulsive differential equations, we obtain the globally asymptotically stable condition of a population-extinction solution; we also present the permanent condition for the investigated system. The results indicate that the nontransient and transient impulsive harvesting rate play important roles in system permanence. Finally, numerical analyses are carried out to illustrate the results. Our results provide effective methods for biological resource management in a polluted environment.


Introduction
In recent years, the aggravation of environmental pollution not only affects human life style, but also poses a serious threat to the long-term survival of the species. The European Environmental Protection Agency released a new "Health and Environmental Assessment Report" on September 8, 2020, which said that 13 percent of the deaths in 28 European countries are related to environmental pollution. Therefore, the study of the population models in polluted environments is becoming more important. To date, some work has been carried out to study the population models in polluted environments [1][2][3][4][5][6][7]. Many workers have adopted a mathematical modeling approach to study the influence of environmental pollution on the survival of biological populations [8][9][10]. Most of the previous studies assumed that the input of the toxicant was continuous. The toxicants, however, are often emitted to the environment in regular pulses [11]. For example, the spraying of agricultural chemicals can be regarded as time-pulse discharge, though the discharge of the toxin is transient, the influence of the toxin will be long lasting.
Currently, the population system with a stage structure has become another focus of many studies [12][13][14][15]. Cai [16] presented a stage-structured single-species model with pulse input in a polluted environment and revealed that a long mature period of the population in a polluted environment can cause it to go extinct. Kang [17] proposed and studied an age-structured population with nonlocal diffusion. Jiao [18] investigated a stagestructured single-population model with nontransient and transient impulsive effects.

The model
In real life, when facing pollutants from the environment, a mature population and an immature population have different reactions. Considering the population with different stage structures has more practical significance. Inspired by the above discussions, we consider a new stage-structured population model with transient and nontransient impulsive effects in a polluted environment: where it is assumed that system (2.1) consists of two lakes that are connected by underground rivers. Environmental toxins will be dispersed between the two lakes due to weather conditions, such as rainy season or flood outbreaks. x(t), y(t) represent the densities of the immature and mature populations, which depend on drinking the water from lake 1, at time t, respectively. c o (t) represents the average concentration of toxins in the organism of the immature and mature populations at time t. c e1 (t) represents the concentration of environmental toxins in lake 1 at time t. c e2 (t) represents the concentration of environmental toxins in lake 2 at time t. c 1 > 0 represents the rate of immature population x turning into mature population y on (nτ , (n + l)τ ]. d 1 > 0 represents the natural death rate of population x on (nτ , (n + l)τ ]. d 2 > 0 represents the natural death rate of population y on (nτ , (n + l)τ ]. β 1 > 0 and β 2 > 0 represent the mortality coefficient of the immature population and the mature population due to the influence of toxins, respectively. f > 0 represents the uptake rate of toxin from lake 1 per unit biomass. g > 0 represents the toxin-consumption coefficient of the population by means of excretion and so on. m > 0 represents the toxin-consumption coefficient in the population by means of biochemical reactions in the body. h 1 > 0 represents the consumption coefficient of the environmental toxins with lake 1 as the water source is affected by processes such as chemical hydrolysis, volatilization, microbial degradation and photosynthesis on ((n + l)τ , (n + 1)τ ]. h 2 > 0 represents the consumption coefficient of environmental toxins with lake 2 as the water source affected by processes such as chemical hydrolysis, volatilization, microbial degradation and photosynthesis on ((n + l)τ , (n + 1)τ ]. 0 < u 1 < 1 represents the transient impulsive harvesting rate of population x at time t = (n + l)τ . 0 < u 2 < 1 represents the transient impulsive harvesting rate of population y at time t = (n + l)τ . c 2 > 0 represents the rate of immature population x turning into mature population y on ((n + l)τ , (n + 1)τ ]. d 3 > 0 represents the natural death rate of population x on ((n + l)τ , (n + 1)τ ]. d 4 > 0 represents the natural death rate of population y on ((n + l)τ , (n + 1)τ ]. The population is birth pulse as intrinsic rate of nature increase and density dependence rate of the population, which are denoted by a > 0, b > 0, respectively. E 1 > 0 represents the nontransient impulsive harvesting rate of the immature population x on ((n + l)τ , (n + 1)τ ]. E 2 > 0 represents the nontransient impulsive harvesting rate of the mature population y on ((n + l)τ , (n + 1)τ ]. 0 < d < 1 is the dispersal rate between the two lakes. v 1 > 0 represents the concentration of toxins that input into lake 1 due to environmental changes at time t = (n + 1)τ . v 2 > 0 represents the concentration of toxins that input into lake 2 due to environmental changes at time t = (n + 1)τ . τ is the period of the population-impulsive harvesting or pulse-input period of toxins.

The dynamics
Denoting U(t) = (x(t), y(t), c o (t), c e1 (t), c e2 (t)) T , the solution of system (2.1), is a piecewise continuous U : is continuous on (nτ , (n + l)τ ] and ((n + l)τ , (n + 1)τ ]. According to Ref. [19], the global existence and uniqueness of the solution of system (2.1) is guaranteed by the smoothness properties of f 1 , which denotes the mapping defined by the right-side of system (2.1). The subsystem of system (2.1) is ⎫ ⎬ ⎭ t ∈ (n + l)τ , (n + 1)τ , Considering the first and second equations and the fifth and sixth equations of system (3.1), we can obtain the analytic solution of system (3.1) between pluses as Considering the third and fourth equations and the seventh and eighth equations of system (3.1), we obtain the stroboscopic map of system (3.1) The system (3.3) has two fixed points as F 1 (0, 0) and Proof For convenience, we make a notation as (x n , y n ) = (x(nτ + ), y(nτ + )). The linear form of (3.3) can be written as Obviously, the near dynamics of F 1 (0, 0) and F 2 (x * , y * ) of (3.3) are determined by the linear system (3.5). The stabilities of the two fixed points of (3.3) are determined by the eigenvalues of J 1 less than 1. We can determine the eigenvalue of J 1 less than 1, if J 1 satisfies the Jury criteria [20] Calculating From the Jury criteria, F 1 (0, 0) is locally stable. Then, it is globally asymptotically stable.
From the Jury criteria, F 2 (x * , y * ) is locally stable. Then, it is globally asymptotically stable.
This completes the proof.
According to Theorem 1, and similar to reference [18], the following lemma can be easily proved.

Numerical simulations
Using numerical simulations, we analyze the influences of E 1 and u 1 on system (2.1). If it is assumed that x(t) = 1, y(t) = 0.

Discussion
In this paper, we propose a new stage-structured population model with impulsive effects in a polluted environment. The condition for the globally asymptotic stability of the triviality periodic solution (0, 0) of system (2.1) is obtained, and the permanent condition of system (2.1) is also obtained. It can be seen from the analyses that the nontransient harvesting rate and transient impulsive harvesting rate play important roles in system (2.1).