Next Article in Journal
Mapping Urban Spatial Structure Based on POI (Point of Interest) Data: A Case Study of the Central City of Lanzhou, China
Next Article in Special Issue
Combining Satellite Remote Sensing and Climate Data in Species Distribution Models to Improve the Conservation of Iberian White Oaks (Quercus L.)
Previous Article in Journal
Open Geospatial Software and Data: A Review of the Current State and A Perspective into the Future
Previous Article in Special Issue
Urban Ecological Corridor Network Construction: An Integration of the Least Cost Path Model and the InVEST Model
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

GIS-Based Assessment of Habitat Networks for Conservation Planning in Kas-Kekova Protected Area (Turkey)

1
Ministry of Environment and Urbanization, General Directorate for Protection of Natural Assets, Cankaya, 06530 Ankara, Turkey
2
Department of Geomatics Engineering, Konya Technical University, 42075 Konya, Turkey
3
Department of Geomatics Engineering, Hacettepe University, 06800 Ankara, Turkey
*
Author to whom correspondence should be addressed.
ISPRS Int. J. Geo-Inf. 2020, 9(2), 91; https://0-doi-org.brum.beds.ac.uk/10.3390/ijgi9020091
Submission received: 23 December 2019 / Revised: 19 January 2020 / Accepted: 27 January 2020 / Published: 1 February 2020
(This article belongs to the Special Issue Application of GIS for Biodiversity Research)

Abstract

:
The determination of protected area (PA) boundaries and the level of restrictions is very important for sustainable conservation, and such decisions must involve biodiversity data and assessment. In a PA, the consensus of the government and the locals is crucial to ensure protection–use balance. The PA restrictions constrain legal human activities, and the boundary determination and the restrictions should be based on various scientific analyses to achieve consensus. In this study, a GIS-based approach is proposed to utilize the biodiversity data for efficient conservation and land use planning in Kas-Kekova PA, which is among the most important PAs in Turkey. Spatial analysis methods, i.e., kernel density estimation, natural breaks classification and integrated density index, were performed for the assessment of the habitat networks using georeferenced biodiversity datasets, and the results were evaluated with respect to the actual land use data and the land ownership pattern. The developed spatial analysis approach is efficient to produce the conservation base maps required for regional land use planning, for defining sustainable conservation strategies, and to provide a widely accepted base for land use planning and biodiversity monitoring in the PA; although careful investigations and expert opinions are still required for data deficient areas.

1. Introduction

In the Earth Summit on Development carried out in Rio de Janeiro in 1992, it was accepted that the loss of biodiversity is an important issue that can be prevented by coordinated global efforts. The Convention on Biological Diversity (CBD) was adopted for “sustainable development” in the same year [1]. The Convention mainly focused on the conservation of biological diversity. Despite all efforts, the loss of biodiversity continues, and both climate change and environmental pollution play significant roles in this process [2,3,4]. Therefore, integrated global efforts should be increased to protect natural habitats and to preserve them for future generations.
Many species are in danger of extinction on Earth. The importance of nature protection is emphasized by many national and international organizations. As a leading global environmental authority, the United Nations Environment Program (UNEP) includes the conservation and sustainable use of biodiversity in the global environmental agenda. One of the global biodiversity conservation goals, Aichi Target 11, defined in the CBD strategic plan 2011–2020, has the aim that “By 2020, at least 17% of terrestrial and inland water areas and 10% of coastal and marine areas, especially areas of particular importance for biodiversity and ecosystem services, are conserved through effectively and equitably managed, ecologically representative and well-connected systems of protected areas and other effective area-based conservation measures, and integrated into the wider landscapes and seascapes.” [5]. According to [6], there were a total of 202,467 terrestrial and inland water protected areas (PAs) recorded in the World Database on Protected Areas (WDPA) by 2016, covering 14.7% (~20 million km2) of the world’s extent of these ecosystems. Although the value falls short of Aichi Target 11, the complexity of the problem and the amount of the geospatial data involved in the processes suggest efficient utilization of geographical information systems (GIS) and related technologies.
Designation of a conservation area as PA by governments does not stop the loss of biodiversity by itself [5]. There are many issues that should be organized in a PA to implement a successful plan [7,8]. In addition, strict regulations based on ecological networks must be determined [9,10,11]. Essentially, effective protection can only be performed at a local scale. Environmental planning helps to determine the suitable land use by clarifying the biological, physical and social systems of a PA [12]. The first step to understand ecological processes is to identify the ecological patterns. For this purpose, the spatial distribution and densities of flora and fauna elements can be analyzed using spatial statistics [13]. For example, plant species with sparse distributions can be expressed with discrete points [14]. As output of spatial analysis, habitat patterns can be obtained and assessed for defining the regional ecological network. Such a network can be used for improving the habitat connections and proposed as a solution to the land fragmentation problems while ensuring the protection of threatened species and biodiversity [15,16]. As part of an ecological network, ecological corridors and their preservation are also important for sustainable protection. The determination of ecological corridors has important objectives, such as creating a scattered habitat, compensating for habitat losses, and ensuring the protection of endangered species and populations [15].
Several methods were proposed in the literature to identify ecological networks using spatial analysis methods. Opdam et al. [17] proposed a spatial indexing approach to analyze the conservation potential of landscapes for different species. McHugh and Thompson [18] studied an ecological network model for conservation planning using eco-profiles of species. Ferretti and Pomarico proposed an effective tool for spatial planning by integrating GIS with multi-criteria analysis [19]. Gurrutxaga et al. [20] proposed a least-cost path analysis for the determination of ecological corridors based on core habitat zones. Guo and Liu [21] and Santos et al. [22] also proposed an ecological network determination approach based on least-cost analysis at regional scales. Chang et al. [23] proposed a green infrastructure planning approach-based ecological connectivity assessment in Longgang District of Shenzhen in China. Xie et al. [24] proposed an integrated index which evaluates the biodiversity conservation and the disaster protection at the same time. Hong et al. [25] also developed an index system based on ecological sensitivity analysis. The kernel density estimation (KDE) method has been employed very often in the literature to assess the spatial patterns of habitats (e.g., [26,27,28,29]); although its use for evaluating other discrete point pattern analysis applications are also available (e.g., [30,31]). However, the determination of ecological networks is a complex problem and an active research field.
In this study, a GIS-based approach for the determination of habitat networks in Kas-Kekova PA, which is located in the south of Turkey, was investigated. The PA contains diverse endemic flora and threatened fauna species. In addition, there are archeological PAs inside the region. The PA is under threat from biodiversity deterioration or disappearance due to the anthropogenic pressures such as tourism and rapid urbanization. In the area, the use of resources by humans should be managed through sustainable planning and, at the same time, their interest in protecting the biodiversity and the environment should be ensured. An integrated planning approach can assist the development and implementation of conservation principles and support the needs of locals. Therefore, the main goals of the present study are; (a) to develop an easy-to-use approach for producing the base maps to present and involve the biodiversity data for defining land use restrictions in Kas-Kekova PA; (b) to integrate the threatened biodiversity elements (discretely collected flora and fauna data) in land use planning; (c) and to propose additional conservation principles on top of the legal principles defined by the Ministry of Environment and Urbanization (MoEU) based on the spatial analysis results of the biodiversity data, the actual land use and the property ownership pattern. In addition, utilizing a GIS platform to store and present diverse data can help to increase the public and government awareness and to support the decision making.
The article is structured as follows. The main characteristics of Kas-Kekova PA are described in Section 2. The spatial analysis methods and the datasets are provided in Section 3. In Section 4, the results are provided and integrated with the land ownership and land use data, and conservation strategies for land use planning are proposed. The main findings of the study are discussed in the final section.

