Fragmented Riverine Habitats in Taiwan Have Spatio-Temporal Consequences , ReDistributing Caprimulgus affinis into Urban Areas Leading to a Human – Wildlife Conflict

Caprimulgus affinis is an endemic species commonly known as the Taiwan nighthawk, forest nighthawk, or South Asian nighthawk. In recent years, the C. affinis population has gradually moved outward from river habitats into the metropolitan areas of Taiwan. Because male birds’ booming sounds at night can reach up to 90 dB or higher and they can be intermittently tweeting for more than 10 hours, they often disturb the sleep of urban residents and can even cause nightmares. In this study, we analyzed long-term survey data to assess the distribution history of C. affinis in Taiwan. By using 1738 entries of observation data collected from 1999 to 2014, a model of C. affinis distribution was constructed, and a geostatistical method was used to improve the accuracy of the model estimate. In addition, the Mann–Kendall trend test was applied to predict future C. affinis distribution. Based on the results, four variables were selected to construct the C. affinis distribution matrix. Urban population, coastal proximity, and distance to upstream river location represented positive driving forces, whereas city elevation was a negative driving force. C. affinis is currently distributed on the plains of central and southern Taiwan and in eastern Taiwan. The C. affinis emergence trend diagram was plotted using the time–space trend diagram, which showed yearly increases in the C. affinis populations in urban and settled areas including central, western, northeastern, and southeastern Taiwan and yearly decreases in the populations in rural areas including northern and southern Taiwan, especially in fragmented riverine habitats. Regression kriging can correctly describe the distribution of the entire C. affinis population, which leads to the correct understanding of the biological corridor of C. affinis in their migration through the graphical contours in GIS. The analytic model in this study contributes to the establishment of the time–space trend diagram.


Introduction
In recent years, the high-decibel booming calls of Caprimulgus affinis Horsfield, 1821 are now being heard at night during spring and summer in metropolitan areas of Taiwan.These birds are distributed around cities and substantially impact the lives of urban residents, who have often resorted to flashlights, firecrackers, and other measures to drive away the birds.However, these measures have been ineffective.Therefore, C. affinis has become a new hazard for urban residents.In addition, because airports in Taiwan are located near cities, increased bird strike incidents with aircraft in recent years have caused aviation safety problems.
Caprimulgus affinis is one of ten subspecies with an average prevalence in Taiwan that is also called the Savanna nightjar, South Asian nighthawk or "breakwater earner".Body length, weight, and wing span are approximately 20-26 cm, 54-110 g, and 64 cm, respectively.First described in Tainan City in Taiwan by Robert Swinhoe in 1863 [1], C. affinis is widely distributed in Taiwan, India, Indonesia, central China, and Pakistan, besides other areas [1].In its native habitat, the yellowish-brown plumage provides superior camouflage from predators, resembling a rock when in its preferred sitting position on the ground.Having a short beak and wide rictus with a huge oral cavity, Savanna nightjars are often described as "hideous" creatures that make "disturbing booming sounds".Due to their unflattering appearance, humans rarely deliberately capture these birds [1][2][3].Nocturnal in nature, they feed on flying insects-their preferred food-in the evening.
It appears that the male birds are a nuisance.Their courtship behavior to attract females is established by the male birds' vocalization in mating strategies [4][5][6][7].The male birds also defend their territories using their squeaky sounds in general breeding behaviors.The female birds, therefore, seldom contribute to the noise [4], doing so only while they defend their solitary territories during the hatching and brooding phase.Possessing the same territorial breeding behaviors as female birds, every male remains in its own territory guarding and tending to ground nests, which are in low-elevation riverbed gravel, rocky soil, or open spaces.It is during the breeding season, which extends from January to July, that males produce their shrill booming sounds or continual caw calls while perched or in flight, respectively.In comparison, in the urban environment the males perch on the rooftops of buildings while making these well-recognized booming calls.A male bird can produce intermittent booming and bubbling trills of approximately 2 s in duration.These nocturnal sounds can reach up to 4.2 kHz, 90 decibels (dB) or higher, and can last for more than half an hour [5][6][7].However, male birds' booming sounds can go on for more than 10 hours as detected from our overnight audio recordings.
There are no significant differences between the nightjars currently living in the city and those from rural areas.That is both have "high-decibel tweets" according to Sharps et al. and Lowe et al. [8,9].Caprimulgus affinis originally occupied low-altitude open grasslands, rocky soil, wheat fields, sugar cane plantations, and around streams and rivers in forests.Because these locations are far from cities and towns, the breeding calls of the birds were not a disturbance to humans.Although recent surveys have shown C. affinis populations to be growing, reproductive success rates continue to be low in natural environments due to environmental impacts (e.g., flooding) and predation (e.g., cats, dogs, rats, snakes, and birds such as crested goshawks, Accipiter trivirgatus, and collared scops-owls, Otus lettia glabripes).In contrast, "urbanized" C. affinis birds prefer to lay eggs on the flat roofs of penthouses in sparsely populated areas and particularly on the roofs of schools or larger institutions.Given the lack of natural predators, the urban environment provides safer breeding grounds and increased fecundity.Other than disturbances from humans and rodents, the absence of threats from natural predators enhances breeding success rates and contributes to higher reproduction rates for the species [10][11][12].In addition, the phenomenon of nightjars' "urbanization" has been observed, and the birds have begun to enter metropolitan areas [10][11][12][13][14][15][16].This phenomenon of "urbanization" is occurring in other parts of the world where C. affinis exists [17].
Although riverbeds are the preferred habitats [18], in some cases, the concentration of nightjars in cities has been found to be higher than in the surrounding riverbeds and riverine sites [19][20][21][22][23][24].River fragmentation from natural and anthropogenic influences, such as flooding, dredging, channelization, and dams has influenced water flow regimes throughout the country [25][26][27][28].The resulting disturbed riparian and aquatic habitats may lead to shifts in avian distribution for breeding and foraging sites [27].Avian communities are sensitive to spatial and phenological mismatches with aquatic prey following flow disturbances in rivers [28].For instance, disturbance of gravel and sand along the river edge has destroyed natural habitat of C. affinis.However, the effects of river fragmentation on C. affinis distribution are still poorly understood and, to our knowledge, there have been no related reports in Taiwan or even in Asia on this topic.Therefore, this is the first study performed to evaluate trends in C. affinis' distribution history to determine whether this species has changed habitats due to disturbances.The two purposes of this study were to: (1) determine current and past C. affinis distributions and (2) predict the likelihood of C. affinis sightings in major cities in the future by using past data.

