Next Article in Journal
The Role of Vegetation in Mitigating Urban Land Surface Temperatures: A Case Study of Munich, Germany during the Warm Season
Previous Article in Journal
A Hybrid Inexact Optimization Method for Land-Use Allocation in Association with Environmental/Ecological Requirements at a Watershed Level
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Habitat Mapping of the Leopard Cat (Prionailurus bengalensis) in South Korea Using GIS

1
Korea Adaptation Center for Climate Change, Korea Environment Institute (KEI), 370 Sicheong-daero, Sejong-si 339-007, Korea
2
Department of Landscape Architecture, Dankook University, 119 Dandae-ro, Dongnam-gu, Cheonan-si, Chungnam 330-714, Korea
3
Geological Research Division, Korea Institute of Geoscience and Mineral Resources (KIGAM), 124 Gwahang-no, Yuseong-gu, Daejeon 305-350, Korea
4
Korea University of Science and Technology, 217 Gajeong-ro, Yuseong-gu, Daejeon 305-350, Korea
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Sustainability 2015, 7(4), 4668-4688; https://0-doi-org.brum.beds.ac.uk/10.3390/su7044668
Submission received: 22 January 2015 / Revised: 9 April 2015 / Accepted: 9 April 2015 / Published: 21 April 2015
(This article belongs to the Section Environmental Sustainability and Applications)

Abstract

:
The purpose of this study was to create maps of potentially sustainable leopard cat (Prionailurus bengalensis) habitats for all of South Korea. The leopard cat, which is on the International Union for Conservation of Nature (IUCN) Red List, is the only member of the Felidae family in Korea. To create habitat potential maps, we selected various environmental factors potentially affecting the species’ distribution from a spatial database derived from geographic information system (GIS) data: elevation, slope, distance from a forest stand, road, or drainage, timber type, age, and land cover. We analyzed the spatial relationships between the distribution of the leopard cat and the environmental factors using a frequency ratio model and a logistic regression model. We then overlaid these relationships to produce a habitat potential map with a species potential index (SPI) value. Of the total number of known leopard cat locations, we used 50% for mapping and the remaining 50% for model validation. Our models were relatively successful and showed a high level of accuracy during model validation with existing locations (frequency ratio model 82.15%; logistic regression model 81.48%). The maps can be used to manage and monitor the habitat of mammal species and top predators.

1. Introduction

Long-term persistence of biodiversity is the ultimate goal of most conservation plans (Minor [1]). Maintaining biodiversity is important because identifying the species critical to the maintenance of ecosystem stability or those that will be useful human resources in the future is not always easy (Burton [2]). The most effective way to conserve biodiversity globally is to focus on the protection of high biodiversity areas (Myers [3]). One method of identifying these areas is to model top predator habitats, which are often biodiversity hot spots (Schmitz [4], Sergio [5], Gavashelishvili [6]).
Carnivores are major predators and scavengers in terrestrial ecosystems. Understanding their status is important for understanding ecosystem integrity in regions of high human disturbance. The leopard cat (Prionailurus bengalensis) is a small, wild cat native to North and East Asia. Since 2002, the International Union for Conservation of Nature (IUCN) has listed it as a species of least concern. Although it is widely distributed, it is also threatened by habitat loss and hunting in parts of its range. Of the 12 leopard cat subspecies, which differ widely in appearance (Sanderson [7]) only one exists in South Korea.
GIS is a useful tool for determining the spatial relationships of an event and its controlling factors. GIS has been successfully used for habitat mapping of an event distribution based on probability and statistical models (Kim [8], Lee [9]), Analytic Hierarchy Process (AHP) decision models (Matsuura [10]), fuzzy relations (Choi [11]), and artificial neural networks (Song [12]). More recently, many studies have employed GIS to produce habitat maps for various species. Poplar‐Jeffers [13] and Ottaviani [14] used a GIS-based model to quantify and indicate the habitat of mammals. Newton‐Cross [15], Tien Bui [16], and Huck [17] mapped the distribution of badgers using a logistic regression model. In the case of bats, some studies have applied statistical models to analyze their habitat distributions Jaberg [18] and Greaves [19], Clement [20]. Northrup [21] used logistic regression models in a geographic information system to map the probability of bear-human conflict and the relative probability of grizzly bear habitat selection based on global positioning system radiotelemetry data. Kuemmerle [22] applied a probabilistic model to habitat mapping of the European bison (Bison bonasus), Speed [23] used remote sensing data for habitat mapping of deer, and Gavashelishvili [6] analyzed leopard habitat in central Asia using a logistic regression model. Studies of the spatial relationships between species habitats and various ecological environments have also been attempted (Kocev [24], Meixler [25], Walker [26], Walters [27]). However, none of these studies has analyzed the habitat of the leopard cat in South Korea. Therefore, the purpose of our study was to identify relationships between the distribution of the leopard cat and various environmental factors using probability and statistical models. Furthermore, we used the resulting leopard cat habitat range to map the species in South Korea. We used frequency ratios and logistic regressions as probabilistic and statistical models, respectively.
South Korea is located in North Asia (37°16′32″N, 127°03′05″E), in the southern half of the Korean Peninsula (Figure 1). South Korea is geomorphologically stable and contains three major mountain ranges: the Taebaek Mountains, the Sobaek ranges, and the Jiri Massif. Furthermore, no active volcanoes exist and no strong earthquakes have occurred. It has no extensive plains and its lowlands, which make up approximately 30% of its land area, are the product of mountain erosion. Uplands and mountains comprise the remaining area.
Figure 1. Study area with leopard cat occurrences points.
Figure 1. Study area with leopard cat occurrences points.
Sustainability 07 04668 g001

