J. Indones. Math. Soc. Vol. No. , pp. 261Ae269. ANALYSIS OF FOREST EXISTENCE AGAINST POPULATION OF Eos histrio IN NORTH SULAWESI Nastitie1 and Gantina Rachmaputri2 Teaching of Mathematics. Faculty of Mathematics and Natural Sciences. Bandung Institute of Technology, aku. nastitie@gmail. Algebra Research Group. Faculty of Mathematics and Natural Sciences. Bandung Institute of Technology, gantina. r@math. Abstract. Red List International Union for the Conservation of Nature and Natural Resource (IUCN) states there are more than 8000 animals are critically endangered and it will be increase. One of the Indonesias endemic animal that is included in the IUCN Red List is Eos histrio or Burung Nuri Talaud which lives in North Nusa Islands. North Sulawesi. Their natural habitat is a forest with tall trees as a place to Based on data from the Central Statistics Agency (BPS) the area of forest in the region is decreasing every year, causing Eos histrio to lose his home and causing a decrease in the population. Based on data from the IUCN, the population of Eos histrio in 1999 was at 8230 21400, while in 2016 it was in the range of 5500 14000. The number will continue to decrease due to the large amount of deforestation and the lack of protection for animals. By using the Pontryagin Maximum method, it is found that one way to protect the population of Eos histrio is to reduce the rate of deforestation in North Sulawesi by at least 10% so as to increase the population of Eos histrio within 4 year. Key words and Phrases:Eos histrio, deforestation. Pontryagin Maximum method INTRODUCTION The International Union for the Conservation of Nature and Natural Resource (IUCN) Red List reported that there are over 8. 000 critically endangered species and it seems likely to climb. One of the species included in the IUCN . Red List is Eos histrio, also known as the Red-and-blue lory, a bird endemic to Indonesia that inhabits Pulau Karakelang Island. Nusa Utara Islands. North Sulawesi. Adi Widyanto . said that Karakelang is a conservation area in North Sulawesi Province with an area of 24669 hectares as stipulated by the Decree of the Minister of Forestry and Plantation Number 97/Kpts-II/2000. More than just a lowland rain forest and 2020 Mathematics Subject Classification: 49-11 Received: 02-07-2021, accepted: 29-09-2021. Nastitie and G. Rachmaputri mountains with karst soil, which is a unique habitat, the area also houses Talaud Archipelagos typical biodiversity, including the Red-and-blue lory. According to Adi Widyanto . , the decline rate of vegetation cover on the Talaud Islands, which includes the three main islands of Karakelang. Kabaruan, and Salibabu, was 4457 hectares in 2001 2017. On the other hand, the forest cover increase rate was 1320 hectares in 2001 - 2012. The forest loss locations are found not just in plantation regions, but also in forest areas. This data confirms the reports stated that there are still forest encroachment activities to open a farm around the forest area. From Diah Irawati . we know that Eos histrioAos roosting trees criterion are branch-free bole length, canopy diameter, canopy density and the number of trees around the roost tree on diameter 20 - 40 cm and up to 41 cm. Eos histrios native environment is a forest with tall trees and focused in the primary forest that is a closed classification forest with low habitat disturbance. In this case. IUCN . reported that there are 8230 21400 Eos histrio in 1999, however, in 2016, it was only around 5500 14000 full-grown Eos histrio. Due to the many threats in the conservation area, such as illegal trading of the bird and severe deforestation encroachment for farming on the main forest and commercial logging, the bird population is expected to continue to decline. Based on Adi Widiyanto . that focus on community empowerment on forest protection near Karakelang wildlife sanctuary and Jon Riley . that focus on population size of Eos histrio, hence we will analyze the effect of forest existence on Eos histrio population growth on Karakelang Island. North Sulawesi. MAIN RESULTS The mathematical model that used to analysis of forest carrying capacity on Eos histrio is Lotka-Volterra . rey-predato. model from A. Panvilov . where prey is assumed to be the population of Eos histrio while predators are assumed to be the rate of deforestation, as follows: Figure 1. Flowchart of Lotka-Volterra = aR Oe bR2 Oe cRN = dRN Oe eN. Analysis of Forest Existence Againts Population of Eos histrio in North Sulawesi From IUCN . Adi Widyanto . Jon Riley . , and Birdlife International . , we will find parameter values as follows: Because . is non-linear model, so we Table 1. Description and parameter values Parameter Description Value Eos histrioAos population Deforestation Increased Eos histrios population Competition between Eos histrio Number of bird who died due to deforestation Forest growth rate Rate of deforestation Ae Ae need to linearization that equation using Jacobi matrix around their equilibrium Fisrt, we will find the equilibrium point from equation . , as follows: dRN Oe eN = 0 N . R Oe . = 0. Hence the equilibrium point from . Using equilibrium point from . , we will find equilibrium point for . , as follows : = aR Oe bR2 Oe cRN = 0. If N = 0, then aR Oe bR2 = 0 N =0OR= When R = de then R. Oe bR) = 0 R=0OR= . Oe b( )2 Oe N = 0 ad Oe be . Thus there are three equilibrium points can be found from equation . , . , and . , which are . , . , ( ab , . , ( de , adOebe cd ) and if we substitute the parameter values in Table 1, we get the results are . , . Nastitie and G. Rachmaputri STABILITY ANALYSIS : Perform the linearization process by forming the Jacobian matrix as followss as in A. PanvilovAos Book . OCR f OCR g a Oe 2bR Oe cN OecR OCN f OCN g dR Oe e Determine the stability around the equilibrium point by determining the eigen values using equation . For the equilibrium point . , . then |J Oe I| = 0 was aOe Oee Oe 1 = a O 2 = Oee. For the equilibrium point . 74, . then |J Oe I| = 0 was obtained : Oe0. 49 Oe Oe0. 72 Oe 1 = Oe0. 49 O 2 = 4. For the equilibrium point . , . then |J Oe I| = 0 was obtained : Oe0. 01788 Oe Oe0. Oe0. 0001 Oe 1 = Oe0. 00899 Oe 0. 29i O 2 = Oe0. From the above calculation the equilibrium point . is a stable equilibrium point based on A. Panvilov . and de Bour . CONTROLLABILITY ANALYSIS : The objective function of this research is to maximize the population growth of Eos histrio by minimizing the rate of deforestation on Karakelang Island. The controllability of models . will be reviewed around the equilibrium point . by forming a control matrix. OCe f OCe g Oee Next step is find the matrix Mc = (B|JB|J 2 B|. |J nOe1 B) based on SubionoAos book . Mc = Oee e2 Oe edR Because rank of matrix Mc = J then the system is controlled . DETERMINE THE OBJECTIVE FUNCTION : In this research, the objective function is to control forest damage so as to maximize the growth of the Eos histrio population on Karakelang Island. The objective function of the system is Z tf H. = R C1 e . , dt. t OO . 0 , tf ]. R > 0. OPTIMAL CONTROL SOLUTION : The optimal control problem for the Analysis of Forest Existence Againts Population of Eos histrio in North Sulawesi system of equations . and the objective function above using the Pontryagin Maximum method with the following steps based on Naidu . Step 1: Determine the Hamiltonian function. H(R. N, e, ) = R C1 e . R Oe bR2 Oe cRN ) 2 . RN Oe eN ). Step 2: Minimize H to e. OCH OCe C1 e. Oe 2 N = 0 Step 3 : Determine state and costate The formed of state equation is e. = OCH = aR Oe bR2 Oe cRN OC1 OCH = dRN Oe eN. OC2 And the formed of costate equation is Oe OCH = Oe. Oe 2bR Oe cN )). OCR SIMULATION AND ANALYSIS : This section will be discussed about the results of numerical simulations before and after being subjected to control of the forest carrying capacity model on the population of Eos histrio where the control variable is deforestation rate in conservation areas on Karakelang Island and using the Runge Kutta method of order 4 based from Oruh . The initial values that used in the simulation was the population density of Eos histrio per km2 in primary forest is 15. 1 individuals while in secondary habitat it is 4. 3 individuals per km2 and the rate of forest destruction is 44 hectares based on Adi Widyanto . The numerical simulation will be given a control weight of 10% means that the rate of forest destruction will be reduced by 10% and the observation period is 4 years with the following results. In Figure 2, the red line shows the uncontrolled deforestation rate during 4-year observation period. After the control was adjusted to the deforestation rate during the 4-year observation period, a blue line shows the fluctuating damage rate, but the level is not as severe as when it has not been It means that the damage rate, although still relatively high, can still be controlled. It has an impact on the Eos histrio population in the same observation period of four years as shown in Figure 3. The red line shows the Eos histrio population when the deforestation rate is not controlled. It can be seen it influences the population growth even to its lowest level of 0. Meanwhile, if the deforestation rate Nastitie and G. Rachmaputri Figure 2. Comparison graph of the rate deforestation before and after given control. Figure 3. Comparison graph of the Eos histrio before and after given control. is controlled, the population of the bird is fluctuating. It correlates with the blue line in Figure 2. If the deforestation rate is controlled, the Eos histrio population could increase. On the other hand, if deforestation is uncontrolled, the Eos histrio population would be threatened. Figure 4 shows deforestation rate has to be controlled by 10% for 4 years to optimally maintain the Eos histrio population. The value is constant so that the desired objective function can be achieved. It is indicated by the result of the objective function after the control value of 78. 2973 was added, which is higher than Analysis of Forest Existence Againts Population of Eos histrio in North Sulawesi Figure 4. Comparison of control that given to model for 4 years. Figure 5. Comparison of control that given to model for at least 10 years. the objective function before the addition of the control value that is 0. Figure 4 show how if we compare deforestation rate has to be controlled by 5% and 10% for 4-years. The result is no different between them, although the value of objective function decrease into 67. But, if we look for a long time, example in 10 years, there will have different in graphics control. In figure 5. there is graphic with controlled 5% for 10 years and the result is not constant. But, in figure 5. there is a graphic with controlled 10% for 15 years and more constant. So do in figure 6 there is a graphic with controlled 15% for 15 years and the result is constant. Therefore, deforestation rate has to be controlled at least by 10% to optimally maintain the Eos histrio population. Nastitie and G. Rachmaputri Figure 6. Deforestation rate controlled by 15% for 15 years. CONCLUDING REMARKS The existence of forest have an impact to Eos histrios population. When the rate of deforestation can be controlled, the population of Eos histrio has the possibility to increase, and vice versa. With a forest destruction suppression rate at least 10%, the population of Eos histrio has the opportunity to increase when compared to no efforts to suppress the deforestation. REFERENCES