Journal of Natural Resources and Environmental Management http://dx. org/10. 29244/jpsl. RESEARCH ARTICLE The Habitat Suitability Modelling of Rhinoceros Hornbills (Buceros rhinocero. in Java Island. Indonesia Taufik Setiawana. Lilik Budi Prasetyob. Yeni A. Mulyanib. Jarulisc Tropical Biodiversity Conservation Study Program. Graduate School. IPB University. IPB Darmaga Campus. Bogor, 16680. Indonesia Department of Forest Resources Conservation and Ecotourism. IPB University. IPB Darmaga Campus. Bogor, 16680. Indonesia Department of Biology. Faculty of Mathematics and Natural Sciences. University of Bengkulu. Kandang Limun. Bengkulu, 38122. Indonesia Article History Received 06 September 2023 Revised 04 October 2023 Accepted 05 December 2023 Keywords environmental variables. Java Island habitat suitability, presence data. Rhinoceros hornbills ABSTRACT Rhinoceros hornbills (Buceros rhinocero. are a bird species belonging to the Bucerotidae family, which is vulnerable based on the IUCN red list of species. This is due to habitat fragmentation, which reduced the Rhinoceros hornbill habitat on Java. Efforts and strategies are needed to maintain Rhinoceros hornbill habitats. Information on the suitability of the Rhinoceros hornbill habitat on Java Island is required to develop a Rhinoceros hornbill conservation strategy. This study aimed to determine a habitat suitability model that produces the highest accuracy, analyze hornbill habitat suitability, and identify environmental variables that affect the existence of rhinoceros hornbills. Habitat suitability models were processed using three algorithms: random forest, support vector regression, and MaxEnt. The data used to model habitat suitability were presence and environmental variables. The model was evaluated using various accuracy measures, namely overall accuracy, sensitivity . , specificity . Area Under Curve (AUC), and kappa coefficient. The results of model processing showed that the random forest algorithm produced the highest average accuracy of 0. The most important environmental variables for the habitat suitability model were the distance from the road . 62%), distance from the forest . 73%), and land cover . 47%). The habitat suitability model was divided into three classes: low suitability, covering 75,048 km 2 . 94%). medium suitability, covering 52,911 km2 . 44%). and high suitability, covering 6,213 km2 . 63%). The results of the habitat suitability model showed that the habitat suitability class was the smallest in the area. Introduction Rhinoceros hornbills (Buceros rhinocero. are bird species belonging to the Bucerotidae Family. Rhinoceros hornbills are characterized by their unique shape and orange-head protection . According to Mackinnon et al. , hornbills are spread across several islands of Indonesia, including Sumatra. Java, a nd Kalimantan. According to the International Union for Conservation of Nature (IUCN), rhinoceros hornbills are categorized as vulnerable (VU). Conservation efforts are needed to protect and preserve the Rhinoceros hornbill so that the population does not decrease. Java Island is one of the Rhinoceros hornbill distribution islands in Indonesia. Java Island is one of the islands in Indonesia that has experienced massive forest loss and habitat fragmentation . Habitat fragmentation on Java Island has led to the isolation of many wildlife populations and the threat of extinction . According to the Ministry of Environment and Forestry, forest area data in Java Island was 3,212. 7 ha. Forest loss on Java Island and habitat fragmentation are feared to impact habitat loss that threatens the existence and sustainability of rhinoceros hornbills on Java Island. Based on these conditions, effective and efficient conservation strategies are needed. Complete information on hornbill habitat suitability models is required to provide an overview of hornbill habitats and can be considered when developing conservation strategies. Habitat suitability models can help determine locations Corresponding Author: Taufik Setiawan taufiksetiawan437@gmail. Graduate School. IPB University. IPB Darmaga Campus. Bogor. Indonesia. Tropical Biodiversity Conservation Study Program. A 2024 Setiawan et al. This is an open-access article distributed under the terms of the Creative Commons Attribution (CC BY) license, allowing unrestricted use, distribution, and reproduction in any medium, provided proper credit is given to the original authors. Think twice before printing this journal paper. Save paper, trees, and Earth! that need to be maintained or improved. Modeling hornbill habitat suitability on Java Island is essential to provide information on the level of hornbill habitat suitability on Java Island. The existence of habitat suitability model information can be considered by policymakers when developing hornbill conservation plans on Java Island. Based on previous studies, many algorithms have been used to model the suitability of a species habitat. This research used three different algorithms, namely Random Forest (RF). Maximum Entropy (MaxEn. , and Support Vector Regression (SVR). These three algorithms have been widely used to model the distribution of habitat suitability for both flora and fauna. Research on hornbill distribution modeling has been conducted only by Jarulis . Hidayat and Febriani . , who used the MaxEnt model. SVR and RF models have never been used to model hornbill habitat suitability. Therefore, this research is critical for comparing and determining the best model among MaxEnt. SVR, and RF. The determination of the best model among the three algorithms was based on the accuracy value generated by each algorithm. According to Cahyana . , the difference in the level of accuracy produced by each method can also be considered separately when choosing the model used. This research will provide information on the best algorithm for modeling Rhinoceros hornbill distribution. In addition, this research is also expected to be considered by policymakers in hornbill conservation efforts, especially on Java. Material and Method Research Location and Time The scope of the research location analyzed in this study was Java Island. The process of collecting secondary data on Rhinoceros HornbillAos presence on Java Island and data on environmental variables was conducted in January 2023. The hornbill presence data in this study were in the range of 2000Ae2023. Figure 1. Map of research location. Data Collection Habitat suitability models used hornbill presence distribution data and data on environmental variables that affect hornbill distribution. The presence distribution data for hornbills used in this study were derived from secondary sources, including the scientific literature and websites. The presence of data collected in this study is in the range of 2000Ae2023 because the availability of data before 2000 is still relatively small. The scientific literature data were divided into two groups: books and journals. The book data comprised 32 This journal is A Setiawan et al. JPSL, 14. | 254 points of presence, consisting of textbooks, theses, dissertations, and research reports. Journal data totaled 95 attendance points from journal articles, proceedings, and seminar/conference/symposium/workshop Secondary data from books and journals collected were in publications between 2000Ae2023 because the availability of hornbill presence data from books and journal sources is still relatively small. In addition, data on the distribution of presence were obtained from the websites. The website used in this study is the Global Biodiversity Information Facility (GBIF). The data from the website obtained in this study comprised a total of 1,158 attendance points. This website provides complete data on the existence of fauna species and other organisms . The initial presence data collected before the filtering process included 1,285 points of presence for one region in Indonesia. Environmental factor data consisted of altitude, slope, temperature, rainfall. Normalized Difference Vegetation Index (NDVI), land cover, distance from the road, distance from the river, distance from the forest, distance from the farmland, distance from the plantation, distance from the settlement, and distance from food trees. Elevation and slope data were generated from DEM data processing downloaded from Temperature data was downloaded from w. Rainfall data was downloaded from w. The NDVI and closure data of Landsat 8 images were processed and downloaded to Google Earth Engine. Distance from the road and river data using RBI data downloaded from Distance from the forest, distance from the farmland, distance from the plantation, and distance from the settlement was obtained from the processing of land cover data. Meanwhile, distance from the food trees data was obtained from the GBIF website. The environmental factor data obtained were then converted into raster data in the GeoTIFF format. Data Analysis This study uses hornbill presence data and environmental factors. Presence data and environmental factors are analyzed using 3 types of algorithms, namely MaxEnt. RF, and SVR. The selection of these algorithms is due to machine learning algorithms that are often used and produce good accuracy. The data analysis process used in this study is shown in the flowchart (Figure . Figure 2. Research flowchart. Data Pre-Processing Data pre-processing is the initial stage before the data are used for modeling. Before being used as a model processing variable, a resampling process was performed on the variable to equalize its spatial resolution. Rhinoceros hornbill presence data were filtered previously. This process is essential for modeling the distribution of species habitats because it improves the quality of the data and accuracy of the model produced . The filtering process was carried out with various considerations, namely, the presence points have the exact coordinates, eliminating duplicate presence data, and eliminating presence data outside Java Island. The number of Rhinoceros hornbill presence data collected for the Java Island area was 172 points, and then filtering was performed by eliminating duplicate presence data. thus, the presence data used in the model processing was 63 points. Using Pseudo-absence Data Presence and absence data are needed to model hornbill habitat suitability on Java. It is difficult to obtain animal absence data, which can be anticipated using pseudo-absence data. The amounts of pseudo-absence and presence data were the same. Furthermore. Barbet-Massin . recommended using the same pseudoabsence points as presence points. The estimation of pseudo-absence points in this study was based on the presence of forests in the study landscape, because Rhinoceros hornbills are highly dependent on the presence of forests. Therefore, pseudo-absence points were randomly selected from the areas outside the http://dx. org/10. 29244/jpsl. JPSL, 14. | 255 Data Partitioning The presence and absence of data was then divided into training and test data. Training data are used to build models, and testing data are used to evaluate model performance . Division was performed using the K-fold cross-validation method. Data division or partitioning was conducted using the dismo package in the R-Statistic application with a number K of 5 folds. Selection of Environmental Variables The initial stage that needs to be done in modeling is the selection of the environmental variables. The goal was to select the variables used in the model-processing process. One way to choose environmental variables is to conduct a multicollinearity test for all environmental variables. Multicollinearity is a condition of a linear relationship or high correlation between each independent variable. The multicollinearity test is used to test the linearity between one variable and another, such that the variables used in building the model do not affect each other . A multicollinearity test was performed using the olsrr package in R-Statistics. This test was performed to obtain the Variance Inflation Factor (VIF). If the VIF value obtained is more than ten (>. , then one of the variables must be eliminated or reduced . Model Processing Model processing in this study used the MaxEnt, random forest, and support vector regression algorithms. All data from the three types of machine algorithms were processed using the R Studio software. MaxEnt is a machine-learning technique that can be used to predict the presence of data information only . MaxEnt was applied using R-Studio using the dismo package . Random forest (RF) is a machine learning technique with a bootstrap-based classification and regression tree method . Random forest was implemented using R Statistics with the random forest package . SVR is part of the Support Vector Machine (SVM) SVR is the implementation of an SVM for regression. The regression output is a natural or continuous number . Support Vector Regression was applied using the R-Statistic through package e1071 with the function SVM . Model Evaluation The results of model processing were then evaluated to determine the best model that produced the highest The model was evaluated using various accuracy measures: sensitivity, specificity. Area Under Curve (AUC), and kappa coefficient. Overall accuracy describes the success of attendance and absenteeism being correctly classified in the resulting model . Sensitivity is the chance of attendance being correctly classified and specificity is the chance of attendance being correctly classified . AUC is the area under the Receiver Operating Curve (ROC) and is a standard method for identifying the prediction accuracy of distribution models . The kappa statistic measures the model performance using the information in a mixture matrix . that ranges from Ae1 to 1, where 1 indicates good results and a value of zero or less than 0 indicates poor performance . To classify potential presence, a probability threshold above 0. was used . Hornbill Habitat Suitability Class Range The hornbill habitat suitability is divided into three suitability classes. The three classes are low suitability class, medium suitability class, and high suitability class. This map provides a habitat suitability map for hornbills and determines which areas are most suitable for conservation and protection. The classes of the habitat suitability intervals are listed in Table 1. Table 1. Habitat suitability category of rhinoceros hornbills . Range Suitability class Minimum < S O . ean Ae Std. Low suitability . ean Ae St. < S O . ean St. Medium suitability S > . ean St. High suitability Note: Std = standard deviation. S = score. mean = average. This journal is A Setiawan et al. JPSL, 14. | 256 Results and Discussion Selection of Environmental Variables The selection of environmental variables must be performed in the early stages of modelling. This selection aims to select the variables that are used in the modelling process. One way to select environmental variables is to conduct a multicollinearity test for all environmental variables. A multicollinearity test determined the correlation between each environmental variable used in the model processing. Multicollinearity between these variables affects the accuracy of the model . The results of the multicollinearity test between the environmental variables are presented in Table 2. Based on the results of the multicollinearity test, we found multicollinearity between the variables of altitude and temperature. This can be seen from the VIF values of the temperature and altitude variables. The VIF value for the temperature and altitude variables was more than 10 (>. If multicollinearity exists between two variables, then one of the environmental variables must be selected for use in model processing . The altitude variable was chosen because the higher a region, the lower . its surface temperature. The lower the area, the higher . the surface temperature of the hotter . In addition, several previous studies on multicollinearity tests have shown that altitude is widely chosen to model species distribution, such as in elephants . , dare monkeys . , and litter frogs . Table 2. Multicollinearity test results. Environmental Variables VIF Altitude Slope Temperature Rainfall NDVI Land cover Distance from road Distance from river Distance from forest Distance from agricultural area Distance from plantation Distance from settlement Distance from food tree Note: the red color in the table indicates a correlation between variables. Model Evaluation The resulting model was evaluated by determining the overall accuracy, sensitivity, specificity. AUC values, and kappa coefficient. A comparison of the model performance is presented in Table 3. Table 3 presents the evaluation results of the Rhinoceros hornbill distribution model for each algorithm. Table 3. Rhinoceros hornbill distribution model accuracy test. Algorithm Kappa AUC Maximum Entropy Random Forest Support Vector Regression Note: OA= overall accuracy. Sn= sensitivity. Sp= specificity. AUC= area under curve. The models that produce high accuracy are colored in red. Based on the accuracy values generated, the random forest algorithm produced high accuracy from four accuracy measures, namely overall accuracy, sensitivity, kappa, and AUC, with consecutive values of 0. 06, 0. 83 A 0. 11, 0. 78 A 0. 12, and 0. 97 A 0. 02, respectively. The MaxEnt algorithm exhibits specificity accuracy with a value of 0. 98 A 0. This accuracy comparison showed differences in the sensitivity and specificity http://dx. org/10. 29244/jpsl. JPSL, 14. | 257 between the three algorithms used in this study. The sensitivity value was used to determine the chance of attendance that was correctly classified in the model results. In contrast, the specificity value was used to determine the opportunity of absence being correctly classified in the model results. The random forest model produced the highest accuracy in terms of sensitivity, while the MaxEnt model produced the highest accuracy in terms of specificity. There were differences in the accuracy values produced by each of these models, so that the best model selection could be seen from the average value generated from the overall accuracy, sensitivity, specificity, kappa, and AUC values. The average calculation results show that random forest produces the highest average value of 0. MaxEnt and support vector regression produced an average accuracy of 0. 70 and 0. The random forest algorithm produces a higher model accuracy than the MaxEnt and support vector regression algorithms. According to Wijaya et al. , a random forest is an ensemble method consisting of several decision trees. The random forest algorithm is also a development of the decision-tree The development in a random forest uses several decision trees that are trained and broken down in each tree and then randomly subsets . The random forest algorithm produces high accuracy because this algorithm is an algorithm that combines and averages the results of several decision trees, so the random forest algorithm is the algorithm used in the subsequent analysis. Importance Variable The Variable Importance Method (VIM) was used to determine the environmental variables that are important to the RF algorithm. The VIM method can be used to determine the significance of environmental variables in model processing. The method measures the extent to which an environmental variable can reduce model accuracy . Furthermore, the calculation can use the mean decreased accuracy (%IncMSE) . The greater the %IncMSE value, the greater the error value caused by randomizing an environmental variable's permutation, such that the importance value of the environmental variable increases . The percentage increase in the Mean Square Error (%IncMSE) can identify environmental variables that are important for modeling using the random forest algorithm. %IncMSE expresses the increase in mean square error when a variable is randomly permuted (MSEpermute. Mathematically, this formula is as follows: %IMSE = ( ycAycIyaycEyceycycoycycyceycc Oe ycAycIy. ycAycIya x 100% . The random forest algorithm produces environmental variables that are important for the habitat suitability The model results showed that distance from road resulted in the highest %IncMSE of 16. followed by distance from forest . 73%), land cover . 47%), distance from settlement . 26%), rainfall . 04%), distance from plantation . 11%), altitude . 61%). NDVI . 61%), slope . 59%), distance from farmland . 54%), distance from food trees . %), and distance from the river . 41%). Table 4. Important environmental variables on habitat suitability. Environmental variables %IncMSE Distance from road Distance from forest Land cover Distance from settlement Rainfall Distance from plantation Altitude NDVI Slope Distance from farmland Distance from food trees Distance from river Based on Table 4, the distance from the road is an essential variable for modeling because if the distance from the road is omitted or not used, it will significantly reduce the accuracy value of the %IncMSE by 16. Conversely, the distance from the river variable produced the lowest %IncMSE value, indicating that if the distance from the river variable is not used in the model processing, the accuracy will be reduced by 1. This journal is A Setiawan et al. JPSL, 14. | 258 However, distance from the river is still considered in the model processing because the distance from the river influences the presence of rhinoceros hornbills ecologically. According to Hidayat and Febriani . , rivers are essential for animal survival in fulfilling water requirements. They are used as energy sources and for animal livelihoods, including Rhinoceros hornbills. Effect of Environmental Variables on Prediction of Rhinoceros hornbill Presence The effect of environmental variables on the presence of Rhinoceros hornbill can be observed from the response curve results obtained from the model processing results. The response curve predicts the presence of rhinoceros hornbills in each environmental variable. According to Phillips et al. , the closer the value is to 1, the higher is the chance of animal presence. Topographical Factors Based on the response curve of the altitude variable obtained from the model processing results, the possibility of hornbill presence was between 500 and 2,000 masl. This is in contrast to Hadiprakarsa and Kinnaird . , who explained that the distribution of Rhinoceros hornbill species on the island of Sumatra is found at an altitude of 0Ae2,200 masl. This condition differs from the research of Jarulis . , which states that, on average, all hornbill species in Indonesia prefer lower altitudes, namely, <1,500 masl. The effects of altitude on the presence of hornbills are more diverse among species. According to Bibby et al. , altitude is one of the most important factors for bird distribution. In addition, altitude also affects the vegetation used by birds as both nest trees and food trees. The fruit of ficus trees from the Moraceae family is a tree much favored by Rhinoceros hornbills. Ficus trees were found in the tropics at <1,500 masl. However, a few species of ficus trees can grow above 3,000 masl . The slope is a topographic factor that affects hornbill distribution. The land slope response curve results showed that the probability of hornbill presence increased with increasing degrees of land slope (Figure . According to Hadiprakarsa . , slope is one of the essential parameters for predicting hornbill presence. The possibility of hornbill presence was mainly found on land with a slope of more than 10A. This indicates that the more the slope level increases in an area, the more likely is the presence of Rhinoceros hornbills will be The Bucerotidae family lives in hilly regions with an undulating topography and steep slopes . Kubangun et al. stated that flat topographic areas generally experience more land change than steep topographic areas. Therefore, flat topographic regions are less suitable habitats for Rhinoceros hornbills. Figure 3. Land elevation and slope response curves. Climatic Factors The results of the response curve for rainfall showed that the probability of Rhinoceros hornbill presence increased in the range 2,000 to 4,500 mm (Figure . The possibility of hornbill presence decreased when rainfall was below 4,500 mm. Rhinoceros hornbill habitat is generally found in areas with moderate rainfall. In contrast. Kinnaird and OAoBrien et al. stated that hornbill species are usually located in areas with more than 3,000 mm/year of rainfall. Rainfall affects the availability of water for the hornbills. High rainfall can produce water used by hornbills to meet their needs, such as for drinking . http://dx. org/10. 29244/jpsl. JPSL, 14. | 259 Figure 4. Rainfall response curve. Biophysical Factors In species distribution models, land-cover variables can influence the presence of species . The land cover classes in the study included settlements . , forests . , water bodies . , open land . , shrubs . , agricultural land . , and plantations . The response curve results showed that Rhinoceros hornbill was mainly found in the forest. Hornbills use forest habitats dominated by trees to find food, perform activities, and breed. Trainor et al. revealed that forests are the most essential habitat for the existence of birds because forests provide food, water, and shelter for birds in living life. The NDVI is an index value that represents the greenness of an area. Generally. NDVI values range from Ae1 NDVI values close to 1 indicate excellent vegetation greenness, whereas NDVI values less than 0 indicate poor innocence. The range of NDVI values produced in this study is Ae0. 23 to 0. The response curve results showed that the probability of Rhinoceros hornbill presence was in the NDVI range of 0. 25 to 0. According to Departemen Kehutanan . NDVI values are grouped into three classes: Ae1Ae0. 32 for low density, 0. 32Ae 42 for medium density, and 0. 42Ae1 for dense canopy density. Based on the NDVI value class, it was concluded that the hornbill's chance of existence is in the medium- and thick-density classes. According to Rahman . , dense vegetation in forested areas is essential for species survival, including Rhinoceros Hornbill survival depends on forests dominated by large and tall trees. Hornbills use these trees for breeding, perching, and nesting purposes. The response curve results for distance from the river showed that the probability of presence increased with distance from the river. Nevertheless, distance from the river remains essential for the life of rhinoceros hornbills because rivers provide water that can be an energy supply and a source of livelihood for trees used as food and nest trees for rhino hornbills. Previous research conducted by Hidayat and Febriani . stated that the distance from the river is essential in the presence of hornbills because the river is directly adjacent to the forest, thus showing the high need for water for the survival of Rhinoceros hornbills. Rivers are one of the sources of energy and a source of livelihood for animals, including Rhinoceros hornbills. Food trees are an essential factor for the life of rhinoceros hornbills. Rhinoceros hornbills require food to fulfill their energy requirements. Fruits of Ficus spp. are the leading food source and are also favored by According to Anggriawan et al. , the distribution of Ficus spp. is positively correlated with the hornbill distribution. The response curve results showed that the greater the distance from the food tree, the smaller the chance of hornbill presence (Figure . Therefore, food trees are essential for the survival of rhinoceros hornbills on Java. Human Activity Factors The environmental variable that can influence the distribution of rhinoceros hornbills is disturbance caused by human activities. Human activity factors are based on five types of environmental variables: distance from roads, distance from forests, distance from agricultural land, distance from plantations, and distance from The response curves of the five variables demonstrated different trends (Figure . The curves of distance from roads, distance from farmland, and distance from settlements showed an increase in the This journal is A Setiawan et al. JPSL, 14. | 260 value of environmental variables that aligned with the increased likelihood of rhino hornbill presence. contrast, the distance from forests and the distance from plantation curves showed a decrease in the value of environmental variables aligned with a reduction in the probability of rhinoceros hornbills. Figure 5. Response curves of NDVI, land cover, distance from the river, and distance from food trees. http://dx. org/10. 29244/jpsl. JPSL, 14. | 261 Figure 6. Response curves of distance from roads, distance from forest, distance from farmland, distance from plantation, and distance from settlement. The distance from the road response curve in Figure 6 reveals that the greater the distance from the road, the greater the probability of hornbill presence. According to Hidayat et al. , the closer the road is to the presence of hornbills, the more threatened their existence. Easy access to hornbill habitats can facilitate hunting activities carried out by humans against hornbills. The results of the distance from the forest response curve showed that the greater the distance from the forest, the less likely the presence of hornbills. According to Kumara . , hornbills generally live in forested areas far from human activities. In addition, according to research by Jarulis . , hornbill encounters were found at a distance from the forest. This is because Rhinoceros hornbills do not resemble habitats disturbed by human presence or activity. Based on the results of the response curve generated by the model, the greater the distance from the plantation, the lower the probability of Rhinoceros hornbillAos presence. Therefore. Rhinoceros hornbills can avoid human disturbance. According to Syechalad and Muhammad . , plantations are agricultural businesses that become an economic activity by cultivating crops that are controlled separately. In addition to being managed individually, plantations are governed by companies. Some hornbill species can also be found in plantation areas . Distance from settlements is a human disturbance factor that affects the distribution of Rhinoceros hornbills. Rhinoceros hornbills are hornbill species that are sensitive to human presence, causing them to live away from human existence. The response curve generated from the model processing showed that the farther away from the settlement, the greater the chance of hornbill presence. According to Hidayat et al. , the closer to the settlement area, the higher the threat to hornbills due to later human access to hornbill habitats. This journal is A Setiawan et al. JPSL, 14. | 262 Rhinoceros hornbill Habitat Suitability Random forest processing resulted in raster data with the lowest pixel value of 0. 0034 and highest pixel value The standard deviation of the resulting data was 0. 20, and the mean was 0. Based on the results of data processing that has been carried out, the division of hoses for the Rhinoceros hornbill habitat suitability class is presented in the following table 5. Based on the results of the suitability class area calculation and the suitability class map in Figure 7, the results show that the high suitability class has the smallest area. The area of the low habitat suitability class is 75,048 km2 . 94%) with two presence points, then the medium suitability class with an area of 52,911. 77 km2 . 44%) with 15 presence points, and the high suitability class with an area of 6,213. 76 km2 . 6%) with 46 presence points. Table 5. Hornbill habitat suitability classes on Java Island. Range Suitability class Occurrence point Area . Percentage (%) 0034 Ae 0. Low suitability 048,00 29 Ae 0. Medium suitability 911,77 70 Ae 0. High suitability 213,76 Figure 7. Map of Rhinoceros hornbill habitat suitability on Java Island. The results of modeling hornbill habitat suitability on Java Island showed that suitable habitat for Rhinoceros hornbills on Java Island has the most minor proportion of area at 4. The low area of suitable habitat for Rhinoceros hornbills indicates a severe threat to the preservation of Rhinoceros hornbills. Jarulis . revealed that the danger of hornbills in Indonesia is high. Therefore, it is necessary to increase hornbill conservation efforts to reduce and prevent the rate of hornbill extinction, especially Rhinoceros hornbills on Java Island. Therefore, this must be a concern for all parties, including the central government, local governments, and associated stakeholders, to preserve the Rhinoceros hornbill habitat on Java Island. Thus, policies such as making areas with high suitability levels into priority zones in habitat development and Rhinoceros hornbill protection are needed. In addition, there needs to be a plan and action to restore forest and non-forest areas to become habitats for Rhinoceros hornbills on Java Island. Pramana et al. stated that to conserve animals, there must be collaborative conservation efforts between various parties, including the government, area managers, and http://dx. org/10. 29244/jpsl. JPSL, 14. | 263 local communities. Therefore, cooperation from various related parties, such as the Ministry of Forestry and Environment, the relevant Regional Forestry Service, and related stakeholders who have concerns about hornbill conservation efforts must be carried out so that these plans and actions can be carried out effectively and efficiently. Conclusion The model processing results revealed that random forest produced a higher average accuracy of 0. Therefore, the random forest was the best algorithm that could be used in processing habitat suitability The environmental variables that were important to the Rhinoceros hornbill habitat suitability model on Java Island were distance from roads . 62%), distance from forests . 73%), and land cover . 47%). High habitat suitability resulted in a small area of 6,213. 76 km2 . 6%). The small proportion of high habitat suitability indicates that this is a severe concern for policymakers and related stakeholders. Therefore, there is a need to plan to preserve the habitat and population of rhinoceros hornbills on Java. A critical step in conserving Rhinoceros hornbill habitats and populations is to protect habitats in areas with high suitability levels and to restore hornbill habitats in areas with low suitability levels. The suggestion obtained from this study is the need for additional Rhinoceros hornbill distribution points to improve the resulting model. Additional distribution points were obtained through direct survey activities in the field. However, actions will require considerable costs and human resources, considering the vast area of Java. In addition, the algorithm used in this study can be applied to model hornbill habitat suitability on islands other than Java. Thus, it can determine the algorithm that produces the highest accuracy and also know the level of hornbill suitability other than on Java Island. References