Effects of greenhouse gases and hypoxia on the population of aquatic species: a fractional mathematical model.

Study of ecosystems has always been an interesting topic in the view of real-world dynamics. In this paper, we propose a fractional-order nonlinear mathematical model to describe the prelude of deteriorating quality of water cause of greenhouse gases on the population of aquatic animals. In the proposed system, we recall that greenhouse gases raise the temperature of water, and because of this reason, the dissolved oxygen level goes down, and also the rate of circulation of disintegrated oxygen by the aquatic animals rises, which causes a decrement in the density of aquatic species. We use a generalized form of the Caputo fractional derivative to describe the dynamics of the proposed problem. We also investigate equilibrium points of the given fractional-order model and discuss the asymptotic stability of the equilibria of the proposed autonomous model. We recall some important results to prove the existence of a unique solution of the model. For finding the numerical solution of the established fractional-order system, we apply a generalized predictor-corrector technique in the sense of proposed derivative and also justify the stability of the method. To express the novelty of the simulated results, we perform a number of graphs at various fractional-order cases. The given study is fully novel and useful for understanding the proposed real-world phenomena.

Introduction

In the study of greenhouse effects, we know that in the day the sun warms up the atmosphere of earth. But when the Earth supercools at the night, then the presented heat is radiated again into the environment. In the duration of this process, the heat is exploited by the greenhouse gases in the environment of earth. This process makes the layer of the earth thermal, which causes the possibility of living being’s survival on earth. However, because of the increment in the level of greenhouse gases, the earth’s temperature has raised simultaneously. This has caused a number of drastic impacts. In the list of reasons of greenhouse effect, deforestation, burning of fossil fuels, farming, industrial waste, and landfills play a major role. The major effects of increased greenhouse gases are depletion of ozone layer, global warming, air and smog pollution, water bodies acidification, etc. Since the starting of the industrial revolution, the concentration of carbon dioxide, chlorofluorocarbon (CFC), nitrous oxide, and methane have enhanced in the environment, and there is firm witness that the venomous impacts of greenhouse gases on our ecological systems have been taken account as a outcome of human bustles.

Aquatic life simply means to stay in surface water, and water in this paragon is specified as a marine habitat. Living beings that live in the water either permanently or momentarily are called aquatic animals and plants, and these compose the beings in water aquatic life. It is well known that the increment in the temperature of water causes the reduction in the concentration of mixed oxygen of the aquatic environment and also rises the requirement of mixed oxygen for the aquatic animals. Invertebrates, fish, and other aquatic species rely upon the amount of oxygen decomposed in the water, and in the absence of it, they may not live. A small changes in concentration of mixed oxygen can effect the conformation of aquatic society [1]. So the rate of survival of the aquatic density (Fig. 1) goes down under hypoxia, and the oxygen necessary for their living raises with growth in temperature [2, 3]. Hence, because of the combined influences of reduced concentration of mixed oxygen and enhanced demand of oxygen by the animals, the warming of species bodies rises the death rate of species population [3, 4]. To define these dynamics, a number of models have been proposed, but only few models [5, 6] have been given to simulate the effects of dissolved oxygen and temperature on the population of aquatic species.

Figure 1

Some aquatic animals

