<figure><img src="image_1.png" /></figure>

HAL
open science

# Acoustic localization at large scales: a promising method for grey wolf monitoring

Morgane Papin, Julian Pichenot, François Guérold, Estelle Germain

► To cite this version:

Morgane Papin, Julian Pichenot, François Guérold, Estelle Germain. Acoustic localization at large scales: a promising method for grey wolf monitoring. Frontiers in Zoology, 2018, 15, pp.11. 10.1186/s12983-018-0260-2 . hal-03163531

HAL Id: hal-03163531

https://hal.univ-lorraine.fr/hal-03163531

Submitted on 10 Jul 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

# Acoustic localization at large scales: a promising method for grey wolf monitoring

Morgane Papin<fnref n="1" />,² ID, Julian Pichenot³, François Guérold² and Estelle Germain<fnref n="1" />*

## Abstract

**Background:** The grey wolf (*Canis lupus*) is naturally recolonizing its former habitats in Europe where it was extirpated during the previous two centuries. The management of this protected species is often controversial and its monitoring is a challenge for conservation purposes. However, this elusive carnivore can disperse over long distances in various natural contexts, making its monitoring difficult. Moreover, methods used for collecting signs of presence are usually time-consuming and/or costly. Currently, new acoustic recording tools are contributing to the development of passive acoustic methods as alternative approaches for detecting, monitoring, or identifying species that produce sounds in nature, such as the grey wolf. In the present study, we conducted field experiments to investigate the possibility of using a low-density microphone array to localize wolves at a large scale in two contrasting natural environments in north-eastern France. For scientific and social reasons, the experiments were based on a synthetic sound with similar acoustic properties to howls. This sound was broadcast at several sites. Then, localization estimates and the accuracy were calculated. Finally, linear mixed-effects models were used to identify the factors that influenced the localization accuracy.

**Results:** Among 354 nocturnal broadcasts in total, 269 were recorded by at least one autonomous recorder, thereby demonstrating the potential of this tool. Besides, 59 broadcasts were recorded by at least four microphones and used for acoustic localization. The broadcast sites were localized with an overall mean accuracy of 315 ± 617 (standard deviation) m. After setting a threshold for the temporal error value associated with the estimated coordinates, some unreliable values were excluded and the mean accuracy decreased to 167 ± 308 m. The number of broadcasts recorded was higher in the lowland environment, but the localization accuracy was similar in both environments, although it varied significantly among different nights in each study area.

**Conclusions:** Our results confirm the potential of using acoustic methods to localize wolves with high accuracy, in different natural environments and at large spatial scales. Passive acoustic methods are suitable for monitoring the dynamics of grey wolf recolonization and so, will contribute to enhance conservation and management plans.

**Keywords:** acoustic monitoring, autonomous recorders, *Canis lupus*, field research, localization estimation, microphone array, wolf howl

## Background

Passive acoustic monitoring is being used increasingly to study species that produce sounds in their natural environments (e.g. vocalizations and stridulations) [1]. The current protocols based on passive acoustics methods allow the study of elusive and/or nocturnal species that live in harsh environments (e.g. dangerous access, thick vegetation or limited visibility) [2–4]. These protocols

are focused on species detection [5], density estimation [6, 7], territory use [8], and localization [9, 10]. They are not technically limited to a time period, non-invasive and so, avoid interference with animal behavior in contrast to other monitoring methods (e.g. direct capture or the intrusive presence of observers in the field) [2, 11]. Passive acoustics may also help to reduce the time and human resources required in the field [12, 13]. These main features of passive acoustics suggest that this interesting approach could be employed for monitoring elusive species that require conservation or management plans, such as the grey wolf (*Canis lupus*).

<footnote>* Correspondence: estelle.germain@croc-asso.org</footnote>

<footnote n="1">¹Centre de Recherche et d'Observation sur les Carnivores (CROC), 4 rue de la banie, 57590, Lucy, France</footnote>

<footnote>Full list of author information is available at the end of the article</footnote>
During the two last centuries, the grey wolf was extir-
pated in many areas throughout Europe and North Amer-
ica [14]. In Europe, the species is now legally protected by
the Bern Convention (1979) and the Habitats Directive
(1992). As a consequence, wolves have been recolonizing
their former areas in recent decades [14, 15]. However, con-
flicts emerge with humans where their ranges overlap with
human settlement and agriculture mainly due to the preda-
tion on livestock [16, 17]. Thus, understanding and moni-
toring the expansion of the grey wolf's range is important
for preventing or mitigating conflicts as well as for conser-
vation and management purposes. However, the monitor-
ing of wolves is still challenging in the field because it is a
wide-ranging habitat generalist, which lives at low densities
and is often secretive and elusive [18, 19]. Moreover, the
conventional methods used for detecting the presence of
grey wolves and estimating their number and population
dynamics can be very time-consuming and costly.

Studying howls may be a powerful approach for moni-
toring grey wolf populations, especially in the summer
and during the mating season when howls are produced
widely [20–22]. For instance, wolf howls can allow scien-
tific and wildlife managers to identify a pack due to their
acoustic structure [23, 24]. In addition, several studies per-
formed in captivity have shown that wolf has individual
vocal signature [25–28]. Other studies have highlighted
the potential use of bioacoustics for detecting wolves [29]
as well as for counting them [28, 30–33] or detecting
reproduction events [13]. The results of these studies sup-
port the possibility of using acoustics for monitoring
wolves in the wild. However, to our knowledge, very few
studies have employed passive acoustics for monitoring
wolves (e.g. [29]) and none for localizing them.

