
==== Front
Sci Rep
Sci Rep
Scientific Reports
2045-2322
Nature Publishing Group UK London

39223196
71182
10.1038/s41598-024-71182-z
Article
Global stability and bifurcations in a mathematical model for the waste plastic management in the ocean
Parsamanesh Mahmood mparsamanesh@tvu.ac.ir

1
Izadi Mohammad 2
1 https://ror.org/00854zy02 grid.510424.6 0000 0004 7662 387X Department of Mathematics, Technical and Vocational University, Tehran, Iran
2 https://ror.org/04zn42r77 grid.412503.1 0000 0000 9826 9569 Department of Applied Mathematics, Shahid Bahonar University of Kerman, Kerman, Iran
2 9 2024
2 9 2024
2024
14 2032823 6 2024
26 8 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by/4.0/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
The use of plastic is very widespread in the world and the spread of plastic waste has also reached the oceans. Observing marine debris is a serious threat to the management system of this pollution. Because it takes years to recycle the current wastes, while their amount increases every day. The importance of mathematical models for plastic waste management is that it provides a framework for understanding the dynamics of this waste in the ocean and helps to identify effective strategies for its management. A mathematical model consisting of three compartments plastic waste, marine debris, and recycle is studied in the form of a system of ordinary differential equations. After describing the formulation of the model, some properties of the model are given. Then the equilibria of the model and the basic reproduction number are obtained by the next generation matrix method. In addition, the global stability of the model are proved at the equilibria. The bifurcations of the model and sensitivity analysis are also used for better understanding of the dynamics of the model. Finally, the numerical simulations of discussed models are given and the model is examined in several aspects. It is proven that the solutions of the system are positive if initial values are positive. It is shown that there are two equilibria E0 and E∗ and if BR<1, it is proven that E0 is globally stable, while when BR>1, the equilibrium E∗ exists and it is globally stable. Also, at BR=1 the model exhibits a forward bifurcation. The sensitivity analysis of BR concludes that the rates of waste to marine, new waste, and the recycle rate have most effect on the amount of marine debris.

Keywords

Waste plastic management
Compartmental model
Basic reproduction number
Bifurcation
Subject terms

Applied mathematics
Biological physics
issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

The use of plastic is very widespread in the world and the spread of plastic waste has also reached the oceans. Observing marine debris is a serious threat to the management system of this pollution. Because it takes years to recycle the current wastes, while their amount increases every day. The importance of mathematical models for plastic waste management is that it provides a framework for understanding the dynamics of this waste in the ocean and helps to identify effective strategies for its management. In this paper, we study the asymptotic behavior of waste plastic process management (WPM) in the ocean using a three-compartment mathematical model consisting of plastic waste, marine debris, and recycled materials. This model can be described by the following nonlinear ordinary differential equation system1.1 W′(t)=λ-γW(t)-βM(t)W(t)+μR(t),M′(t)=βM(t)W(t)-αM(t),R′(t)=γW(t)+αM(t)-μ+θR(t).

In this system, the time t (with a unit such as one day) is considered positive (t>0) and the initial values of the variables are assumed to be non-negative values as follows W(0)=W0,M(0)=M0,R(0)=R0. In this model, by W(t) we denote the amount of waste plastic material, M(t) is the marine debris, and R(t) represents the process of recycling. The descriptions of the coefficients of the WPM system are as follows. The parameter α is the marine debris rate to recycle and the waste rate to become marine denoted by β. Moreover, γ is the waste rate to be recycled directly without utilizing marine debris and λ describes the new waste rate to be reproduced. The parameter μ is the recycled waste rate to be reproduced as new waste and finally by θ we denote the recycled waste rate to be lost.

Using the properties of the Levenberg–Marquardt backpropagation (LBMBP), authors in1 designed an artificial neural network method for solving this model computationally. In2, the model was solved using two schemes based on a modification of the Morgan–Voyce (MV) functions by applying directly matrix collocation procedure and by using the quasi-linearization together with the modified MV collocation method. As mentioned, in previous studies, this model has been investigated numerically, but so far, its dynamics have not been comprehensively discussed in any work. In the present study, we intend to analyze this model with mathematical tools, although we also test the obtained results by simulation and numerical examples.

Compartmental models have been used in mathematical modeling to study the epidemics and ecological systems3–6. Also, there are some three-compartmental epidemic models to study diseases, including SIR or SIS models with vaccination or quarantine7–11, which may be useful for studying this model due to the nature of the cycle of propagation and recycling of plastic waste in the ocean. But according to our knowledge, the system studied in none of the works are completely compatible with the system related to this phenomenon in this paper.