In the matter of the above discussion, in this study, we prepare a fractional-order mathematical system to simulate the joint influences of low mixed oxygen density, exalted water temperature, and raised oxygen demand on the extinction or survival of aquatic population. Fractional derivatives are one of the most effective tools for simulations and have been proposed in many different ways (for instance, see [7-9]). Fractional-order models have been widely used to define a number of real-world problems because their memory effects make these models more visible in the literature. Recently, a number of fractional-order models have been prepared by researchers. In this regard, in [10-18] the authors have proposed a number of fractional-order mathematical models to describe the dynamics of Covid-19 epidemic. In [19, 20] the authors have simulated the fractional-order dynamics of well-known lassa hemorrhagic fever. The applications of fractional derivatives in ecology can be seen in [21]. Regarding some more specific areas, nonclassical derivatives have been successfully used to derive the structure of tuberculosis [22], malaria [23], mosaic disease [24], Nipah epidemic [25], canine distemper virus [26], and huanglongbing transmission [27]. In [28] the authors used a fractional-order time-delay mathematical model to describe the process of oncolytic virotherapy. A study on analytic solution for oxygen diffusion from capillary to tissues via fractional derivatives is proposed in [29]. Also, an application of a new generalized Caputo derivative to define the famous love story of Layla and Majnun is given in [30]. So the literature of fractional-order calculus is increasing exponentially day by day. Also, a number of true and false results come on the various fractional derivatives. Recently, in [31] the authors have proved that in the case of evolution equations in terms of the Caputo–Fabrizio and Atangana–Baleanu fractional derivatives, intrinsic discontinuities occur. The geometry of fractional-order derivatives is still not well-defined, but their applications in different scientific fields make them more visible to the literature. Some important studies related to the properties of fractional derivatives, special functions, and different types of inequalities can be learned from [32-35]. Nonstandard Chebyshev collocation and finite difference schemes for solving fractional diffusion equations are proposed in [36]. Some novel analysis on the fractional differential equations for the generalized Mittag-Leffler function are discussed in [37]. A study on the analytical solutions of the fractional-order equations with uncertainty is proposed in [38]. Alderremy et al. [39] have discussed some novel models of the multispace-fractional Gardner equation. A study on spectral collocation method for solving smoking model is proposed in [25]. In [40] the authors have proposed a study on Darcy–Brinkman–Forchheimer model for nanobioconvection stratified MHD flow through an elastic surface. In [41] a reduced differential transform scheme for simulating nonlinear biomathematics models is given. In [42] a study on dynamical features and signal flow graph of nonlinear noninteger order smoking mathematical model is explored. In [43], some numerical methods for a model of relativistic electrons arising in the laser thermonuclear fusion are investigated. The manuscript is designed as follows: In Sect. 2, firstly, we remind some important definitions and results. In Sect. 3, we give a complete description of the proposed fractional-order nonlinear model, where we define the significance and importance of every small part of the model. Then in Sect. 4, we give a complete mathematical analysis related to the solution existence, derivation, and stability. To show the correctness of our results, in Sect. 5, we present the necessary graphs at various fractional-order values and parameter weights. At the end, a conclusion gives a comfortable end to the paper.

Preliminaries

Firstly, we remind some important definitions and results.

Definition 1

([44])

The new definition of the Caputo-type fractional derivative of order (called a new generalized Caputo) for the function is given by where , , and .

Lemma 1

([45])

For and a nonnegative integer ϱ, there exist positive constants and , dependent only on b, such that and

Lemma 2

([45])

Let for and for , let , and , where r is a positive integer. Let for . If then where is a positive constant not dependent on r and h.

Model dynamics

Now we propose a fractional-order mathematical model to study the proposed dynamics. In [1] the authors have already given an idea on the proposed topic by using an integer-order model. We propose a fractional-order model because it is well known that the memory effects, which cannot be studied in the classical case, can be easily observed by fractional-order derivatives. It is very important that when we propose a fractional-order model, it should have the same time dimension on both sides of the system. Taking care of all these aspects, we define the novel fractional-order model as follows: where is the new generalized Caputo-type fractional-order operator of order σ. In this model, and In this model, we have five different classes, in which class N shows the logistically crescent aquatic species density whose rate of growth is taken as a function of temperature and mixed oxygen, U justifies the dissolved oxygen concentrations, T defines the water temperature average of the species, C expresses the greenhouse gases accumulative concentrations, and Z justifies the concentration of ozone. Also, the term expresses the specific rate of growth of the species, which is in fact an exponentially decreasing function of T for and increasing function of U. The function denotes the quantity of dissolved oxygen demanded by species population, which rises with temperature increase.