2. Main Characteristics of the Kas-Kekova Protected Area

Kas-Kekova Region is located on the coast of Antalya Province, Turkey, between the towns of Kas in the west and Kale in the east (Figure 1). The study area covers a total area of approximately 260 km2 (~100 km2 terrestrial and ~160 km2 marine parts). Three different PAs have been declared in the region by two governmental organizations (MoEU and Ministry of Culture and Tourism) based on different laws (Articles 644 and 2863 in [32,33]). The study area was declared as a Highest (First) Degree Natural Protection Site by MoEU in 1989, then as a PA in 1990 under the UNEP and the Barcelona Convention [34]. The PA also contains archaeological ruins, mostly from Lycians, and coastal landscapes that are subject to protection (Figure 1).
The Kas-Kekova PA Biological Diversity Project was carried out in the study area between 2008 and 2010 by the MoEU [36]. Within the project, European Union Nature Information System (EUNIS) habitat classes, plant species, mammalian, amphibian, reptile and bird species in the PA were determined according to the EU Habitats Directive [37] and EU Birds Directive [38]. The extensive field works held by specialists throughout the project (four consecutive spring and fall seasons in two years) ensured the representative data collection for biodiversity. The PA also has high accessibility by foot and, thus, the whole area could be wandered by the ecology experts in the project time frame. In order to cover the whole PA, the region was divided into 35 square grids with 2.5 × 2.5 km grid size and each grid was numbered [36]. During the eight field campaigns between October 2008 and June 2010, the habitat classes were assessed according to EUNIS. In addition, the field work results were supported with the satellite imagery for accurate delineation of vegetation class boundaries. The endemic and threatened flora were identified during the field works and their locations were measured with ca. 1 m positional accuracy using hand-held GPS (global positioning system) devices. Although mammals and reptiles were observed by the experts, their locations have a lower spatial accuracy since the measured locations were mostly the positions of the observer with an approximated correction. Still, a spatial accuracy of better than 10 m can be expected. The birds’ locations have bigger measurement ambiguity, since the identifications were performed by the observers and by voice recordings, and the positions mostly represent the locations of the recorders and observes. Thus, for the birds, a spatial accuracy of better than 50 m can be expected.
A total of 51 families, 187 genera and 272 flora species were identified in the terrestrial habitats of the Kas-Kekova PA. Among those, 26 species are endemic to Turkey. The species were also assessed according to the International Union for Conservation of Nature (IUCN) Red List Categories and Criteria [39,40], which is an important indicator of the sustainability of the world’s biodiversity. The IUCN categories of threat are listed as; extinct (EX), extinct in the wild (EW), critically endangered (CR), endangered (EN), vulnerable (VU), near threatened (NT), and least concern (LC). The species in the PA fall into the categories of EN, VU, NT, and LC. Those flora and fauna species which were classified as CR, EN, VU and NT (Table 1) were employed in the spatial analysis algorithms in this study, since they are threatened and their habitats are subject to protection. The reptiles and amphibian data were grouped and mentioned as reptiles in the following sections, since few amphibians were identified in a very small area in the PA.
Kas-Kekova PA also has four major plant communities [36,41], which are pine forests (Aetheorhizo Bulbosae–Pinetum Brutiae), maquis (Quercus aucheri–Oleetum europaeae), phrygana (Alysso-Genistetum acanthocladae) and halophyte (Salicornietum ramosissimae). The maquis community is the main plant association in the PA [36]. This community has a rich and stable biodiversity, includes endemic species, and provides the climax phase and the floristic composition in the PA. The phrygana community was developed as a result of the burning of the maize vegetation in the region and is represented in smaller parts of the PA. The halophyte community develops only in the salt marsh and a very small halophyte plant region exists in the PA, which occurred by the tidal movements coming from the sea. The forest community is mainly composed of red pine and was determined as the union of Aetheorhizo bulbosae and Pinetum brutiae although the Pinus brutia is very rare. These communities were integrated into the land use map, which was generated according to the EUNIS classifications.
There are also anthropogenic affects in Kas-Kekova and several rural settlements in and around the PA exist (Figure 2). The human population inside the PA is ca. 2000 (as of year 2013). Since the PA has a stunning coastal landscape, the tourism activities in the summer season are high and double the population. The daily tourism activities are mostly in the bays and may expand to spring seasons as well. The main income sources of the locals are agriculture, tourism and fishing [36]. Since the 1990s, touristic accommodations and roads in the PA have increased rapidly (Figure 3). Due to improper infrastructure planning, marinas and hotels prompt illegal construction, which causes the uncontrolled use of natural resources [36]. Although the PA has a master land use plan with a scale of 1:25.000 from 1991, this plan was not based on biodiversity data and does not comply with the human needs either, despite the fact that 20% of the cadastral parcels were registered under private ownership. Since the economic expectations are more important than environmental protection for locals [42], the conservation strategies must be developed to obtain a protection–use balance.

3. Materials and Methods

In order to integrate the threatened biodiversity elements into the final land use plans, the spatial habitat patterns of these elements (i.e., endemic and threatened flora and fauna species) should be revealed. This task can be achieved with spatial analysis methods although several complexities should be considered, such as; (i) the biodiversity is high (i.e., many species types exist) and performing extensive habitat analysis for every species type in the PA is extremely difficult; (ii) only a small number of observations for the threatened species is available and the habitat pattern may not be revealed adequately; and (iii) there are large variations between the numbers of observations for different species types. Various spatial analysis approaches can be utilized for habitat pattern determination purposes, but special care must be taken when selecting the method and its parameters. For rarely observed and threatened fauna species, a home range estimation method can be used but it may require exhaustive field observations over a long term. If the number of observations is sufficient for demonstrating the species spatial distributions, a statistical density estimation method can be employed by considering the distribution pattern (e.g., random or clustered). In this study, an easy-to-use methodology was developed for the spatial habitat pattern determination to be used as a base for land use planning in the PA, where the different elements were grouped under four categories (i.e., mammals, reptiles, birds and endemic flora) and every observation in each group was represented in the final map properly (e.g., for the threatened mammals) (Figure 4).
First, the flora and fauna point data collected in the fieldwork [36] were analyzed to determine the patterns of plant and animal ecology and to produce ecological sensitivity maps [43,44]. After the initial investigations on the biodiversity reports (i.e., endemic species, habitat types, plant community, threatened categories, location etc.), a geodatabase was populated with spatial and attribute data. The flora and fauna points were filtered for endemism and/or threatened category (EN, VU, NT) attributes and divided into four groups; namely, flora, mammals, reptiles (including amphibians), and birds. After the visual investigations on the filtered point locations, overlaid with the EUNIS habitat map and the topography, a spatial methodology was developed. The KDE method was used for defining the point patterns and an optimal bandwidth determination approach was developed according to the spatial statistics of the points. The natural breaks classification was used to obtain the spatial configurations of each group. The integrated density index was used to combine the patterns of all four groups and thus to obtain a habitat network. The final results were analyzed together with the EUNIS land use classes, land ownership data (cadastral data), and the roads. Several conservation strategies for the PA management are proposed based on the integrated analysis results. The spatial analyses were performed using ArcGIS 10.2 software from ESRI, Redlands, CA, U.S.A. [45].