In the present work, we first give some basic properties of the model in Section “The description of the model”. After showing that the solutions of the system are always non-negative, we obtain two equilibrium points for the system along with their existence conditions. Then, by applying the next generation method, we determine the basic reproduction number BR for the model. The investigation of the long-term behavior of the model as a dynamical system is done through the analysis of the WPM system (1) stability at the equilibrium points in Section “Stability of the model”. Also, we state and prove the conditions in terms of this quantity under which the system is stable at each equilibrium point. The bifurcations of the model and the sensitivity analysis of the parameters are discussed in Sections “Bifurcations of the model” and “Sensitivity analysis”, respectively. Numerical experiences and simulation of the solutions of the system are also carried out via some examples in Section “Numerical experiments”. Finally, the results of the paper are summarized as conclusions.

The description of the model

In this section we obtain some basic properties of the model such as the invariant set of solutions, the equilibria of the model and the basic reproduction number.

Non-negativity of solutions

The following lemma states that the solutions of system (1) remain non-negative with positive initial values.

Lemma 1

The set Ω={(W,M,R):W≥0,M≥0,R≥0} is a positive invariant for system (1).

Proof

Assume W(0)>0, M(0)>0 and R(0)>0. We letτ=sup{t>0:W(t)≥0,M(t)≥0,andR(t)≥0},

and so τ>0 since all variables are continuously differentiable. Now, if τ=+∞ then the positivity of solutions holds, but if 0<τ<+∞ then one of the variables is zero at τ and for t>τ is negative. Assume that W(τ)=0 and W(t)<0 for t>τ. From the differential equation corresponding to W we findW′(τ)=λ-γW(τ)-βM(τ)W(τ)+μR(τ)=λ+μR(τ)≥0,

Therefore W is non-decreasing at τ and W(t)≥0 for t>τ. This is a contradiction with the previous assumption for τ<+∞ and thus for all t we have W(t)≥0. We can use a similar argument for variables M and R and conclude M(0)≥0 and R(0)≥0 for all t>0. Therefore the solutions of the system (1) are always non-negative and Ω={(W,M,R):W≥0,M≥0,R≥0} is positive invariant for solutions of system (1). □

Equilibria of the model

Solving the following system of equations gives the equilibrium points (in the format of (M¯,W¯,R¯)) for model (1):2 λ-γW¯-βM¯W¯+μR¯=0βM¯W¯-αM¯=0γW¯+αM¯-μ+θR¯=0.

Two equilibrium points are obtained; the equilibriumE0=(M0,W0,R0)=0,(μ+θ)λθγ,λθ

if M¯=0, and the equilibriumE∗=(M∗,W∗,R∗)=(μ+θ)λθα-γβ,αβ,λθ.

when M¯>0.

Basic reproduction number

First, we note that M∗ in equilibrium point E∗ can be written asM∗=(μ+θ)λθα-γβ=(μ+θ)λβ-γαθαθβ.

Thus, the equilibrium point E∗ exists if and only ifM∗>0⇔γβ(βλ(μ+θ)γαθ-1)>0.

Therefore if we let BR=βλ(μ+θ)αγθ, then M∗=γβ(BR-1)>0 if and only if BR>1 and thus E∗ exists if and only if BR>1. On the other hand, considering that the components of E0 are always non-negative, we have stated the following lemma.

Lemma 2

For model (1), the equilibrium point E0 always exists, while the equilibrium point E∗ exists under the condition that BR>1.

We call BR as the basic reproduction number of the WPM model. The basic reproduction number for a compartmental model is actually defined as the number of secondary reproductions in a completely homogeneous population by a member of a compartment that causes the spread of undesirable properties in other compartments12. In the preceding discussions, this quantity was implicitly given by using the condition that the amount of marine debris (M∗) is positive at the equilibrium point E∗. Moreover, it can be also obtained by the next generation matrix method as the spectral radius of the Jacobian matrix of the model at E0 as follows13.

In this method we rewrite the second equation in (1) corresponding to the marine debris (compartment M) asM′(t)=βM(t)W(t)-αM(t)=F-V,

where F=βMW and V=αM. Here, F refers to the terms in the equation that generate new amounts of marine debris (M), and V consists of terms that represent the transmissions between compartment M and other compartments.