The term is defined for the consideration that if the water temperature level is high, close to the optimum temperature, then the natural loaded dissolved oxygen concentration reduces. The significance of all other parameters is completely given in Table 1. The more deep texture of the given model in classical sense can be learned from [1].

Table 1

Description of model parameters

g_0	Intrinsic growth rate
β_10	Dissolved oxygen’s minimum natural concentration needed by the aquatic species
β_11	Increment rate in the mass of mixed oxygen demanded for the species per unit rise in the level of temperature above the suitable temperature
T_opt	Optimal water temperature for the aquatic species maximum rate of growth
γ_0	Carrying capacity of the environment
D_s_0	Dissolved oxygen’s natural saturated concentration at T=T_opt
A_0	Ejection rate of greenhouse gases cause of anthropogenic bustles
w	Increment rate in the temperature of water cause of greenhouse gases
ζ_1	Coefficient of heat transfer of surface
O_c	Physic manufacture of concentration of ozone per unit time in the environment
Λ_1	Natural deterioration rate of concentration of ozone
Λ	Deterioration rate of concentration of ozone cause of greenhouse gases
γ_1	Coefficient of reaeration at the reference temperature
δ_2	Deterioration rate of mixed oxygen because of breathing by the species
ζ	Deterioration rate of mixed oxygen because of a rise in the temperature above the suitable temperature
β	A constant that succumbs upon the tincture state of the water body
γ	Variations rate in the water temperature because of changes in the ozone concentration level associated with its threshold value
Z_0	Threshold of concentration of ozone below which temperature will rise
T_10	Temperature of the environment
δ_1	Depletion rate of greenhouse gases
T_0	Context temperature (associated with the turbulence degree in the water, in which turn succumbs on the depth and speed of the river)
b	Constant that incarnates the toxic influence of divergence of T from T_opt and divergence of T from T_max
T_max	Maximum temperature of water at which growth can occur
N(0)	Initial population of N
T(0)	Initial population of T
C(0)	Initial population of C
Z(0)	Initial population of Z
U(0)	Initial population of U

The equilibria of the given fractional-order mathematical model can be obtained by solving the following system: Equation (5) gives Equation (6) gives Equation (4) gives Here we have two different types of equilibrium points.

1. Boundary equilibrium point :

(no species population), . Here are given by (8), (9), (10), respectively.

A boundary equilibrium point Ē exists if and .

2. Interior equilibrium point , where (species population exists) and , provided that . Here are given by Equations (8), (9), (10), respectively, and is the positive root of the quadratic equation where , , . When the given conditions are taken account, the quadratic equation (11) has at least one positive root if , and . Now we derive the following nonautonomous system after solving the given model (2) for C: since .

Hence the fractional-order nonautonomous model (12)–(13) can be specified in the following equivalent fractional-order autonomous model: The equilibrium points of the dynamic system (14)–(15) are calculated by the following group of equations: 1. Boundary equilibrium point : The existence of boundary equilibrium point provides 2. Interior equilibrium point :

(aquatic population exists) and , provided that where is a positive root of the quadratic equation where , , .

If the given constraints are satisfied, then the quadratic equation specified by (18) has at least one positive root if

Lemma 3

For the fractional-order mathematical system (14)–(15), is locally asymptotically stable if and is an unstable saddle point if .

Proof

After the linearization, taking the Laplace transform of both sides of system (14)–(15), the Jacobian matrix for system (14)–(15) simulated at is given by where . The eigenvalues associated with the matrix are .

The eigenvalue is negative if and positive if .

The other eigenvalue is negative. Hence the required results are obtained. □

Lemma 4

The given equilibrium point of the fractional-order system (14)–(15) is always locally asymptotically stable.

The Jacobian matrix of system (14)–(15) with respect to is where , .

The behavior of the eigenvalues is estimated by using Hurwitz’s criteria in the quadratic equation Using Hurwitz’s criteria, we observe that the eigenvalues of the matrix are negative if . Thus we get that is locally asymptotically stable under the restriction . □