3.1. Spatial Datasets

Kas-Kekova PA has fourteen EUNIS habitat classes (Table 2) defined according to EUNIS Habitat Classification [46]. The class boundaries were delineated manually using satellite images and inspected by field work [36]. The formats, sources, and the accuracies of datasets evaluated in the study are given in Table 3. Figure 5 shows the EUNIS Habitat Map of the PA and the locations of threatened fauna and endemic flora. The region is covered largely by forests (82.4%). In Figure 6, the cadastral parcels and their ownerships are demonstrated. The ownerships are categorized as; (i) private ownership; (ii) private ownership status obtained after deforestation in the area, (iii) parcels owned by the state, and (iv) state forest areas that are not subject for private ownership.
The topographical characteristics of the PA were analyzed in terms of slope, aspect and elevation values and investigated with respect to mammals, reptiles and endemic flora species. The slope values in the PA range between 0–33 degrees (Figure 7 and Figure 8) and the elevation values range between 0–1170 m. Although topography may form barrier to biodiversity elements, considering the spatial distributions of the mammals (i.e., Capra aegagrus Erxleben) and the reptiles locations identified in the PA, the topography was not found to be an effective factor or barrier (Figure 7 and Figure 8). Therefore, it was not considered in the ecological network assessment. In addition, no clear indications of spatial correlation with the topographical characteristics (slope, aspect and elevation ranges) and the flora data were observed as shown in Figure 9.

3.2. Spatial Analysis Methods

The KDE method has often been used in the literature for the density estimation of point patterns (e.g., [26,27,28,29,30,31]). It is a statistical method developed by Silverman [47] and acts as a smoothing filter for discrete points (Equation (1)). The KDE method is based on the value of i at each sampling point (x, y) in a kernel. The density contribution of each sampling point was calculated using the kernel functions. The density value changes depending on the search radius range (bandwidth) and grid cell distance in each sample point. The method is already implemented in ArcGIS 10.2 from ESRI, Redlands, CA, U.S.A. and uses a quartic kernel, as proposed by Silverman [47]. In the search radius range, the sampling points closer to the grid cell center have a higher density contribution value. A cumulative density value for all kernels in the entire working area are computed using the method.
D e n s i t y = 1 r a d i u s 2 i = 1 n { 3 π p o p i [ 1 d i s t i r a d i u s ] 2 } 2 ; For   d i s t i < r a d i u s
where;
i = 1,…,n are the input points. Only include points in the sum if they are within the radius distance of the (x, y) location.
popi: population field value of point I
disti: the distance between point i and the (x, y) location.
The parameter selection for KDE (e.g., kernel type and radius) was investigated in several studies (e.g., [48,49,50,51,52]) and yet there is no standard approach for this purpose. In this study, a spatial statistical analysis, average nearest neighbor (ANN) [53] implemented in ArcGIS 10.2, was used to analyze the distributions of the points and an approach based on these results was developed for the radius selection. The method computed the average nearest neighbor ratio and expected mean distances (DE) and observed mean distances (DO), as shown in Equation (2).
D o = i = 1 n d i n ; D E = 0.5 n / A ;   A N N = D o D E
where;
D O is the observed mean distance between each feature and its nearest neighbor.
D E is the expected mean distance for the features given in a random pattern.
di equals the distance between feature i and its nearest neighbor.
n corresponds to the total number of features.
A is the area of a minimum enclosing rectangle around all features, or it is a user-specified area value [52].
If the ANN index is less than 1, the distribution is clustered. If the index is greater than 1, the points are dispersed. The results of the ANN analysis for each group are provided in Table 4. The flora points are clustered and the DE represents the mean sizes of the clusters. Due to the random distribution of bird and reptile points, the DO values are considered to be representative. On the other hand, there are only five observations (points) for mammals and they are dispersed. Since there is only one species type (Capra aegagrus Erxleben) in this group, its home range was obtained from [54] as 150–400 ha. It should be noted that identifying home ranges precisely is an exhaustive process and a simple and easy-to-use approach is needed for planning purposes. In this study, a 1000 m kernel radius was selected for all groups based on the DO and DE values, and the home range obtained from [54]. In addition, during the initial investigations on kernel radius, it was observed that using a smaller radius results in a disconnected and highly fragmented ecological network. On the other hand, using a much larger kernel radius would cause over-smoothing of the data. A 1000 m radius was found to be visually optimal and it also results in a ca. 3.14 ha home range for the mammals. The following strategy can be derived from the current implementation based on the ANN analysis:
If distribution pattern is clustered → check DE
If distribution pattern is random → check DO
If distribution pattern is dispersed → check home range for fauna (if possible) or EUNIS habitat class for flora.
The natural breaks classification is used to reclassify the output KDE values of each group into five categories. The method implemented in ArcGIS 10.2 uses the Jenks optimization algorithm [55] that minimizes the squared deviations of the class means. After this step, the classes were considered as habitat zones, as shown in Table 5. The four habitat classes were combined using an integrated density index approach to obtain a joint habitat pattern. The approach used here compares the four class values (flora, mammals, reptiles and birds) for every grid point and takes the maximum class value as output. Although other integration approaches (e.g., averaging the class values, calculating the sums of all class values for every pixel and reclassifying into five classes again, etc.) could be used; the proposed approach allows us to represent the importance level of species when having a small number of samples (e.g., mammals) better and therefore is preferred.
The results of different bandwidth tests with 500 m and 800 m for endemic flora and a single class for fauna data (instead of three different classes for mammals, reptiles and birds) are presented in Figure 10. These radii cause disconnected and fragmented ecological network. On the other hand, as can be seen from Figure 10b, the core zones of threatened mammals (denoted with red triangles) cannot be represented in the final map when all fauna data are grouped together. The KDE analysis results of all four groups classified with natural breaks are provided in Figure 11. The same anomaly occurs if other integration approaches (e.g., averaging the class values, etc.) were applied and therefore using maximum class value allows the optimal representation of the core zones, buffer zones, etc., as presented in the following section.

4. Results

4.1. The PA Network

The integrated results clipped using the coastline and the PA boundary are presented in Figure 12. Unlike using discrete points for demonstrating the habitats, the results indicate a continuous pattern as expected, and show that the PA must be protected against habitat fragmentation and losses caused by it [56,57]. The five levels were interpreted as core habitat zones, buffer patches, ecological corridors, stepping zones and isolated elements for endemic and threatened species (Table 5). It should be noted that the collected data are still likely incomplete, due to the time span of the field campaigns, although the species were precisely identified and their locations were accurately measured [36]. Undersampling is an expected limitation of biodiversity determination projects due to the data collection methods (e.g., often manually by human interpretation) and the discrete nature of the measurements (discrete time intervals). The results mainly point out the areas for strict protection and the remaining parts should be analyzed by experts, if possible, and interpreted carefully by planners.

4.2. Interactions between the PA Network and the Human Use