Study Area
Taiwan has a population of about 23.5 million people and the annual population growth rate is 1.22%.Taiwan is approximately 36,000 square kilometers in size, with 32% of the whole island having mean sea level elevations higher than 1000 meters (Figure 1).Distribution of the species based on elevation is as follows: 70.5% of the population occurs at elevations between 0 and 50 m, 27.0% occurs between 50 and 500 m, and 2.5% occurs between 500 and 600 m.The average annual precipitation in Taiwan is about 2150 mm, but rainfall in northern Taiwan is about 20% higher than that in the southern Taiwan.Most of this precipitation occurs from typhoons during the wet season that extends from May to October.However, the spatial and temporal distribution of precipitation is extremely uneven, leading to great differences in river flows during the dry and wet seasons.Although riverbeds are the preferred habitats [18], in some cases, the concentration of nightjars in cities has been found to be higher than in the surrounding riverbeds and riverine sites [19][20][21][22][23][24].River fragmentation from natural and anthropogenic influences, such as flooding, dredging, channelization, and dams has influenced water flow regimes throughout the country [25][26][27][28].The resulting disturbed riparian and aquatic habitats may lead to shifts in avian distribution for breeding and foraging sites [27].Avian communities are sensitive to spatial and phenological mismatches with aquatic prey following flow disturbances in rivers [28].For instance, disturbance of gravel and sand along the river edge has destroyed natural habitat of C. affinis.However, the effects of river fragmentation on C. affinis distribution are still poorly understood and, to our knowledge, there have been no related reports in Taiwan or even in Asia on this topic.Therefore, this is the first study performed to evaluate trends in C. affinis' distribution history to determine whether this species has changed habitats due to disturbances.The two purposes of this study were to: (1) determine current and past C. affinis distributions and (2) predict the likelihood of C. affinis sightings in major cities in the future by using past data.

Study Area
Taiwan has a population of about 23.5 million people and the annual population growth rate is 1.22%.Taiwan is approximately 36,000 square kilometers in size, with 32% of the whole island having mean sea level elevations higher than 1000 meters (Figure 1).Distribution of the species based on elevation is as follows: 70.5% of the population occurs at elevations between 0 and 50 m, 27.0% occurs between 50 and 500 m, and 2.5% occurs between 500 and 600 m.The average annual precipitation in Taiwan is about 2150 mm, but rainfall in northern Taiwan is about 20% higher than that in the southern Taiwan.Most of this precipitation occurs from typhoons during the wet season that extends from May to October.However, the spatial and temporal distribution of precipitation is extremely uneven, leading to great differences in river flows during the dry and wet seasons.