In the present study, we conducted field experiments to
investigate the possibility of using a low-density micro-
phone array to localize wolves at a large scale in two areas
located in the colonization front of the species in north-
eastern France [34–37]. For scientific and social reasons,
the experiments were based on a synthetic sound with
similar acoustic properties to wolf howls. As these areas
were characterized by two contrasting environmental con-
texts (mid-mountain and lowland), the synthetic sound
was broadcast at several sites defined according to a strati-
fied sampling technique based on topography and land-
use. We calculated localization estimates and the accur-
acy. Finally, we identified the parameters and biases that
influenced the localization accuracy.

**Methods**

**Study areas**

The study was conducted in two different areas located in
the colonization front of grey wolf in north-eastern France
(Fig. 1). The first study area was located in a mid-
mountain environment in the Massif des Vosges (VM),

where the presence of a wolf pack (at least two individ-
uals) has been attested since 2011 [34, 35, 37]. This area is
covered by mainly herbaceous vegetation (22%), shrub
(51%), and coniferous forest (27%), and the altitude ranges
from 518 to 1305 m above sea level (mean: 930 m).

The second study area was located in the Côtes de Meuse (CM) at altitude ranging from 247 to 381 m above sea level (mean: 329 m), where the presence of the grey wolf was observed in 2012 [36]. The area is covered mainly by deciduous forest (90%) and open land with herbaceous vegetation accounts for only 10% of the area.

The grey wolf howls throughout the year but the pe-
riods with the most frequent howling activity are the
breeding season (January to April: [20]) and the months
following the birth of pups (August to October: [21]).
Thus, this study was conducted during August 2015 in
VM and August 2016 in CM. These periods also coinc-
dited with good conditions for access to the study areas.

Sampling methods and microphone arrays

Twenty autonomous recorders were placed on a systematic grid with an area of 30 km² (6 × 5 km; Fig. 2) at a regular spacing of 1 km, conducing to a relatively low recorder density (0.67 recorders per km²) for both study areas. The automatic recording units employed were Wildlife Acoustics Song Meters (model: SM3; Wildlife Acoustics Inc., Concord, MA, USA) with two built-in omnidirectional microphones (SM3-A1, bandpass: 20–20,000 Hz, frequency response: 20–20,000 Hz ± 10 dB). All of the recorders were associated with a global positioning system (GPS) unit (Garmin International Inc., Olathe, KS, USA) to synchronize their clock time automatically with high precision. The recorders collected 40 acoustic information channels in stereo using 16-bit .wav files at a sampling rate of 16,000 Hz. The recorders were programmed to operate from 8:55 PM to 8:54 AM and to generate 59-min files separated by a break of 1 min (ensure time synchronization). The gain was set to 24 dB for each channel.

The recorders were fixed to tree trunks at a height of
2.88 ± 0.49 m (mean ± standard deviation [SD]). Their
locations were measured with a Trimble GPS (model:
Juno 5B EGPS, real-time accuracy: 2–4 m; Trimble
Navigation Limited, Sunnyvale, CA, USA).

Broadcast sites and periods

In each study area, 60 broadcast sites were randomly dis-
tributed by stratified sampling according to the topog-
raphy using the "Topographic Position Index" [38–40] and
land-use using the "Corine Land Cover" code (European
Union – SoeS, Corine Land Cover 2006) with QGIS soft-
ware (version 2.8.1: [41]). The sites located away from
roads were then moved to the closest road to allow access
with a vehicle. The spatial sampling of the broadcast sites
in terms of the distances to the autonomous recorders
<figure><img src="image_4.png" /><figcaption>**Fig. 1** Locations of the two study areas in north-eastern France</figcaption></figure>

was similar in the two study areas (z-test: z = -1.8180, α = 0.05).

The sound was broadcast during three consecutive
nights from 9 PM to 6 AM, where each night was di-
vided into three periods (dusk: 9 PM to 12 AM;
night: 12 AM to 3 AM; dawn: 3 AM to 6 AM). For
each night, a different itinerary was used so that each
broadcast site was visited once during the three dif-
ferent periods. All of the broadcast sites locations
were measured using the Trimble GPS.

<figure><img src="image_5.png" /><figcaption>**Fig. 2** Sampling method employed in each study area. Red circles: autonomous recorders (20). Black stars: dataloggers (10). Yellow triangles: broadcast sites (60)</figcaption></figure>

Synthetic sound and broadcast equipment

As the study of large carnivores is a sensitive subject
[17, 42], we chose to use a synthetic sound with similar
acoustic properties to wolf howls rather than using real
howls. This sound also permitted to exclude the effects
of wolves' individual acoustic characteristics [25–28, 43].
It was created with the Seewave package [44] in R soft-
ware (version 3.1.2). The sound comprised mixed pure
tones of 7 s with fundamental frequencies ranging from
300 to 1000 Hz, which was accompanied by four har-
monics that covered a wide range of the frequencies that
can be found in wolf howls [20, 28, 45].

The sound was broadcast from four directional loud-
speakers (model: MSH 30/BT, bandpass: 90–20,000 Hz,
output: 50 W at 8 Ω; Work Pro CA, Valencia, Spain) con-
nected to a mixing amplifier (model: PA 90/2 USB, fre-
quency response: 80–18,000 Hz ± 3 dB; output: 30 W RMS;
Work Pro CA) and a 12 V battery. The loudspeakers were
attached to a car roof. During each broadcast, a digital
sound level meter was employed to control the intensity
level at 1 m (model: FI 70SD, bandpass: 31.5–8000 Hz, fre-
quency response: 8000 Hz ± 5.6 dB, settings: fast response,
A-weighting; Distrame S.A, Sainte-Savine, France).