2. Data

2.1. Leopard Cat Survey

From 1997 to 2005, field experts from the National Institute of Environmental Research, universities, and research institutes conducted surveys on species occurrences and spatial distributions of wildlife via direct observations, surveys of local residents, and field signs such as tracks, feces, and footprints. This study is part of the second national environment survey in South Korea. The leopard cat, or traces of the species, was found at 630 points in the study area. Some of the data were excluded because they identified the locations based on indirect research methods such as a local resident survey based on residents’ hearing the animals rather than a direct observation of the species. This resulted in the removal of 201 points, with a further 429 points utilized for modeling and validation.
The leopard cat is one of the top predators in the mountainous regions of Korea, but few habitat models have been developed at the regional level for this animal. The species is designated as an Endangered Species Type II by the Wildlife Conservation Act of Korea. Recent habitat loss and fragmentation are critical factors threatening leopard cat populations in mountainous regions. Despite the leopard cat’s wide distribution, little is known about its ecology or behavior in the wild.
Since 1986, Korea has undertaken national environmental research consisting of nine sectors: (1) geography; (2) vegetation; (3) flora; (4) mammals; (5) birds; (6) herptiles; (7) freshwater fishes; (8) land insects; and (9) benthic macroinvertebrates. Under the supervision of the National Institute of Environmental Research (NIER), the first and second national environmental research efforts occurred between 1986 and 1990 and between 1997 and 2005, respectively. Field experts from the National Institute of Environmental Research, universities, and research institutes have conducted surveys on species occurrence and the spatial distributions of wildlife. The third survey has been in progress since 2006, and is expected to conclude in 2015. In the case of mammals, research has been conducted annually from February to October. Literature reviews, and geographical and vegetation statuses were used to determine the most appropriate locations for identifying local mammals. When species identification is difficult in the field, researchers resort to capturing the animals after obtaining necessary permissions. After the investigation, all captured mammals are released at the locations in which they were captured. Since direct observation of mammals is difficult, interviews have been used to complement fieldwork. Through interviews with local mammal experts, museum workers, hunting license issuers, and local residents, these efforts have identified apparent species locations, seasons, and population sizes. Since leopard cats are distributed throughout forested areas, research has focused not only on locations related to agricultural lands, but also on ridges, considering the total forested land. Since direct observation of leopard cats is difficult, trace investigations, including that of excrement, has been necessary. In addition, data have been acquired through interviewing local residents. Locations in which leopard cats have been observed or detected via traces were recorded by a global positioning system (GPS) or map, and used for constructing geographic information system (GIS) data.

2.2. Controlling Factors

