Translator Disclaimer
15 October 2021 Effectiveness of signs of activity as relative abundance indices for wild boar
Daishi Higashide, Takeo Kuriyama, Shun Takagi, Yoshihiro Nakashima, Keita Fukasawa, Gota Yajima, Minoru Kasada, Mayumi Yokoyama
Author Affiliations +
Abstract

Population growth in wild boars and feral pigs Sus scrofa has negative environmental and economic implications worldwide. Accordingly, it is necessary to monitor population trends for appropriate management. Despite the potential for bias, relative abundance indices based on signs of activity have the potential to be practical, low-cost monitoring tools for data collection at a local scale and over large areas. However, few studies have examined the effectiveness of specific signs as measures of relative wild boar abundance in comparison with reliable density estimates. In this study, we examined whether three activity signs (digging marks, rubbing marks and fecal pellet groups) could be used as relative abundance indices for wild boar. In particular, we conducted transect surveys for signs of activity as well as camera trap surveys from September 2017 to January 2018 at six sites in Hyogo Prefecture and eight sites in Chiba Prefecture, Japan. We modeled the relations between the number of activity signs and wild boar densities or abundances estimated from camera trap data in a hierarchical framework. Fecal pellet counts is a well-established method for estimating the abundance of herbivores, including wild boar; however, we found few fecal pellet groups in our study area, and the counts were not related to wild boar densities. Instead, we found that digging marks were strongly associated with estimates of wild boar density. Although fluctuations in the number of digging marks may be affected by factors other than boar density, including environmental conditions and seasonality, our results suggest that digging marks are an effective relative abundance index for evaluations of the spatial patterns of wild boar in Japan.

Introduction

The wild boar Sus scrofa is one of the most widely distributed terrestrial mammals in the world, and its distribution has expanded both naturally and artificially by management and intentional introductions (Oliver and Leus 2008). Additionally, wild boar populations have increased over the past few decades in various parts of the world (e.g. Europe: Massei et al. 2015; United States: Snow et al. 2017). Wild boar populations, depending on the distribution and density, may have diverse ecosystem and economic impacts (Barrios-Garcia and Ballari 2012). For example, crop damage (Schley et al. 2008), the spread of diseases (Gortázar et al. 2007), and changes in the composition of plant and animal communities (Massei and Genov 2004) have been reported.

In Japan, wild boar populations have increased and the species range has expanded over the past few decades (MOE 2015). The Japanese government estimated that wildlife damage to crops in 2016 resulted in losses of 17.2 billion yen, of which 5.1 billion yen was caused by wild boars (MAFF 2018). In 2013, the Japanese government set a target of reducing populations of free-ranging wild boars and sika deer Cervus nippon in half by 2023 (MOE and MAFF 2013). To achieve this goal, each prefecture has conducted culling based on the Specified Wildlife Conservation and Management Plans for wild boar and sika deer. Furthermore, in September 2018, classical swine fever (CSF) reemerged in Japan after 26 years, affecting both domestic pigs S. scrofa domesticus and wild boars (Postel et al. 2019). The delivery of bait vaccines is being promoted as a control measure to limit the spread of CSF in wildlife (Isoda et al. 2020).

In most cases, the resources (funds, materials and labor) available for these measures are limited; thus, efficient strategies are needed to meet management goals. In particular, spatial culling quotas according to local densities are effective in wildlife management, including wild boar management (McMahon et al. 2010, Pepin et al. 2017). Therefore, reliable estimates of population abundance and trends are important for the adaptive management of wildlife (Nichols and Williams 2006).

Various methods have been proposed to estimate the absolute density or relative abundance of wild boars (Engeman et al. 2013, Enetwild Consortium et al. 2018). For example, nocturnal distance sampling with thermal imagers (Franzetti et al. 2012, Focardi et al. 2020) and capture–mark–recapture with ear tags (Hebeisen et al. 2008) or DNA (Ebert et al. 2012) are used to estimate the absolute density of wild boar. Recently, methods for the estimation of the absolute density of mammals from camera trap data without individual identification have been devised and used, such as the random encounter model (REM; Rowcliffe et al. 2008), spatially explicit N-mixture model (Chandler and Royle 2013), distance sampling (CTDS; Howe et al. 2017), and random encounter and staying time model (REST; Nakashima et al. 2018), and some of these have been applied to wild boar (Palencia et al. 2021). However, methods for estimating absolute density are often not suitable for monitoring over large spatial areas (e.g. prefectures, regions or states) or over long time periods because they are labor-intensive and expensive. In these cases, relative abundance indices (proxies) are widely used owing to relatively low costs and labor requirements (Schwarz and Seber 1999, Stephens et al. 2015). This approach assumes that the proxy is directly related to animal density and the relation is constant over space and/or time. An increasing dependence on proxies without the verification of prerequisites has been criticized (Pollock et al. 2002). However, some studies have shown that proxies are linearly associated with true abundances (Nimmo et al. 2015). Additionally, Iijima et al. (2013) suggested that evaluations of population dynamics based on multiple proxies and the number of captured animals are effective for estimating deer population sizes. In any case, relationships between proxies and population density should be verified to improve their use in wildlife management.

Hunting data are commonly used to evaluate relative abundance (Acevedo et al. 2007, Enetwild Consortium et al. 2018). However, these data have some limitations; for example, they can be collected only when/where hunting occurs and they may be skewed by hunter-related bias (Acevedo et al. 2007, Enetwild Consortium et al. 2018). Alternatively, relative abundance indices based on signs of mammalian activity are effective monitoring tools. Feces (pellets, dung and scat) are typically used to estimate the relative abundance of wild boar (Vicente et al. 2004, Acevedo et al. 2007, Plhal et al. 2014, Ferretti et al. 2016), and Abaigar et al. (1994) used digging (rooting) and tracks. Ideally, indices of relative abundance should have a linear relationship with the true abundance and this relationship should not change spatially or temporally (Enetwild Consortium et al. 2018). However, few studies have examined the effectiveness of signs of activity as measures of relative wild boar abundance by comparisons with reliable estimates of density (Anderson and Stone 1994).