Now supposing f=∂F∂M(E0)=βW0=βλ(μ+θ)γθ and v=∂V∂M(E0)=α, then the basic reproduction number of the model is obtained by3 BR=ρ(fv-1)=βλ(μ+θ)αγθ.

Stability of the model

In the paper2, the authors studied the local asymptotic stability of the model. By examining the eigenvalues of the Jacobian matrix at the equilibrium points of the model and using the Routh-Hurwitz criterion14, they determined the conditions under which the eigenvalues of the Jacobian matrix have a negative real part and showed that when BR<1, the system is stable at equilibrium point E0 and if BR>1, this point is unstable and equilibrium point E∗ is stable.

Now, in the following we investigate the global stability of the model at the equilibria of the system (1).

Theorem 1

The equilibrium E0 of model (1) is globally asymptotically stable when BR≤1.

Proof

For this purpose, we use the basic concept of global stability that states under the given conditions, all solutions of the model converge to the equilibrium point E0 starting from any point15.

From third equation in system (1) we haveR(t)=e-(μ+θ)t∫0tγW(s)+αM(s)e(μ+θ)sds+R(0)=R(0)e-(μ+θ)t+∫0tγW(s)+αM(s)e-(μ+θ)(t-s)ds.

The first equation in (1) implies W′(t)≤λ-γW(t)+μR(t) and thus we getW(t)≤e-γt∫0tλ+μR(s)eγsds+W(0)=W(0)e-γt+∫0tλ+μR(s)e-γ(t-s)ds=W(0)e-γt+λ∫0te-γ(t-s)ds+μR(0)∫0te-(μ+θ)e-γ(t-s)ds+μγ∫0t∫0sW(v)e-(μ+θ)(s-v)dve-γ(t-s)ds+μα∫0t∫0sM(r)e-(μ+θ)(s-r)dre-γ(t-s)ds=W(0)e-γt+λγ1-e-γt-μR(0)μ+γ-θe-(μ+θ)t-e-γt-μγμ+γ-θ∫0tW(s)e-(μ+θ)(t-s)ds-∫0tW(s)e-γ(t-s)ds-e-γt-μαμ+γ-θ∫0tM(s)e-(μ+θ)(t-s)ds-∫0tM(s)e-γ(t-s)ds-e-γt.

Therefore we haveW(t)≤W(0)e-γt+λγ1-e-γt-μR(0)μ+γ-θe-(μ+θ)t-e-γt-μμ+γ-θ∫0tγW(s)+αM(s)e-(μ+θ)(t-s)ds-∫0tγW(s)+αM(s)e-γ(t-s)ds-(γ+α)e-γt.

Now, from the second equation in system (1) we obtainM′(t)=βM(t)W(t)-αM(t)≤βMλγ+W(0)-λγ+μR(0)μ+γ-θ+μ(γ+α)μ+γ-θe-γt-μR(0)μ+γ-θe-(μ+θ)t-μμ+γ-θ∫0tγW(s)+αM(s)e-(μ+θ)(t-s)ds+μμ+γ-θ∫0tγW(s)+αM(s)e-γ(t-s)ds-αβ.

We consider two cases; if μ+γ>θ, thenM′(t)≤βMλγ+W(0)-λγ+μR(0)μ+γ-θ+μ(γ+α)μ+γ-θe-γt+μμ+γ-θ∫0tγW(s)+αM(s)e-γ(t-s)ds-αβ,

and thus there exists τ1>0 such that for t>τ1M′(t)≤βMλγ-αβ,

since θμ+θ<1 and BR≤1 implies λγ<αβ. Therefore limt→∞M(t)=0 and in the followig limiting system4 W′(t)=λ-γW(t)+μR(t),R′(t)=γW(t)-μ+θR(t),

we find that limt→∞W(t)=λ(μ+θ)γθ=W0 and limt→∞R(t)=λθ=R0. Therefore all trajectories of (4) converge to (W0,R0) and the equilibrium E0 is globally asymptotically stable.

But, if μ+γ<θ we haveM′(t)≤βMλγ+W(0)-λγ+μR(0)μ+γ-θ+μ(γ+α)μ+γ-θe-γt-μμ+γ-θ∫0tγW(s)+αM(s)e-(μ+θ)(t-s)ds-αβ,

and with similar arguments we find again limt→∞M(t)=0 and thus the desired result is obtained. □