Meteorological context

The nights were selected according to the optimal me-
teorological conditions for acoustic experimentation, i.e.
very low wind speed and no rainfall. The wind speed
was measured for 1 min at each broadcast site with an
anemometer (model: WS9500; La Crosse Technology,
Geispolsheim, France) and it was always less than 2 m.s⁻¹. In addition, 10 weather dataloggers (model: DT-174B; Center for Educational Measurement Inc., Makati, Philippines) were installed below 10 recorders to record the air temperature every 2 min (Fig. 2). The temperature data acquired by all of the dataloggers were averaged per night period and per night. They were used subsequently to calculate the speed of sound during the nocturnal broadcast period, which was required for localization estimation.

## Analysis of recordings and localization estimates

The two channels in all of the recordings were analyzed with Raven Pro software to detect the synthetic sound (version 1.5: [46]; Spectrogram view preset: Hann, 1024 samples, 90% overlap). We used the Sound Finder package in R software for localization estimation (see [47]). This free tool has a higher accuracy than other software [47]. To estimate the localization of a sound, the algorithm in this package requires the time of arrival (TOA) of sound to at least four microphones, the temperature (mean temperature in the study area during the night period), and the coordinates of the microphones. Among the two microphones on each recorder that recorded the sound (ideally four different recorders), we chose that with the best signal-to-noise ratio. When the sound was recorded by only three different recorders, the second microphone on the recorder with the best signal-to-noise ratio was used to obtain a total of four microphones. As the signal-to-noise ratio was too low to use cross-correlation or automatic detection algorithms, the TOA were measured manually based on the spectrogram view (Fig. 3). The

TOA measures were repeated three times and then averaged.

Sound Finder was used to estimate the coordinates of the broadcast sites as well as the temporal error values. The temporal error is defined as the root-mean-squared error of the combined discrepancies between the theoretical and observed delays in the TOA for each pair of microphones [47]. It was used to evaluate the reliability of the localization estimates, where perfect localization had a temporal error of 0 ms.

The distance between the estimated localization (coordinates given by Sound Finder) and the actual broadcast site position (coordinates given by the GPS) corresponded to the localization accuracy. It was calculated using the distance matrix tool in QGIS software.

## Statistical analysis

All of the statistical analyses were conducted with R software (version 3.1.2: [48]) and results were considered to be statistically significant when *P* ≤ 0.05. All of the values were reported as the mean ± SD.

Linear mixed-effects models (lmer function in the lme4 package: [49]) were used to identify parameters that influenced the localization accuracy (“locendority”). All combinations of the fixed effects and their interactions were used to construct the models. The four fixed effects (see Fig. 4) comprised the microphones area (“areamic” in m²), the distance between the microphones area centroid and the broadcast site (“dist” in m), the broadcast period (dusk, night, or dawn: “period”), and the broadcast site position compared with the microphones area (in or out:

<figure><img src="image_6.png" /><figcaption>**Fig. 3** Time of arrival (TOA) measures based on a spectrogram obtained using Raven Pro software. Spectrogram view preset: Hann, 1024 samples, 90% overlap, time window length of 21.5 s, frequency range from 0 to 2000 Hz, greyscale color. Red crosses: pointers placed at the start positions detected in the broadcast sound signal</figcaption></figure>
<figure><img src="image_7.png" /><figcaption>**Fig. 4** Representation of fixed effects when the broadcast site was positioned in the microphones area</figcaption></figure>

"inout"). Given our data structure, random effect was built as a nested random effect because the data came from two study areas (VM or CM: "array") on three nights (N1, N2, or N3: "night") and in three broadcast periods.

The best model was selected according to the lowest Akaike's information criterion (AIC) [50]. The significance of fixed effects was tested using one-way analysis of variance (anova function in the MASS package: [51]) and random effect with the restricted likelihood ratio test (exactRLRT function in the RLRsim package: [52]).

# Results

## Broadcasts

In VM, two sites were excluded from the study because they were too dangerous to access during field nights. Thus, the synthetic sound was broadcast 174 times in VM and 180 times in CM with a total of 354 broadcasts. The broadcast sound amplitude remained constant during the three nights with a mean sound intensity level at 1 m of 115.04 ± 3.07 dBA in VM and 116.53 ± 3.59 dBA in CM. These values are close to the natural amplitude of wolf howls ([29], MP *unpublished observations*).

## Effectiveness of the recorders

All of the autonomous recorders were functional so the effectiveness of the experiments was 100%, with nearly

1200 h of acoustic recordings. According to visual and audio inspections of the recordings, 269 broadcasts were recorded by at least one autonomous recorder. In total, 101 broadcasts were recorded by one recorder (56 in VM and 45 in CM), 85 by two (36 in VM and 49 in CM), 55 by three (25 in VM and 30 in CM), 21 by four (three in VM and 18 in CM), and seven by five only in CM.

## TOA measures

The distances separating the broadcast sites and autonomous recorders ranged from 67 to 3595 m in CM and from 144 to 2751 m in VM. For several recordings, measures of the TOA were impossible to achieve because the signal-to-noise ratio was very low or the synthetic sound was only partially recorded and/or conspicuous. Thus, some recordings could not be included in the analysis. Finally, 59 broadcasts (17%), i.e. 14 in VM (8%) and 45 in CM (25%), recorded by at least four microphones were used for acoustic localization.

