Optimum Control for Spread of Pollutants through Forest Resources

Pollution has become the most critical factor spread by forest resources through wood-based and non-wood based industries. In other words, pollution is omnipresent. In this paper, the major pollutants caused due to wood and non-wood based industries are discussed which are the primary resources of the forest in spreading the pollution. In order to study the impact of industrialization and associated pollution on forest resources, the system of non-linear ordinary differential equations is formulated. The controls are advised on both types of industries to reduce the pollution.


Introduction
The forest resources mean a large area covered by trees. It means the various types of vegetation automatically growing on forest land where forest is considered as if it grew trees in the past, or will grow trees in the future. Wood-based industries are the branch of production and employment based on the fabrication, processing and preparation of products from raw materials and merchandises of wood and wood-based pulp products. Non-wood based industries are the branch of production and employment based on fabrication, processing and preparation of products from water, energy, chemicals etc. The term "pollution" is a substance into the environment which has harmful or poisonous effects on living beings. "Pollutants" are the components of pollution. It is observed that them merged in air as well as absorbed by forest region can be harmful for nature and can kill human organisms, essential microbes etc.
The forest resources are significant for human and some organisms. But after industrial revolution in 18 th century, industries are also growing very speedy [1].
This growth of wood and non-wood based industries has reduced the density of forest region. Damodar Valley, Nowamundi, Saranda are example of reduced forest resources [2]. Once Damodar Valley was covering 65% of forest area, now-a-days it is surrounded by only 0.05% [3]. In past, few years, the temperature of the environment is increasing due to the emission of pollutants which is examined by scientists and ecologists. This gives opposite impact on humans and environment [4] [5]. Absorption of pollutants by the plants is harmful and which affected the growth of forest resources [6] [7] [8] [9]. This motivated to formulate the system in which the effect of industrialization on the forest resources is analyzed. Some researchers have studied the mathematical model for the effects of industrialization and pollution on forest resources. [10] studied the models for the effect of toxicant in single-species and predator-prey system. [11] analysed the modeling effect of an intermediate toxic product formed by uptake of a toxicant on plant biomass. [12] introduced the effects of industrialization and pollution on resource biomass with the help of a mathematical model. [13] have performed the modeling effects of industrialization, population and pollution on renewable resources. [14] deliberated modeling effects of primary and secondary toxicants on renewable resources.
In this paper, a mathematical model is formulated with hereditary transmission of SIRS model in Section 2. The stability analysis of the transmission model is derived in Section 3. Sensitivity analysis is carried out in Section 4. Optimal control for the forest resources is discussed in Section 5. In Section 6, the model validated with numerical simulation and analysis.

Mathematical Modeling
In the society, there are different types of industries and pollution. Everybody in the society plays a role to decrease the pollution. Therefore, in the proposed model, five discrete compartments viz. the density of forest resources (F), the density of wood based industries (W), the density of non-wood based industries (I), the pollutants through wood based industries (P W ) and the pollutants through non-wood based industries (P I ) are considered. u 1 is the rate which decreases wood based industries to control the usage of forest resources. u 2 and u 3 are the control rates which decreases pollutants due to wood and non-wood based industries, respectively.
The notations and parametric values for the dynamical model are exhibited in Table 1.
Using these notations and assumptions which are required for formulating the mathematical model, the transmission diagram of forest resources is shown in  The loss of pollutants generated by wood based industries due to forest resources 0.02 The loss of pollutants generated by non-wood based industries due to forest resources 0.01 The depletion rate of forest resources caused by the pollutants generated through wood based industries 0.5 The depletion rate of forest resources caused by the pollutants generated through non-wood based industries 0.5 The growth rate of pollutants generated by wood based industries 0.1 The growth rate of pollutants generated by non-wood based industries 0.7 industries with associated pollutants is described as follows: Equations (1) to (5) is described as system (1) This gives, Therefore, the feasible region for system (1) is Thus, the equilibrium state of the system (1) is where ( ) Here, v is non-singular matrix, so the basic reproduction number 0 R is where ( )

Equilibrium
The equilibrium for the local and global stability of the forest transmission model are discussed here.

Local Stability
The forest resources equilibrium is locally asymptotically stable if all the eigenvalues of the matrix have positive real values [15]. The Jacobian matrix for sys- Using the parametric values given in the Table 1, Hence, system (1) is locally stable.

Global Stability
The forest resources transmission model is globally stable is ( ) Therefore, system (1) is also globally stable.

Sensitivity Analysis
In this section, the sensitivity analysis for all parameters are discussed in Table   2.
The normalised sensitivity index of the parameters is computed by using the following formula: where α denotes the model parameter.
The rate of compactness degree of forest resources, the constant rate of resources, migration of wood based industries to forest region, the depletion rate of forest resources due to wood based industries, the growth rate of wood based industries due to forest resources, the rate of competitive effects of I on W, the loss of pollutants generated by wood based industries due to forest resources, the loss of pollutants generated by non-wood based industries due to forest resources and the growth rate of pollutants generated by wood based industries have positive effect on 0 R which means they are helping us to save forest resources. Other parameters have negative impact on model.

Optimal Control
The objective of the model is to minimize the number of pollutants through wood and non-wood based industries to revive forest resources. The control functions are united to achieve the objective. The objective function for the mathematical model of forest resources in system (1)  , , w w w are weight constants for control variables 1 2 3 , , u u u respectively.
As, the weight parameters 1 2 , w w and 3 w are constants of forest resources control ( ) The partial derivative of the Lagrangian function with respect to each variable of the compartment gives the adjoint equation variables (  )   1  2  3  4  5 , , , , i A λ λ λ λ λ = corresponding to the system (1) which is as follows: The necessary condition for Lagrangian function L to be optimal for controls are ( ) To find the values of Thus, the required optimal control condition is computed as ( ) In next section the optimal control is calculated numerically to support the analytical results.

Numerical Simulation
Using the data given in Table 1 and Table 2, the sensitivity on model parameters is carried out.
The Figure 2 specifies that as the depletion rate of forest resources due to wood based industries increases the density of forest resources decreases.
From Figure 3, it is observed that increase in growth rate of wood based industries due to forest resources the density of forest resources decreases.
The Figure 4 indicates that with the loss of pollutants generated by non-wood based industries the forest resources increases.   From Figure 5 it is concluded that as the constant rate of resources provided to non-wood based industries increases the forest resources decreases.
From Figure 6, one can see that forest should be controlled 12% in 50 years, wood based industries should be controlled 33% in 68 years and non-wood  With control forest resources degradation reduces at a lower rate compare to no effective majors are taken up as shown in Figure 7. Figure 8 suggest that wood based industries can be controlled with effective majors at a lower rate compare to no control over it and when control is applied to wood based industries it decreases by 7%. Similar observation is from Figure   9 for non-wood based industries. In fact, it decreases by 4%.

Conclusions
In this paper, a mathematical model is formulated to study the spread of the forest transmission with wood and non-wood based industries. An optimal control for spread of the pollutants through the forest resources to study the effects of Wood, Non-wood based industries and the pollution emitted through them on the density of forest resources. The wood-based industries reduce the density of forest resources directly by harvesting as well as indirectly by pollutants. But non-wood based industries reduce the density of forest resources only indirectly by pollutants. Therefore, by more industries the forest resources are affected and may be wiped out.
The stability of forest resources model discussed with numerical data. The basic reproduction number is computed as 0.4960, which shows that controls on construction of wood and non-wood based industries will be beneficial to reduce the pollution. This suggested growing more and more forest resources and putting up less number of industries per human usage.