The next theorem states the conditions under which the equilibrium E∗ is globally stable. For this purpose we use the Lyapunov’s direct method16 which has also been employed by many authors17–21

Theorem 2

The equilibrium E∗ of model (1) is globally asymptotically stable if BR>1 and α=γ.

Proof

Using combinations of composite quadratic and common quadratic terms17, we consider function L:(M,W,R)∈Ω:M,W,R>0→R as followsL(M,W,R)=θ2(R-R∗)2+α2(W-W∗)+(M-M∗)+(R-R∗)2.

Function L is C1 on the interior of defined domain set, and we also see L≥0, the equilibrium point E∗ is the global minimum of L, and moreover L(E∗)=0. By differentiating L we getdLdt=θ(R-R∗)R′+α((W-W∗)+(M-M∗)+(R-R∗)W′+M′+R′=θ(R-R∗)γW+αM-(μ+θ)R+α(W-W∗)+(M-M∗)+(R-R∗)(λ-θR)

From system (2) we haveγW∗+αM∗-μ+θR∗=0λ=θR∗.

Thus we can writedLdt=θ(R-R∗)γ(W-W∗)+α(M-M∗)-(μ+θ)(R-R∗)+α(W-W∗)+(M-M∗)+(R-R∗)-θ(R-R∗)=(γθ-αθ)(W-W∗)(R-R∗)-(μ+θ)θ+αθ(R-R∗)2,

and thus since γ=α, we havedLdt=-(μ+θ)θ+αθ(R-R∗)2≤0.

Also, dLdt=0 if and only if R=R∗ (which implies from (1) M=M∗ and W=W∗) and(M,W,R)∈Ω:dLdt=0=E∗.

Therefore the equilibrium E∗ is globally asymptotically stable according to LaSalle’s invariant principle16. □

Bifurcations of the model

We now consider the behavior of the model when BR=1. For this purpose we choose β as the bifurcation parameter and from BR=1 we findβ∗=αγθλ(μ+θ).

To identify the bifurcations that model (1) exhibits, we apply center manifold theory22. According to this method we need to carry out the following change of model variables. Let x1=W, x2=M and x3=R by using vector X=(x1,x2,x3)T, the waste plastic management model can be rewritten in the form X′=F(X) with F=(f1,f2,f3)T as follows5 x1′=λ-γx2-βx1x2+μx3,x2′=βx1x2-αx2,x3′=γx1+αx2-μ+θx3.

The jacobian matrix is obtained asJ=-γ-βx2-βx1μβx2βx1-α0γα-(μ+θ),

and at equilibrium E0 for β∗ (when BR=1) isJβ∗0=J(E0,β∗)=-γ-αμ000γα-(μ+θ),

because β∗λ(μ+θ)γθ=α. The eigenvalues of Jβ∗0 are λ1=0 and those for the sub-matrix6 Ψ=-γμγ-(μ+θ).

For this matrix we obtaintrace(Ψ)=-(μ+θ+γ)<0anddet(Ψ)=γ(μ+θ)-μγ=γθ>0,

and thus its eigenvalues have negative real part14.

The right eigenvector u=u1,u2,u3T of matrix Jβ∗0 corresponding to the zero eigenvalue is obtained by solving the system Jβ∗0u=0 and we get u1=-α, u2=γ and u3=0. Besides, the components of the left eigenvector v=v1,v2,v3 corresponding to the zero eigenvalue can be found from vu˙=1 as are gotten as v1=0 and v2=1γ ad v3=0. It can easily be checked that vJβ∗0=0. According to the center manifold theory we have to calculate the two following coefficientsa=∑k,i,j=13vkuiuj∂2fk∂xi∂xj,b=∑k,i=13vkui∂2fk∂xi∂β.

All second order partial derivatives of fk,(k=1,2,3) with respect to xi,(i=1,2,3) and β are calculated at (E0,β∗) with E0=(x1∗,x2∗,x3∗)=λ(μ+θ)γθ,0,λθ. Since v1, v3 and u3 are zero, we only need to calculate the following second order partial derivatives:∂2f2∂x2∂x1(E0,β∗)=∂∂x2βx2|(E0,β∗)=β∗,∂2f1∂x1∂β(E0,β∗)=∂∂x1x1x2|(E0,β∗)=0,∂2f1∂x2∂β(E0,β∗)=∂∂x2x1x2|(E0,β∗)=λ(μ+θ)γθ.

Thusa=v2u1u2∂2f1∂x1∂x2(E0,β∗)=-αβ∗<0,b=v2u1∂2f1∂x1∂β+v2u2∂2f1∂x2∂β=0+λ(μ+θ)γθ>0.

Therefore, according to Theorem 4.1 in22 since a<0 and b>0 we conclude that there is a forward bifurcation (transcritical bifurcation) at β∗ (when BR=1) for the waste plastic management model.

Sensitivity analysis

To find out how sensitive a function is to changes in the variables in its formula, the method of sensitivity analysis is used. This method uses a quantity called the normalized forward sensitivity index as a measure of the sensitivity of each variable. Since the basic reproduction number has a significant impact on the behavior of the model, we calculate the sensitivity indices of BR for its variables, which are β, λ, μ, α, γ, and θ.

The normalized forward sensitivity index of variable BR for a variable v, is defined by SIv0=vBR×∂BR∂v23,24.

If SIvBR>0 then the variable v has positive impact on BR and the value of variable BR increases by increasing the value of v. While SIvBR<0 shows the reverse impact of v on BR; increasing the value of v implies decreasing in the value of R0. Also, the magnitude of SIvBR shows the proportion of changes in BR with respect to v.

The normalized forward sensitivity indices for BR are caculating as follows:SIβBR=1,SIλBR=1,SIμBR=μμ+θ<1,SIαBR=-1,SIγBR=-1,SIθBR=-μμ+θ<0.

We see that SIvBR>0 for v=β,λ and μ, but SIvBR<0 for v=α,γ and θ. Thus, any increase (or decrease) in values of β,λ and μ has direct impact on the value of R0. However, any increase (or decrease) in values of v=α,γ and θ has reverse impact on the value of BR.

Numerical experiments

We consider the values β=0.15,γ=0.41,α=0.65,μ=0.4,λ=0.36,θ=0.15 for the parameters in model (1). For these values we have BR=0.743<1 and according to Theorem 1 the equilibrium E0=(0,3.2195,2.4000) is stable. Figure 1 illustrates the solutions of the model with these parameter values and 20 different initial values for sub-populations. Moreover, the left picture in Fig. 3 shows the solutions of the system with initial values as W0=1.5, M0=2, and R0=1.Fig. 1 Solutions of model for different initial values for sub-populations. The parameter values are β=0.15,γ=0.41,α=0.65,μ=0.4,λ=0.36,θ=0.15 which imply BR=0.743<1.

Now, we change the parameter values to β=0.4,γ=0.21,α=0.5,μ=0.4,λ=0.66,θ=0.2. Here, we have BR=7.5429>1 and by Theorem 2 the equilibrium E∗=(1.25,3.44,3.30) is stable. In Fig. 2 the solutions of the system (1) has been depicted for 20 initial values for sub-populations. In the case of that the initial values are supposed as W0=1.5, M0=2, and R0=1, the right picture in Fig. 3 shows the solutions of the model.Fig. 2 Solutions of model for different initial values for sub-populations. The parameter values are β=0.4,γ=0.21,α=0.5,μ=0.4,λ=0.66,θ=0.2 which imply BR=7.5429>1.

Fig. 3 Solutions of the model for two cases BR<1 (the left picture) and BR>1 (the right picture) with their theoretical solutions when initial values are same.

For parameter values γ=0.41,α=0.65,μ=0.4,λ=0.36,θ=0.15 (as in the first example), if we solve the equation BR=1 for parameter β as the bifurcation parameter, we get β∗=0.2019. Thus for this value the model exhibit a forward bifurcation as it can be seen in bifurcation diagram presented in Fig. 4.Fig. 4 Bifurcation diagram of the model for parameter β∈[0,2].

Now, we investigate the impact of the parameters of the model on the dynamics of the waste management system by using sensitivity analysis of the basic reproduction number BR with respect to each parameter as it was explained in Section “Sensitivity analysis”. Let us consider the parameter values in Table 1 for parameters in the model introduced by system (1) and their corresponding sensitivity indices with respect to BR as a differentiable function.Table 1 Parameters in the waste management system (1) and their sensitivity index.

Parameter	Value	Sensitivity Index	
β	0.75	+1	
λ	0.66	+1	
μ	0.4	+0.8889	
α	0.5	-1	
γ	0.21	-1	
θ	0.05	-0.8889	

Fig. 5 Normalized sensitivity indices for model parameters.

The normalized sensitivity indices for parameters have been shown as a chart in Fig. 5. According to the Table 1 we find that the sensitivity indices for parameters β, λ and μ have positive value and their impact on BR is direct. While, the parameters α, γ and θ have negative sensitivity indices and thus they have reverse impact on BR. Thus for example, an increase (decrease) in values of λ and μ by %10 yields to a %10 and %8.889 increase (decrease) in BR, respectively. For instance, for parameter values in Table 1 we have BR=42.4286 and if we decrease λ by %10, the basic reproduction number also decreases by %10 and becomes BR=38.1857. On the other hands, a %10 increase (decrease) for example in γ and θ causes a %10 and %8.889 decrease (increase) in BR. For example, if the value of θ is increased by %20 (to θ=0.06), then the value of BR decreases by %17.778 and becomes BR=34.8763. Therefore, according to Table 1 we find that either decreasing the rates of waste to marine (β) and new waste (λ) or increasing the recycle rate (α and γ), have the most impact in reducing the value of BR and as a result most impact on controlling the amount of marine debris. The impact of parameters β and α on the model have also been depicted in Figs. 6 and 7, respectively. The parameter values are assumed as β=0.75,γ=0.21,α=0.5,μ=0.4,λ=0.66,θ=0.05 and initial values are W0=1.5, M0=2, and R0=1. In Fig. 6 the values of β change in interval [0.1, 2.1] and it is seen that by increasing β ,the final solution corresponding to the amount of the marine debris (M) leaves the zero and will take positive values. Indeed, the stability of the model changes from E0 to E∗. With the same terms and for values of α in [0.1, 0.9], the solutions of the model have been shown in Fig. 7. By increasing the value of α, we observe that the values of M finally take zero value. This shows that by increasing α the equilibrium point E∗ becomes unstable and E0 becomes E∗ stable.Fig. 6 Solutions of the model for γ=0.21,α=0.5,μ=0.4,λ=0.66,θ=0.05 and β∈[0.1,2.1].

Fig. 7 Solutions of the model for β=0.75,γ=0.21,μ=0.4,λ=0.66,θ=0.05 and α∈[0.1,0.9].

Conclusions

In this paper, we studied the waste plastic management (WPM) system in the ocean through a mathematical three-compartmental model. The basic reproduction number BR, and two equilibria of the model were found, in addition to positivity of solutions of model. The dynamics of the model was determined in terms of threshold BR; if BR<1, it was proved that the equilibrium E0 is globally stable, while the equilibrium E∗ exists and it is stable when BR>1. Also, it was shown that the model exhibit a forward (transcritical) bifurcation at BR=1. The sensitivity of the model has been analyzed by calculating normalized forward sensitivity index for each parameter for BR and it was concluded that decreasing the rates of waste to marine (β) and new waste (λ) or increasing the recycle rate (α and γ), are most effective for controlling the amount of marine debris. Finally, the theoretical results were discussed also numerically for different parameter values and various initial values for sub-populations via several examples, solutions of model and bifurcation diagram.

The global stability of the equilibrium point E∗ has been proved for the case that the marine debris recycling rate (α) and the direct recycling rate of waste materials (γ) are equal. Constructing a more appropriate Lyapunov function that does not impose such an additional condition on stability can be the subject of future studies. Investigating the impact of seasonality on the system behavior may also complement the present study, since the waste rate to become marine (β) can be considered as a periodic function.

Author contributions

M.P.: Conceptualization, Writing- Original draft preparation, Software, Investigation, Supervision. M.I.: Visualization, Investigation, Writing- Reviewing and Editing.

Funding

The authors declare that no funds, grants, or other support were received during the preparation of this manuscript.

Data availability

All data generated or analyzed during this study are included in this article.

Competing interests

The authors declare no competing interests.

Publisher's note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

1. Al Nuwairan M Sabir Z Asif Zahoor Raja M Aldhafeeri A An advance artificial neural network scheme to examine the waste plastic management in the ocean AIP Adv. 2022 12 4 045211 10.1063/5.0085737
Al Nuwairan, M., Sabir, Z., Asif Zahoor Raja, M. & Aldhafeeri, A. An advance artificial neural network scheme to examine the waste plastic management in the ocean. AIP Adv. 12(4), 045211 (2022).10.1063/5.0085737
2. Izadi M Parsamanesh M Adel W Numerical and stability investigations of the waste plastic management model in the ocean system Mathematics 2022 10 23 4601 10.3390/math10234601
Izadi, M., Parsamanesh, M. & Adel, W. Numerical and stability investigations of the waste plastic management model in the ocean system. Mathematics 10(23), 4601 (2022).10.3390/math10234601
3. Thirthar AA Jawad Sh Shah K Abdeljawad Th How does media coverage affect a COVID-19 pandemic model with direct and indirect transmission? J. Math. Comput. Sci. 2024 35 2 169 181 10.22436/jmcs.035.02.04
Thirthar, A. A., Jawad, Sh., Shah, K. & Abdeljawad, Th. How does media coverage affect a COVID-19 pandemic model with direct and indirect transmission?. J. Math. Comput. Sci. 35(2), 169–181 (2024).10.22436/jmcs.035.02.04
4. Mondal B Thirthar AA Sk N Alqudah MA Abdeljawad T Complex dynamics in a two species system with Crowley–Martin response function: Role of cooperation, additional food and seasonal perturbations Math. Comput. Simul. 2024 221 415 434 10.1016/j.matcom.2024.03.015
Mondal, B., Thirthar, A. A., Sk, N., Alqudah, M. A. & Abdeljawad, T. Complex dynamics in a two species system with Crowley–Martin response function: Role of cooperation, additional food and seasonal perturbations. Math. Comput. Simul. 221, 415–434 (2024).10.1016/j.matcom.2024.03.015
5. Thirthar AA Jawad S Majeed SJ Nisar KS Impact of wind flow and global warming in the dynamics of prey-predator model Results Control Optim. 2024 15 100424 10.1016/j.rico.2024.100424
Thirthar, A. A., Jawad, S., Majeed, S. J. & Nisar, K. S. Impact of wind flow and global warming in the dynamics of prey-predator model. Results Control Optim. 15, 100424 (2024).10.1016/j.rico.2024.100424
6. Thirthar AA A mathematical modelling of a plant-herbivore community with additional effects of food on the environment Iraqi J. Sci 2023 64 7 3551 5366
Thirthar, A. A. A mathematical modelling of a plant-herbivore community with additional effects of food on the environment. Iraqi J. Sci 64(7), 3551–5366 (2023).
7. Lu Z Chi X Chen L The effect of constant and pulse vaccination on SIR epidemic model with horizontal and vertical transmission Math. Comput. Model. 2002 36 9–10 1039 1057 10.1016/S0895-7177(02)00257-1
Lu, Z., Chi, X. & Chen, L. The effect of constant and pulse vaccination on SIR epidemic model with horizontal and vertical transmission. Math. Comput. Model. 36(9–10), 1039–1057 (2002).10.1016/S0895-7177(02)00257-1
8. Brauer F Backward bifurcations in simple vaccination models J. Math. Anal. Appl. 2004 298 2 418 431 10.1016/j.jmaa.2004.05.045
Brauer, F. Backward bifurcations in simple vaccination models. J. Math. Anal. Appl. 298(2), 418–431 (2004).10.1016/j.jmaa.2004.05.045
9. Moneim I Greenhalgh D Threshold and stability results for an SIRS epidemic model with a general periodic vaccination strategy J. Biol. Syst. 2005 13 02 131 150 10.1142/S0218339005001446
Moneim, I. & Greenhalgh, D. Threshold and stability results for an SIRS epidemic model with a general periodic vaccination strategy. J. Biol. Syst. 13(02), 131–150 (2005).10.1142/S0218339005001446
10. Parsamanesh M Global dynamics of an sivs epidemic model with bilinear incidence rate Ital. J. Pure Appl. Math. 2018 40 544 557
Parsamanesh, M. Global dynamics of an sivs epidemic model with bilinear incidence rate. Ital. J. Pure Appl. Math. 40, 544–557 (2018).
11. Parsamanesh M Farnoosh R On the global stability of the endemic state in an epidemic model with vaccination Math. Sci. 2018 12 313 320 10.1007/s40096-018-0271-3
Parsamanesh, M. & Farnoosh, R. On the global stability of the endemic state in an epidemic model with vaccination. Math. Sci. 12, 313–320 (2018).10.1007/s40096-018-0271-3
12. Hethcote HW The mathematics of infectious diseases SIAM Rev. 2000 42 4 599 653 10.1137/S0036144500371907
Hethcote, H. W. The mathematics of infectious diseases. SIAM Rev. 42(4), 599–653 (2000).10.1137/S0036144500371907
13. Van den Driessche P Watmough J Reproduction numbers and subthreshold endemic equilibria for compartmental models of disease transmission Math. Biosci. 2002 180 1–2 29 48 10.1016/S0025-5564(02)00108-6 12387915
Van den Driessche, P. & Watmough, J. Reproduction numbers and subthreshold endemic equilibria for compartmental models of disease transmission. Math. Biosci. 180(1–2), 29–48 (2002).12387915 10.1016/S0025-5564(02)00108-6
14. Ortega JM Matrix theory: A second course 2013 Springer
Ortega, J. M. Matrix theory: A second course (Springer, 2013).
15. Shub M Global stability of dynamical systems 2013 Springer
Shub, M. Global stability of dynamical systems (Springer, 2013).
16. La Salle, J. & Lefschetz, S. Stability by Liapunov’s direct method with applications (Elsevier, 2012).
17. Vargas-De-León C Constructions of Lyapunov functions for classic SIS, SIR and SIRS epidemic models with variable population size Foro-Red-Mat. Revista Electrónica de Contenido Matemático 2009 26 5 1 12
Vargas-De-León, C. Constructions of Lyapunov functions for classic SIS, SIR and SIRS epidemic models with variable population size. Foro-Red-Mat. Revista Electrónica de Contenido Matemático 26(5), 1–12 (2009).
18. Vargas-De-Leon C On the global stability of SIS, SIR and SIRS epidemic models with standard incidence Chaos Solitons Fractals 2011 44 12 1106 1110 10.1016/j.chaos.2011.09.002
Vargas-De-Leon, C. On the global stability of SIS, SIR and SIRS epidemic models with standard incidence. Chaos Solitons Fractals 44(12), 1106–1110 (2011).10.1016/j.chaos.2011.09.002
19. Parsamanesh M Global stability analysis of a veisv model for network worm attack Univ. Politeh. Buchar. Sci. Bull.-Ser. A-Appl. Math. Phys. 2017 79 4 179 188
Parsamanesh, M. Global stability analysis of a veisv model for network worm attack. Univ. Politeh. Buchar. Sci. Bull.-Ser. A-Appl. Math. Phys. 79(4), 179–188 (2017).
20. Parsamanesh M Erfanian M Global dynamics of an epidemic model with standard incidence rate and vaccination strategy Chaos Solitons Fractals 2018 117 192 199 10.1016/j.chaos.2018.10.022
Parsamanesh, M. & Erfanian, M. Global dynamics of an epidemic model with standard incidence rate and vaccination strategy. Chaos Solitons Fractals 117, 192–199 (2018).10.1016/j.chaos.2018.10.022
21. Cangiotti N Capolli M Sensi M Sottil S A survey on lyapunov functions for epidemic compartmental models Boll. dell’Unione Mat. Ital. 2023 17 241 257 10.1007/s40574-023-00368-6
Cangiotti, N., Capolli, M., Sensi, M. & Sottil, S. A survey on lyapunov functions for epidemic compartmental models. Boll. dell’Unione Mat. Ital. 17, 241–257 (2023).10.1007/s40574-023-00368-6
22. Castillo-Chavez C Song B Dynamical models of tuberculosis and their applications Math. Biosci. Eng. 2004 1 2 361 404 10.3934/mbe.2004.1.361 20369977
Castillo-Chavez, C. & Song, B. Dynamical models of tuberculosis and their applications. Math. Biosci. Eng. 1(2), 361–404 (2004).20369977 10.3934/mbe.2004.1.361
23. Chitnis N Hyman JM Cushing JM Determining important parameters in the spread of malaria through the sensitivity analysis of a mathematical model Bull. Math. Biol. 2008 70 5 1272 10.1007/s11538-008-9299-0 18293044
Chitnis, N., Hyman, J. M. & Cushing, J. M. Determining important parameters in the spread of malaria through the sensitivity analysis of a mathematical model. Bull. Math. Biol. 70(5), 1272 (2008).18293044 10.1007/s11538-008-9299-0
24. Abidemi A Abd Aziz M Ahmad R Vaccination and vector control effect on dengue virus transmission dynamics: Modelling and simulation Chaos Solitons Fractals 2020 133 109648 10.1016/j.chaos.2020.109648
Abidemi, A., Abd Aziz, M. & Ahmad, R. Vaccination and vector control effect on dengue virus transmission dynamics: Modelling and simulation. Chaos Solitons Fractals 133, 109648 (2020).10.1016/j.chaos.2020.109648