The distribution of the leopard cat is the result of the interaction of complex factors. The selection of these factors and the preparation of corresponding thematic data layers are crucial components of any model for leopard cat habitat potential mapping. The important factors for mammalian habitat include ground elevation, slope steepness, aspect, timber distribution, land cover, land use, and human activity (Table 1). These factors were collected from available maps and field investigations. A digital elevation model (DEM) was prepared through the digitization of contours at 5-m intervals from the topographical maps. Using the DEM, the slope gradient and slope aspect were calculated. The forest map is a series of polygons with a scale of 1:25,000 that is published by the Korea Forest Research Institute (KFRI). The land-cover types were identified from a panchromatic SPOT-5 (Système Probatoire d’Observation de la Terre 5) image taken in November 2007.
The land-cover map was also a series of polygons with a scale of 1:5000 that was published by the Korea Ministry of Environment. The diameter, type, density, and age of the timber were obtained from the forest maps. The maps relevant to leopard cat occurrence were constructed in a vector format spatial database using the ArcGIS (ESRI) software package. To calculate the frequency ratio for the class or type of each factor, the scale factors were divided into 10 classes with equal area using ArcGIS. Therefore, the range of each class was automatically determined based on equal areas. Nine factors, both calculated and extracted (Table 1) from the maps, were converted to a 30 × 30-m grid format (ArcGIS GRID type). As a result, the dimensions of the study area grid were 7224 rows by 6792 columns, and the total number of cells was 12,307,439 (except those with no data). Then, the factors were converted to ASCII data for use with the statistics program. All of the factors were placed into one of the 10 classes. Each of the analyzed factors (Table 2) were made by utilizing the data of Table 1.
Table 1. Data layer related to leopard cat of study area.
Table 1. Data layer related to leopard cat of study area.
CategoryFactorsData TypeScale
HabitatLeopard catPoint-
Forest map aTimber type
Timber age
Polygon1:25,000
Land CoverDistance from road (m)
Distance from water (m)
Distance from forest (m)
Land cover
Polygon1:25,000
Topographic map bGround elevation (m)
Slope gradient (°)
Slope aspect
GRID1:5000
a The forest map produced by Korea Forest Service (KFS; the http://www.forest.go.kr); b Topographical factors were extracted from digital topographic map by National Geographic Information Institute (NGII; http://www.ngii.go.kr).

3. Methods

The general progression of leopard cat potential habitat mapping is illustrated in Figure 2. The process of mapping the habitat potentials of two leopard cat communities in South Korea included six major steps: (1) a field survey to determine the occurrence of the leopard cat; (2) the determination and construction of a database of the controlling factors, in combination with a GIS analysis; (3) the construction of a spatial database based on the two leopard cat communities and nine factors influencing their distribution; (4) the division of leopard cat individuals into a training set (50%) to analyze habitat potential using models and a test set (50%) to validate the predicted potential habitat map; (5) data processing using the frequency ratio and logistic regression models; and (6) validation of the leopard cat potential habitat maps using the known distributions of leopard cat that were not used in the analysis.
Figure 2. Flowchart of study.
Figure 2. Flowchart of study.
Sustainability 07 04668 g002
For the application and validation of habitat potential models for the leopard cat, known locations were identified via interviews and field surveys in 2005 (Figure 1). We used these locations as the dependent variable, and nine factors believed to influence leopard cat habitat were set as independent variables: elevation, slope, aspect, timber distribution, land cover, land use, and human activity. Using known locations and calculated or database-extracted factors for model training, we conducted a habitat analysis using a frequency ratio model. A logistic regression model was also used for leopard cat habitat analysis. The frequency ratio model is a simple and basic technique that can be used to explain spatial relationships between known locations and potential habitat-influencing factors. For the application of these models, a statistical package was used in the GIS program. Finally, the resulting models were validated using known leopard cat locations that were not used to train the model.
As stated earlier, the frequency ratio and logistic regression were used as a probabilistic and statistical model, respectively. The frequency ratio, as a probability model, can be easily represented as the frequency ratio of each factor. The frequency ratio is the probability of occurrence of a certain attribute (Bonham-Carter [28]). The frequency ratio is the ratio of the area in which an event in a class or type for a given factor occurs divided by the overall study area. In Equation (1), P(P) denotes the area ratio for the class or type for a given number of unit cells containing a percentage of the pixels in the domain for the class, and P(O) denotes the percentage of occurrence in the total event. The frequency ratio of each factors type or class, C, is then expressed by:
F R C = P C ( O ) P C ( P )
Logistic regression enables investigation of multivariate regression relations between one dependent and several independent variables. Logistic regression is limited in that the calculation process cannot be traced because it repeats calculations to find the optimized regression equation for determining the possibility that the dependent variable will occur. However, logistic regression does have the following advantages: (1) the assumption of a normal distribution is not applied for independent variables because the relationship between the dependent variable and independent variables is identified as a non-linear relationship; (2) It is able to explain complex phenomena because a range of data types can be used for the independent variable, including discrete, gradational, nominal, and continuous types. Thus, the method is suitable for analyzing complex spatial relationships in a quantitative manner; (3) The result of a logistic regression analysis includes individual values for each factor related to the habitat. These factor values can be used in similar studies targeting other regions.
By determining the frequency ratio, the area ratio for leopard cat habitat was calculated for the range or type of each factor, and the area ratio of the range or type of each factor to the total area was calculated. Finally, the frequency ratios for the range or type of each factor were calculated by dividing the distribution area ratio by the area ratio. The frequency ratio was assigned to each factor’s class. The frequency ratio of the habitat potential was created using the overlay functions in the GIS, which were used to merge different factors that were assigned the ratio.
To apply the logistic regression model for analysis of leopard cat habitat potential mapping, the dependent variable was binary, representing the presence or absence of habitat. Quantitatively, the relationship between the occurrence and its dependency on several variables can be expressed as:
P ( T = 1 | X ) = Λ ( X )
with the logistic function Λ, and introducing the actual logistic regression model explicitly as:
logit ( P ( T = 1 | X 0 X n ) ) = z P ( T = 1 | X 0 X n ) = Λ ( z )
where P is the probability of an event’s occurring, and e is the natural logarithm. In the present situation, P is the estimated probability of a habitat based on intrinsic properties only, which we term “susceptibility.” The probability varies from 0 to 1 on an S-shaped curve, and z is the linear combination. It follows that logistic regression involves fitting the data to an equation of the form:
z = b0x0 + b1x1 + b2x2 + … + bnxn
where b0 is the intercept of the model, bi (i = 0, 1, 2, …, n) represents the slope coefficients of the logistic regression model, and xi (i = 0, 1, 2, …, n) are independent variables (Dai [29]). The linear model that is formed is then a logistic regression for the presence or absence of leopard cat habitat (present conditions) on the independent variables (pre-failure conditions).
Using these formulae, a habitat potential map was constructed. The logistic regression analysis was performed by dividing the study area into grid squares of 30 × 30 m. Data for the 11 factors were converted to an ASCII format for use in the statistical package.
Although a “best-fit” equation is found in logistic regression using the same least-squares method as linear regression, the principles on which it does so are rather different. Instead of using least-square deviations criteria for the best fit, logistic regression employs a maximum likelihood method, which maximizes the probability of obtaining the observed results given the fitted regression coefficients. A consequence of this is that the goodness of fit and overall significance statistics used in logistic regression are different from those used in linear regression. Here, logistic regression was used to calculate and map the probability of habitat potential, and logistic regression values for the study area were applied.
Log likelihood is a key concept for understanding the tests used in logistic regression. Normally, overall significance is determined by a chi-square test, which is derived from the likelihood of observing the actual data under the assumption that the model that has been fitted is accurate. Table 2 and Table 3 contain the base model results for the frequency ratio and logistic regression analysis.

4. Results

4.1. Factors that Influence Leopard Cat Distributions

We evaluated the spatial data using the frequency ratio model to reveal correlations between the distribution of the leopard cat and various environmental factors (Table 1) in the study area. A positive correlation designates higher habitat potential, while a negative correlation indicates lower habitat potential.
Relationships between the distribution of the leopard cat and topography-related environmental factors derived from the digital elevation model (DEM; Table 2) are as follows. All of the topography-related factors (elevation, slope, and aspect) were positively correlated with the distribution of the leopard cat, indicating that at higher elevations and steeper slopes, habitat potential increases for this species. Higher elevations and steeper slopes may provide a safe habitat from competition, including that from humans. Accordingly, lower elevations and gentler slopes could produce the opposite result. The frequency ratio model results indicated that elevation, slope-related factors, and timber age were positively correlated with leopard cat locations. Areas closer to water and forest and farther from roads had higher habitat potential. With respect to timber type, oak forests showed the most positive correlation, and broad-leaved forest classes presented the most positive correlation among land cover types. In general, grasslands and edge areas formed adjacent to water and forests are known to be productive. These results indicate that the leopard cat uses habitats containing both safe topological features and rich food sources.
Table 2. Frequency ratio values between leopard cat and related factors.
Table 2. Frequency ratio values between leopard cat and related factors.
FactorClassNo. of Leopard Cat% of Leopard CatNo. of Pixels in Domain% of Pixels in DomainFrequency Ratio
Timber TypeNo data136.054,804,10139.030.15
Farmland20.9363,4870.521.80
Larch73.26583,5204.740.69
Pinus rigida fores52.33420,9833.420.68
Non-stocked forest land00.0017,8620.150.00
Chestnut artificial forest00.0057,5920.470.00
Cut-over area00.00440.000.00
Pinus densiflora Forests3717.211,961,70715.941.08
Pinus densiflora artificial forest00.0017,9960.150.00
Pinus koraiensis forest00.00245,0091.990.00
Left-over area31.4075,3770.612.28
Oak forest31.4025,7840.216.66
Oak artificial forest00.00350.000.00
Grassland00.0059,6280.480.00
Conifer mixed forest00.009280.010.00
Mixed forest of soft and hardwood4520.931,971,62816.021.31
Poplar forest73.26154,4321.252.59
Water00.00100.000.00
Broadleaved forest9343.261,838,13514.942.90
Denuded land00.001190.000.00
Bamboo stand00.0082760.070.00
Other00.007860.010.00
Timber ageNon forest area209.304,807,56839.060.24
1st age83.72348,1702.831.31
2nd age115.121,119,3939.100.56
3rd age7333.953,093,56625.141.35
4th age6027.912,229,05218.111.54
5th age2411.16558,1524.542.46
6th age198.84151,5381.237.18
Distance from forest (m)020093.027,599,55661.751.51
0–376156.983,758,17530.540.23
376–75200.00590,0224.790.00
752–112900.00173,5951.410.00
1129–150500.0056,3640.460.00
1505–188100.0047,6620.390.00
1881–263300.0039,6950.320.00
2633–526600.0017,8030.140.00
5266–978000.0011,9950.100.00
9780–96,29900.0012,5720.100.00
Distance from road (m)010.47998,1658.110.06
0–8924219.537,167,23558.230.34
892–17837333.952,503,82220.341.67
1783–26755626.05878,5847.143.65
2675–35662612.09321,6282.614.63
3566–5346125.58204,6221.663.36
5346–802452.3376,2070.623.76
8024–12,48100.0062,5330.510.00
12,481–21,39600.0050,7920.410.00
21,396–228,22500.0043,8510.360.00
Distance from Water (m)0–3004922.796,341,83951.530.44
300–600 13261.404,277,69634.761.77
600–900 3315.351,307,76110.631.44
900–1200 10.47228,6511.860.25
1200–1500 00.0075,7390.620.00
1500–1800 00.0047,2200.380.00
1800–2200 00.0073380.060.00
2200–2900 00.0031960.030.00
2900–4100 00.0097330.080.00
4100–231,60000.0082660.070.00
Land CoverNo Data00.0025,3660.210.00
Other00.0060.000.00
Residential Area10.47368,2592.990.16
Manufacturing Area00.0087,1770.710.00
Commercial Area00.0051,0800.420.00
Recreational Area00.0072660.060.00
Trafficker Area10.47150,2791.220.38
Public Area00.0063,6230.520.00
Agricultural Area00.00300.000.00
Paddy41.861,673,34213.600.14
Field73.261,086,5078.830.37
Growing in Plastic Greenhouse00.0055,5900.450.00
Orchard00.00166,2131.350.00
Other Plantations00.0044,7300.360.00
Broadleaved Forest8740.472,121,55617.242.35
Coniferous Forest5626.053,398,50227.610.94
Mixed Forest5726.512,038,54116.561.6
Forest Area00.00120.000.00
Natural Grassland00.0064,1090.520.00
Artificial Pasture00.0014,7780.120.00
Other Grassland10.47201,1261.630.29
Wetland00.0040.000.00
Inland Wetland10.4784,9900.690.68
Coastal Wetland00.0043,7140.360.00
Bare Land00.00200.000.00
Natural Bare Land00.0010,6300.090.00
Other Bare Land00.00207,7581.690.00
Beach00.0010.000.00
Inland Water00.00271,5652.210.00
Marine Water00.0070,6650.570.00
Ground Elevation (m)1–2010.471,477,42812.000.04
20–6031.401,528,01712.420.11
60–10062.791,252,32210.180.27
100–149104.651,162,5049.450.49
149–200156.981,215,8079.880.71
200–26494.191,145,7979.310.45
264–342156.981,133,7059.210.76
342–4564219.531,135,6289.232.12
456–6404119.071,135,5809.232.07
640–19407333.951,120,6519.113.73
Slope gradient (°)000.001,598,77912.990.00
0–2.762.791,221,8229.930.28
2.7–5.3125.581,227,4359.970.56
5.3–8.0167.441,187,7459.650.77
8.0–11.02612.091,175,8799.551.27
11.0–13.82813.021,166,9709.481.37
13.8–17.03817.671,227,3539.971.77
17.0–20.0209.31,188,7809.660.96
20.0–24.53013.951,167,2779.481.47
24.5–56.63918.141,145,3999.311.95
Slope aspectFlat00.001,598,77912.990.00
N219.771,212,6649.850.99
NE3616.741,358,78911.041.52
E2913.491,395,19111.341.19
SE3114.421,367,29211.111.30
S167.441,235,55710.040.74
SW2411.161,389,11211.290.99
W2712.561,396,31811.351.11
NW3114.421,353,73711.001.31

4.2. Habitat Potential Mapping

The frequency model was used to derive and calculate correlation ratings between the leopard cat distribution and each factor influencing habitat. Each factor’s rating was assigned as the relationship between leopard cat distribution and each factor’s type or range (Table 2). The ratio of the number of cells where the leopard cat was not founded to the number of cells where the leopard cat was founded is shown in Table 2. The habitat potential index (HPIFR), Equation (5), was calculated by a summation of each factor ratio value, as shown in Table 2 (Lee [30]):
HPIFR = ΣFR
where FR is the rating of each factor type or range. A FR of 1 indicates that the class has a density of habitat area proportional to the size of the class in the map. If the value is greater than 1, then there is a high correlation, and a value of less than 1 means a lower correlation.
Table 3. Logistic regression coefficient between leopard cat and related factors.
Table 3. Logistic regression coefficient between leopard cat and related factors.
FactorClassLogistic Regression CoefficientSignificance Level
Timber TypeNo data0.0000000.000
Farmland2.147514
Larch2.502028
Pinus rigida fores0.948505
Non-stocked forest land2.462551
Chestnut artificial forest4.025706
Cut-over area1.895348
Pinus densiflora forest−18.043562
Pinus densiflora artificial forest3.942215
Pinus koraiensis forest−16.933439
Left-over area0.000000
Oak forest2.252224
Oak artificial forest2.532997
Grassland−22.151046
Conifer mixed forest2.291914
Mixed forest of soft and hardwood0.000000
Poplar forest0.000000
Water0.000000
Broadleaved forest0.000000
Denuded land0.000000
Bamboo stand0.000000
Other0.000000
Timber ageNon forest area−0.3057750.002
1st age0.874957
2nd age−0.260902
3rd age−0.081416
4th age0.087840
5th age−0.655312
6th age0.000000
Distance from forest (m)-0.004420.935
Distance from road (m)-−0.000260.480
Distance from Water (m)-−0.0001330.001
Land CoverNo Data0.0000000.722
Other0.000000
Residential Area1.252586
Manufacturing Area2.305031
Commercial Area2.030913
Recreational Area2.078151
Trafficker Area19.626055
Public Area20.999587
Agricultural Area0.000000
Paddy18.571218
Field18.817534
Growing in Plastic Greenhouse0.000000
Orchard1.801892
Other Plantations3.061980
Broadleaved Forest18.682290
Coniferous Forest19.251522
Mixed Forest19.127733
Forest Area0.000000
Natural Grassland17.532816
Artificial Pasture−0.795913
Other Grassland17.843009
Wetland0.000000
Inland Wetland21.357972
Coastal Wetland0.000000
Bare Land0.000000
Natural Bare Land0.442838
Other Bare Land18.492102
Beach0.000000
Inland Water19.498249
Marine Water0.000000
Ground Elevation (m)-0.0038150.000
Slope gradient (°)-−0.0132090.130
Slope aspect-−0.0016900.004
The spatial databases of each variable were converted to ASCII files using ArcGIS for use in the statistical package SPSS 20. Using this approach, logistic multiple regression coefficients (B), standard errors of slope coefficients (S.E), the Wals tests (Wals), the significance levels (Sig.), and the exponentiated slope coefficients (Exp(B)) of the related variables were calculated (Table 3). The coefficients were estimated using the maximum-likelihood model. Because the relationship between the independent variables and the probability was nonlinear in the logistic multiple regression model, an iterative algorithm was necessary for parameter estimation (Oh [31]). Coefficients denote the meaning of the influences of related factors or classes to habitat potential. A negative value means that the factor or class has a negative effect on the occupancy of the leopard cat at the study site, such as slope gradient, timber type (Pinus densiflora and Pinus koraiensis forest, and grassland), road distance, and distance from water.
Z = TIMBER_TYPE + TIMBER_AGE + (−0.0003 × DIST_ROAD) + (−0.0001 × DIST_WATER) + (0.004 × DIST_FOREST) + LAND_COVER + (0.004 × DEM) + (−0.013 × SLOPE) + (−0.002 × ASPECT) − 21.853
Here, DEM is the intertidal ground elevation value, SLOPE is the slope gradient value, ASPECT is the slope aspect value, DIST_ROAD is the distance from a road, DIST_WATER is the distance from water, and DIST_FOREST is the distance from forest area. TIMBER_TYPE, TIMBER_AGE, and LAND_COVER are the values of each categorical factor in Table 3, and Z is a prediction parameter. Using the logistic regression coefficient (Table 3), the probability of a species was computed and mapped as the habitat potential index (HPI).
Leopard cat habitat potential maps were quantitatively developed using the HPI values. These were calculated using the logistic regression and frequency ratio models for the interpretation (Figure 3). Figure 3a is result of applying frequency ratio (Table 2) and Figure 3b is result of applying logistic regression (Table 3). The index was composed of five classes based on area for easy visual interpretation. Index ranges of very high, high, moderate, low, and very low in 5%, 10%, 15%, 20%, and 50% of the study area, respectively, were used. The classification was useful to visually delineate the predicted habitat potential areas.
Figure 3. Leopard cat habitat potential map (a) Frequency Ratio; (b) Logistic Regression.
Figure 3. Leopard cat habitat potential map (a) Frequency Ratio; (b) Logistic Regression.
Sustainability 07 04668 g003aSustainability 07 04668 g003b

4.3. Validation

A leopard cat habitat potential map should effectively predict future leopard cat potential habitat areas. This could be validated using new potential locations as the cats become distributed. In the study, many locations of leopard cats were detected from survey data. These locations were divided into a training set to analyze the habitat potential using the frequency ratio and logistic regression models, and a validation set to validate the predicted habitat potential map. The leopard cat habitat potential analysis result was validated using a validation set that was not used for training the model. Validation was performed by comparing the known leopard cat distribution locations with the habitat potential maps.
A rate curve was created, and the area under the curve was calculated. The rate explains how well the model and factors predict leopard cat distribution. The area under the curve qualitatively assesses the accuracy of the prediction. To obtain the relative ranks for each prediction pattern, the calculated index values of all cells in the study area were sorted in descending order. The ordered cell values were then divided into 100 classes with accumulated 1% intervals. The rate validation result appears as a line in Figure 4. For example, in the case of the logistic regression model, an index rank above 10% of the HPI could explain 33% of all the leopard locations. To obtain quantitative results, the areas under the curve were recalculated as if the total area were 1.0, which would mean perfect prediction accuracy. Using this method, the area under a curve can be used to assess the prediction accuracy qualitatively. In the case of the frequency ratio, the area under the curve was 0.821, and the prediction accuracy was 82.15%. In the case of the logistic regression, the area under the curve was 0.8148, and the prediction accuracy was 81.48%, as shown in Figure 4.
Figure 4. Success rate curves showing the cumulative percentage of each species occurrence (y-axis) for the descending ordered species potential index (SPI) rank (x-axis).
Figure 4. Success rate curves showing the cumulative percentage of each species occurrence (y-axis) for the descending ordered species potential index (SPI) rank (x-axis).
Sustainability 07 04668 g004

5. Conclusions and Discussion

In this study, we applied the frequency ratio and logistic regression models to habitat potential mapping for the leopard cat. The first step was to select the nine most important variables potentially affecting leopard cat habitat. We then mapped habitat potential using frequency ratio and logistic regression models representing the relationships between leopard cat distribution and environmental variables. We assembled factors associated with habitat potential in a spatial database and created habitat potential maps for the leopard cat. Finally, we validated the maps using location data that had not been used for model training. We arrived at the following conclusions:
(1)
The results of the frequency ratio model indicated that elevation, slope-related factors, and timber age had a positive correlation with locations used by the leopard cat. Areas closer to water and forest and farther from roads, oak forests, and broad-leaved forest classes showed the most positive correlations.
(2)
According to the logistic regression coefficients, the factors of slope gradient, timber type (Pinus densiflora and Pinus koraiensis forests, and Grassland), the distance from roads and distance from water were negatively correlated with the locations used by the leopard cat. In contrast, the factors of ground elevation and distance from a forest had a positive effect on leopard cat habitat potential. Some factors contrasted with the results of the frequency ratio, i.e., slope gradient and distance from water.
(3)
Generally, the maps resulting from the frequency ratio and logistic regression models had similar spatial distribution patterns. The central south and northeastern parts of the inland area of South Korea and the central part of Jeju Island were expected to have high and very high potential. The results of this study can be used in future studies of predator reintroduction on Jeju Island. In particular, the reintroduction of the leopard cat is being considered because it can play the role of top predator on Jeju Island. This study indicated high availability of potential habitat. These areas have high elevation, steep slopes, and forest, and they are hilly or mountainous. Such areas of high and very high potential should be given priority during land-use or wildlife management planning. The western and eastern coastal parts of the site were shown to have low and very low potential in all of the habitat potential maps. Almost all areas in this region are low-lying, with coastal and non-forest habitat.
(4)
Using the frequency ratio and logistic regression models, we created leopard cat habitat potential maps. Half of known leopard cat locations were used as training data and the remaining half was used to validate the maps. The resulting frequency ratio and logistic regression models were 82.15 and 81.48% accurate, respectively. Therefore, the results had an overall agreement of more than 80%, which we regarded as satisfactory.
Some limitations exist to detecting exact leopard cat locations, and the locations used in this study were based on surveys, not exact figures. Inaccurate location data can lead to difficulties in spatial analysis.
The frequency ratio model is somewhat simplistic, but the process of input, calculation, and output can be easily understood. Moreover, large amounts of data can be quickly and easily processed in the GIS environment. The spatial database can be used in other studies. The logistic regression model requires the conversion of the data to ASCII or other formats for use in the statistical package and subsequent re-conversion to incorporate them into the GIS database. It is hard to process the large amounts of data in the statistical package. However, the correlations of leopard cat locations with other factors can be analyzed qualitatively in the statistical package. The frequency ratio model had better accuracy than the logistic regression model used in this study, and the use of all factors produced better results. In the case of a similar statistical model (discriminant analysis), the factors must have a normal distribution, and in the case of multi-regression analysis, the factors must be numerical. However, for logistical regression, the dependent variable must be input as 0 or 1; therefore, the model applies well to habitat potential analysis.
Remote sensing technology and GIS provide ways to introduce information from various data sources into the decision-making process and aid in the handling and manipulation of classified remote sensing data (Adinarayana [32]). Using GIS enables quantitative assessment of the consequences of heterogeneity in ecological systems over a broad range of spatial and temporal scales. The integration of several surface features that indicate mammal habitat potential is an important aspect of ecological management studies.
This study identified factors that may be associated with leopard cat habitat, and our methods and results can also be applied to habitat potential mapping of other mammalian species. Moreover, the resulting habitat potential map can be used as basic data for establishing plans to manage mammalian species, such as locating monitoring sites. However, more case studies and models are needed to generalize factors associated with mammalian habitats.

Acknowledgments

This research was supported by the Basic Research Project of the Korea Institute of Geoscience and Mineral Resources (KIGAM) funded by the Ministry of Science, ICT & Future Planning, and the Korea Environment Institute funded by the Development Project of Environmental Technology for Climate Change by the Korea Environmental Industry & Technology Institute.

Author Contributions

Moung-Jin Lee suggested the idea. In addition, he collected data and processed input data. Wonkyong Song interpreted the input data and this resulted in the habitat aspect. Saro Lee organized the paperwork and applied the frequency ratio and logistic models. All of authors contributed to the writing of each part.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Minor, E.S.; Lookingbill, T.R. A multiscale network analysis of protected-area connectivity for mammals in the united states. Conserv. Biol. 2010, 24, 1549–1558. [Google Scholar] [CrossRef] [PubMed]
  2. Burton, P.; Balisky, A.; Coward, L.; Kneeshaw, D.; Cumming, S. The value of managing for biodiversity. For. Chron. 1992, 68, 225–237. [Google Scholar] [CrossRef]
  3. Myers, N.; Mittermeier, R.A.; Mittermeier, C.G.; da Fonseca, G.A.; Kent, J. Biodiversity hotspots for conservation priorities. Nature 2000, 403, 853–858. [Google Scholar] [CrossRef] [PubMed]
  4. Schmitz, O.J. Top predator control of plant biodiversity and productivity in an old-field ecosystem. Ecol. Lett. 2003, 6, 156–163. [Google Scholar] [CrossRef]
  5. Sergio, F.; Newton, I.; Marchesi, L. Conservation: Top predators and biodiversity. Nature 2005, 436, 192–192. [Google Scholar] [CrossRef] [PubMed]
  6. Gavashelishvili, A.; Lukarevskiy, V. Modelling the habitat requirements of leopard panthera pardus in west and central asia. J. Appl. Ecol. 2008, 45, 579–588. [Google Scholar] [CrossRef]
  7. Sanderson, J.; Sunarto, S.; Wilting, A.; Driscoll, C.; Lorica, R.; Ross, J.; Hearn, A.; Mujkherjee, S.; Khan, J.; Habib, B. Prionailurus bengalensis. In IUCN 2011; IUCN Red List of Threatened Species. Version 2011.2; IUCN: Gland, Switzerland, 2008. [Google Scholar]
  8. Kim, K.D.; Lee, S.; Oh, H.J.; Choi, J.K.; Won, J.S. Assessment of ground subsidence hazard near an abandoned underground coal mine using gis. Environ. Geol. 2006, 50, 1183–1191. [Google Scholar] [CrossRef]
  9. Lee, S.; Sambath, T. Landslide susceptibility mapping in the damrei romel area, Cambodia using frequency ratio and logistic regression models. Environ. Geol. 2006, 50, 847–855. [Google Scholar] [CrossRef]
  10. Matsuura, T.; Sugimura, K.; Miyamoto, A.; Tanaka, N. Knowledge-Based Estimation of Edible Fern Harvesting Sites in Mountainous Communities of Northeastern Japan. Sustainability 2014, 6, 175–192. [Google Scholar] [CrossRef]
  11. Choi, J.K.; Kim, K.D.; Lee, S.; Won, J.S. Application of a fuzzy operator to susceptibility estimations of coal mine subsidence in taebaek city, Korea. Environ. Earth Sci. 2010, 59, 1009–1022. [Google Scholar] [CrossRef]
  12. Song, J.; Wang, X.; Liao, Y.; Zhen, J.; Ishwaran, N.; Guo, H.; Yang, R.; Liu, C.; Chang, C.; Zong, X. An Improved Neural Network for Regional Giant Panda Habitat Suitability Mapping: A Case Study in Ya’an Prefecture. Sustainability 2014, 6, 4059–4076. [Google Scholar] [CrossRef]
  13. Poplar-Jeffers, I.O.; Petty, J.T.; Anderson, J.T.; Kite, S.J.; Strager, M.P.; Fortney, R.H. Culvert replacement and stream habitat restoration: Implications from brook trout management in an appalachian watershed, USA. Restor. Ecol. 2009, 17, 404–413. [Google Scholar] [CrossRef]
  14. Ottaviani, D.; Panzacchi, M.; Jona Lasinio, G.; Genovesi, P.; Boitani, L. Modelling semi-aquatic vertebrates’ distribution at the drainage basin scale: The case of the otter Lutra lutra in italy. Ecol. Model. 2009, 220, 111–121. [Google Scholar] [CrossRef]
  15. Newton-Cross, G.; White, P.C.; Harris, S. Modelling the distribution of badgers meles meles: Comparing predictions from field-based and remotely derived habitat data. Mamm. Rev. 2007, 37, 54–70. [Google Scholar] [CrossRef]
  16. Tien Bui, D.; Lofman, O.; Revhaung, I.; Dick, O. Landslide susceptibility analysis in the Hoa Binh province of Vietnam using statistical index and logistic regression. Nat. Hazards 2011, 59, 1413–1444. [Google Scholar] [CrossRef]
  17. Huck, M.; Davison, J.; Roper, T.J. Predicting european badger meles meles sett distribution in urban environments. Wildl. Biol. 2008, 14, 188–198. [Google Scholar] [CrossRef]
  18. Jaberg, C.; Guisan, A. Modelling the distribution of bats in relation to landscape structure in a temperate mountain environment. J. Appl. Ecol. 2001, 38, 1169–1181. [Google Scholar] [CrossRef]
  19. Greaves, G.J.; Mathieu, R.; Seddon, P.J. Predictive modelling and ground validation of the spatial distribution of the new zealand long-tailed bat (Chalinolobus tuberculatus). Biol. Conserv. 2006, 132, 211–221. [Google Scholar] [CrossRef]
  20. Clement, M.J.; Castleberry, S.B. Estimating density of a forest-dwelling bat: A predictive model for rafinesque’s big-eared bat. Popul. Ecol. 2013, 55, 205–215. [Google Scholar] [CrossRef]
  21. Northrup, J.; Stenhouse, G.; Boyce, M. Agricultural lands as ecological traps for grizzly bears. Anim. Conserv. 2012, 15, 369–377. [Google Scholar] [CrossRef]
  22. Kuemmerle, T.; Perzanowski, K.; Chaskovskyy, O.; Ostapowicz, K.; Halada, L.; Bashta, A.-T.; Kruhlov, I.; Hostert, P.; Waller, D.M.; Radeloff, V.C. European bison habitat in the carpathian mountains. Biol. Conserv. 2010, 143, 908–916. [Google Scholar] [CrossRef]
  23. Speed, J.D.; Woodin, S.J.; Tommervik, H.; Tamstorf, M.P.; van der Wal, R. Predicting habitat utilization and extent of ecosystem disturbance by an increasing herbivore population. Ecosystems 2009, 12, 349–359. [Google Scholar] [CrossRef]
  24. Kocev, D.; Džeroski, S.; White, M.D.; Newell, G.R.; Griffioen, P. Using single-and multi-target regression trees and ensembles to model a compound index of vegetation condition. Ecol. Model. 2009, 220, 1159–1168. [Google Scholar] [CrossRef]
  25. Meixler, M.S.; Bain, M.B.; Todd Walter, M. Predicting barrier passage and habitat suitability for migratory fish species. Ecol. Model. 2009, 220, 2782–2791. [Google Scholar] [CrossRef]
  26. Walker, B.K.; Jordan, L.K.; Spieler, R.E. Relationship of reef fish assemblages and topographic complexity on southeastern florida coral reef habitats. J. Coast. Res. 2009, 53, 39–48. [Google Scholar] [CrossRef]
  27. Walters, D.; Roy, A.; Leigh, D. Environmental indicators of macroinvertebrate and fish assemblage integrity in urbanizing watersheds. Ecol. Indic. 2009, 9, 1222–1233. [Google Scholar] [CrossRef]
  28. Bonham-Carter, G. Geographic Information Systems for Geoscientists: Modelling with GIS; Elsevier: Amsterdam, The Netherlands, 1994. [Google Scholar]
  29. Dai, F.; Lee, C. Landslide characteristics and slope instability modeling using GIS, Lantau Island, Hong Kong. Geomorphology 2002, 42, 213–228. [Google Scholar] [CrossRef]
  30. Lee, S.; Min, K. Statistical analysis of landslide susceptibility at Yongin, Korea. Environ. Geol. 2001, 40, 1095–1113. [Google Scholar] [CrossRef]
  31. Oh, H.J.; Lee, S.; Chotikasathien, W.; Kim, C.H.; Kwon, J.H. Predictive landslide susceptibility mapping using spatial information in the pechabun area of Thailand. Environ. Geol. 2009, 57, 641–651. [Google Scholar] [CrossRef]
  32. Adinarayana, J.; Krishna, N.R. Integration of multi-seasonal remotely-sensed images for improved landuse classification of a hilly watershed using geographical information systems. Int. J. Remote Sens. 1996, 17, 1679–1688. [Google Scholar] [CrossRef]

Share and Cite

MDPI and ACS Style

Lee, M.-J.; Song, W.; Lee, S. Habitat Mapping of the Leopard Cat (Prionailurus bengalensis) in South Korea Using GIS. Sustainability 2015, 7, 4668-4688. https://0-doi-org.brum.beds.ac.uk/10.3390/su7044668

AMA Style

Lee M-J, Song W, Lee S. Habitat Mapping of the Leopard Cat (Prionailurus bengalensis) in South Korea Using GIS. Sustainability. 2015; 7(4):4668-4688. https://0-doi-org.brum.beds.ac.uk/10.3390/su7044668

Chicago/Turabian Style

Lee, Moung-Jin, Wonkyong Song, and Saro Lee. 2015. "Habitat Mapping of the Leopard Cat (Prionailurus bengalensis) in South Korea Using GIS" Sustainability 7, no. 4: 4668-4688. https://0-doi-org.brum.beds.ac.uk/10.3390/su7044668

Article Metrics

Back to TopTop