The property ownership rights should be taken into account in a PA for effective protection. The private property areas are of particular interest, since protection should be carried out with the mutual agreement of the landowners. The protection–use balance can be achieved with the voluntary cooperation of the stakeholders [58,59,60,61,62] and the conservation strategies should be developed by social and financial analysis of the land property ownership and land use patterns. The distributions of different ownership types (i.e., private properties, private properties obtained after deforestation, forests under state ownership, and non-forest state properties), as shown in Figure 6, demonstrate the serious fragmentation of land registered under private ownership, both private parcels and newly registered parcels after deforestation. A major conservation strategy can be expropriation by the government in strict protection areas. The private ownership percentage (~20%) given in Table 6 also indicates the need for expropriation for nature conservation. For this purpose, a preliminary analysis was made using the privately owned parcels, the habitat network obtained from the developed methodology, and the EUNIS classes. The analysis extracts the privately owned parcels from the dataset that fall inside the core habitat zones and the buffer patches (red polygons in Figure 13). However, parcel areas that fall into the EUNIS classes of residential areas (villages), cemeteries, agricultural constructions and mixed crops of market gardens and horticulture were excluded from the analysis. According to the national regulations [33] for the use permissions in unplanned areas, the borders of these four classes were expanded by 100 m prior to the analysis. Although the results given in Figure 13 show the areas that need to be expropriated, the remaining parcels should still be examined by experts for a final decision, due to the complex nature of the biodiversity data.
Infrastructure facilities in rural areas provide economic development but may destroy natural assets [63] and threaten the biodiversity of PAs [64,65,66,67]. Roads affect the PAs by their width, traffic level and speed [68]. Amphibians, reptiles, birds and mammals are mostly affected by the roads [69,70]. The impact of roads on the ecological network is in the range of 100–600 m [65] and according to different road types (provincial, district and rural roads), the impact area can be determined as 500, 100 and 50 m, respectively [71,72]. The roads in the PA were buffered using different radii (500, 100 and 50 m) based on their type (province, district, rural) and visualized together with the determined ecological network for defining the conservation strategies and to support the decision making for the construction of new roads. Figure 14 shows the road impact zones (buffers) and settlement areas overlaid with the determined ecological network. From the Figure 14, it can be seen that they have impacts on the core habitats and the buffer patches.

4.3. Conservation Strategies

The strategies reported by biodiversity researchers [36] mainly include observations on biodiversity elements and actual land use within different parts of the PA with qualitative geographical definitions (e.g., west of Kekova Island, west of the PA, etc.). Although these recommendations describe the biodiversity and land use status, the PA management strategies and regulations are mainly defined by the MoEU. The strategies and regulations can be summarized as:
  • A PA should be determined based on the assessments of ecological, geological, geomorphological, hydrological and landscape characteristics in the area.
  • The PA strategies, which should support the decision makers, must be based on qualitative and quantitative assessment results.
  • The biological/ecological importance of a PA should be exhibited with the existence of critical species and the actual status and pattern of their habitats. The main aim of biodiversity studies carried out under the direction of the MoEU is to determine the biodiversity (flora and fauna) and habitat characteristics. The ecological network should also be demonstrated with such studies.
  • The whole ecosystem should be assessed together with the spatiotemporal distributions of critical species and their habitats, and also their relations.
  • Thus, planning and protection can be based on ecological studies and scientific data.
Here, these principles were expanded for Kas-Kekova PA based on the outcomes of this study. After the determination of the ecological characteristics and the PA habitat network, conservation–use principles for the protection of endangered species and sustainable land use could be proposed. First, joint efforts were carried out with the government authorities responsible for the development of the Kas-Kekova Conservation Plans. For sustainable land use, the following conservation strategies can be proposed here:
  • The ecological network should be integrated into the land use plans. Strict protection measures, such as no construction, must be taken in the core habitat zones. Restoration of the core habitats should be improved.
  • Buffer patches (zones) are helpful to maintain the core zones and they contain important habitats. Buffer patches ensure distinctive habitat isolation and must be defined as strictly PA as well.
  • Infrastructure facilities in rural areas provide economic development, but easily lead to the destruction of natural assets [58,59,60,61,62]. To provide effective land management in the PA, the expansion of the settlements should be prevented.
  • Land consolidation should be carried out in the PA to support the protection of the habitat and obtain an optimal land use pattern. Existing orchards can be supported with ecological agriculture and used as connecting elements for the other zones, especially for endangered species and threatened habitats.
  • The Government should expropriate the private property areas that are located in the core habitats and the buffer patches. If the expropriation process cannot be executed rapidly, their owners can be supported by the social system to be able to prohibit any economic activity (e.g., tourism, farming, etc.) or the farmers can be trained and encouraged for ecological agriculture activities to protect and enhance the environment on their farmland.
  • The roads in the PA should be analyzed based on their type (province, district, rural) and compared with the determined ecological network for defining the conservation strategies and supporting the decision making for the construction of new roads and the reconstruction of existing ones.

5. Conclusions and Future Works

A quick and feasible method based on KDE and the integrated density index to define the habitat patterns, buffer zones and ecological corridors of threatened species in Kas-Kekova PA, Turkey was investigated in the present study. According to the results of spatial analysis and interpretation, a protection strategy was developed and proposed as a comprehensive approach for effective land management. From the results, it can be concluded that the KDE is powerful to derive a connected habitat network from discrete points of threatened species locations, although the selection of a kernel radius that may be suitable for diverse species types can be challenging. Here, ANN was employed to analyze the point patterns (random, clustered, disperse) and a kernel radius selection approach is proposed based on the results. However, the ANN method may not be suitable with few samples (e.g., for endangered mammals) and more expertise and field studies are needed for the acquisition of more knowledge about the threatened species and thus about the kernel radius selection. The integrated density index approach was used for combining the KDE results of different point groups (i.e., mammals, reptiles, birds and flora) and using the maximum input value for each output index value avoids over-smoothing and enables the representation of the importance value (e.g., core zones) of species with small numbers of samples. The output of the methodology is a base map which demonstrates the ecological network of threatened flora and fauna that can connect the fragmented habitats and be used to analyze the habitat patterns.
Nature conservation activities are often limited by the property ownership rights and the actual human use (e.g., roads) in PAs. As stated in the conservation strategies, expropriation and land consolidation tools should be developed by the local authority to improve this situation. The biodiversity studies must be integrated into the regional land use plans in Pas, since the land use plans are the most important means for sustainable nature conservation. For effective land management in a PA, a spatial environmental plan must be prepared by giving the priority to the biodiversity reports. If the priorities are not clarified with biodiversity maps, urban planners could easily get away from the protection purpose or cannot determine the land use boundaries in spatial environmental plans.
Despite the importance of biodiversity determination projects in PAs, they require a high level of expertise, extensive manual labor, and field works which should be spread over many seasons and years. On the other hand, with the recent advances in spatial data collection, interpretation and analysis methods—such as volunteer geographical information (VGI), citizen science (CitSci), remote sensing, new sensor technologies, and GIS as data storage, analysis, sharing and visualization platform—the task becomes easier. The increase in the studies on biodiversity and land use monitoring with remote sensing techniques indicates that PA planning and management will also be supported better with the help of geospatial technologies. Future work in Kas-Kekova PA and other PAs in Turkey should involve the development of a system-based monitoring approach complementary to the biodiversity determination projects.
Kas-Kekova PA is also a marine protected area. In this study, only the terrestrial habitats of the PA have been assessed. In future works, the analyses should be expanded to the marine part of the PA and a combined set of conservation strategies should be developed for the whole region, and the terrestrial and marine PAs should be evaluated together.

Author Contributions

Conceptualization, methodology, D.T. and S.K.; formal analysis, investigation, data curation, writing—original draft preparation, D.T.; validation, writing—review and editing, S.K.; supervision, S.I. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