In this study, we examined the effectiveness of various signs of activity as relative abundance indices to evaluate spatial differences in local densities of wild boars. Among signs of activity derived from wild boars (Allwin et al. 2015), we focused on digging marks, rubbing marks and fecal pellet groups, which are likely to be present in forested, warm and temperate areas. We conducted transect surveys and camera trap surveys simultaneously at 14 sites with different wild boar densities in two prefectures in Japan. We evaluated the effectiveness of the three signs as relative abundance indices of wild boars by determining relationships with local densities and abundances estimated from camera trap data.

Material and methods

Study area

The field study was conducted at six sites (H1–H6) in Hyogo Prefecture and at eight sites (C1–C8) in Chiba prefecture, Japan (Fig. 1). The Hyogo prefectural government employed a rectangular grid system (2′30″ in latitude, 3′45″ in longitude, with cells of about 4.6 × 5.5 km) established by the Ministry of the Environment of Japan for the spatial unit to evaluate the wildlife population status. The Chiba prefectural government divided 13 municipalities into 71 units for wildlife management. Therefore, we adopted these management units for the selection of study sites. Site elevations ranged between 10 and 950 m a.s.l. in Hyogo and between 40 and 320 m a.s.l. in Chiba. The mean annual temperatures during the last 10 years (2009–2018) at Toyooka (Hyogo; 3.4 m a.s.l.), Himeji (Hyogo; 38.2 m a.s.l.) and Katsuura (Chiba; 11.9 m a.s.l.) meteorological stations (Fig. 1) were 14.8, 15.6 and 16.2°C, respectively, and the mean annual precipitation estimates were 2167, 1459 and 2068 mm, respectively (JMA 2019). The mean annual maximum snow depth is 492 mm at the Toyooka meteorological station, and the other two stations do not measure snow-related parameters owing to negligible snowfall (JMA 2019).

About 67% of Hyogo Prefecture is forested area (MAFF 2019) mainly composed of secondary deciduous broadleaved forests dominated by Konara oak Quercus serrata and Chinese cork oak Q. variabilis, secondary evergreen coniferous forests dominated by Japanese red pine Pinus densiflora, and artificial evergreen coniferous forests of Japanese cedar Cryptomeria japonica or Hinoki cypress Chamaecyparis obtusa. Only 30% of Chiba Prefecture is covered with forest (MAFF 2019). These areas are composed mainly of secondary deciduous broad-leaved forests dominated by Konara oak and evergreen coniferous forests of Japanese cedar, Hinoki cypress and Japanese red pine. Also, secondary evergreen broad-leaved forests dominated by Itajii Castanopsis sieboldii are seen in the southeastern area. The forest conditions at each site are shown in the Supporting information.

Between 1960 and 2000, about 4000–6000 wild boars year–1 were harvested or culled in Hyogo Prefecture; recently, these hunting bags have increased, reaching 18 000 wild boars year–1 in 2010 (Yokoyama 2014). Agricultural damage caused by wild boars has become a serious issue throughout Hyogo Prefecture (Kuriyama et al. 2018). Wild boars in Chiba Prefecture were believed to be extinct in the mid-1970s. The current population likely originated from artificially introduced boars in the mid-1980s with subsequent expansion, which caused agricultural damage in the southern part of the prefecture (Asada 2011). The annual counts of harvested or culled wild boar were less than 1000 until 2000 in Chiba Prefecture, but these hunting bags have recently increased, reaching 22 000 wild boars in 2015 (Chiba Prefecture 2017).

Figure 1.

Study sites (black circles) for signs of activity on transects and camera trap surveys conducted in Hyogo and Chiba Prefectures with locations of meteorological stations (white squares). Forested area (dark gray shaded area) is drawn based on data provided by the Geographical Survey Institute of Japan.

img-z3-1_wlb.00869.jpg

Selection of the study period

The farrowing season of wild boar typically lasts from April to July in Japan (Tsuji et al. 2013). Furthermore, wild boar have larger litter sizes compared with other similarly sized ungulates (Carranza 1996) and adult females have high fertility rates (Bieber and Ruf 2005, Tsuji et al. 2013). Therefore, the wild boar population is expected to increase steadily and drastically during this period.

Wild boar populations decrease through harvesting or culling. The hunting season is from 15 November to 15 March in Hyogo and from 15 November to 15 February in Chiba. However, in recent years, culling has been carried out year-round to protect crops and has accounted for the majority of captures (MOE 2019). It is therefore difficult to select a study period that avoids the risk of human-derived mortality in Japan.

The rate of disappearance of fecal pellet groups, often used as a relative abundance index for ungulates, changes seasonally due to changes in temperature, rain and dung beetles (Massei et al. 1998). Although data for wild boar are lacking, pellet groups of sika deer disappear rapidly in the summer (during June–August) under the Japanese climate (Koike et al. 2013). We selected autumn to winter as the study period, when variation in the wild boar density is expected to be low and fecal pellets are likely to remain over long periods. We conducted a survey for signs of activity along a transect once at each site from October to December 2017 and a camera trap survey for one or two months, including this date (Supporting information).

Signs of activity along the transect

Prior to the field survey, we established one transect of >5 km within each site using Quantum GIS (ver. 2.18.19, < www.qgis.org/>). Transects were set along the forested mountain ridge to minimize variation in topographical conditions among transects (sites) and to ensure that it was possible to walk along the route, avoiding steep terrain. Transects were separated by at least 3 km. Based on the reported average monthly home range size of wild boar of 50–400 ha (Caley 1997, Massei et al. 1997, Kay et al. 2017, Koh et al. 2018), individuals were not likely to traverse multiple transects during the study period.

