Journal of Natural Resources and Environmental Management http://dx. org/10. 29244/jpsl. RESEARCH ARTICLE Spatial Distribution of Local Sustainability Index in the Ciletuh-Palabuhanratu UNESCO Global Geopark Using Spatial Autocorrelation Siti Salsiaha. Izuru Saizenb. Endriatmo Soetartoc. Andrea Emma Pravitasaria,d Department of Soil Science and Land Resource. Faculty of Agriculture. IPB University. IPB Dramaga Campus. Bogor, 16680. Indonesia Graduate School of Global Environmental Studies. Kyoto University. Sakyo-ku. Kyoto, 606-8501. Japan Department of Communication Science and Community Development. Faculty of Human Ecology. IPB University. IPB Dramaga Campus. Bogor, 16680. Indonesia Center for Regional System Analysis. Planning, and Development (CRESTPENT). IPB University. Bogor, 16127. Indonesia Article History Received 15 June 2023 Revised 21 January 2024 Accepted 24 July 2024 ABSTRACT Keywords local moran index, local sustainability index, spatial The Ciletuh-Palabuhanratu UNESCO Global Geopark (CPUGG) is a geopark area that is also a tourist destination in Sukabumi. West Java. Indonesia. CPUGG is experiencing rapid regional development, which can be seen from the increase in tourist visitors, changes in land use from agricultural land to built-up areas, the rapid growth of infrastructure and facilities, and the emergence of various economic activities. This study aims to determine local sustainability performance based on the spatial distribution pattern of the local sustainability index (LSI) in CPUGG using a spatial autocorrelation approach. The results show that the CPUGG development is not evenly distributed The spatial distribution pattern of the local sustainability index in economic, social, and environmental are clustered. Each clustering shows that each village influences the surrounding So, increasing the economic, social, and environmental sustainability index in several villages with a low sustainability index is necessary. Good cooperation is needed between the government, society, and several stakeholders to improve economic, social, and environmental development in The Ciletuh-Palabuhanratu UNESCO Global Geopark. West Java. Indonesia. Introduction The Ciletuh-Palabuhanratu UNESCO Global Geopark (CPUGG) is a tourist destination in Sukabumi District. West Java. Indonesia. Geographically. CPUGG is advantageous because it is located in a district close to the metropolitan cities of Jakarta. Bogor. Depok, and Bekasi (Jabodebek in bahas. This supports the rapid development of the Sukabumi District area. Road access, which continues to be built, supports the growth of foreign and domestic visitors. In this case, the high number of visitors may cause environmental issues . addition, rapid regional development in the CPUGG area can be observed from changes in land use from agricultural land to built-up areas, the rapid development of infrastructure and facilities, and the emergence of various economic activities. If left unchecked, this will result in sustainable development that is difficult to In line with the results of Fauzi and Oktavianus . , economic development tends to sacrifice the environment if development implementation is not comparable across dimensions. Over the last few decades, humanity's overall development has negatively impacted the environment and socio-economic instability, endangering the survival of the Earth and future generations . Therefore, tourist destinations are expected to implement a long-term sustainable tourism strategy by considering economic, social, and environmental factors . According to data from the Regional Development Agency (BAPPEDA) Sukabumi regency . , there are several strategic issues regarding CPUGG management: an increase in the number of visitors and its correlation with ecosystem sustainability, unclear management systems, and a lack of development policy Corresponding Author: Siti Salsiah sitisalsiah97@gmail. Agriculture. IPB University. IPB Dramaga Campus. Bogor. Indonesia. Department of Soil Science and Land Resource. Faculty of A 2025 Salsiah 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! The government must overcome this regional transformation phenomenon in the CPUGG area. Geopark development to benefit geo-tourism must follow sustainable development standards, as directed at the 30th International Geological Congress in Beijing in 1966 . There are several definitions of sustainable One well-known definition was provided by the Brundtland Commission in 1987 . There are three pillars of sustainable development: economy, society, and ecology . Ae. Conceptually, the achievement of sustainable development is the magnificence between those primarily concerned with nature and the environment, those who value economic development, and those dedicated to improving human living conditions . The concept of the development approach is known as sustainable development, which is a development concept that can meet the needs of the current generation without compromising future generations . An assessment of sustainable development must be determined as a measure of sustainability. Research on sustainable development has been frequently discussed at the regional, national, and global levels . but has rarely been discussed at the local level . Specifically, the spatial dependency between locales is greater locally than at the national or global level. The sustainability performance of a site is influenced by its surroundings . The law on village autonomy in Indonesia requires the creation of development indicators at the village level . These indicators serve as measurement tools and recommendations for local government authorities while conducting development projects. This indicates that an analysis of sustainable development at a local level is required. The government should recognize the significance of village development. Consequently, village development must be well planned and based on the findings of a thorough study of all the potential as well as the issues encountered by the village to meet the true requirements of the village community. Against this background, this study aims to determine the local sustainability performance based on the spatial distribution pattern of the local sustainability index (LSI) in CPUGG using a spatial autocorrelation approach. Although several studies have been conducted to measure sustainable tourism in Indonesia, in contrast to this study, the indicators focus on a sustainable development approach. Materials and Methods Study Area This study was conducted in the CPUGG area. Sukabumi Regency. West Java Province. Indonesia, which consists of eight sub-districts and 74 villages with borders with Banten Province in the North. Jampang Kulon. Lengkong. Bantargadung, and Cikidang in the East, the Indian Ocean in the South and West, as shown in Figure 1. Figure 1. Study area. http://dx. org/10. 29244/jpsl. JPSL, 15. | 235 The total area of CPUGG is approximately 126,000 ha, which is 30. 3% of the Sukabumi District. The population of the CPUGG area in 2020 reached 502,019, with an average density of 412 people/km 2. Ciletuh Geopark was designated as a national geopark by the Komite Nasional Geopark Indonesia (KNGI) and the Indonesian National Committee for UNESCO on December 22, 2015. The geopark area was expanded, and its name was changed to Ciletuh-Palabuhanratu, inaugurated by the same committee on June 21, 2016, and designated as a UNESCO Global Geopark on April 17, 2018, with the geopark theme "Subduction Fossil. Plato Jampang. Magmatic Path ShiftAy. The CPUGG offers several advantages regarding geology, biology, and cultural The geological diversity includes subduction fossils. Jampang plateau, and shifting magmatic Biodiversity was found in the Cikongga Wildlife Reserve. Cibanteng Nature Reserve. Turtle conservation, and Halimun Salak National Park. The cultural diversity in this area extends from the megalithic era to the modern era, including megalithic sites in Cengkuk. Kasepuhan Banten Kidul, and historical heritage buildings preserved today. In addition, there is a diversity of intangible cultures, such as dances, songs, and Data Collection In this study, we used secondary data. There are two secondary data sources: village potential data . sourced from the Badan Pusat Statistik (BPS) and CPUGG land-use and administration maps sourced from BAPPEDA. Village potential data were processed to develop indicators . to measure LSI. Land use and administration maps were then processed to identify the spatial distribution patterns of the LSI. Data Analysis The LSI is a local-level development sustainability index created by the author in 2016 for the research region in the Greater Jakarta Metropolitan region (Jabodetabek in Bahas. , with a unit of analysis per village. After evaluation in the Jabodetabek metropolitan region, this index is now being tested in several other places with distinct features and units of analysis, employing a variety of modified variables or indicators. The LSI was calculated using factor analysis (FA) to select the most notable variables / indicators . Ae. Table 1 presents the variables used. FA is a statistical technique used to describe the variability of observed correlated variables in terms of a potentially small number of unobserved variables called factors . The FA model used in this study is expressed as follows: LSIki = Ocycuyco yaycoyco. ycIycoycoycn yco=1 . where LSIki = LSI for k to n dimensions in village i to n. k = dimension . = 1: economy. k = 2: social. k = 3: Ekm = eigen value for the k to n dimension in the m to n factor. Skmi is the factor score for the k to n dimension . actor score from the results of the STATISTICA7 calculation, factor m to n in village i to n. = 1,2,3, . , n. To standardize the LSI value (LSIki . ) on a scale of 0Ae100, we used the following formulation: LSIki. = (LSIki-LSIki. )* yaycIyaycoycn . Oe yaycIyaycoycn . where LSIki. is the LSI standardization for k to n dimensions in villages i to n. LSIki = LSI for the k to n dimensions in the i to n villages. LSIki. is the value minimum LSI for k to n dimensions in i to n villages. and LSIki. = value the maximum LSI for the k to n dimensions in villages i to n . Table 1. List of variables of LSI. Code of variables Economy . Percentage of households using electricity (%) Number of Base Transceiver Station (BTS) tower Number of communication service operator Distance to the closest market . Distance to the closest central business district (CBD) . Number of industries per 1,000 population Number of markets, mini markets, shops per 1,000 population Number of hotels, hostels, motels, and inns per 1,000 population Number of restaurants or food stalls X10 Percentage of developed or urban land per total area (%) This journal is A Salsiah et al. JPSL, 15. | 236 X11 Distance to the closest bank . Social . X12 Number of formal education facilities . indergarten to universit. per 1,000 population X13 Number of informal education facilities . indergarten to universit. per 1,000 population X14 Distance to the closest education facility . indergarten to universit. X15 X16 Number of health facilities . ospitals, clinics, health centers, pharmacie. per 1,000 Number of medical workers per 1,000 population X17 Number of places of worship per 1,000 population X18 Number of non-profit organizations per 1,000 population X19 Distance to the closest health facility . ospital, clinic, health center, pharmac. Environment . X20 Percentage of households living along the river . iparian are. (%) X21 Percentage of households living in the slum area (%) X22 Percentage of water bodies per total area (%) X23 Percentage of forest area per total area (%) X24 Percentage of plantation area per total area (%) X25 Percentage of paddy fields per total area (%) X26 Percentage of other agricultural land per total area (%) Source: . The results of the LSI were subsequently analyzed for spatial autocorrelation, both globally and locally, using the GeoDa Software. The global spatial autocorrelation estimation for spatial information can be calculated using global Moran's index method. The Moran Index is a measure of the correlation between observations in an area and other neighboring areas. A positive spatial autocorrelation indicates similar values from locations that are close together and tend to cluster. An autocorrelation negative spatial distribution indicates that adjacent locations have different values and tend to spread . The range of values of the Moran Index is Ae1 O I O 1, and a value of Ae1 O I < 0 indicates a negative spatial correlation, whereas a value of 0 < I O 1 indicates the presence of positive spatial autocorrelation. The MoranAos Index Value is zero . , indicating no autocorrelation and that it tends to be random. The local Moran's index method, also known as the local indicator of spatial autocorrelation (LISA) approach, can be used to evaluate local spatial autocorrelation in spatial data . A positive spatial autocorrelation occurs when a high value correlates with a high neighbor or a low value correlates with a low neighbor. Negative spatial autocorrelation occurs when high values correlate with low neighbors or low values correlate with high neighbors . Identifying the relationship or spatial association of LSI is also important because there is a spatial interdependency between locations at the local level . ,20,. The relationships between a local observation location and another observation location were determined using LISA. There are several ways to calculate spatial autocorrelation using the local MoranAos index . yaycn = ycycn Oeyc Ocyc ycOycnyc ycyc Oeyc ycycycc Where: Ii : LISA Index yj : value at j-locations wij : spatial weighted matrix yi : value at i-locations std : standard deviations value Results The FA produced the LSI to demonstrate sustainability performance in each village and dimension. There were 26 variables in the 74 villages grouped by three dimensions: economy (LSI. , social (LSI. , and environment (LSI. The distributions of LSI1. LSI2, and LS3 are shown in Figure 2. The color gradation shows that a dark color indicates a village with a high sustainability index value. Black circles on the maps focus on villages in the Palabuhanratu Sub-district. The Palabuhanratu Sub-district is the capital of the Sukabumi Regency and the leading tourist center in the CPUGG. http://dx. org/10. 29244/jpsl. JPSL, 15. | 237 Figure 2. Spatial distribution maps. LSI1, . LSI2, . LSI3. The villages with high LSI1 values are located in the Palabuhanratu Sub-district . n the black circle are. The Palabuhanratu Sub-district is the capital of the Sukabumi Regency. This area is a strategic and developing economic area compared with other regions. The number of infrastructure and economic facilities in the area is greater than in other areas because the Palabuhanratu Sub-district is one of the main points of the CPUGG tourism area. Compared to LSI2 . ocial dimensio. , the villages surrounding the Palabuhanratu Sub-district have a low LSI value, with the fewest formal education institutions. Villages with high LSI2 values were located outside the district capital. Education is a concern because it has a positive effect on sustainable tourism. LSI3 is LSI in the environmental dimension. The villages in the Palabuhanratu Sub-district have a high LSI3 value, with visible stretches of rice fields and farmland. Based on the results, the villages with the highest LSI3 values were generally concentrated in villages with large forests, plantations, and farmland. In this study, we employed spatial autocorrelation analysis using the global and local Moran index / LISA to analyze the spatial distribution pattern of LSI in the CPUGG area. The results of global MoranAos index analysis are shown in Figure 3. Based on the analysis. LSI1. LSI2, and LSI3 are 0. 476, 0. 279, and 0. 269 of the Global MoranAos I, respectively, which shows that the spatial distribution patterns of LSI1. LSI2, and LSI3 are clustered. Figure 3. Global Moran Index of . LSI1, . LSI2, and . LSI3. This journal is A Salsiah et al. JPSL, 15. | 238 In this study, we used GeoDA software to identify the spatial distribution pattern of the LSI using the Local Moran Index. The spatial distribution pattern of LSI appears when the spatial association (SA) of its data is positive (High-High or Low-Low typ. or negative (High-Low or Low-High typ. A positive SA occurs when a high value correlates with a high neighbor, or a low value correlates with a low neighbor. The results are shown in Figure 4. Not significant . High- High Low- Low . Low- High High- Low Figure 4. Local Moran Index . LSI1, . LSI2, and . LSI3. The LISA analysis results in Figure 4 show that the LSI values in the economic dimension have a positive spatial autocorrelation . igh-hig. , specifically in the village around the district capital of the Sukabumi Regency (Palabuhanratu Sub-distric. Palabuhanratu is the leading tourist centre in CPUGG, with the highest concentration of tourism and hospitality facilities and infrastructure. This HH-type value indicates that the LSI1 value in the village is high and surrounded by villages with high LSI1 values. Compared to LSI2 . ocial dimensio. , the villages surrounding the Palabuhanratu Sub-district have a positive spatial autocorrelation . ow-lo. This LL-type value indicates that the LSI2 value in the village is low and surrounded by villages with low LSI2 values. Villages in this area have a lower number of social facilities compared to other villages. The LSI3 values in the villages surrounding the Palabuhanratu Sub-district had a positive spatial autocorrelation . igh-hig. This HH-type value indicates that the LSI3 value in the village is high and surrounded by villages with high LSI3 values. LH-type village have more potential to increase their sustainability index compared to villages with HL-type. Villages with a low sustainability index can be encouraged and influenced by villages with a high sustainability index around them. Discussion Based on the Local Sustainability Index analysis, the villages surrounding the Palabuhanratu Sub-District have a high economic and environmental sustainability index but a low social sustainability index. The high number of economic facilities in the main points of the CPUGG tourism area is essential because it can accelerate economic circulation, increase and equalize the actors' income, create job opportunities, and increase business opportunities or job diversification . Unfortunately, the number of formal education surrounding the Palabuhanratu Sub-district is low. So, increasing the number of educational facilities in the area is essential. However, education is a concern because it positively affects sustainable tourism. Education is necessary because it can produce high-quality human resources and has good knowledge and understanding of sustainable tourism . The high quality of the local community has become an essential component of sustainable tourism because local communities play a crucial role in developing tourist villages as they depend on their level of acceptance and support . http://dx. org/10. 29244/jpsl. JPSL, 15. | 239 This study employed spatial autocorrelation analysis using the global and local Moran index / LISA. The results show that the spatial distribution pattern of the local sustainability index in economic, social, and environmental are clustered. Each clustering shows that each village influences the surrounding villages. The results align with the spatial distribution of LSI discussed previously . However, using the local Moran index analysis, we can find how the spatial association of LSI in the CPUGG is more specific and more clearly defined because the sustainability index in a particular location is affected by the sustainability condition in its surrounding areas . Overall, the CPUGG development is not evenly distributed yet. So, increasing the economic, social, and environmental sustainability index in several villages with a low sustainability index is Good cooperation is needed between the government, society, and several stakeholders to create sustainable CPUGG development for the present and future. An interesting finding in this research is that the environmental sustainability index in the surrounding Palabuhanratu Sub-district is high. Economic growth in the region does not interfere with environmental The extensive forests, plantations, and farmland are beneficial because they provide landscape services essential for ecosystem survival . A sustainable environment can be seen in the region's preserved vegetation while the economic and social sectors are growing. Sustainability is the concept of development that improves the well-being of people . eople and interest. without harming the environment . and supports the lives of future generations . This phenomenon differs from previous research . , which shows a trade-off between economic and environmental dimensions. The statement in their research states that economic development will be able to increase the growth of a region very well, but this also has negative impacts, one of which is the high demand for land. However, this research shows that economic progress in the Palabuhanratu area, a tourist center, still maintains the environment The analysis of local sustainability index in this study can be used not only for tourist areas. However, it can also be used to see the local sustainability index in other areas, both urban and rural because the methods and variables used are very adaptive, so they can be adjusted to certain regional / local conditions. This study calculates the local sustainability index based on village potential data without considering its status or it would be better if, in the future, the calculation of the number of economic and social facilities is carried out, taking into account the feasibility conditions of the facilities. So that the study results will be more precise and government recommendations or policies will be more appropriate. Conclusions This study detected the sustainability performance in the CPUGG area at the local and village levels. Thus, we assessed the differences in each village's economic, social, and environmental conditions in the CPUGG area. The LSI in the CPUGG area was still not evenly distributed based on the results of the global and local Moran index analyses. The spatial autocorrelation results show that based on the global Moran index, the spatial distribution patterns of LSI1. LSI2, and LSI3 are clustered. The local Moran index shows that spatial interdependency of local sustainability performance exists between local-scale spatial units and that the sustainability state of the surrounding areas influences the LSI of a village. An interesting finding of this research is that the villages around Palabuhanratu, which have a high economic index, have a high environmental index value compared to other villages. This shows that economic progress in the Palabuhanratu area, a tourist center, is carried out while still paying attention to environmental aspects. Author Contributions SS: Conception, and Design of study. Acquisition of data. Analysis and/or interpretation of data. Drafting the manuscript. Critical review - revision. IS: Analysis - Interpretation of data. Critical review - revision. ES: Analysis - Interpretation of data. Critical review - revision. and AEP: Conception. Design of study. Analysis Interpretation of data. Critical review - revision. Conflicts of Interest There are no conflicts to declare. This journal is A Salsiah et al. JPSL, 15. | 240 Acknowledgements The authors are sincerely grateful to IPB University and Kyoto University for supporting the research in the Double MasterAos Degree Program and would like to thank the BAPPEDA Sukabumi regency for providing the References