Now for the deformation of fractional-order system (2), we convert it to an equivalent compact form in the case of singular kernels as follows: Here are the proposed kernels with respect to the given classes , respectively.

Fractional-order analysis on the proposed model

Analysis of the existence and uniqueness of the solution

Proving the existence of the solution for fractional-order systems is always a sensitive part because not all fractional differential equations have their proof of the existence of a solution. In this area a number of works have been done, and lots of researchers work. Here, before deriving the solution of the proposed model, we first prove that the given fractional-order model has a unique solution. We give the results only for the class , and the results are as for the other model classes. So we recall the model equation for N, and the relative Volterra integral equation

Theorem 1

([46] (Existence))

Let , and . Define , and let the mapping be continuous. Further, define and Then the IVP (21a)–(21b) has a solution .

Lemma 5

([46])

By considering the result of Theorem 1a function solves the IVP (21a)–(21b) if and only if it solves the Volterra integral equation (22).

Theorem 2

([46] (Uniqueness))

Let , and . Further, let and . For the set given in Theorem 1, let be a continuous function that satisfies the Lipschitz condition with respect to the second variable, that is, with a constant independent of , and . Then the IVP (21a)–(21b) has a unique solution .

Numerical solution of the proposed model with application of the generalized predictor–corrector technique

In the last few years, a number of fractional-order numerical schemes have been proposed by the scientists to solve various types of dynamical models. Very recently, the authors of [47] have proposed a new numerical method in the generalized Caputo derivative sense. Here we solve the proposed model with the help of generalized P-C scheme for the solution of the IVP (21a)–(21b) by following the methodology proposed in [44]. Also, we will analyze the stability of the given scheme. In that way, we first recall the above given Volterra integral equation (22), which gives Now with supposing that a unique solution exists for the function on the interval , we divide the adopted interval into N unequal subparts using the mesh points where . Now let us try to analyze the approximations , to get a numerical solution of the given IVP. Suppose that we have already derived the approximations and want to derive approximations by means of the integral equation By substitution we get that is, Now, to simulate the right-side of Eq. (28), applying the trapezoidal quadrature rule with respect to the weight function and shifting the function by its piecewise linear interpolant with nodes , we get So, fitting the above-proposed approximations in Eq. (28), we establish the corrector term for : where The final task for our solution is changing the quantity on the right-hand side of formula (30) with the predictor value , which can be calculated by applying the one-step Adams–Bashforth technique to the integral equation (27). In this case, by changing the mapping by the quantity at each integral in Eq. (28) we get So our P-C method for deriving the approximations is totally evaluated by the formula where , and the predicted value can be simulated as mentioned in Eq. (32) with the terms estimated according to (31).

Therefore the derivation for the approximate solution of the proposed system (2) is derived successfully and defined by the following equations: where

Method stability

Theorem 3

Let satisfy the Lipschitz condition, and let be approximate solutions of the derived P-C method (34) and (35), respectively. Then the proposed numerical algorithm (34)–(35) is conditionally stable.

Proof

Let , and be perturbations of , and , respectively. Then the given below perturbation equations are estimated with the help of Eqs. (34) and (35). where , Using the Lipschitz condition, we obtain where . Also, from Eq. (3.18) in [45] we derive where . Substituting from Eq. (39) into Eq. (38) results in where , is a positive constant depending on σ (by Lemma 1), and h is assumed to be small enough. Using Lemma 2, we have . which concludes the proof. □

Experimental simulations

After finishing all necessary theoretical analysis, we start to perform some experimental calculations to show the correctness of our results. We use Mathematica software for performing the number of graphs. For the case of interior equilibrium point , we use the following set of parameter values: . Here we observe that for the case of fractional order (when the model behaves like an integer-order system), the authors of [1] have calculated the value of the interior equilibrium point and then, in this case, have specified the constraints for the solution boundedness, equilibrium point stability, and positivity of the solution. Our target is to explore the dynamics of all model classes with respect to the interior equilibrium points at different fractional-order values σ.