We trekked once along the transect at each site, counted signs of wild boar activity (digging marks, rubbing marks and fecal pellet groups) within 1 m on both sides of the transect (total width, 2 m), and recorded counts in sections of about 100 m. Digging (rooting) marks were pits formed on the ground by the foraging behavior of wild boar, and a pit with distinct boundaries (mounds) was counted as an independent mark (Supporting information). Signs formed only on the surface of the ground, such as disturbances of fallen leaves, may also be derived from wild boar but were excluded owing to uncertainty. When it was difficult to distinguish the independence of digging marks, they were recorded as a single large-scale mark. However, because large-scale digging marks accounted for only about 3.1% of the total digging marks and were likely derived from the foraging behavior of the same individual (or herd), we used both marks for analyses without distinction. Rubbing marks formed by wild boars scratching their bodies were identified by mud attached to the trunk or by the wear of bark (Supporting information). We counted the trees with these features. We defined one fecal pellet group as a pile of feces believed to be produced from one defecation from a single wild boar. To determine whether adjacent groups were independent, the level of freshness (drying or weathering) and size of pellets were evaluated. These signs of wild boar activity are relatively easy to distinguish from those of other mammals.

Camera trap survey

We placed infrared-triggered cameras at 30 locations during a one-month period at six sites (H1–H6) in Hyogo Prefecture as well as at 20 locations for about one month at one site (C1) or 10 locations for about two months (with a location change after one month) at seven sites (C2–C8) in Chiba Prefecture. We estimated the densities of wild boar from camera data by the REST model (Nakashima et al. 2018). The estimated densities by the REST model are relatively unbiased even when survey effort is low; however, the precision is improved by increasing survey effort (the number of locations and survey duration). However, increasing the survey period by 20 days or more has little effect on precision (Nakashima et al. 2018). Therefore, we set the survey period to one month at each location. Each camera was randomly located within an area of 20 m from the transect route or the nearby forested ridgeline to unify the topographical conditions of the transect survey.

Three infrared-triggered camera models were used, the Bushnell Trophy Cam HD or Bushnell Trophy Cam Aggressor No Glow in Hyogo Prefecture and Browning Strike Force HD Pro in Chiba Prefecture. Cameras were set to video mode (HD: 30 s, AG: 10 s, SF: 20 s) with the minimum wait time (≤1 s) to prevent missed shooting. We defined the equilateral–triangular focal areas for each of the three camera models considering its field of view and detection area and analyzed only the records obtained in the areas (Supporting information). Although the focal areas differed among camera models, it is clear from the definition of the REST model described later that there is no effect on the expected values for densities. To clearly record the stay of wild boar in the focal area, cameras were mounted on the tree trunk about 20–30 cm above the ground and marking stakes were installed as landmarks at each apex of the triangular focal area. Then, we adjusted the direction and angle of the camera in reference to the position of the marking stakes while checking the field of view on a monitor (built into SF) or with a tablet device connected via USB (HD and AG). After the adjustment, we removed the marking stakes to avoid effects on the behavior of wild boars.

We played all videos captured by cameras and recorded the number of wild boars entering the focal area, the length of staying time and the research period for each camera. The research period was the total number of effective camera trap days, excluding the days of our field survey and periods when the camera did not work properly due to battery depletion, accidental animal attacks on the cameras and so on. When the same individual stayed in the focal area continuously across multiple videos, we collectively regarded the stays as one observation record.

Density and abundance estimation from camera trap data

Among various methods, we used the REST model (Nakashima et al. 2018) and Poisson–Poisson N-mixture model (PP model; Kéry and Royle 2015) to estimate densities or abundances of wild boar at each site without individual identification from camera trap data.

The REST model requires some assumptions, including the certain detection of animals entering the focal area (Nakashima et al. 2018); however, it is an efficient and realistic method for estimating the density of ground-dwelling mammals lacking individually recognizable markings (Palencia et al. 2021). In the model, the relationship between animal density D and camera trap data is described by the following equation:

e01_wlb.00869.gif

where s is the focal area of the camera trap defined prior to the survey, H is the research period, E(Y) is the expected number of animal encounters and E(T) is the expected staying time per encounter. Nakashima et al. (2018) pointed out the necessity of estimating the proportion of time that animals are active (a) to correct the total research period H because inactive animals are inevitably less detectable by cameras. Thus, we estimated a of the wild boars in each prefecture using the activity package in R described by Rowcliffe et al. (2014).

Prior to population density estimation, we chose probability distributions of observed staying time and the number of observed animal encounters using model selection criteria. The observed staying time (T) for each encounter was regarded as a non-negative, continuous random variable following an arbitrary temporal probability distribution, such as the exponential, gamma, Weibull and log-normal distributions. Similarly, the number of observed animal encounters (Y) during the sampling period (H) can be regarded as a non-negative, discrete random variable following an arbitrary probability distribution, such as Poisson, negative binomial, zero-inflated Poisson and zero-inflated negative binomial distributions. We estimated parameters for probability distributions of T and Y for each site using Bayesian methods. However, we estimated T for each site by treating inter-site variation as a random effect with mean zero and variance σ2 because it was not possible to estimate the parameter independently for some sites with low wild boar densities and few observed passes. The posterior distribution was estimated using the Markov chain Monte Carlo (MCMC) method by simulating posterior samples of parameters and their variances. A gamma distribution with shape (0.001) and rate (0.001) parameters was used as the prior distribution for staying time T, dispersion parameter for the number of animal encounters Y and boar density D. We selected the gamma distribution for T (Supporting information) and zero-inflated negative binomial distribution for Y (Supporting information) based on the WAIC (Watanabe–Akaike information criterion, Watanabe 2010).