The authors thank to the Ministry of Environment and Urbanization, General Directorate for Protection of Natural Assets for provision the flora, fauna and basemap data for the Kas-Kekova protected area. The authors would also like to thank anonymous reviewers for their constructive comments and suggestions, which greatly helped to improve the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. United Nations Secretariat, C. B. D. Convention on Biological Diversity. 1992. Available online: https://www.cbd.int/doc/legal/cbd-en.pdf (accessed on 29 January 2020).
  2. Bosso, L.; Ancillotto, L.; Smeraldo, S.; D’Arco, S.; Migliozzi, A.; Conti, P.; Russo, D. Loss of Potential Bat Habitat Following a Severe Wildfire: A Model-Based Rapid Assessment. Int. J. Wildland Fire 2018, 27, 756. [Google Scholar] [CrossRef] [Green Version]
  3. Marshall, L.; Biesmeijer, J.C.; Rasmont, P.; Vereecken, N.J.; Dvorak, L.; Fitzpatrick, U.; Francis, F.; Neumayer, J.; Odeqaard, F.; Paukkunen, J.P.T.; et al. The interplay of climate and land use change affects the distribution of EU bumblebees. Glob. Chang. Biol. 2018, 24, 101–116. [Google Scholar] [CrossRef] [PubMed]
  4. Skogen, K.; Helland, H.; Kaltenborn, B. Concern about Climate Change, Biodiversity Loss, Habitat Degradation and Landscape Change: Embedded in Different Packages of Environmental Concern? J. Nat. Conserv. 2018, 44, 12–20. [Google Scholar] [CrossRef]
  5. CBD, UNEP. Strategic Plan for Biodiversity 2011–2020 and the Aichi Targets. Available online: https://www.cbd.int/kb/record/decision/12268 (accessed on 23 December 2019).
  6. UNEP-WCMC. IUCN Protected Planet Report 2016. Available online: https://wdpa.s3.amazonaws.com/Protected_Planet_Reports/2445%20Global%20Protected%20Planet%202016_WEB.pdf (accessed on 23 December 2019).
  7. Stanciu, E.; Steindlegger, G. RAPPAM. Rapid Assessment and Prioritization of Protected Area Management—Methodology Implementation in Romania: Key Findings and Results. Available online: https://www.protectedplanet.net/c/protected-areas-management-effectiveness-pame/reports (accessed on 23 December 2019).
  8. Tezel, D.; Buyukdemircioglu, M.; Kocaman, S. Accurate Assessment of Protected Area Boundaries for Land Use Planning Using 3D GIS. Geocarto Int. 2019, 1–14. [Google Scholar] [CrossRef]
  9. Ferrier, S.; Funk, V. Mapping Spatial Pattern in Biodiversity for Regional Conservation Planning: Where to from Here? Syst. Biol. 2002, 51, 331–363. [Google Scholar] [CrossRef]
  10. Pressey, R.L. Priority conservation areas: Towards an operational definition for regional assessments. In National Parks and Protected Areas: Selection, Delimitation and Management; Pigram, J.J., Sundell, R.C., Eds.; Centre for Water Policy Research, University of New England: Armidale, UK, 1997; pp. 337–357. [Google Scholar]
  11. Reid, W.V. Beyond protected areas: Changing perceptions of ecological management objectives. In Biodiversity in Managed Landscapes: Theory and Practice; Szaro, R.C., Johnston, D.W., Eds.; Oxford University Press: Oxford, UK, 1996; pp. 442–453. [Google Scholar]
  12. Steiner, F.; Brooks, K. Ecological planning: A review. Environ. Manag. 1981, 5, 495–505. [Google Scholar] [CrossRef]
  13. Fortin, M.J.; Dale, M.R.T.; ver Hoef, J. Spatial Analysis in Ecology. Available online: http://geog.ucsb.edu/~chris/readings/Spatial.Analysis.in.Ecology.Encyclopedia.Environmetrics.pdf (accessed on 23 December 2019).
  14. Bailey, T.C.; Gatrell, A.C. Interactive Spatial Data Analysis; Longman Scientific & Technical: Essex, UK, 1995; p. 413. [Google Scholar]
  15. Berges, L.; Roche, P.; Avon, C. Establishment of a National ecological network to conserve biodiversity (pros and cons of ecological corridor). Sci. Eaux Territ. 2010, 3, 34–39. [Google Scholar]
  16. Jongman, R. Ecological Networks Are an Issue for All of US. J. Landsc. Ecol. 2008, 1, 7–13. [Google Scholar] [CrossRef] [Green Version]
  17. Opdam, P.; Verboom, J.; Pouwels, R. Landscape Cohesion: An Index for the Conservation Potential of Landscapes for Biodiversity. Landsc. Ecol. 2003, 18, 113–126. [Google Scholar] [CrossRef]
  18. McHugh, N.; Thompson, S. A Rapid Ecological Network Assessment Tool and Its Use in Locating Habitat Extension Areas in a Changing Landscape. J. Nat. Conserv. 2011, 19, 236–244. [Google Scholar] [CrossRef]
  19. Ferretti, V.; Pomarico, S. Ecological Land Suitability Analysis through Spatial Indicators: An Application of the Analytic Network Process Technique and Ordered Weighted Average Approach. Ecol. Indic. 2013, 34, 507–519. [Google Scholar] [CrossRef]
  20. Gurrutxaga, M.; Lozano, P.J.; Del Barrio, G. GIS-Based Approach for Incorporating the Connectivity of Ecological Networks into Regional Planning. J. Nat. Conserv. 2010, 18, 318–326. [Google Scholar] [CrossRef]
  21. Guo, Y.; Liu, Y. Connecting regional landscapes by ecological networks in the Greater Pearl River Delta. Landsc. Ecol. Eng. 2017, 13, 265–278. [Google Scholar] [CrossRef]
  22. Santos, J.S.; Leite, C.C.C.; Viana, J.C.C.; Dos Santos, A.R.; Fernandes, M.M.; Abreu, V.D.S.; Nascimento, T.P.D.; Dos Santos, L.S.; Fernandes, M.R.D.M.; Da Silva, G.F.; et al. Delimitation of Ecological Corridors in the Brazilian Atlantic Forest. Ecol. Indic. 2018, 88, 414–424. [Google Scholar] [CrossRef]
  23. Chang, Q.; Li, X.; Huang, X.; Wu, J. A GIS-Based Green Infrastructure Planning for Sustainable Urban Land Use and Spatial Development. Procedia Environ. Sci. 2012, 12, 491–498. [Google Scholar] [CrossRef] [Green Version]
  24. Xie, H.; Yao, G.; Liu, G. Spatial Evaluation of the Ecological Importance Based on GIS for Environmental Management: A Case Study in Xingguo County of China. Ecol. Indic. 2015, 51, 3–12. [Google Scholar] [CrossRef]
  25. Hong, W.; Guo, R.; Su, M.; Tang, H.; Chen, L.; Hu, W. Sensitivity Evaluation and Land-Use Control of Urban Ecological Corridors: A Case Study of Shenzhen, China. Land Use Policy 2017, 62, 316–325. [Google Scholar] [CrossRef]
  26. Downs, J.A.; Horner, M.W. Network-Based Kernel Density Estimation for Home Range Analysis. Available online: https://0-link-springer-com.brum.beds.ac.uk/chapter/10.1007/978-3-642-12156-2_3#citeas (accessed on 23 December 2019).
  27. Downs, J.A. Network analysis of animal space-use patterns. Ph.D. Thesis, Florida State University, Tallahassee, FL, USA, 2008. [Google Scholar]
  28. Biondi, E.; Casavecchia, S.; Pesaresi, S.; Zivkovic, L. Natura 2000 and the Pan-European Ecological Network: A new methodology for data integration. Biodivers. Conserv. 2012, 21, 1741–1754. [Google Scholar] [CrossRef]
  29. Li, P.; Lv, Y.; Zhang, C.; Yun, W.; Yang, J.; Zhu, D. Analysis and Planning of Ecological Networks Based on Kernel Density Estimations for the Beijing-Tianjin-Hebei Region in Northern China. Sustainability 2016, 8, 1094. [Google Scholar] [CrossRef] [Green Version]
  30. Erdogan, S.; Yilmaz, I.; Baybura, T.; Gullu, M.; Yılmaz, I. Geographical Information Systems Aided Traffic Accident Analysis System Case Study: City of Afyonkarahisar. Accid. Anal. Prev. 2008, 40, 174–181. [Google Scholar] [CrossRef] [PubMed]
  31. Cai, X.; Wu, Z.; Cheng, J. Using Kernel Density Estimation to Assess the Spatial Pattern of Road Density and Its Impact on Landscape Fragmentation. Int. J. Geogr. Inf. Sci. 2013, 27, 222–230. [Google Scholar] [CrossRef]
  32. Law 2863, Kültür ve Tabiat Varlıklarını Koruma Kanunu. 1983; 18113.21: 7; Mevzuatı Geliştirme Yayın Genel Müdürlüğü Mevzuat Bilgi Sistemi, e.mevzuat. Available online: https://www.mevzuat.gov.tr/Metin1.Aspx?MevzuatKod=1.5.2863&MevzuatIliski=0&sourceXmlSearch=2863&Tur=1&Tertip=5&No=2863 (accessed on 9 November 2019).
  33. Law 644, Çevre ve Şehircilik Bakanlığının Teşkilat Ve Görevleri Hakkında Kanun Hükmünde Kararname. 2011; 27984; Mevzuatı Geliştirme Yayın Genel Müdürlüğü Mevzuat Bilgi Sistemi, e.mevzuat. Available online: https://www.mevzuat.gov.tr/Metin1.Aspx?MevzuatKod=4.5.644&MevzuatIliski=0&sourceXmlSearch=644&Tur=4&Tertip=5&No=644 (accessed on 9 November 2019).
  34. UNEP. Mediterranean Action Plan. Convention for the Protection of the Mediterranean Sea against Pollution and Its Related Protocols. Available online: http://wedocs.unep.org/bitstream/handle/20.500.11822/520/mts170.pdf?sequence=2 (accessed on 23 December 2019).
  35. Ministry of Environment and Urbanization. Kas-Kekova Protected Area Boundary. Available online: www.csb.gov.tr (accessed on 9 November 2019).
  36. MoEU. Kaş-Kekova Özel Çevre Koruma Bölgesi Biyolojik Çeşitliliğin Tespiti Projesi. Available online: https://tvk.csb.gov.tr/kas-kekova-ozel-cevre-koruma-bolgesi-karasal-biyolojik-cesitlililigin-tespiti-projesi-proje (accessed on 9 November 2019).
  37. Directive, Habitats. Council Directive 92/43/EEC of 21 May 1992 on the conservation of natural habitats and of wild fauna and flora. Off. J. Eur. Union 1992, 206, 7–50.
  38. Directive, EU Birds. Council Directive of 2 April 1979 on the conservation of wild birds (79/409/EEC). Off. J. Eur. Communities 1979, L103, 25-04.
  39. IUCN, International Union for Conservation of Nature. Available online: www.iucn.org (accessed on 9 November 2019).
  40. International Union for Conservation of Nature. 2006 IUCN Red List of Threatened Species. 2006. Available online: http://www.iucnredlist.org/ (accessed on 29 January 2020).
  41. Aytaç, Z.; Duman, H. A New Species and 2 New Records from Turkey. Turk. J. Bot. 2013, 37, 1055–1060. [Google Scholar] [CrossRef]
  42. Kay, R.; Alder, J. Coastal Planning and Management; First published 1999 by E & FN Spon; Second edition published 2005 by Taylor & Francis: Canada and USA; CRC Press: Boca Raton, FL, USA, 2017. [Google Scholar]
  43. Tezel, D.; Inam, S. Determination of Plant Sensitivity in Protected Areas For Land Management: Kas-Kekova Region. In Proceedings of the FIG Congress 2018, Embracing Our Smart World Where the Continents Connect: Enhancing the Geospatial Maturity of Societies, Istanbul, Turkey, 6–11 May 2018. [Google Scholar]
  44. Tezel, D. Kaş-Kekova (Antalya) Bölgesi Flora ve Fauna Duyarlılığı Üzerine Ekolojik Koridorların Mekansal veri Analizi ile Belirlenmesi. Ph.D. Thesis, Selcuk University, Graduate School of Natural Sciences, Konya, Turkey, August 2018; 164p. Available online: http://acikerisimarsiv.selcuk.edu.tr:8080/xmlui/handle/123456789/14341 (accessed on 29 January 2020).
  45. ESRI. Available online: www.esri.com (accessed on 29 January 2020).
  46. Davies, C.E.; Moss, D.; Hill, M.O. EUNIS Habitat Classification Revised 2004. Report to: European Environment Agency-European Topic Centre on Nature Protection and Biodiversity. 2004, pp. 127–143. Available online: https://www.eea.europa.eu/data-and-maps/data/eunis-habitat-classification/documentation/eunis-2004-report.pdf (accessed on 23 December 2019).
  47. Silverman, B.W. Density Estimator for Statistics and Data Analysis. In Statistics and Applied Probability; Chapman and Hall: New York, NY, USA, 1986. [Google Scholar]
  48. Worton, B.J. Kernel Methods for Estimating the Utilization Distribution in Home-Range Studies. Ecology 1989, 70, 164–168. [Google Scholar] [CrossRef]
  49. Tremaine, G. Home Range Estimation. Int. Encycl. Primatol. 2016, 1–4. [Google Scholar]
  50. Fleming, C.H.; Fagan, W.F.; Mueller, T.; Olson, K.A.; Leimgruber, P.; Calabrese, J.M. Rigorous Home Range Estimation With Movement Data: A New Autocorrelated Kernel Density Estimator. Ecology 2015, 96, 1182–1188. [Google Scholar] [CrossRef] [Green Version]
  51. Shi, X. Selection of Bandwidth Type and Adjustment Side in Kernel Density Estimation over Inhomogeneous Backgrounds. Int. J. Geogr. Inf. Sci. 2010, 24, 643–660. [Google Scholar] [CrossRef]
  52. Kuter, S.; Usul, N.; Kuter, N. Bandwidth Determination for Kernel Density Analysis of Wildfire Events at Forest Sub-District Scale. Ecol. Model. 2011, 222, 3033–3040. [Google Scholar] [CrossRef]
  53. Average Nearest Neighbor Analysis. Available online: https://pro.arcgis.com/en/pro-app/tool-reference/spatial-statistics/h-how-average-nearest-neighbor-distance-spatial-st.htm (accessed on 23 December 2019).
  54. Gundogdu, E.; Ogurlu, I. The distribution of Wild Goat Capra aegagrus Erxleben 1877 and population characteristics in Isparta, Turkey. J. Anim. Veter. Advan. 2009, 8, 2318–2324. [Google Scholar]
  55. Jenks, G.F. The data model concept in statistical mapping. Int. Yearb. Cartogr. 1967, 7, 186–190. [Google Scholar]
  56. Fahrig, L. Effects of habitat fragmentation on biodiversity. Annu. Rev. Ecol. Evol. Syst. 2003, 34, 487–515. [Google Scholar] [CrossRef] [Green Version]
  57. Debinski, D.M.; Holt, R.D. A Survey and Overview of Habitat Fragmentation Experiments. Conserv. Biol. 2000, 14, 342–355. [Google Scholar] [CrossRef]
  58. Eklund, J.F.; Cabeza-Jaimejuan, M.D.M. Quality of governance and effectiveness of protected areas: Crucial concepts for conservation planning. Ann. N. Y. Acad. Sci. 2017, 1399, 27–41. [Google Scholar] [CrossRef] [PubMed]
  59. Franks, P. Social Assessment of Protected Areas: Early Experience and Results of a Participatory, Rapid Approach; International Institute for Environment and Development: London, UK, 2014. [Google Scholar]
  60. Lockwood, M. Good Governance for Terrestrial Protected Areas: A Framework, Principles and Performance Outcomes. J. Environ. Manag. 2010, 91, 754–766. [Google Scholar] [CrossRef]
  61. Reed, M.S. Stakeholder Participation for Environmental Management: A Literature Review. Biol. Conserv. 2008, 141, 2417–2431. [Google Scholar] [CrossRef]
  62. Ward, C.; Stringer, L.C.; Holmes, G. Protected Area Co-Management and Perceived Livelihood Impacts. J. Environ. Manag. 2018, 228, 1–12. [Google Scholar] [CrossRef]
  63. Chomitz, K.M.; Gray, D.A. Roads, Land Use, and Deforestation: A Spatial Model Applied to Belize. World Bank Econ. Rev. 1996, 10, 487–512. [Google Scholar] [CrossRef] [Green Version]
  64. Ashley, E.P.; Robinson, J.T. Road mortality of amphibians, reptiles and other wildlife on the Long Point Causeway, Lake Erie, Ontario. Can. Field Nat. 1996, 110, 403–412. [Google Scholar]
  65. Forman, R.T.T.; Alexander, L.E. Roads and Their Major Ecological Effects. Annu. Rev. Ecol. Syst. 1998, 29, 207–231. [Google Scholar] [CrossRef] [Green Version]
  66. Laurance, W.F.; Peletier-Jellema, A.; Geenen, B.; Koster, H.; Verweij, P.; Van Dijck, P.; Lovejoy, T.E.; Schleicher, J.; Van Kuijk, M. Reducing the Global Environmental Impacts of Rapid Infrastructure Expansion. Curr. Biol. 2015, 25, R259–R262. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  67. Arriagada, R.A.; Echeverria, C.M.; Moya, D.E. Creating protected areas on public lands: Is there room for additional conservation? PLoS ONE 2016, 11, e0148094. [Google Scholar]
  68. Caro, T.; Dobson, A.; Marshall, A.J.; Peres, C.A. Compromise Solutions between Conservation and Road Building in the Tropics. Curr. Biol. 2014, 24, R722–R725. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  69. Fahrig, L.; Pedlar, J.H.; Pope, S.E.; Taylor, P.D.; Wegner, J.F. Effect of Road Traffic on Amphibian Density. Biol. Conserv. 1995, 73, 177–182. [Google Scholar] [CrossRef] [Green Version]
  70. Hodson, N.L. A survey of road mortality in mammals (and including data for the grass snake and common frog). J. Zool. 1966, 148, 576–579. [Google Scholar] [CrossRef]
  71. Mao, X.; Meng, J.; Xiang, Y. Cellular automata-based model for developing land use ecological security patterns in semi-arid areas: A case study of Ordos, Inner Mongolia, China. Environ. Earth Sci. 2013, 70, 269–279. [Google Scholar] [CrossRef]
  72. Shi, H.; Shi, T.; Yang, Z.; Wang, Z.; Han, F.; Wang, C. Effect of Roads on Ecological Corridors Used for Wildlife Movement in a Natural Heritage Site. Sustain. 2018, 10, 2725. [Google Scholar] [CrossRef] [Green Version]