Methods
To sample the entire island nation and account for C. affinis having various degrees of mobility, we used stratified random sampling methods [25] suitable for C. affinis niches from sea level up to 1000 m elevation.We conducted monthly surveys at 335 locations across Taiwan between January 1999 to December 2014 (Figure 1) and developed a database for the 16 years of field surveys that were conducted by 234 field birdwatchers, who were members of a citizen science training program, which we managed along with the birdwatchers from the Chinese Wild Bird Federation, Taiwan, R.O.C.More than 45,000 observation hours went into the field survey component of this study.To ensure quality control and data reliability, the records in the database were proofread and verified by the authors.The database included 60 riparian areas (Appendix A) for all river names and their total lengths (TLs = 3696.5km, 2310.3 miles) of riparian sections along these rivers in the monitored areas.Using the observation data, coordinate calibration and cross-matching were performed.Corresponding data for every urban and rural settlement (n = 335) in Taiwan were compiled to establish a record of C. affinis sightings from 1999 to 2014.Surveys for C. affinis were conducted from evening to the following morning because this nocturnal species roosts and sleeps during the day.Sampling data were collected and no birds were harmed or disturbed.All survey data are maintained in an open database with free access upon request.All available data (i.e., 1738 entries of observation data) were mapped to 335 urban and rural settlements.No specific permissions were required for field permits in this study and no prior approval was required by the Council of Agriculture, Executive Yuan, R.O.C.
Caprimulgus affinis distributions in urban and rural settlements were set as analysis locations in the logistic regression model as well as the following variables: urban and rural settlement populations, coastal proximity, city elevation, urban and rural settlement location, distance to main river location (calculating the shortest distance from main river), and distance to upstream river location (calculating the shortest distance from the source of the main river).Logistic regression was conducted to analyze C. affinis distribution and the aforementioned variables.The distance from the edge between an urban and rural settlement to a peripheral location was defined as the short and direct distance from a core point of the urban and rural settlement to a line (i.e., distance to coastline defined as the shortest distance to nearest ocean coast, distance to main river location, and distance to upstream river location).A stepwise regression was performed using SPSS software.During the analyses, the recursive loop continued to calculate the weight of each variable and the insufficiently heavy variables were eliminated to retain the more influential variables; the significance of the variables was also calculated.
The Mann-Kendall trend test [29,30] was used to explore any changes in the distribution of C. affinis from 1999 to 2014.Indicator kriging estimates the probability that the concentration of a pollutant exceeds a specific threshold value at a given location [29][30][31].The random function is then modeled as a combination of the trend and a random variable.Regression kriging combines these two approaches; regression is used to fit the explanatory variation and simple kriging, with the expected value of zero used to fit the residuals (i.e., unexplained variation) [29].In addition, regression kriging is mathematically equivalent to universal kriging and kriging with external drift [32,33].
The Mann-Kendall trend test estimated only space distribution by using the observation data of C. affinis and the kriging method was used to estimate the incidence of C. affinis sightings.To enhance estimation accuracy, the logistic regression results were analyzed using the regression kriging method.Thus, in addition to the analysis of the spatial distribution of C. affinis, other assessment factors could be added to the estimation analysis (i.e., urban population, coastal proximity, city elevation, distance to main river location, and distance to upstream river location).In the present study, kriging and indicator kriging were performed using Geostatistical Software Library and related software (GSLIB) [31].
We, therefore, designed a first situational category and a second situational category to describe the following situations.The estimation results of the first situational category and the regression kriging method underwent correctness analysis [25,[33][34][35].We assessed whether the estimation model would be more accurate from regression kriging associated with indicator kriging method than that of logistic regression after the driving forces analysis.
A receiver operating characteristic curve (ROC curve) was performed in this model.This graph shows the performance of a classification model at all classification thresholds.This curve plots true positive rate (TPR) and false positive rate (FPR) at different classification thresholds.Lowering the classification threshold classifies more items as positive, thus increasing both false positives and true positives.To compute the points on a ROC curve, we could evaluate a logistic regression model many times with different classification thresholds.Therefore, the area under the curve (AUC under the ROC curve) is an efficient, sorting-based algorithm that can provide this information for us.That is, AUC measures the entire two-dimensional area underneath the entire ROC curve (think integral calculus) from (0,0) to (1,1).AUC ranges in value from 0 to 1.A model whose predictions are 100% wrong has an AUC of 0.0; one whose predictions are 100% correct has an AUC of 1.0.

Driving Forces of Diversity Detected Using Logistic Regression
Caprimulgus affinis observations were conducted in 335 urban and rural settlements, of which 213 had no record of C. affinis sightings and 122 did (Figure 1).The data are primarily from the southern and central parts of western Taiwan and only a few data records were from northern and eastern Taiwan.The analysis steps and results are shown in Table 1.Four variables were selected to construct the C. affinis distribution matrix: urban population, coastal proximity (the distance to the nearest ocean coast), and distance to upstream river location (calculating the shortest distance from the source of the main river) represented positive driving forces; and city elevation represented a negative driving force.The main-river distance parameter (calculating the shortest distance from a main river) was not included as a variable in the equation because it was not significant.The results of the driving-force analysis are consistent with the biological characteristics of C. affinis, which indicates the species prefers urban habitats rather than the ocean coast.The main river distance parameter was non-significant and was not included in the LR equation.b,c This is the Wald chi-square that tests the null hypothesis that the constant equals zero.This hypothesis is rejected because the p value (listed in the column called "Sig.") is smaller than the critical p value of 0.05 (or 0.01, or 0.001).Hence, we conclude that the constant is not zero.* p < 0.05, ** p < 0.01, *** p < 0.001.

Estimation of C. affinis Space Distribution Using Indicator Kriging Method
Using the index model, the sightings of C. affinis underwent situational categorization as follows: First situational category: If C. affinis is sighted, then this is coded as 1, otherwise it is coded as 0 (single-threshold method).
Second situational category: The frequency of C. affinis sightings from 1999 to 2014 was calculated.For example, the frequency of one sighting was calculated as 1/16 = 0.0625 and that of eight sightings was calculated as 8/16 = 0.5 (multi-threshold method).
The semi-variogram results of the incidence of C. affinis sightings are shown in Table 2.The first and second situational categories were suitable for the index model.The C 0 of the first situational category should be smaller than that of the second situational category.The impact range of the first situational category was approximately nine times smaller than that of the second situational category.This finding indicated a more suitable space structure for the first situational category and showed regional location was the main influencing factor.By contrast, the second situational category provided a greater estimation result and more suitable outcomes for long distances.The estimation results of the two categories are shown in Figure 2. The basis of the first situational category (Figure 2a) is all of the locations of C. affinis sightings, which shows the species is spatially clustered and distributed mostly in the southern and central areas of Taiwan with few occurrences in northern Taiwan.However, the basis of the second situational category (Figure 2b) is the frequency of sightings.The diagram shows central Taiwan possesses the most C. affinis sightings, which agrees with the current phenomenon of C. affinis disturbing people's sleep and the current location of this species.(single-threshold method).
Second situational category: The frequency of C. affinis sightings from 1999 to 2014 was calculated.For example, the frequency of one sighting was calculated as 1/16 = 0.0625 and that of eight sightings was calculated as 8/16 = 0.5 (multi-threshold method).
The semi-variogram results of the incidence of C. affinis sightings are shown in Table 2.The first and second situational categories were suitable for the index model.The C0 of the first situational category should be smaller than that of the second situational category.The impact range of the first situational category was approximately nine times smaller than that of the second situational category.This finding indicated a more suitable space structure for the first situational category and showed regional location was the main influencing factor.By contrast, the second situational category provided a greater estimation result and more suitable outcomes for long distances.The estimation results of the two categories are shown in Figure 2. The basis of the first situational category (Figure 2a) is all of the locations of C. affinis sightings, which shows the species is spatially clustered and distributed mostly in the southern and central areas of Taiwan with few occurrences in northern Taiwan.However, the basis of the second situational category (Figure 2b) is the frequency of sightings.The diagram shows central Taiwan possesses the most C. affinis sightings, which agrees with the current phenomenon of C. affinis disturbing people's sleep and the current location of this species.

Estimation of C. affinis Space Distribution Using Regression Kriging Method
The estimation results are shown in Figure 3, which reveals several phenomena.First, the sightings mostly occurred in central, southern, and eastern (i.e., Lanyang Plains and Hualien County) Taiwan.Second, the line of equal probability for the incidence of C. affinis sightings is partially horizontal in the space.This finding indicates the range of bird activity is less limited and confirmed that a biological corridor exists for C. affinis.Third, although urban population was selected as a parameter for determining the C. affinis distribution matrix and represented a positive driving force, C. affinis was rarely observed in Taipei City (see Figure 1).Thus, C. affinis migration to northern Taiwan did not appear to be strong.
C. affinis was rarely observed in Taipei City (see Figure 1).Thus, C. affinis migration to northern Taiwan did not appear to be strong.
The results show the estimation result of the first situational category had an AUC of 0.688.The estimation result of the regression kriging method had an AUC of 0.765.Therefore, the regression kriging method yielded slightly higher results.The ROC curve is shown in Figure 4.The results show the estimation result of the first situational category had an AUC of 0.688.The estimation result of the regression kriging method had an AUC of 0.765.Therefore, the regression kriging method yielded slightly higher results.The ROC curve is shown in Figure 4.

Time-Space Distribution of C. affinis Sightings
This study used observation data of C. affinis (continuity data collected over 16 years and shown in Figure 1) to calculate the incidence trend of C. affinis sightings in urban and rural settlements (n = 335).Using the Mann-Kendall trend analytic method [29,30,37,38], the τ and p values were obtained.The τ value indicates whether the 16-year data show an upward or downward trend.The p value indicates whether the test data volume is sufficient to show the trend.The trends histogram of C. affinis sightings from 1999 to 2014 revealed C. affinis sightings in towns and cities showed an upward trend (Figure 5).The data for 26 locations in Taiwan showed a significantly upward trend (τ > 0 and p value < 0.05) [31][32][33][34], and the sightings were mostly concentrated in central Taiwan.It can be seen from the reduction of the overall trends in the south and the increase in the trend in the central region that C. affinis had moved north from the southern parts of Taiwan and these trends might be related

Time-Space Distribution of C. affinis Sightings
This study used observation data of C. affinis (continuity data collected over 16 years and shown in Figure 1) to calculate the incidence trend of C. affinis sightings in urban and rural settlements (n = 335).Using the Mann-Kendall trend analytic method [29,30,36,37], the τ and p values were obtained.The τ value indicates whether the 16-year data show an upward or downward trend.The p value indicates whether the test data volume is sufficient to show the trend.The trends histogram of C. affinis sightings from 1999 to 2014 revealed C. affinis sightings in towns and cities showed an upward trend (Figure 5).The data for 26 locations in Taiwan showed a significantly upward trend (τ > 0 and p value < 0.05) [31][32][33][34], and the sightings were mostly concentrated in central Taiwan.It can be seen from the reduction of the overall trends in the south and the increase in the trend in the central region that C. affinis had moved north from the southern parts of Taiwan and these trends might be related to the population density of people.Since the onset of the 21st century, C. affinis numbers have increased in the central/western regions of Taiwan, predominantly in settlements.

Time-Space Distribution of C. affinis Sightings
This study used observation data of C. affinis (continuity data collected over 16 years and shown in Figure 1) to calculate the incidence trend of C. affinis sightings in urban and rural settlements (n = 335).Using the Mann-Kendall trend analytic method [29,30,37,38], the τ and p values were obtained.The τ value indicates whether the 16-year data show an upward or downward trend.The p value indicates whether the test data volume is sufficient to show the trend.The trends histogram of C. affinis sightings from 1999 to 2014 revealed C. affinis sightings in towns and cities showed an upward trend (Figure 5).The data for 26 locations in Taiwan showed a significantly upward trend (τ > 0 and p value < 0.05) [31][32][33][34], and the sightings were mostly concentrated in central Taiwan.It can be seen from the reduction of the overall trends in the south and the increase in the trend in the central region that C. affinis had moved north from the southern parts of Taiwan and these trends might be related to the population density of people.Since the onset of the 21st century, C. affinis numbers have increased in the central/western regions of Taiwan, predominantly in urban settlements.Although the main habitats of C. affinis are the Kaohsiung and Pingtung riverbeds in southern Taiwan [38,39], it is well documented that the species has migrated to western and eastern Taiwan [40][41][42].While urban river construction projects have impacted the birds' original habitats [25], suitable habitat appears to have been inadvertently created for C. affinis (or selected/deemed appropriate by the species) in urban areas.The large number of artificial urban river construction projects could be an indirect cause of C. affinis migration, but whether habitat change, degradation or elimination, change in breeding habitat, or reduced food supplies (i.e., moths) caused C. affinis, one of the important nocturnally insectivorous birds [43,44], to migrate from the countryside to metropolitan areas is not known and requires further investigation.[41][42][43].While urban river construction projects have impacted the birds' original habitats [25], suitable habitat appears to have been inadvertently created for C. affinis (or selected/deemed appropriate by the species) in urban areas.The large number of artificial urban river construction projects could be an indirect cause of C. affinis migration, but whether habitat change, degradation or elimination, change in breeding habitat, or reduced food supplies (i.e., moths) caused C. affinis, one of the important nocturnally insectivorous birds [44,45], to migrate from the countryside to metropolitan areas is not known and requires further investigation.

Conclusions
In this study, long-term (16 years) observational field survey data were used to determine the current and past distribution of C. affinis in Taiwan and analyze distribution changes for the species (n = 1738 observations of C. affinis) in 335 urban and rural settlements.Based on the results, three driving forces were determined to be important: urban population and far coastal proximity

Figure 1 .
Figure 1.(a) Study area and location of Caprimulgus affinis sightings; (b) C. affinis sightings superimposed on human density values.

Figure 2 .
Figure 2. Probability of C. affinis sightings.(a) First situational category and (b) Second situational category.

Figure 2 .
Figure 2. Probability of C. affinis sightings.(a) First situational category and (b) Second situational category.

Figure 3 .
Figure 3. Estimation of the probability of C. affinis sighting distribution using regression kriging.

Figure 3 .
Figure 3. Estimation of the probability of C. affinis sighting distribution using regression kriging.

Figure 5 .
Figure 5. Distribution of the τ value, which shows the upper probability of C. affinis sightings in 1999-2014.

Figure 5 .
Figure 5. Distribution of the τ value, which shows the upper probability of C. affinis sightings in 1999-2014.The incidence trends of C. affinis sightings in Taiwan are presented in Figure6.Larger τ values indicate increased C. affinis sightings, while negative τ values show reducing sighting trends.Figure6brepresents whether the sighting trend of C. affinis in Taiwan is statistically significant.The smaller the p value, the better the confidence of the τ value.Caprimulgus affinis sightings in central, northeastern, and southeastern Taiwan show upward/increased sighting trends, which is consistent with the observed urbanization of C. affinis in these regions in recent years.Caprimulgus affinis sightings in northern and southern Taiwan showed decreased sighting trends, which is consistent with observations and historical records.At present, the number of C. affinis in the central region is large, widely distributed, and has an upward trend (p < 0.05).It is statistically significant and suggests the number of C. affinis in the central region is not likely to decrease in the future.The p value obtained from the Mann-Kendall trend analytic method indicates the significance of the data[31,32].Although the main habitats of C. affinis are the Kaohsiung and Pingtung riverbeds in southern Taiwan[38,39], it is well documented that the species has migrated to western and eastern Taiwan[40][41][42].While urban river construction projects have impacted the birds' original habitats[25], suitable habitat appears to have been inadvertently created for C. affinis (or selected/deemed appropriate by the species) in urban areas.The large number of artificial urban river construction projects could be an indirect cause of C. affinis migration, but whether habitat change, degradation or elimination, change in breeding habitat, or reduced food supplies (i.e., moths) caused C. affinis, one of the important nocturnally insectivorous birds[43,44], to migrate from the countryside to metropolitan areas is not known and requires further investigation.

Figure 6 .
Figure 6.(a) τ map of the probability of C. affinis sightings and (b) p map of the probability of C. affinis sightings.

Figure 6 .
Figure 6.(a) τ map of the probability of C. affinis sightings and (b) p map of the probability of C. affinis sightings.

Table 1 .
Variables used in the logistic regression (LR).
RSS: Residual sum of square; C 0 : Nugget.In the early development of geostatistics, the term 'nugget effect' was coined for the apparent discontinuity at the beginning of many semi-variogram graphs; C 0 + C (sill): The value at which the model first flattens out; C/C 0 + C is the nugget coefficient.
RSS: Residual sum of square; C0: Nugget.In the early development of geostatistics, the term 'nugget effect' was coined for the apparent discontinuity at the beginning of many semi-variogram graphs; C0 + C (sill): The value at which the model first flattens out; C/C0 + C is the nugget coefficient.