The N-mixture model is widely used to estimate animal abundance (when the home range includes the survey location) based on spatially and temporally replicated counts (Kéry and Royle 2015). Its derivative the PP model can effectively explain the structure of camera trap data because it can accommodate multiple counts for the same individual within a single sampling occasion by assuming that the number of detections per individual follows a Poisson distribution (Kéry and Royle 2015, Nakashima 2020). The number of animals for each site i in location j, Nij is assumed to follow a Poisson distribution:

e02_wlb.00869.gif

where λi is the expected number of animals per sampling location at site i. The observation process is described by the following equation:

e03_wlb.00869.gif

where yijt is the number of individuals detected at site i location j during t sampling events (day), and µi is the expected number of times in which an individual is detected at each location at site i during a sampling event. We estimated the parameters λi and µi for each site using MCMC simulations. Uniform distributions from 0 to 100 and from 0 to 1 were used as the prior values for λi and µi, respectively.

Relationships between signs of activity and wild boar densities

We modeled the relations between the observed number of activity signs ASij (digging marks, rubbing marks and fecal pellet groups) found on transect i section j and wild boar densities or abundances Di at site i in a hierarchical framework as follows:

e04_wlb.00869.gif

where pij is the probability parameter and r is the size parameter for the negative binomial distribution. We treated Di as a latent variable associated with the REST and PP models. In the model, the length of transect i section j, TLij (i.e. observation effort), is used as an offset term. We estimated the coefficients α and β using MCMC simulations. The prior distributions of α and β were uniform distributions from –100 to 100 and that of r was the uniform distribution from 0 to 100.

MCMC simulations were performed using JAGS (ver. 4.3.0, < www.mcmc-jags.sourceforge.net>, accessed 11 October 2017) called from R (ver. 3.3.2, < www.r-project.org>, accessed 3 November 2016) with the runjags package (Denwood 2016). We ran three MCMC chains in parallel for 20 000 iterations, following a burn-in period of 20 000 iterations for each chain, with thinning to every 20th sample. The medians, 50% and 95% credible limits, and variances of wild boar density, abundance and coefficient estimates were calculated as the posterior summary. Model convergence was examined using the potential scale reduction factor (PSRF; Gelman and Rubin 1992) and by verifying no divergent transitions. Convergence is approximately reached if PSRF is close to one.

Then, we evaluated the model by comparing the median of the posterior distribution of signs of activity (as expected values) and observed values. We used the following three indices (Moriasi et al. 2007) for model evaluation: 1) coefficient of determination, which describes the degree of collinearity between observed and expected values; 2) Nash–Sutcliffe efficiency (NSE), which describes the relative magnitude of the residual variance compared with the explained variance; and 3) percent bias (PBIAS), which estimates the model's tendency to overpredict (PBIAS > 0) or underpredict values (PBIAS < 0).

Results

We found, on average, 29.3 (range, 1.2–66.5) digging marks km–1, 9.6 (0.0–51.0) rubbing marks km–1 and 0.3 (0.0–1.2) fecal pellet groups km–1 along the transects (Table 1). We obtained 523 and 2388 video files of wild boar and 225 and 606 independent staying records within the focal area during 6042 and 4985 effective camera trap days for Hyogo and Chiba Prefectures, respectively. Wild boars showed crepuscular and nocturnal activity patterns in both prefectures (Supporting information), and the proportions of time spent active were 0.35 (95% CI, 0.27–0.44) and 0.46 (0.38–0.52) in Hyogo and Chiba, respectively. The medians of the posterior distributions of boar density determined by the REST model varied from 0.7 to 114.3 wild boars km–2 among sites and the abundance determined by the PP model varied from 0.02 to 0.70 among sites (Table 2). All MCMC simulations converged (PSRF close to one).

Table 1.

Counts of signs of wild boar activity and the length of transects at each site.

img-z6-2_wlb.00869.gif

For the relations between signs of activity and wild boar densities, the median values (95% CI) of parameters α and β estimated by MCMC simulations were 1.70 (0.77–2.37) and 0.59 (0.37–0.89) for digging marks (Fig. 2a) and –1.04 (–2.53 to 0.47) and –0.11 (–0.61 to 0.52) for fecal pellet groups (Fig. 2c). However, the MCMC sampling distribution for the relationship between rubbing marks and wild boar densities did not converge well. The expected values of digging marks from MCMC simulations were positively correlated with the observed values (R2 = 0.82, p < 0.001), and NSE was high (0.82) and PBIAS was close to zero (2.3) (Table 3). The expected values of fecal pellet groups from MCMC simulations were not correlated with the observed values (Table 3).

Similarly, α and β estimated from the relationship with the abundance determined by the PP model were 4.77 (4.23–5.48) and 0.97 (0.62–1.42) for digging marks (Fig. 3a), 5.02 (4.26–6.09) and 2.09 (1.56–2.85) for rubbing marks (Fig. 3b), and –0.35 (–2.15 to 2.14) and 0.60 (–0.47 to 1.99) for fecal pellet groups (Fig. 3c). The expected values of digging marks from MCMC simulations were positively correlated with the observed values (R2=0.79, p < 0.001), and NSE was high (0.78) and PBIAS was close to zero (3.5) (Table 4). The expected values of rubbing marks and fecal pellet groups from MCMC simulations were weakly or not correlated with the observed values (Table 4).

Table 2.

Estimated wild boar densities obtained from the REST model, abundances obtained from PP model and camera-trapping details for each study site.

img-z6-10_wlb.00869.gif

Discussion

Of the three signs of wild boar activity, we found that only digging marks were strongly associated with density estimates obtained using the REST model (Fig. 2a) and with abundance estimates obtained by the PP model (Fig. 3a). Furthermore, the model showed very good performance based on assessment indices (Moriasi et al. 2007). Based on these results, the number of digging marks reflected the wild boar density and abundance in our study area and is reliable, to some extent, as a relative abundance index to examine spatial patterns, consistent with a previous report by Anderson and Stone (1994). However, the median β value was 0.59 in the regression against the wild boar density by the REST model (that is, the number of digging marks is saturated at high densities), whereas it was 0.97 in the regression against the abundance by the PP model (that is, close to linear). Therefore, it should be noted that the relation with wild boar density may vary. Our results conflict with those of previous studies suggesting that there is no relationship between digging marks and wild boar density and digging marks are therefore not a suitable predication index (Hone et al. 1988, Massei et al. 2018). These conflicting results may be related to a difference in wild boar densities, which were 0.89 wild boars km–2 (Hone et al. 1988) and 0.71–6.99 wild boars km–2 (Massei et al. 2018) in previous studies, compared with 0.7–114.3 wild boars km–2 at our study sites. According to Welander (2000), digging is affected by various factors, including year, season, habitat type, soil category and availability of natural food resources above and below ground. Accordingly, factors other than wild boar density can result in substantial variation in the number of digging marks. Therefore, as in this study, if densities differ on the order of dozens of wild boars km–2 or more, it may be possible to roughly evaluate abundance using the number of digging marks; however, this may not be applicable when the variation in density is low. Additionally, given that the rooting frequency is related to the production of hard masts (Bruinderink and Hazebroek 1996), it should be noted that the relationship between the wild boar density and the number of digging marks is affected by annual fluctuations and spatial synchronization.

Figure 2.

Relationship between wild boar densities (50% credible intervals, gray dashes; 95%, black bars) estimated by the REST model from camera trap data and counts per kilometer for digging marks (a), rubbing marks (b) and fecal pellet groups (c) at 14 sites in Hyogo (black dots) and Chiba (white dots) Prefectures. Predictions by the hierarchal model are shown as median values (black dashes) and 95% credible intervals (gray shaded areas).

img-z7-1_wlb.00869.jpg

Table 3.

Model evaluation indices for MCMC simulations of the relationships between signs of activity and wild boar densities by the REST model.

img-z7-5_wlb.00869.gif

Rubbing behaviors in wild boars function to scratch the skin and remove mud, hair and ectoparasites (Allwin et al. 2015). Our prediction that the number of rubbing marks reflects the density of wild boars was not supported by our results. There are several potential explanations for this observation. First, rubbing in wild boars is often associated with wallowing behavior (Allwin et al. 2015); however, the distribution of wallows is spatially nonuniform. In addition, other than marks left by wild boars covered in mud, signs of rubbing are difficult to detect. Second, rubbing marks detected as bark wear form more easily for some tree species than others, and once formed, remain for a long time. Finally, although we counted trees with rubbing marks in this study, tree species and tree density may differ among sites. Because it is difficult to consider all of these factors, we conclude that rubbing marks are not a reliable relative abundance index for wild boars.

Fecal counts are a standard and effective indicator of herbivore abundance (e.g. rabbit Oryctolagus cuniculus, Wood 1988; kangaroo Macropodidae, Vernes 1999; red deer Cervus elaphus, Forsyth et al. 2007). They are also used for estimating the relative abundance of wild boars (Acevedo et al. 2007, Plhal et al. 2014) and are related to observed density (Hone 2012). Although the survey period was determined to maximize fecal pellet groups remaining under the Japanese climate, few fecal pellet groups were found, and these counts were not related to wild boar densities. Previous studies have reported averages of 4.5–16.3 fecal pellet groups/200 m2 in the Czech Republic (Plhal et al. 2014) and 0.9 fecal pellets m–2 (18 pellets/200 m2) in Australia (Hone and Martin 1998); we found, on average, only about 0.03 fecal pellet groups/200 m2 for 14 sites. Although differences in boar densities can lead to differences in numbers of fecal pellet groups, it is unlikely that about 100 times more wild boars inhabited these areas than our study area. Hone and Martin (1998) reported that the higher the temperature and precipitation, the faster the disappearance of fecal pellets of wild boars. Therefore, the differences among studies might be explained in part by the higher temperatures and precipitation in our study area than in the previous study areas (Supporting information). Additionally, dense understory vegetation obstructs visibility and prevents the detection of pellets (Zabek et al. 2016). However, the effectiveness of pellet counts as a relative abundance index for sika deer has been confirmed in Japan, including our study area, and the method is used nationwide (Hamasaki et al. 2007). Therefore, the low pellet detection in our study area may be explained by properties specific to wild boar as well as environmental conditions. Wild boars defecate less frequently than other ungulates and have site preferences for defecation (Ebert et al. 2009, Plhal et al. 2014). In addition, various factors, including dung beetles, digging behavior of wild boars, moisture and composition of feces, might affect the rate of feces disappearance (Massei et al. 1998). Kosugi and Sakamoto (1994) reported that dung beetles (Onthophagus ater and O. lenzii) strongly favor feces of wild boar over those of sika deer and Japanese hare Lepus brachyurus in Japan. These characteristics of wild boars and their feces may explain the scarcity of fecal pellet groups on forested mountain ridges, at our study sites, and/or in temperate Japan.

Figure 3.

Relationship between wild boar abundances (50% credible intervals, gray dashes; 95%, black bars) estimated by the PP model from camera trap data and counts per kilometer for digging marks (a), rubbing marks (b) and fecal pellet groups (c) at 14 sites in Hyogo (black dots) and Chiba (white dots) Prefectures. Predictions by the hierarchal model are shown as median values (black dashes) and 95% credible intervals (gray shaded areas).

img-z8-1_wlb.00869.jpg