## Localization estimation

Localization estimates were calculated for the 59 broadcast sites (14 in VM and 45 in CM). The mean localization accuracy was about 315 ± 617 m and the mean temporal error was 685.57 ± 2049.73 ms (N = 59; Table 1). All of the usable broadcast sites in VM were
<table><caption>**Table 1** Localization accuracy and temporal error values estimated with Sound Finder. (a) With all localization estimations. (b) According to the 200 ms error reliability threshold</caption><thead><tr><th rowspan="2"></th><th colspan="3">Localization accuracy (m)</th><th colspan="3">Temporal error (ms)</th><th rowspan="2">N</th></tr><tr><th>Mean ± SD</th><th>Min</th><th>Max</th><th>Mean ± SD</th><th>Min</th><th>Max</th></tr></thead><tbody><tr><td colspan="8">(a)</td></tr><tr><td>VM</td><td>442 ± 597</td><td>28</td><td>1987</td><td>903.37 ± 2321.32</td><td>1.60</td><td>7865.69</td><td>14</td></tr><tr><td>CM</td><td>276 ± 625</td><td>1</td><td>2983</td><td>617.81 ± 1981.50</td><td>0.21</td><td>9389.91</td><td>45</td></tr><tr><td>Total</td><td>315 ± 617</td><td>1</td><td>2983</td><td>685.57 ± 2049.73</td><td>0.21</td><td>9389.91</td><td>59</td></tr><tr><td colspan="8">(b)</td></tr><tr><td>VM</td><td>316 ± 540</td><td>28</td><td>1987</td><td>29.67 ± 32.08</td><td>1.60</td><td>118.00</td><td>12</td></tr><tr><td>CM</td><td>123 ± 185</td><td>1</td><td>937</td><td>42.03 ± 53.29</td><td>0.21</td><td>189.84</td><td>41</td></tr><tr><td>Total</td><td>167 ± 308</td><td>1</td><td>1987</td><td>39.23 ± 49.29</td><td>0.21</td><td>198.84</td><td>53</td></tr></tbody></table>

located out of the microphones area whereas in CM, 28 were “out” and 17 were “in”. The mean distance between the microphones area centroid and the broadcast site was about 656.55 ± 422.09 m. The mean microphones area was 746,823 ± 342,362 m² (N = 59).

There was a positive correlation between the localization accuracy and the temporal error value (Pearson's correlation coefficient, *r* = 0.83, *P* < 0.001; Fig. 5), which indicated that the localization accuracy decreased when the temporal error increased. Based on this relationship, we identified a threshold in the temporal error above which the estimates were unreliable. After setting this reliability threshold to 200 ms, six inaccurate data (two in VM and four in CM) with high localization accuracy values (ranging between 534 and 3083 m) were excluded (see Fig. 5). The mean localization accuracy was then 167 ± 308 m (N = 53; Table 1). Considering the error threshold, all of the remaining data had accuracies less than 400 m, except three values having high localization accuracy and low temporal error values. These three aberrant data were excluded from the dataset used in the analysis of parameters influencing the localization accuracy.

### Parameters that influenced the localization accuracy

Eleven linear mixed-effects models were built in order to identify the parameters that influenced the localization accuracy (Table 2). Among the four fixed effects, only the broadcast site position relative to the microphones area (in or out) was not tested because both conditions were not present in VM (all sites positions were out). The mixed-model with the fixed effect “dist” had the lowest AIC (i.e. “m2”: AIC = 591.35). Thus, the distance between the microphones area centroid and the broadcast site significantly affected the localization accuracy ($\chi^2$ (1) = 11.27, *P* < 0.001). In particular, the localization accuracy was lower when the broadcast site

was far from the microphones area centroid (estimate ± SE: 0.14 ± 0.04 m).

Considering the random effect in the selected model, the localization accuracy did not vary between the two study areas (restricted likelihood ratio test [RLRT] = 0.06, *P* > 0.05) and there was no effect of period (RLRT = 0.05, *P* > 0.05). However, the localization accuracy varied significantly among the different nights inside each study area (RLRT = 4.53, *P* < 0.05; Fig. 6).

### Discussion

Since its natural return to France from the Italian population [53], the grey wolf first recolonized mountainous areas (French Alps) and its range is currently expanding west and northward into mid-mountain and lowland environments [19, 34–37]. Documenting and updating presence and localization of wolves is crucially important for managing this protected species and for preempting potential conflicts with human activities, especially livestock attacks. Thus, in this study, we investigated a new, non-invasive, and large-scale acoustic method for localizing wolves.

### Acoustic localization estimates

In our study, 76% of the broadcasts were recorded by at least one recorder, thereby demonstrating the potential for using a low-density microphone array to detect howls over large areas (30 km²) with contrasting environmental contexts. The 59 broadcasts recorded by at least four microphones were used to estimate localizations. Although accuracies did not differ significantly between the two study areas, we observed a difference in the sample size of the broadcasts recorded and used for localization estimation, particularly in the lowland environment (45 in CM) compared with the mid-mountain environment (14 in VM).

After considering the relationship between the localization accuracy and the temporal error value, we defined a reliability threshold for the temporal error. We set this threshold to 200 ms. Then, most of the inaccurate values were excluded conducing to a mean localization accuracy of less than 200 m. This value may be considered a poor localization estimate when compared with most studies of acoustic localization [2, 47, 54, 55]. However, these previous studies were conducted in much smaller study areas. Thus, considering the distance between the autonomous recorders in our experiments, a localization accuracy of 200 m appears to be consistent.

In addition, we showed that some parameters could influence the localization accuracy and so, should be considered to optimize future protocols. First, the localization accuracy varied among different nights in each study area and at the same broadcast site (relicate). The wind speed was negligible during the
<figure><img src="image_8.png" /><figcaption>**Fig. 5** Relationship between the localization accuracy and temporal error in both study areas. Pearson's correlation coefficient: *r* = 0.83, *P* < 0.001, *N* = 59. Vertical dotted line: 200 ms reliability threshold for the temporal error. Blue triangles: VM. Red circles: CM</figcaption></figure>

experiments, but variations in other meteorological conditions among different nights may explain the differences in accuracy (e.g. air temperature or wind direction). Indeed, the meteorological conditions are known to have strong effects on sound propagation and signal detection, and thus on the localization accuracy [56]. We also showed that the distance between the microphones area centroid and the broadcast site had a significant effect on the localization accuracy. The localization accuracy was lower when the broadcast site was far from the microphones area centroid, as shown in previous studies [2, 54, 57].

<table><caption>**Table 2** Linear mixed-effects models used to identify parameters that influenced the localization accuracy. Models were built with fixed effects alone, summed ("+") and with interactions ("*"). The nested random effect was similar for all of the models (see "m0"). Akaike's information criterion (AIC) was calculated to select the best statistical model, i.e. the "m2" model</caption><thead><tr><th colspan="4">Models</th></tr><tr><th>Name</th><th>Fixed effects</th><th>Random effect</th><th>AIC</th></tr></thead><tbody><tr><td>m0</td><td>1</td><td>(1|array) + (1|array:night) + (1|array:night:period)</td><td>600.62</td></tr><tr><td>m1</td><td>period</td><td></td><td>604.35</td></tr><tr><td>m2</td><td>dist</td><td></td><td>591.35</td></tr><tr><td>m3</td><td>areamic</td><td></td><td>602.62</td></tr><tr><td>m4</td><td>period+dist</td><td></td><td>594.98</td></tr><tr><td>m4.int</td><td>period*dist</td><td></td><td>591.70</td></tr><tr><td>m5</td><td>period+areamic</td><td></td><td>606.31</td></tr><tr><td>m5.int</td><td>period*areamic</td><td></td><td>604.10</td></tr><tr><td>m6</td><td>dist+areamic</td><td></td><td>593.26</td></tr><tr><td>m6.int</td><td>dist*areamic</td><td></td><td>594.79</td></tr><tr><td>mc</td><td>period+dist+areamic</td><td></td><td>596.95</td></tr></tbody></table>

## Recommendations and perspectives for grey wolf monitoring

According to our results, some recommendations may be made regarding the development of effective acoustic methods for grey wolf monitoring. The measures of the TOA were performed manually and this was a time-consuming task. Automatic and autonomous methods for detecting wolf howls in recordings and then for localizing them, such as methods based on temporal cross-correlation, could improve the results and save time. However, these methods are still very complicated [58, 59]. Moreover, amplitude and frequency modulations in wolf howls may make difficult to parameterize a unique automatic detector that could be trusted without human verification.

As shown in the present study and previous investigations (e.g. [2, 54, 57]), the distance between the sound source and the microphones area centroid influenced the localization accuracy. Similarly, during field recordings, large distances between the study species and recorders may also influence the localization accuracy because of a low signal-to-noise ratio (as found in our study). Thus, the selection of the recording sites should be optimized according to the ecology and behavior of wolves but also based on local expert knowledge in order to increase the likelihood of collecting acoustic data.

The structure and the composition of the landscape, such as the topography and vegetation (e.g. composition and stand density), could also influence the localization estimations, and thus they should be considered when defining protocols based on acoustic methods. This may partly explain the difference in the sample sizes for the broadcasts used in the lowland and mid-mountain environments. Previous studies also demonstrated that the optimal placement of recorders is important for ensuring maximum cover of the study area [29, 54].
<figure><img src="image_9.png" /><figcaption>**Fig. 6** Variations in the localization accuracy on the three different nights in both study areas. Random effect of linear mixed-effect model: RLRT = 4.53, *P* < 0.05</figcaption></figure>

These recommendations highlight the necessity to find
a compromise between the distance that separates the
microphones, the area covered by the microphone array,
the areas where vocalizations or sounds are produced,
and the desired localization accuracy [54]. Considering
our results and soundscape parameters, it would be in-
teresting to model the sound detection space of the au-
tonomous recorders in order to place them optimally in
the field and to improve the localization accuracy.

Finally, the acoustic localization protocol may concern
much more wolves living in pack rather than dispersers
or lone wolves (less frequent howls; [20]). This potential
limit could be balanced by combining autonomous
recorders with howling playback method to elicit wolves
to howl [60]. This would be even more recommended in
the colonization fronts (like in north-eastern France) for
monitoring wolf dispersion but also for detecting new
pack installation [61].

Conclusions

Currently, monitoring of the distribution and demo-
graphic dynamics of the grey wolf in France is based on
the standardized collection of presence signs by a network
of 3500 trained volunteers [62]. Different methods are
used such as opportunistic survey (scat, hair, saliva, etc.),
non-invasive genetics analysis, intensive snow-tracking

during the winter, and wolf howling in summer to detect
breeding events [62]. However, the potential use of acous-
tic and autonomous recorders has not been considered for
localizing individuals as well as specific areas such as
rendezvous sites in contrasting environments. Thus, the
development of a localization protocol based on passive
acoustic methods could help scientists and decision-
makers to collect new data to understand and monitor
wolf recolonization. Importantly, these data could help
prevent or mitigate conflicts with human activities as well
as being used for conservation and management purposes.

Today, more than ever, large scale studies for monitor-
ing elusive species are necessary and remain challenging
[63]. Localization protocols based on our results and
recommendations could be applied to species producing
long-distance acoustic signals, even in large territories
and contrasting environments. This kind of protocols
will considerably help to monitor the conservation status
of many elusive species in the long term.

Acknowledgments

We are very grateful to CROC members M. Clasquin, M. Marc, P. Germain, Dr. A. Charbonnel, and C. Botta for their help with fieldwork. We are also grateful to the National Game and Wildlife Agency (ONCFS) for authorization to access the wolves database, and to D. Wilson and J. Brzustowski for their advice regarding the Sound Finder package. We would like to thank the funding partners: the European Union within the framework of the Operational Program FEDER-FSE "Lorraine et Massif des Vosges 2014-2020",
the Commissariat à l'Aménagement du Massif des Vosges for the FNADT (Fonds National d'Aménagement et de Développement du Territoire), the DREAL Grand Est (Direction Régionale pour l'Environnement, l'Aménagement et le Logement), the Région Grand Est, the ANRT (Agence Nationale de la Recherche et de la Technologie, CIFRE award), the Zoo d'Amnéville, and the Parc Animalier de Sainte Croix. We also thank the Fondation Le Pal Nature for its complementary financial support. This study could not have been conducted without authorizations and agreements from municipalities, environmental managers (ONCFS, ONF, natural reserves, natural parks, etc.), and owners in the study areas. We also thank the two anonymous reviewers for their careful reading and their perceptive comments and suggestions.

## Funding

This work has been supported by the European Union within the framework of the Operational Program FEDER-FSE "Lorraine et Massif des Vosges 2014-2020", the Commissariat à l'Aménagement du Massif des Vosges for the FNADT (Fonds National d'Aménagement et de Développement du Territoire), the DREAL Grand Est (Direction Régionale pour l'Environnement, l'Aménagement et le Logement), the Région Grand Est, the ANRT (Agence Nationale de la Recherche et de la Technologie, CIFRE award), the Zoo d'Amnéville, and the Parc Animalier de Sainte Croix. It also has been supported by the Fondation Le Pal Nature. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

## Availability of data and materials

The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request. They are archived in CROC research center.

## Authors' contributions

Conceived the ideas and designed the methodology: MP JP EG. Collected the data: MP JP EG. Analyzed the data: MP JP. Led the writing of the manuscript: MP JP FG EG. All authors contributed critically to the drafts and gave final approval for publication.

## Ethics approval and consent to participate

Not applicable.

## Consent for publication

Not applicable.

## Competing interests

The authors declare that they have no competing interests.

## Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

## Author details

Received: 15 December 2017 Accepted: 25 March 2018

Published online: 12 April 2018

## References

1. Blumstein DT, Mennill DJ, Clemins P, Girod L, Yao K, Patricelli G, et al. Acoustic monitoring in terrestrial environments using microphone arrays: Applications, technological considerations and prospectus. J Appl Ecol. 2011;48:758-67.

2. Mennill DJ, Battiston M, Wilson DR, Foote JR, Doucet SM. Field test of an affordable, portable, wireless microphone array for spatial monitoring of animal ecology and behaviour. Methods Ecol Evol. 2012;3:704-12.

3. Borker AL, Halbert P, McKown MW, Tershy BR, Croll DA. A comparison of automated and traditional monitoring techniques for marbled murrelets using passive acoustic sensors. Wildl Soc Bull. 2015;39:813-8.

4. Ulloa JS, Gasc A, Gaucher P, Aubin T, Réjou-Méchain M, Sueur J. Screening large audio datasets to determine the time and space distribution of Screaming Piha birds in a tropical forest. Ecol Inform. 2016;31:91-9.

5. Bardeli R, Wolff D, Kurth F, Koch M, Tauchert KH, Frommolt K-H. Detecting bird sounds in a complex acoustic environment and application to bioacoustic monitoring. Pattern Recogn Lett. 2010;31:1524-34.

6. Dawson DK, Efford MG. Bird population density estimated from acoustic signals. J Appl Ecol. 2009;46:1201-9.

7. Marques TA, Thomas L, Martin SW, Mellinger DK, Ward JA, Moretti DJ, et al. Estimating animal population density using passive acoustics. Biol Rev Camb Philos Soc. 2013;88:287-309.

8. Kalan AK, Piel AK, Mundry R, Wittig RM, Boesch C, Kühl HS. Passive acoustic monitoring reveals group ranging and territory use: a case study of wild chimpanzees (Pan troglodytes). Front Zool. 2016;13:34.

9. Mennill DJ, Burt JM, Fristrup KM, Vehrencamp SL. Accuracy of an acoustic location system for monitoring the position of duetting songbirds in tropical forest. J Acoust Soc Am. 2006;119:2832-9.

10. Griffin A, Alexandridis A, Pavlidi D, Mastorakis Y, Mouchtaris A. Localizing multiple audio sources in a wireless acoustic sensor network. Signal Process. 2015;107:54-67.

11. Bradbury JW, Vehrencamp SL. Principles of Animal Communication. Sunderland: Sinauer Associates; 1998.

12. Hartwig S. Individual acoustic identification as a non-invasive conservation tool: an approach to the conservation of the African wild dog Lycaon pictus (Temminck, 1820). Bioacoustics. 2005;15:35-50.

13. Palacios V, López-Bao JV, Llaneza L, Fernández C. Decoding group vocalizations: The acoustic energy distribution of chorus howls is useful to determine wolf reproduction. PLoS One. 2016;11:e0153858.

14. Boitani L. Wolf conservation and recovery. In: Mech LD, Boitani L, editors. Wolves: behavior, ecology, and conservation. Chicago: University of Chicago Press; 2003. p. 317-40.

15. Valière N, Fumagalli L, Gielly L, Miquel C, Lequette B, Poulle M-L, et al. Long-distance wolf recolonization of France and Switzerland inferred from non-invasive genetic sampling over a period of 10 years. Anim Conserv. 2003;6:83-92.

16. Muhly TB, Musiani M. Livestock depredation by wolves and the ranching economy in the Northwestern U.S. Ecol Econ. 2009;68:2439-50.

17. Rigg R, Findo S, Wechselberger M, Gorman ML, Sillero-Zubiri C, Macdonald DW. Mitigating carnivore-livestock conflict in Europe: lessons from Slovakia. Oryx. 2011;45:272-80.

18. Latham MC, Latham ADM, Webb NF, Mccutchen NA, Boutin S. Can occupancy-abundance models be used to monitor wolf abundance? PLoS One. 2014;9:e102982.

19. Louvrier J, Duchamp C, Lauret V, Marboutin E, Cubaynes S, Choquet R, et al. Mapping and explaining wolf recolonization in France using dynamic occupancy models and opportunistic data. Ecography. 2017;40:1-13.

20. Harrington FH, Mech LD. Wolf howling and its role in territory maintenance. Behaviour. 1979;68:207-49.

21. Gazzola A, Avanzinelli E, Mauri L, Scandura M, Apollonio M. Temporal changes of howling in south European wolf packs. Ital J Zool. 2002;69:157-61.

22. Nowak S, Jędrzejewski W, Schmidt K, Theuerkauf J, Mysłajek RW, Jędrzejewska B. Howling activity of free-ranging wolves (Canis lupus) in the Bialowieża Primeval Forest and the Western Beskidy Mountains (Poland). J Ethol. 2007;25:231-7.

23. Passilongo D, Buccianti A, Dessi-Fulgheri F, Gazzola A, Zaccaroni M, Apollonio M. The acoustic structure of wolf howls in some eastern tuscany (Central Italy) free ranging packs. Bioacoustics. 2010;19:159-75.

24. Zaccaroni M, Passilongo D, Buccianti A, Dessi-Fulgheri F, Facchini C, Gazzola A, et al. Group specific vocal signature in free-ranging wolf packs. Ethol Ecol Evol. 2012;24:322-31.

25. Tooze ZJ, Harrington FH, Fentress JC. Individually distinct vocalizations in timber wolves, Canis lupus. Anim Behav. 1990;40:723-30.

26. Palacios V, Font E, Marquez R. Iberian wolf howls: acoustic structure, individual variation, and a comparison with north american populations. J Mammal. 2007;88:606-13.

27. Root-Gutteridge H, Bencsik M, Chebli M, Gentle LK, Terrell-nield C, Bourit A, et al. Improving individual identification in captive Eastern grey wolves (Canis lupus lycaon) using the time course of howl amplitudes. Bioacoustics. 2014;23:1-15.

28. Papin M, Pichenot J, Germain E. La bioacoustique: un outil prometteur pour l'estimation des effectifs de loups gris. In: 11e Rencontres Bourgogne-Nature et du 37e Colloque francophone de Mammalogie, Les Mammifères sauvages - Recolonisation et réémergence, Revue Scientifique Bourgogne. Nature. 2015;21/22:256-65.

<footnote n="1">¹Centre de Recherche et d'Observation sur les Carnivores (CROC), 4 rue de la banie, 57590, Lucy, France. ²Université de Lorraine, CNRS, Laboratoire Interdisciplinaire des Environnements Continentaux, F-57500 Metz, France.</footnote>

<footnote n="3">³Biologiste Écologue Consultant (BEC), 8A rue principale, 57590, Fonteny, France.</footnote>
29. Suter SM, Giordano M, Nietlispach S, Apollonio M, Passilongo D. Non-invasive acoustic detection of wolves. Bioacoustics. 2016;26:237–48.

30. Filibeck U, Nicoli M, Rossi P, Boscagli G. Detection by frequency analyzer of individual wolves howling in a chorus: A preliminary report. Bolletino di Zool. 1982;49:151–4.

31. Harrington FH. Chorus howling by wolves: acoustic structure, pack size and the Beau Geste Effect. Bioacoustics. 1989;2:117–36.

32. Sèbe F, Heintz N, Latini R, Aubin T. Le wolf howling: un outil pour le recensement et la conservation des loups. Possibilités et limites de la méthode. In: Colloque Grands Prédateurs et Pastoralisme; 2004. p. 53–9.

33. Passilongo D, Mattioli L, Bassi E, Szabó L, Apollonio M. Visualizing sound: counting wolves by using a spectral view of the chorus howling. Front Zool. 2015;12:12–22.

34. L'équipe animatrice du réseau. Du nouveau sur le front de colonisation Vosgien et Franc-comtois. Bulletin Loup du réseau. 2012;27:2.

35. Laurent A. Une meute d'au moins 4 loups installée dans les Vosges. Bulletin Loup du réseau. 2014;30:5.

36. Laurent A. Nouveaux indices de présence du loup en Meuse. Bulletin Loup du réseau. 2014;31:2.

37. L'équipe animatrice du réseau. Les données du réseau. Bulletin Loup du réseau. 2017;36:16–27.

38. Reed SE, Boggs JL, Mann JP. SPreAD-GIS: an ArcGIS toolbox for modeling the propagation of engine noise in a wildland setting. San Francisco: The Wilderness Society; 2010.

39. Reed SE, Boggs JL, Mann JP. A GIS tool for modeling anthropogenic noise propagation in natural ecosystems. Environ Model Softw. 2012;37:1–5.

40. De Reu J, Bourgeois J, Bats M, Zwertvaegher A, Gelorini V, de Smedt P, et al. Application of the topographic position index to heterogeneous landscapes. Geomorphology. 2013;186:39–49.

41. QGIS Development Team. Quantum GIS Geographic Information System. Open Source Geospatial Foundation Project: Beaverton, United States; 2014.

42. Mech LD. Where can wolves live and how can we live with them? Biol Conserv. 2017;210:310–7.

43. Theberge JB, Falls JB. Howling as a means of communication in Timber Wolves. Am Zool. 1967;7:331–8.

44. Sueur J, Aubin T, Simonis C. Seewave, a free modular tool for sound analysis and synthesis. Bioacoustics. 2008;18:213–26.

45. Harrington FH, Asa CS. Wolf communication. In: Mech LD, Boitani L, editors. Wolves: behavior, ecology, and conservation. Chicago: University of Chicago Press; 2003. p. 66–103.

46. Charif RA, Waack AM, Strickman LM. Raven Pro 1.4 User's Manual. Ithaca: Cornell Laboratory of Ornithology; 2010.

47. Wilson DR, Battiston M, Brzustowski J, Mennill DJ. Sound Finder: a new software approach for localizing animals recorded with a microphone array. Bioacoustics. 2014;23:99–112.

48. R Development Core Team. R: A language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2014.

49. Bates D, Mächler M, Bolker B, Walker S. Fitting linear mixed-effects models using lme4. J Stat Softw. 2015;67:1–48.

50. Burnham KP, Anderson DR. Model selection and multimodel inference: a practical information-theoretic approach. 2nd ed. New York: SpringerVerlag; 2002.

51. Venables WN, Ripley BD. Modern Applied Statistics with S. 4th ed. New York: Springer; 2002.

52. Scheipl F, Greven S, Küchenhoff H. Size and power of tests for a zero random effect variance or polynomial regression in additive and linear mixed models. Comput Stat Data Anal. 2008;52:3283–99.

53. Fabbri E, Miquel C, Lucchini V, Santini A, Caniglia R, Duchamp C, et al. From the Apennines to the Alps: Colonization genetics of the naturally expanding Italian wolf (Canis lupus) population. Mol Ecol. 2007;16:1661–71.

54. Bower JL, Clark CW. A field test of the accuracy of a passive Acoustic Location System. Bioacoustics. 2005;15:1–14.

55. Muanke PB, Niezrecki C. Manatee position estimation by passive acoustic localization. J Acoust Soc Am. 2007;121:2049–59.

56. Wiley RH, Richards DG. Physical constraints on acoustic communication in the atmosphere: Implications for the evolution of animal vocalizations. Behav Ecol Sociobiol. 1978;3:69–94.

57. Spillmann B, van Noordwijk MA, Willems EP, Mitra Setia T, Wipfli U, Van Schaik CP. Validation of an acoustic location system to monitor Bornean orangutan (Pongo pygmaeus wurmbii) long calls. Am J Primatol. 2015;77:767–76.

58. Wimmer J, Towsey M, Planitz B, Roe P, Williamson I. Scaling acoustic data analysis through collaboration and automation. In: Proceedings - 2010 6th IEEE International Conference on e-Science, eScience 2010; 2010. p. 308–15.

59. Ausband DE, Skrivseth J, Mitchell MS. An automated device for provoking and capturing wildlife calls. Wildl Soc Bull. 2011;35:498–503.

60. Rocha LHS, Ferreira LS, Paula BC, Rodrigues FHG, Sousa-Lima RS. An evaluation of manual and automated methods for detecting sounds of maned wolves (Chrysocyon brachyurus Illiger 1815). Bioacoustics. 2015;24:185–98.

61. Brennan A, Cross PC, Ausband DE, Barbknecht A, Creel S. Testing automated howling devices in a wintertime wolf survey. Wildl Soc Bull. 2013;37:389–93.

62. Duchamp C, Boyer J, Briaudet PE, Leonard Y, Moris P, Bataille A, et al. A dual frame survey to assess time- and space-related changes of the colonizing wolf population in France. Hystrix. 2012;23:14–28.

63. Llaneza L, García EJ, López-Bao JV. Intensity of territorial marking predicts wolf reproduction: Implications for wolf monitoring. PLoS One. 2014;9:e93015.

Submit your next manuscript to BioMed Central and we will help you at every step:

* We accept pre-submission inquiries

* Our selector tool helps you to find the most relevant journal

* We provide round the clock customer support

* Convenient online submission

* Thorough peer review

* Inclusion in PubMed and all major indexing services

* Maximum visibility for your research

Submit your manuscript at
www.biomedcentral.com/submit

BioMed Central