In the set of Fig. 2, we observed the nature of all model classes separately at different fractional-order values σ. In subfigure 2(a) the dynamics of density of aquatic population N is plotted at . Here we observed that at the numerically calculated equilibrium point is satisfied for class N and also at other values of order σ, it changed simultaneously. Similarly, subfigure 2(b) shows the average water temperature of the species (class T), subfigure 2(c) shows the concentration of greenhouse gases (class C), subfigure 2(d) shows the ozone concentration (class Z), and subfigure 2(e) shows the dynamics of dissolved oxygen concentration (class U). The simultaneous changes in the given model classes at particular values of σ can be seen from the set of Fig. 3. Overall, we observed that when the fractional-order σ changes, the dynamics of the model, along with interior equilibrium point changes, justifies the importance of the fractional-order model.

Figure 2

Separate plots of all model classes at various fractional-order values σ for the case of interior equilibrium point

Figure 3

Mixed plots of all model classes at fractional-order values for the case of interior equilibrium point

As investigated above, now we consider the case of boundary equilibrium point . In this case, we consider the following parameter values: . For the given parameter weights, the value of boundary equilibrium point at fractional-order (when the model behaves like an integer-order system given in [1]) is . In that integer-order case, the boundary equilibrium point is linearly asymptotically stable.

For the noninteger-order observations, in the set of Fig. 4, we analyzed the nature of proposed model classes separately at various fractional-order values σ. In subfigure 4(a), the dynamics of density of aquatic population N is plotted at . Here we can see that for , the numerically calculated equilibrium point is satisfied for population N and that at other values of order σ, it changes simultaneously. Following the same way, subfigure 4(b) specifies the average water temperature of the species (class T), subfigure 4(c) demonstrates the concentration of greenhouse gases (class C), subfigure 4(d) shows the ozone concentration (class Z), and subfigure 4(e) shows the dynamics of dissolved oxygen concentration (class U).

Figure 4

Separate plots of all model classes at various fractional-order values σ for the case of boundary equilibrium point

The simultaneous changes in the given model classes at particular value of σ can be analyzed from the set of Fig. 5. Overall, we can see that when the fractional-order σ changes, the dynamics of the model changes along with boundary equilibrium point, which satisfies the role of fractional-order operator.

Figure 5

Mixed plots of all model classes at fractional-order values for the case of boundary equilibrium point

From the above given experimental analysis we see that the fractional-order dynamics with memory effects is much stronger than the integer-order dynamics. Here we have more varieties to understand the structure of the proposed ecosystem dynamics at various fractional-order values along with different values of equilibrium points. The modified Caputo fractional derivative is fully suitable to simulate the novel results with the help of given fractional-order model.

Conclusion

In our study, we have simulated a novel fractional-order mathematical system to study the prelude of deteriorating quality of water because of greenhouse gases on the population of aquatic animals. It has been shown in the given system that greenhouse gases raise the temperature of water, and because of this reason, the dissolved oxygen level goes down, and also the rate of circulation of disintegrated oxygen by the species rises, which causes a decrement in the density of aquatic species. We have used a new generalized Caputo-type fractional-order derivative to simulate the given dynamics. Equilibrium points for the given fractional model have been calculated, and important discussion on the asymptotic stability of the equilibria of a new autonomous system has been evaluated. We have reminded some important results to prove the existence of unique solution for the fractional-order cases. For finding the numerical solution of the given system, we used a generalized predictor–corrector algorithm in the sense of the new generalized Caputo derivative and also justified the stability of the technique. To prove the importance and correctness of the numerically simulated results, we have performed a number of graphs at different fractional-order values. The given derivative and algorithm work very well to understand the dynamics of the given model. From this study the effects of greenhouse gases and hypoxia on the population of aquatic species can be clearly understood with memory effects. For the future scope, the given ecosystem can be further solved by any other fractional-order derivatives. Also, some new mathematical models can be proposed to simulate the structure of given real-world problems.