Populations of wild boars and feral pigs are increasing worldwide and have significant impacts on the environment and economy. In Japan, human–wild boar conflicts, including the appearance of wild boar in urban landscapes and the spread of classical swine fever, are increasing, and prefecture-level administrations require information about population trends to address these issues. However, large-scale surveys for estimating absolute density are generally very expensive. In such cases, signs of activity are alternative indexes to understand spatial abundance. Although fecal counts are the most frequently used relative abundance index for wild boar, our results indicate that it may not be an effective parameter in Japan. Although temporal (year and/or season) fluctuations in the relationship between population density and activity signs should be noted, our results support the effectiveness of digging marks as a relative abundance index of wild boar. Understanding differences in spatial density by extensive surveys using this relative abundance index will guide the proper allocation of effort and resources for effective population management.

Table 4.

Model evaluation indices for MCMC simulations of the relationships between signs of activity and wild boar abundances by the PP model.

img-z8-6_wlb.00869.gif

Acknowledgements –

We are grateful to the students of the Laboratory of Forest Animal Studies, Nihon University and the Laboratory of Biodiversity Science, the University of Tokyo, especially Mr. M. H. Tsutsui and Mr. Y. Yokoyama, for assisting with the camera trap survey in Chiba Prefecture. We also thank the staff of the Wildlife Management Research Center, Hyogo and Hyogo Midori Public Corporation for assisting with the survey in Hyogo Prefecture. We thank Ms. Y. Fukaya of the University of Tokyo for assisting with data processing and Dr. T. Miyashita of the University of Tokyo and Dr. M. Asada for constructive advice.

Funding – This research was supported by the Environment Research and Technology Development Fund (JPMEERF20174004 and JPMEERF20204G01) of the Environmental Restoration and Conservation Agency.

Conflicts of interest – There are no conflicts of interests to declare.

Author contributions

All authors contributed to the sample collection, data analysis or the writing and review of the manuscript.

Data availability statement

Data are available from the Dryad Digital Repository: < http://dx.doi.org/10.5061/dryad.5hqbzkh6h> (Higashide et al. 2021).

References

1.

Abaigar, T. et al. 1994. Habitat preference of wild boar (Sus scrofa L., 1758) in a Mediterranean environment. Indirect evaluation by signs. – Mammalia 58: 201–210. Google Scholar

2.

Acevedo, P. et al. 2007. Estimation of European wild boar relative abundance and aggregation: a novel method in epidemiological risk assessment. – Epidemiol. Infect. 135: 519–527. Google Scholar

3.

Allwin, B. et al. 2015. Wild pig field signs-identity and confirmation. – Int. J. Life Sci. Res. 3: 205–213. Google Scholar

4.

Anderson, S. J. and Stone, C. P. 1994. Indexing sizes of feral pig populations in a variety of Hawaiian natural areas. – Trans. West. Sect. Wildl. Soc. 30: 26–39. Google Scholar

5.

Asada, M. 2011. Distribution and pest control, damage to agricultural production for wild boar in 2009 in Chiba prefecture, Japan. – Chiba Biodivers, Center 3: 49–64, in Japanese. Google Scholar

6.

Barrios-Garcia, M. N. and Ballari, S. A. 2012. Impact of wild boar Sus scrofa in its introduced and native range: a review. – Biol. Invas. 14: 2283–2300. Google Scholar

7.

Bieber, C. and Ruf, T. 2005. Population dynamics in wild boar Sus scrofa: ecology, elasticity of growth rate and implications for the management of pulsed resource consumers. – J. Appl. Ecol. 42: 1203–1213. Google Scholar

8.

Bruinderink, G. G. and Hazebroek, E. 1996. Wild boar (Sus scrofa scrofa L.) rooting and forest regeneration on podzolic soils in the Netherlands. – For. Ecol. Manage. 88: 71–80. Google Scholar

9.

Caley, P. 1997. Movements, activity patterns and habitat use of feral pigs Sus scrofa in a tropical habitat. – Wildl. Res. 24: 77–87. Google Scholar

10.

Carranza, J. 1996. Sexual selection for male body mass and the evolution of litter size in mammals. – Am. Nat. 148: 81–100. Google Scholar

11.

Chandler, R. B. and Royle, J. A. 2013. Spatially explicit models for inference about density in unmarked or partially marked populations. – Ann. Appl. Stat. 7: 936–954. Google Scholar

12.

Chiba Prefecture 2017. Conservation and management plan of specified mammal species in Chiba prefecture: wild boar. – Chiba Prefecture, Japan, in Japanese. Google Scholar

13.

Denwood, M. J. 2016. runjags: an R package providing interface utilities, model templates, parallel computing methods and additional distributions for MCMC models in JAGS. – J. Stat. Softw. 71: 1–25. Google Scholar

14.

Ebert, C. et al. 2009. Is non-invasive genetic population estimation via faeces sampling feasible for abundant mammals with low defecation rates? A pilot study on free ranging wild boar Sus scrofa in South-West Germany. – Acta Silvat. Lign. Hung. 5: 167–177. Google Scholar

15.

Ebert, C. et al. 2012. Estimating wild boar Sus scrofa population size using faecal DNA and capture–recapture modelling. – Wildl. Biol. 18: 142–152. Google Scholar

16.

Enetwild Consortium et al. 2018. Guidance on estimation of wild boar population abundance and density: methods, challenges, possibilities. – EFSA Supp. Publ. 15: 1449E. Google Scholar

17.

Engeman, R. M. et al. 2013. Monitoring wild pig populations: a review of methods. – Environ. Sci. Pollut. Res. 20: 8077–8091. Google Scholar

18.

Ferretti, F. et al. 2016. The use of faeces counts to estimate relative densities of wild boar in a Mediterranean area. – Popul. Ecol. 58: 329–334. Google Scholar

19.

Focardi, S. et al. 2020. Reliable estimates of wild boar populations by nocturnal distance sampling. – Wildl. Biol. 2020: wlb.00694. Google Scholar

20.

Forsyth, D. M. et al. 2007. Modeling the relationship between fecal pellet indices and deer density. – J. Wildl. Manage. 71: 964–970. Google Scholar

21.

Franzetti, B. et al. 2012. Nocturnal line transect sampling of wild boar Sus scrofa in a Mediterranean forest: long-term comparison with capture–mark–resight population estimates. – Eur. J. Wildl. Res. 58: 385–402. Google Scholar

22.

Gelman, A. and Rubin, D. B. 1992. Inference from iterative simulation using multiple sequences. – Stat. Sci. 7: 457–472. Google Scholar

23.

Gortázar, C. et al. 2007. Diseases shared between wildlife and livestock: a European perspective. – Eur. J. Wildl. Res. 53: 241–256. Google Scholar

24.

Hamasaki, S. et al. 2007. Evaluation of three monitoring methods for sika deer population density: block count, pellet group count and sighting per unit effort. – Mamm. Sci. 47: 65–71, in Japanese with English abstract. Google Scholar

25.

Hebeisen, C. et al. 2008. Estimating wild boar abundance and density using capture–resights in Canton of Geneva, Switzerland. – Eur. J. Wildl. Res. 54: 391–401. Google Scholar

26.

Higashide, D. et al. 2021. Data from: Effectiveness of signs of activity as relative abundance indices for wild boar. – Dryad Digital Repository, <  http://dx.doi.org/10.5061/dryad.5hqbzkh6h >. Google Scholar

27.

Hone, J. 1988. Evaluation of methods for ground survey of feral pigs and their sign. – Acta Theriol. 33: 451–465. Google Scholar

28.

Hone, J. 2012. Applied population and community ecology: the case of feral pigs in Australia. – Wiley. Google Scholar

29.

Hone, J. and Martin, W. 1998. A study of dung decay and plot size for surveying feral pigs using dung counts. – Wildl. Res. 25: 255–260. Google Scholar

30.

Howe, E. J. et al. 2017. Distance sampling with camera traps. – Methods Ecol. Evol. 8: 1558–1565. Google Scholar

31.

Iijima, H. et al. 2013. Estimation of deer population dynamics using a Bayesian state-space model with multiple abundance indices. – J. Wildl. Manage. 77: 1038–1047. Google Scholar

32.

Isoda, N. et al. 2020. Dynamics of classical swine fever spread in wild boar in 2018–2019, Japan. – Pathogens 9: 119. Google Scholar

33.

JMA, Japan Meteorological Agency. 2019. Past weather data search. – <  www.data.jma.go.jp/obd/stats/etrn/index.php >, accessed 2 Apr 2019, in Japanese. Google Scholar

34.

Kay, S. L. et al. 2017. Quantifying drivers of wild pig movement across multiple spatial and temporal scales. – Mov. Ecol. 5: 14. Google Scholar

35.

Kéry, M. and Royle, J. A. 2015. Applied hierarchical modeling in ecology: analysis of distribution, abundance and species richness in R and BUGS: Volume 1: Prelude and static models. – Academic Press. Google Scholar

36.

Koh, J. J. M. et al. 2018. Using a spatial mark–resight model to estimate the parameters of a wild pig Sus scrofa population in Singapore. – Raffles B. Zool. 66: 494–505. Google Scholar

37.

Koike, S. et al. 2013. Seasonal changes and altitudinal variation in deer fecal pellet decay. – Eur. J. Wildl. Res. 59: 765–768. Google Scholar

38.

Kosugi, J. and Sakamoto, Y. 1994. Selectivity of dungs by dung beetles. – Jpn. Wildl. Res. Soc. 20: 13–18, in Japanese with English abstract. Google Scholar

39.

Kuriyama, T. et al. 2018. Distribution, density and agricultural damage for wildlife populations in Hyogo Prefecture, Japan. – Wildl. Monogr. Hyogo 10: 9–31, in Japanese with English abstract. Google Scholar

40.

MAFF, Ministry of Agriculture, Forestry and Fisheries Government of Japan 2018. FY2017 Summary of the annual report on food, agriculture and rural areas in Japan. – Ministry of Agriculture, Forestry and Fisheries, Tokyo, Japan. Google Scholar

41.

MAFF, Ministry of Agriculture, Forestry and Fisheries Government of Japan 2019. Ratio of forests in each prefecture in 2017. – <  www.rinya.maff.go.jp/j/keikaku/genkyou/h29/1.html >, accessed 2 Jul 2019, in Japanese. Google Scholar

42.

Massei, G. and Genov, P. V. 2004. The environmental impact of wild boar. – Galemys 16: 135–145. Google Scholar

43.

Massei, G. et al. 1997. Factors influencing home range and activity of wild boar Sus scrofa in a Mediterranean coastal area. – J. Zool. 242: 411–423. Google Scholar

44.

Massei, G. et al. 1998. Fallow deer and wild boar pellet group disappearance in a Mediterranean area. – J. Wildl. Manage. 62: 1086–1094. Google Scholar

45.

Massei, G. et al. 2015. Wild boar populations up, numbers of hunters down? A review of trends and implications for Europe. – Pest Manage. Sci. 71: 492–500. Google Scholar

46.

Massei, G. et al. 2018. Camera traps and activity signs to estimate wild boar density and derive abundance indices. – Pest Manage. Sci. 74: 853–860. Google Scholar

47.

McMahon, C. R. et al. 2010. Spatially explicit spreadsheet modelling for optimising the efficiency of reducing invasive animal density. – Methods Ecol. Evol. 1: 53–68. Google Scholar

48.

MOE and MAFF, Ministry of the Environment and Ministry of Agriculture, Forestry and Fisheries Government of Japan 2013. Radical strengthening measures on culling of birds and mammals. – Ministry of the Environment and Ministry of Agriculture, Forestry and Fisheries, Tokyo, Japan, in Japanese. Google Scholar

49.

MOE, Ministry of the Environment Government of Japan 2015. The wildlife in Japan. – Ministry of the Environment, Tokyo, Japan. Google Scholar

50.

MOE, Ministry of the Environment Government of Japan 2019. Wildlife statistics in 2016. – <  www.env.go.jp/nature/choju/docs/docs2.html >, accessed 2 Jul 2019, in Japanese. Google Scholar

51.

Moriasi, D. N. et al. 2007. Model evaluation guidelines for systematic quantification of accuracy in watershed simulations. – T. ASABE 50: 885–900. Google Scholar

52.

Nakashima, Y. 2020. Potentiality and limitations of N-mixture and Royle–Nichols models to estimate animal abundance based on noninstantaneous point surveys. – Popul. Ecol. 62: 151–157. Google Scholar

53.

Nakashima, Y. et al. 2018. Estimating animal density without individual recognition using information derivable exclusively from camera traps. – J. Appl. Ecol. 55: 735–744. Google Scholar

54.

Nichols, J. D. and Williams, B. K. 2006. Monitoring for conservation. – Trends Ecol. Evol. 21: 668–673. Google Scholar

55.

Nimmo, D. G. et al. 2015. Dingoes can help conserve wildlife and our methods can tell. – J. Appl. Ecol. 52: 281–285. Google Scholar

56.

Oliver, W. and Leus, K. 2008. Sus scrofa in IUCN Red List of threatened species ver. 2019–1. – <  www.iucnredlist.org >, accessed 10 Jun 2019. Google Scholar

57.

Osada, Y. et al. 2015. Exploring the drivers of wildlife population dynamics from insufficient data by Bayesian model averaging. – Popul. Ecol. 57: 485–493. Google Scholar

58.

Palencia, P. et al. 2021. Assessing the camera trap methodologies used to estimate density of unmarked populations. – J. Appl. Ecol. 58: 1583–1592. Google Scholar

59.

Pepin, K. M. et al. 2017. Efficiency of different spatial and temporal strategies for reducing vertebrate pest populations – Ecol. Model. 365: 106–118. Google Scholar

60.

Plhal, R. et al. 2014. An assessment of the applicability of dung count to estimate the wild boar population density in a forest environment. – J. For. Sci. 60: 174–180. Google Scholar

61.

Pollock, K. H. et al. 2002. Large scale wildlife monitoring studies statistical methods for design and analysis. – Environmetrics 13: 105–119. Google Scholar

62.

Postel, A. et al. 2019. Reemergence of classical swine fever, Japan, 2018. – Emerg. Infect. Dis. 25: 1228–1231. Google Scholar

63.

Rowcliffe, J. M. et al. 2008. Estimating animal density using camera traps without the need for individual recognition. – J. Appl. Ecol. 45: 1228–1236. Google Scholar

64.

Rowcliffe, J. M. et al. 2014. Quantifying levels of animal activity using camera trap data. – Methods Ecol. Evol. 5: 1170–1179. Google Scholar

65.

Schley, L. et al. 2008. Patterns of crop damage by wild boar Sus scrofa in Luxembourg over a 10-year period. – Eur. J. Wildl. Res. 54: 589–599. Google Scholar

66.

Schwarz, C. and Seber, G. 1999. Estimating animal abundance: review III. – Stat. Sci. 14: 427–456. Google Scholar

67.

Snow, N. P. et al. 2017. Interpreting and predicting the spread of invasive wild pigs. – J. Appl. Ecol. 54: 2022–2032. Google Scholar

68.

Stephens, P. A. et al. 2015. Management by proxy? The use of indices in applied ecology. – J. Appl. Ecol. 52: 1–6. Google Scholar

69.

Tsuji, T. et al. 2013. Estimation of the fertility rates of Japanese wild boars (Sus scrofa leucomystax) using fetuses and corpora albicans. – Acta Theriol. 58: 315–323. Google Scholar

70.

Vernes, K. 1999. Pellet counts to estimate density of a rainforest kangaroo. – Wildl. Soc. Bull. 27: 991–996. Google Scholar

71.

Vicente, J. et al. 2004. Epidemiological study on porcine circovirus type 2 (PCV2) infection in the European wild boar Sus scrofa. – Vet. Res. 35: 243–253. Google Scholar

72.

Watanabe, S. 2010. Asymptotic equivalence of Bayes cross validation and widely applicable information criterion in singular learning theory. – J. Mach. Learn. Res. 11: 3571–3594. Google Scholar

73.

Welander, J. 2000. Spatial and temporal dynamics of wild boar Sus scrofa rooting in a mosaic landscape. – J. Zool. 252: 263–271. Google Scholar

74.

Wood, D. H. 1988. Estimating rabbit density by counting dung pellets. – Wildl. Res. 15: 665–671. Google Scholar

75.

Yokoyama, M. 2014. Current status and issues of wild boar management in Hyogo Prefecture. – Wildl. Monogr. Hyogo 6: 1–8, in Japanese. Google Scholar

76.

Zabek, M. A. et al. 2016. Estimating distribution and abundance of feral horses Equus caballus in a coniferous plantation in Australia, using line-transect surveys of dung. – Wildl. Res. 43: 604–614. Google Scholar

Appendices

Supplementary information (available online as Appendix wlb-00869 at < www.wildlifebiology.org/appendix/wlb-00869>).

© 2021 The Authors.
Daishi Higashide, Takeo Kuriyama, Shun Takagi, Yoshihiro Nakashima, Keita Fukasawa, Gota Yajima, Minoru Kasada, and Mayumi Yokoyama "Effectiveness of signs of activity as relative abundance indices for wild boar," Wildlife Biology 2021(4), wlb.00869-, (15 October 2021). https://doi.org/10.2981/wlb.00869
Accepted: 15 August 2021; Published: 15 October 2021
JOURNAL ARTICLE
PAGES


SHARE
ARTICLE IMPACT
Back to Top