Figure 1. (a) The study site and the registered protected area (PA) boundaries (upper left coordinates: Lat: 36.250 Lon: 29.630; lower right coordinates: Lat: 36.060 Lon: 29.940) [35]. The red star denotes the location of Simena Castle in Kalekoy, Demre, Turkey. (b) Photo taken from the Simena Castle; which shows the coastal landscapes and archaeological ruins, mostly from Lycians (image credits: MoEU).
Figure 1. (a) The study site and the registered protected area (PA) boundaries (upper left coordinates: Lat: 36.250 Lon: 29.630; lower right coordinates: Lat: 36.060 Lon: 29.940) [35]. The red star denotes the location of Simena Castle in Kalekoy, Demre, Turkey. (b) Photo taken from the Simena Castle; which shows the coastal landscapes and archaeological ruins, mostly from Lycians (image credits: MoEU).
Ijgi 09 00091 g001
Figure 2. The sub-districts of Kas (Bayindir, Belenli, Bogazcik, Kilicli) and Demre (Cevreli, Kapakli, Kale-Ucagiz) municipalities and the corresponding human population (as of year 2013) in Kas-Kekova PA.
Figure 2. The sub-districts of Kas (Bayindir, Belenli, Bogazcik, Kilicli) and Demre (Cevreli, Kapakli, Kale-Ucagiz) municipalities and the corresponding human population (as of year 2013) in Kas-Kekova PA.
Ijgi 09 00091 g002
Figure 3. The changes of settlement in Kas-Kekova PA caused mainly by touristic activity that can be observed in aerial and satellite images from 1992 (a) and 2008 (b). The red polygons denote the rural settlement boundary registered by MoEU.
Figure 3. The changes of settlement in Kas-Kekova PA caused mainly by touristic activity that can be observed in aerial and satellite images from 1992 (a) and 2008 (b). The red polygons denote the rural settlement boundary registered by MoEU.
Ijgi 09 00091 g003
Figure 4. Methodological framework for the Kas-Kekova PA.
Figure 4. Methodological framework for the Kas-Kekova PA.
Ijgi 09 00091 g004
Figure 5. Endemic and/or threatened flora and fauna locations and the EUNIS Habitat classes determined identified and measured throughout the Kas-Kekova PA Biological Diversity Project carried out by MoEU [36].
Figure 5. Endemic and/or threatened flora and fauna locations and the EUNIS Habitat classes determined identified and measured throughout the Kas-Kekova PA Biological Diversity Project carried out by MoEU [36].
Ijgi 09 00091 g005
Figure 6. Land ownership pattern classified according to the land registry type (i.e., forest or cadastral) and owner type (i.e., owned by state or private).
Figure 6. Land ownership pattern classified according to the land registry type (i.e., forest or cadastral) and owner type (i.e., owned by state or private).
Ijgi 09 00091 g006
Figure 7. The distributions of vulnerable mammals observed by field work on the slope (left) and altitude range (right) maps.
Figure 7. The distributions of vulnerable mammals observed by field work on the slope (left) and altitude range (right) maps.
Ijgi 09 00091 g007
Figure 8. The distributions of vulnerable reptiles and amphibians (Table 1) observed during the field works on the slope (left) and altitude range (right) maps.
Figure 8. The distributions of vulnerable reptiles and amphibians (Table 1) observed during the field works on the slope (left) and altitude range (right) maps.
Ijgi 09 00091 g008
Figure 9. The distributions of endemic flora species (Table 1) observed during the field works on the altitude range (a), slope (b) and aspect (c) maps.
Figure 9. The distributions of endemic flora species (Table 1) observed during the field works on the altitude range (a), slope (b) and aspect (c) maps.
Ijgi 09 00091 g009
Figure 10. The results of KDE bandwidth analyses with fauna data (all species types were included in the results) and endemic flora locations. The results were classified according to natural breaks.
Figure 10. The results of KDE bandwidth analyses with fauna data (all species types were included in the results) and endemic flora locations. The results were classified according to natural breaks.
Ijgi 09 00091 g010aIjgi 09 00091 g010b
Figure 11. Results of KDE analysis with 1000 m radius results classified according to natural breaks for: (a) endemic flora locations; (b) reptiles and amphibians; (c) mammals; and (d) birds.
Figure 11. Results of KDE analysis with 1000 m radius results classified according to natural breaks for: (a) endemic flora locations; (b) reptiles and amphibians; (c) mammals; and (d) birds.
Ijgi 09 00091 g011
Figure 12. Habitat network map produced by integrating the results in Figure 11 according to the maximum class value of each grid point. The flora and fauna locations are shown on the map.
Figure 12. Habitat network map produced by integrating the results in Figure 11 according to the maximum class value of each grid point. The flora and fauna locations are shown on the map.
Ijgi 09 00091 g012
Figure 13. Cadastral analysis results and the determined habitat network.
Figure 13. Cadastral analysis results and the determined habitat network.
Ijgi 09 00091 g013
Figure 14. The road impact zones and settlement areas overlaid with the determined ecological network.
Figure 14. The road impact zones and settlement areas overlaid with the determined ecological network.
Ijgi 09 00091 g014
Table 1. List of endemic flora species which had their coordinates collected using global positioning system (GPS) in field works [36].
Table 1. List of endemic flora species which had their coordinates collected using global positioning system (GPS) in field works [36].
Biodiversity ElementsFamilySpeciesIUCN Category
FloraCampanulaceaeCampanula lycicaCritically endangered (CR)
EricaceaeErica boucquetiiCR
UmbelliferaeDaucus conchitaeCR
CaryophyllaceaeSilene delicatulaEndangered (EN)
AristolochiaceaeAristolochia poluniniiEN
BoraginaceaeAlkanna macrosiphonNT
CampanulaceaeCampanula podocarpaNear threatened (NT)
IridaceaeCrocus cancellatusNT
CaryophyllaceaeArenaria pamphylicaVulnerable (VU)
CompositaeAnthemis roseaVU
FagaceaeQuercus aucheriVU
LabiataeBallota glandulosissimaVU
LabiataeSideritis albifloraVU
LabiataeTeucrium montbretiiVU
LiliaceaeAllium sandrasicumVU
LiliaceaeFritillaria elwesiiVU
LabiataeAjuga bombycinaVU
Reptiles & AmphibiansSalamandridaeLyciasalamandra luschaniVU
BufonidaePseudepidalea variabilisVU
TestudinidaeTestudo graecaVU
BirdsStrigidaeAthena noctuaEN
JyngidaeJynx torguillaEN
PicidaeDendrocopus syriacusEN
SylvidaeSylvia rupelliEN
LaridaeLarus audouiniiNT
CoraciidaeCoracias garrulusNT
EmberizidaeEmberiza cineracaeNT
LaridaeLarus audouiniiNT
FalconidaeFalco naumanniVU
MammalsBovidaeCapra aegagrus ErxlebenVU
Table 2. EUNIS Habitat types in Kas-Kekova PA [36].
Table 2. EUNIS Habitat types in Kas-Kekova PA [36].
EUNIS CodeDescription of Habitat TypesHabitat Types
A2.54Low-mid saltmarshesMarine habitats
B1.2Sand beaches above the driftlineCoastal habitats
B3.2Unvegetated rock cliffs, ledges, shores and islets.
B3.3Rock cliffs, ledges and shores, with angiosperms.
F5.2MaquisHeathland, scrub and tundra
F7.3East Mediterranean phrygana
G3.7Lowland to montane Mediterranean pine woodland (excluding black pine Pinus nigra).Woodland, forest and other wooded land
G4. EMixed Mediterranean pine—evergreen oak woodland.
G5.6Early-stage natural and semi-natural woodlands and regrowth.
H3.5PavementsInland unvegetated and sparsely vegetated habitats
I1.2Mixed crops of market gardens and horticultureRegularly or recently cultivated agricultural, horticultural and domestic habitats
J1.2Residential buildings of villages and urban peripheriesConstructed, industrial and other artificial habitats
J2.4Agricultural constructions
J4.7Constructed parts of cemeteries
Table 3. Formats, sources, and accuracies of the datasets evaluated in the study.
Table 3. Formats, sources, and accuracies of the datasets evaluated in the study.
Data Type/NameData Source (Produced by)Format (Vector/Raster)Accuracy
Flora and fauna speciesMoEUPoint1–50 m
EUNIS Habitat MapMoEUPolygon2–2,5 m
RoadsNational databaseLine0.50–1 m
Cadastral parcelsMoEUPolygon0.20 m
Human population (per district) MoEUPointN.A.
Table 4. The results of average nearest neighbor (ANN) analysis applied to the point locations of specie groups.
Table 4. The results of average nearest neighbor (ANN) analysis applied to the point locations of specie groups.
GroupObserved Mean Distances (DO) Expected Mean Distances (DE)Distribution Pattern Radius Selection Strategy
Flora 234 m736 mClusteredDE
Mammals2845 m1235 mDispersedHome range
Birds1037 m1162 mRandomDO
Reptiles1276 m1100 mRandomDO
Table 5. Natural breaks classification of the kernel density estimation (KDE) patterns.
Table 5. Natural breaks classification of the kernel density estimation (KDE) patterns.
Class ValueKernel Density Estimation (KDE)Spatial Configuration
5High density patternCore habitats
4Medium high-density patternBuffer patches
3Medium-density patternEcological corridors
2Medium low-density patternStepping stones
1Low density patternIsolated area
Table 6. The percentages of different ownership types in the PA.
Table 6. The percentages of different ownership types in the PA.
Land Property Ownership TypeArea (ha)Percentage (%)
Private property1.45114,5
Private property obtained after deforestation5455,5
State property (forest)7.38373,9
State property (non-forest)6106,1
Total9.989100

Share and Cite

MDPI and ACS Style

Tezel, D.; Inam, S.; Kocaman, S. GIS-Based Assessment of Habitat Networks for Conservation Planning in Kas-Kekova Protected Area (Turkey). ISPRS Int. J. Geo-Inf. 2020, 9, 91. https://0-doi-org.brum.beds.ac.uk/10.3390/ijgi9020091

AMA Style

Tezel D, Inam S, Kocaman S. GIS-Based Assessment of Habitat Networks for Conservation Planning in Kas-Kekova Protected Area (Turkey). ISPRS International Journal of Geo-Information. 2020; 9(2):91. https://0-doi-org.brum.beds.ac.uk/10.3390/ijgi9020091

Chicago/Turabian Style

Tezel, Dilek, Saban Inam, and Sultan Kocaman. 2020. "GIS-Based Assessment of Habitat Networks for Conservation Planning in Kas-Kekova Protected Area (Turkey)" ISPRS International Journal of Geo-Information 9, no. 2: 91. https://0-doi-org.brum.beds.ac.uk/10.3390/ijgi9020091

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop