Evaluating the Effectiveness of an Ultrasonic Acoustic Deterrent for Reducing Bat Fatalities at Wind Turbines

Large numbers of bats are killed by wind turbines worldwide and minimizing fatalities is critically important to bat conservation and acceptance of wind energy development. We implemented a 2-year study testing the effectiveness of an ultrasonic acoustic deterrent for reducing bat fatalities at a wind energy facility in Pennsylvania. We randomly selected control and treatment turbines that were searched daily in summer and fall 2009 and 2010. Estimates of fatality, corrected for field biases, were compared between treatment and control turbines. In 2009, we estimated 21–51% fewer bats were killed per treatment turbine than per control turbine. In 2010, we determined an approximate 9% inherent difference between treatment and control turbines and when factored into our analysis, variation increased and between 2% more and 64% fewer bats were killed per treatment turbine relative to control turbines. We estimated twice as many hoary bats were killed per control turbine than treatment turbine, and nearly twice as many silver-haired bats in 2009. In 2010, although we estimated nearly twice as many hoary bats and nearly 4 times as many silver-haired bats killed per control turbine than at treatment turbines during the treatment period, these only represented an approximate 20% increase in fatality relative to the pre-treatment period for these species when accounting for inherent differences between turbine sets. Our findings suggest broadband ultrasound broadcasts may reduce bat fatalities by discouraging bats from approaching sound sources. However, effectiveness of ultrasonic deterrents is limited by distance and area ultrasound can be broadcast, in part due to rapid attenuation in humid conditions. We caution that an operational deterrent device is not yet available and further modifications and experimentation are needed. Future efforts must also evaluate cost-effectiveness of deterrents in relation to curtailment strategies to allow a cost-benefit analysis for mitigating bat fatalities. This is an open-access article, free of all copyright, and may be freely reproduced, distributed, transmitted, modified, built upon, or otherwise used by anyone for any lawful purpose. The work is made available under the Creative Commons CC0 public domain dedication.


Introduction
As wind energy production has steadily increased worldwide, bat fatalities have been reported at wind facilities worldwide [1,2,3,4] in a wide range of landscapes. A recent synthesis reported that approximately 650,000 to more than 1,300,000 bats have been estimated to have been killed from 2000-2011 in the U.S. and Canada [5]. Given these fatality rates, accelerating growth of the wind industry [6], and suspected and known population declines in many species of bats [7,8,9], it is imperative to develop and implement solutions to reduce future bat fatalities at wind facilities.
Prior studies have demonstrated that a substantial portion of bat fatalities consistently occur during relatively low-wind conditions over a relatively short period of time during the summer-fall bat migration period [2,4]. Curtailment of turbine operations under these conditions and during this period has been proposed as a possible means of reducing impacts to bats [1,2,10]. Indeed, recent studies in Canada [11] and the U.S. [12] indicate that increasing turbine ''cut-in speed'' (i.e., wind speed at which wind-generated electricity enters the power grid) from the manufactured speed (usually 3.5-4.0 m/s for modern turbines) to between 5.0 and 6.5 m/s resulted in at least a 50% reduction in bat fatalities (and as high as 93%) compared to normally operating turbines [12]. While costs of lost power from curtailment can be factored into the economics and financing and power purchase agreements of new projects, altering turbine operations even on a partial, limited-term basis potentially poses operational and financial difficulties for existing projects, so there is considerable interest in developing other solutions to reduce bat fatalities that do not involve turbine shutdowns. Also, changing turbine cut-in speed may not be effective in other regions that experience bat fatalities although this strategy may ultimately prove sufficiently feasible and economical for reducing bat fatalities. Thus, research on alternative mitigation strategies and their associated costs are warranted.
Studies in Scotland suggest that bat activity may be deterred by electromagnetic signals from small, portable radar units [13]. This study reported that bat activity and foraging effort per unit time were significantly reduced during experimental trials when their radar antenna was fixed to produce a unidirectional signal that maximized exposure of foraging bats to their radar beam. The effectiveness of radar as a potential deterrent has not been tested at an operating wind facility to determine if bat fatalities could be significantly reduced by these means. Moreover, the effective range of electromagnetic signals as well as the number of radar units needed to affect the most airspace near individual turbines would need to be determined to fully evaluate effectiveness and to allow some cost-benefit analysis relative to other potential deterrents or curtailment [11,12].
Echolocating bats produce high frequency vocal signals and perceive their surroundings by listening to features of echoes reflecting from targets in the path of the sound beam [14]. Thus, bats that use echolocation depend heavily on auditory function for orientation, prey capture, communication, and obstacle avoidance. Bats of some species avoid certain territorial social calls emitted by conspecifics [15] and are deterred by ''clicks'' emitted by noxious moths [16]. Because echolocating bats depend upon sensitive ultrasonic hearing, we hypothesized that broadcasting ultrasound from wind turbines may disrupt or ''jam'' their perception of echoes and serve as a deterrent. Such masking of echo perception, or simply broadcasting high intensity sounds at a frequency range to which bats are most sensitive, could create an uncomfortable or disorienting airspace that bats may prefer to avoid.
Few studies have investigated the influence of ultrasound broadcast on bat behavior and activity, particularly under field conditions. Broadband random ultrasonic noise may mask bat echolocation somewhat, but not completely [17]. Ultrasound broadcasts can reduce bat activity, perhaps due to greater difficulty in the bats hearing echoes of insects and thus reduced feeding efficiency [18]. A laboratory test of the response of big brown bats (Eptesicus fuscus) to a prototype eight speaker deterrent device emitting broadband white noise at frequencies ranging from 12.5-112.5 kHz in the laboratory and found that during non-feeding trials, bats landed in a quadrant containing the device significantly less when it was broadcasting broadband noise (J. Spanjer, University of Maryland and E. Arnett, Bat Conservation International, unpublished data). During feeding trials in this experiment, bats never successfully captured a tethered mealworm when the device broadcasted sound but captured mealworms near the device in about 1/3 of trials when it was silent. Field tests of the same acoustic deterrent found that when placed by the edge of a small pond, where nightly bat activity was consistent, nightly activity decreased significantly on nights when the deterrent was activated (J. Szewczak, Humboldt State University and E. Arnett, Bat Conservation International, unpublished data).
Our goal was to improve deterrent devices we previously developed and tested by increasing the effective area of ultrasonic emissions from the nacelle of wind turbines, and to test their effectiveness on reducing bat fatalities. The objectives of this study were 1) to conduct carcass searches and field bias trials (searcher efficiency and carcass removal) to determine rate of bat fatalities at treatment (those with deterrent devices) and control turbines; and 2) compare bat fatality rates at turbines treatment and control turbines to determine effectiveness. We successfully tested our ultrasonic deterrent device at an operating wind facility and offer suggestions for future efforts regarding this potential mitigation strategy to reduce bat fatalities at wind facilities.

Study Area
The study was conducted at the Locust Ridge Wind Project located near the towns of Shenandoah, Mahanoy City, and Brandonville in Columbia and Schuylkill Counties, Pennsylvania, and consisted of two different facilities. The Locust Ridge I (LRI) Wind Farm has 13 Gamesa G87 2.0 MW turbines, each on 80 m monopoles with a rotor diameter of 87 m and a swept area of 5,945 m 2 . There were 51 Gamesa G83 2.0 MW turbines, each on 80 m monopoles with a rotor diameter of 83 m and a swept area of rotor-swept area of 5,411 m 2 , at the Locust Ridge II (LRII) Wind Farm. The facilities lie within the Appalachian mixed mesophytic forests ecoregion and the moist broadleaf forests that cover the plateaus and rolling hills west of the Appalachian Mountains [19,20]. Elevations along ridges where turbines are located range from 530-596 m. All turbines were located along a ridge in deciduous forest, with some species of evergreen trees interspersed. Vegetation across the area included thickets of scrub oak (Quercus berberidifolia) interspersed with chestnut oak (Quercus prinus) and gray birch (Betula populifolia), and mature hardwood forests of red oak (Quercus rubra), red maple (Acer rubrum), yellow birch (Betula alleghaniensis), American beech (Fagus grandifolia) and scrub oak.

Turbine Selection and Deterrent Installation
We randomly selected 15 of the 51 turbines located at LR II to be searched as part of a separate study to determine postconstruction fatality rates and to meet permitting requirements of the Pennsylvania Game Commission's voluntary agreement for wind energy [21]. These 15 turbines were our control turbines for comparing with treatment turbines, those fitted with deterrent devices. In 2009, unforeseen mechanical and safety issues arose at the LRII site and many of these turbines had to be excluded from our potential treatment group due to potential safety hazards. Thus, we included all 13 turbines at LRI as well as the remaining available turbines at LRII (n = 36) when randomly selecting our 10 turbines to be fitted with deterrent devices; 3 turbines were randomly selected from the 13 available at the LRI site and 7 of 36 available at LRII.
We did not assess whether there were any potential inherent differences between the two types of turbines in 2009 and for this year assumed there were no confounding differences in our findings. However, in 2010, we attempted to assess inherent differences between control and treatment turbines by modifying our design and analysis to reflect a Before-After Control-Impact (BACI) design. The same sets of control and treatment turbines were monitored for a period of time prior to implementation of the deterrent treatment (1 May to 26 July 2010), then again during the deterrent implementation period (31 July through 9 October 2010). This design allowed for incorporating initial inherent differences between the two experimental treatment sets prior to implementation of the treatment as a reference for interpreting any differences detected during implementation of the treatment.
The deterrent devices used in our study consisted of a waterproof box (,45645 cm, ,0.9 kg) that housed 16 transducers that emitted continuous broadband ultrasound from 20-100 kHz (manufactured by Deaton Engineering, Georgetown, Texas; Figure S1). We did not test other types of emissions (e.g., short pulses) concurrent with broadband emission because more devices and sample turbines would be required, thus resulting in cost and sample size constraints. Transducers we used had an optimum transmission level at their resonant frequency of 50 kHz transmission and reduced transmit levels at higher and lower frequencies over a broadband range of 20-100 kHz. This frequency range overlaps that of all bats known in the study area. Three factors influence the predicted effective transmitted power at a given distance: 1) the original transmitted power (sound pressure level; SPL); 2) attenuation with distance due to the wave front spreading (inversely proportional to the square of the distance, frequency independent); and 3) attenuation (absorption) in air of the sound wave (dependent on frequency, humidity and distance; Tables S1 and S2). The following discussion describes our estimation to base the target signal level of our deterrent: A typical bat emits calls at about 110 dB sound pressure level (SPL) at 10 cm [22]. During search phase flight a typical North American species of bat emits about 12 calls per second, each about 5 milliseconds in duration [23,24]. Given the speed of sound at 340 m/sec and duration of an open air call, the bat's own call will theoretically mask echoes returning from objects within about 1.5 m (i.e., the bat cannot hear early return echoes while vocalizing). An echo from a target about 1.5 m away will return about 45 dB less than the original 110 dB signal, or at about 65 dB. The bat's next call would mask echoes returning from about 25 m away. By this first order estimation, a bat would theoretically perceive information from returning echoes with amplitudes of #65 dB over a range from about 1.5-25 m. Thus, we estimated that a broadband signal of $65 dB would begin jamming or masking most bat's echo perception from targets beyond about a 1.5 m range.
We attached 8 individual deterrent devices to the nacelle of each of 10 sample turbines. Three devices on each side of the nacelle were evenly spaced and pointed downward with one aimed into the rotor-swept area, one parallel with the monopole, and one aimed toward the back of the nacelle ( Figures S2 and S3). Additionally, two devices were aimed at reflector plates; one that projected emissions into the upper part of the rotor-swept area, and one toward the rear of the nacelle. All devices connected to control boxes that were powered from outlets located in the nacelle and each was set on a timer to operate from K hour before sunset to K hour after sunrise each night of the study.

Delineation of Carcass Search Plots and Habitat Mapping
We delineated a rectangular plot 126 m north-south by 120 m east-west (60 m radius from the turbine mast in any direction; 15,120 m 2 total area) centered on each turbine sampled; this area represents the maximum possible search area for this study. Transects were set 6 m apart within each plot and in an east-west direction, due to the topography and layout of turbines at this facility. However, dense vegetation and the area cleared of forest at this facility was highly varied and, thus, we eliminated unsearchable habitat (e.g., forest) and usually did not search the entire possible maximum area. We used a Trimble global positioning system (GPS) to map the actual area searched at each turbine. The density-weighted area searched was used to standardize results and adjust fatality estimates (see statistical methods). The habitat visibility classes within each plot were also mapped using a GPS unit. We recorded the percent ground cover, height of ground cover (low [,11 cm], medium [11-50 cm], high [.50 cm]), type of habitat (vegetation, brush pile, boulder, etc), and the presence of extreme slope and collapsed these habitat characteristics into visibility classes that reflect their combined influence on carcass detectability (Table S3) [21].

Fatality Searches
We conducted daily searches at 15 control turbines and 10 treatment turbines from 15 August to 10 October 2009 and 1 May to 26 July and 31 July to 9 October 2010. Each searcher completed 5-7 turbine plots each day during the study. Searchers walked at a rate of approximately 10-20 m/min along each transect searching out to 3 m on each side for fatalities. Searches were abandoned only if severe or otherwise unsafe weather (e.g., heavy rain, lightning) conditions were present and searches were resumed that day if weather conditions permitted. Searches commenced at sunrise and all turbines were searched within 8 hr after sunrise.
We recorded date, start time, end time, observer, and weather data for each search at turbines. When a dead bat or bird was found, the searcher placed a flag near the carcass and continued the search. After searching the entire plot, the searcher returned to each carcass and recorded information on date, time found, species, sex and age (where possible), observer name, identification number of carcass, turbine number, perpendicular distance from the transect line to the carcass, distance from turbine, azimuth from turbine, habitat surrounding carcass, condition of carcass (entire, partial, scavenged), and estimated time of death (e.g., #1 day, 2 days, etc.). A field crew leader confirmed all species identifications at the end of each day. Disposable nitrile gloves were used to handle all carcasses to reduce possible human scent bias for carcasses later used in scavenger removal trials. Each carcass was placed into a separate plastic bag and labeled. Fresh carcasses, those determined to have been killed the night immediately before a search, were redistributed at random points on the same day for searcher efficiency and scavenging trials.

Ethics Statement
This study was conducted on private property and authorized by the operator Iberdrola Renewables Locust Ridge Wind LLC. All downed bats were euthanized, even if no physical injury was observed due to the possibility of barotraumas, following requirements by the Pennsylvania Game Commission protocol [21] and using acceptable methods suggested by the American Society for Mammalogists [25]; because sedation or anesthesia was not used in our study, we employed cervical dislocation. Our work did not require approval by an Institutional Animal Care and Use Committee and all aspects of the field work and permission to collect carcasses found during our study was conducted under the auspices of permits issued each year by the state of Pennsylvania, Pennsylvania Game Commission, and the U.S. Fish and Wildlife Service.

Field Bias Trials
Searcher efficiency and removal of carcasses by scavengers (herein referred to as carcass persistence) was quantified to adjust estimates of total bat and bird fatalities for detection bias. We conducted bias trials throughout the entire study period and searchers were never aware which turbines were used or the number of carcasses placed beneath those turbines during trials. Prior to the study's inception, we generated a list of random turbine numbers and random azimuths and distances (m) from turbines for placement of each bat used in bias trials.
We used only fresh killed bats for searcher efficiency and carcass removal trials during the study. At the end of each day's search, a field crew leader gathered all carcasses from searchers and then redistributed fresh bats at predetermined random points within any given turbine plot's searchable area. Data recorded for each trial carcass prior to placement included date of placement, species, turbine number, distance and direction from turbine, and visibility class surrounding the carcass. We attempted to distribute trial bats equally among different visibility classes throughout the study period and succeeded in distributing roughly one-third of all trial bats in each visibility class (easy, moderate, and difficult; difficult and very difficult were combined). We attempted to avoid ''over-seeding'' any one turbine with carcasses by placing no more than 4 carcasses at any one time at a given turbine. Because we used fresh bats for searcher efficiency trials and carcass removal trials simultaneously, we did not mark bats with tape or some other previously used methods [26] that could impart human or other scents on trial bat carcasses. Rather, we used trial bat placement details (i.e. azimuth, distance, sex, species) and signatures from hair and tissue samples (i.e., hair removed between the scapulae and wing punches) to distinguish them from other fatalities landing nearby. Each trial bat was left in place and checked daily by the field crew leader or a searcher not involved with the bias trials at turbines where carcasses were placed. Thus, trial bats were available to be found by searchers on consecutive days during daily searches unless removed by a scavenger. We recorded the day that each bat was found by a searcher, at which time the carcass remained in the scavenger removal trial. If, however, a scavenger removed a carcass before detection it was removed from the searcher efficiency trial and used only in the removal data set. When a bat carcass was found, the searcher determined if a bias trial carcass had been found by looking for markings described above and contacting the crew leader to determine if the location (direction and distance) matched any possible trial bats. All trial bats were left in place for the carcass removal trial. Carcasses were left in place until removed by a scavenger or they decayed and disintegrated to a point beyond recognition. Carcass condition was recorded daily up to 20 days, as present and observable or missing or no longer observable.

Statistical Methods
Carcass persistence/removal. Estimates of the probability that a bat carcass was not removed in the interval between searches were used to adjust carcass counts for removal bias. Removal included scavenging, wind or water, or decomposition beyond recognition. In most fatality monitoring efforts, it is assumed that carcass removal occurs at a constant rate that is not dependent on the time since death; this simplifying assumption allows us to estimate fatality when search intervals exceed one day. The length of time a carcass remains on the study area before it is removed is typically modeled as an exponentially distributed random variable. The probability that a carcass is not removed during an interval of length I can be approximated as the average probability of persisting given its death might have occurred at any time during the interval: where: r r jk is the estimated probability that a carcass in the k th visibility class that died during the interval preceding the j th search will not be removed by scavengers; t t t t jk is the estimated average persistence time of a carcass in the k th visibility class that died during the interval preceding the j th search; I ij is the length of the effective interval preceding the j th search at the i th turbine; Data from 351 and 408 bat carcasses in 2009 and 2010, respectively, were used in our analysis, with carcass persistence time modeled as a function of visibility class. We fit carcass persistence/removal data for bats to an interval-censored parametric failure time model, with carcass persistence time modeled as a function of size and/or visibility class. We used a relatively liberal alpha of 0.15 to identify factors (e.g., carcass size, visibility classes) that influence bias parameter values (i.e., searcher efficiency and carcass persistence) for removal of bat carcasses.
Searcher efficiency. Estimates of the probability that an observer will visually detect a carcass during a search were used to adjust carcass counts for observer bias. Failure of an observer to detect a carcass on the search plot may be due to its size, color, or time since death, as well as conditions in its immediate vicinity (e.g., vegetation density, shade). In most fatality monitoring efforts, because we cannot measure time since death, it is assumed that a carcass' observability is constant over the period of study, which it likely is not. In this study, searches were conducted daily and carcass persistence times were long, providing an opportunity for a searcher to detect a carcass that was missed on a previous search. We used a newly derived estimator [27] that assumes a carcass missed on a previous search will not be observed on a subsequent search (i.e., there are inherent environmental conditions that make the carcass unobservable like heavy foliage, terrain). If this assumption is not met, it can lead to overestimates of fatality. Other estimators [26] assume that a carcass missed on a previous search has the same probability of being observed as it had on the first search (i.e., there is nothing inherent in the environment surrounding the carcass that makes it unobservable), missing it is purely a chance event and that if the carcass is not removed by predators and enough searches are conducted, it will eventually be observed. If this assumption is not met, it can lead to underestimates of fatality. It is likely that neither assumption is appropriate in all cases.
Searcher efficiency trial carcasses were placed on search plots and monitored for 20 days. The day on which a bat carcass was either observed or removed by a scavenger was noted. In these trial data, if a carcass had not been found within the first 8 searches it had essentially no chance of being found. This lends empirical support to the idea that there are some environmental conditions surrounding the carcass that determine its probability of being found. However, several carcasses missed on the first search were found on subsequent searches, lending support to the idea that at least for some carcasses, the probability of missing them is purely a chance event. To allow for some possibility of observing a carcass once having missed it, the set of trial carcasses comprised those found or still observable but not found within the first 8 searches. After accounting for carcasses removed before a searcher had the chance of observing them, we fit data from 139 (2009) and 169 (2010) bat carcasses to a logistic regression model, with odds of observing a carcass given that it persisted, modeled as a function of visibility class. Again, we used a relatively liberal alpha of 0.15 to determine if a significant effect among visibility classes existed. Because we found no bats in the very difficult visibility class, SE was not modeled for this class.
Density of carcasses and proportion of area surveyed. Density of carcasses is known to diminish with increasing distance from the turbine [26], so a simple adjustment to fatality based on area surveyed would likely lead to overestimates, because unsearched areas tend to be farthest from turbines where carcass density is lowest. The calculated function (see below) relating density to distance from a turbine was used to weight each square meter in the plot. The density-weighted fraction of each plot that was actually searched was used as an area adjustment to per-turbine fatality estimates rather than using a simple proportion.
The density of bat carcasses (number of carcasses/m 2 ) was modeled as a function of distance (m) from the turbine. Because searcher efficiency and visibility class are confounded with distance, only fresh bat carcasses found in Easy visibility class were used for this analysis and all non-incidental data from all searched turbines were used, yielding a total of 172 fresh bat carcasses. We assumed that the carcass persistence time and searcher efficiency would be equal for all carcasses within this class and would not change as a function of distance from the turbine. We also assumed that no bat carcasses killed by turbine blades would fall .200 m from the turbine. Carcasses were ''binned'' into 2 m rings extending from the turbine edge out to the theoretical maximum plot distance ( Figure S4). We determined the total area among all search plots that was in the easy visibility class (m 2 ) in each ring and calculated carcass density (number of carcasses/m 2 ) in each ring. Density was modeled as a conditional cubic polynomial function of distance (dist): If distance #50 m, then density = exp (21.77328+0.0346454*dist 20.00271076* dist 2 +0.0000229885* dist 3 ) 2 0.01, else density = 0.009363847*exp (20.05*(distance-50)).
Relative density was derived by dividing the predicted density of each m 2 unit by the total predicted density within 200 m of a turbine, providing a density-weight for each m 2 unit. The density weighted area (DWA) of a plot was calculated as the sum of the density weights for all m 2 units within the searchable area. If no portion of a designated plot was unsearchable, the density weight for the plot would be 1. The physical area surveyed within a plot differed among turbines and ranged from 20-47% of the delineated theoretical maximum search plot, with an average of 31% whereas the weighted density area of plots averaged 62% (range: 44-78%). In addition, using this density weight, we estimated 7.2% of the carcasses killed at a turbine would be found beyond the boundaries of the designated search plot.
Fatality estimates. We adjusted the number of bat fatalities found by searchers by estimates of searcher efficiency and by the proportion of carcasses expected to persist unscavenged during each interval using the following equation: where: f f ijk is the estimated fatality in the k th visibility class that occurred at the i th turbine during the j th search; c ijk is the observed number of carcasses in the k th visibility class at the i th turbine during the j th search; a a i is the density-weighted proportion of the area of the i th turbine that was searched; p p jk is the estimated probability that a carcass in the k th visibility class that is on the ground during the j th search will actually be seen by the observer; r r jk is the probability than an individual bird or bat that died in the k th visibility class during the interval preceding the j th search will not be removed by scavengers; and e e jk is the effective interval adjustment (i.e., the ratio of the length of time before 99% of carcasses can be expected to be removed to the search interval) associated with a carcass in the k th visibility class that died during the interval preceding the j th search.
The value forp p jk was estimated through searcher efficiency trials with estimates given above;r r jk is a function of the average carcass persistence rate and the length of the interval preceding the j th search; andr r jk ,ê e jk andp p jk are assumed not to differ among turbines, but differ with search interval (j) and visibility class (k).
The estimated annual per turbine fatality for bats and birds was calculated using a newly derived estimator [27] and the equation is:f where n i is the number of searches carried out at turbine i, 1 = 1, …, 10, andf f ijk is defined above. The per turbine estimate and confidence limits were multiplied by 64, the total number of turbines, and divided by 0.9279 to adjust for actual densityweighted area searched to give total annual fatality estimates [28]. This estimate assumes that no fatalities occurred during the winter, i.e. prior to April and after November. No closed form solution is yet available for the variance of this estimator, so 95% confidence intervals of this estimate were calculated by bootstrapping [29]. Searcher efficiency was estimated from a bootstrap sample (with replacement) of searcher efficiency data, carcass persistence estimated from a bootstrap sample of carcass persistence data, and these values were applied to the carcass data from a bootstrap sample of turbines to estimate average fatality per turbine. This process was repeated 1000 times. The 2.5 th and 97.5 th quantiles from the 1,000 bootstrapped estimates formed the 95% confidence limits of the estimated fatality [27].
Comparison between treatment and control turbines. In 2009, we compared average fatality at control with treatment turbines for all bats and for each species using one-way analysis of variance with each turbine as the experimental unit and log e transformed estimated total fatalities as the response. In 2010, estimated average bat fatality per turbine at control and treatment turbines, during the treatment phase and the period immediately preceding it (pre-treatment phase) was analyzed using a BACI approach [26,30,31], employing ANOVA repeated measures with the turbine as the experimental unit, repeatedly measured twice. Our approach determined whether the ratio of average perturbine fatality at control turbines (n = 15) to treatment turbines (n = 10) during implementation of the deterrents was significantly greater than it was in the period immediately preceding implementation of the treatments. In both years, the fatality data were log transformed to satisfy assumptions of normality and homogeneity of variance [32].

Results
In 2009, we searched 15 control turbines and 10 treatment turbines each day between 15 August and 10 October, and did not assess inherent variability among turbines. We found 194 carcasses (135 at control, 59 at treatment) of 6 species and two carcasses were not identifiable to species in 2009 (Table 1). During the pretreatment period between 1 May and 26 July 2010, we searched 15 control turbines daily for all but 2 days (16 May and 2 June) and 10 Deterrent turbines daily for all but 4 days (9, 20, 24 25 July 2010) due to heavy rain, or facility maintenance. During the treatment period between 1 August and 15 October, we searched 15 control turbines daily for all but 4 days (26 August; 22, 29, 30 September 2010) and 10 Deterrent turbines daily for all but 3 days (19 August; 9, 30 September 2010) due to heavy rain or facility maintenance. During the pre-treatment period from 1 May to 26 July 2010, we found 59 carcasses comprising 6 species of bats (37 at control, 22 at treatment; Table 2). During the treatment period, we found 223 carcasses comprising 6 species of bats (162 at control, 61 at Deterrent; Table 3). Fatalities were found at all 25 turbines searched and time required to search each plot ranged from 12-100 minutes in both years of the study. Based on data from turbines not equipped with deterrents, the estimated fatality rate for this site ranged from 16-29.

Fatality Estimates in 2009
A total of 278 trial carcasses were used to estimate searcher efficiency in this study. One hundred thirty-nine of the 145 (96%) carcasses in the easy class that persisted .7 days were found by searchers, while 105 of the 123 (85%) carcasses in the moderate class that persisted long enough to be observed were found. Eight of 10 (80%) carcasses in the difficult class were found. A logistic regression model of the odds of detection given persistence as a function of visibility classes was fit to the data and there was strong evidence of a difference in searcher efficiency among the visibility classes (x 2 2 = 10.32, p = 0.006). Data from 351 scavenger removal trial carcasses were fit to an interval-censored parametric failure time model. Average carcass persistence time was found to be strongly related to visibility classes (x 2 2 = 6.58, p = 0.037). Average persistence time was estimated to be 9.4 days (95% CI: 7.7, 11.7 days), 13.9 days (95% CI: 10.8, 18.3 days) and 8.7 days (95% CI: treatment 4.6, 16.1 days) in easy, moderate and difficult visibility classes respectively. Estimates of the probability of a bat carcass persisting for 1 day (r) were 0.948 (95% CI: 0.938, 0.958), 0.964 (95% CI: 0.955, 0.973) and 0.942 (95% CI: 0.900, 0.970), respectively.
The average per-turbine fatality rate at treatment turbines was significantly less than at control turbines (F 1,23 = 14.7, p,0.001).
We estimated twice as many hoary bats ( x x = 2.09, 95% CI = 1.18, 4.04) killed per control turbine than treatment turbine, and nearly twice as many silver-haired bats ( x x1.88, 95% CI = 0.92, 5.14), although the estimated effect was not significant for this species (Tables 3 and 4). Results for other species were highly variable with no statistically significant difference between turbine groups (Tables 4 and 5).

Fatality Estimates in 2010
A total of 169 bat carcasses were used to estimate searcher efficiency in this study. Eighty three of 86 (97%) carcasses in the easy class that persisted .7 days were found by searchers, while 59 of 70 (84%) carcasses in the moderate class that persisted long enough to be observed were found. Eight of 13 (62%) carcasses in the difficult class were found. Because no fatalities were found in the very difficult class, we removed the 6 bats placed in this class from our analysis. A logistic regression model of the odds of detection given persistence was fit to the visibility classes and there was strong evidence of a difference in searcher efficiency among the visibility classes (x 2 2 = 14.59, p = 0.007). Data from 408 scavenger removal trial carcasses were fit to an interval-censored parametric failure time model. Average carcass persistence time was found not to be related to visibility class (x 2 2 = 0.56, p = 0.907), but there was moderate evidence that average persistence time was longer before the treatment period than during the treatment period (x 2 2 = 4.27, p = 0.12). Average persistence time was estimated to be 7.8 days (95% CI: 6.4, 9.6  Little brown 1 2 ---3 Eastern red 6 1 ---7 Bat fatality data from the pre-treatment period were used to evaluate if there were inherent difference between control and treatment turbines. We determined there was marginal evidence that the ratio of control:treatment fatalities was greater during the treatment period than in the pre-treatment period (F 1,23 = 3.9, p = 0.061). During the pre-treatment period, prior to implementation of the deterrents, fatality per control turbine was estimated to be 1.09 times greater than per treatment turbine (95% CI: 0.74-1.61). We determined the initial inherent difference was about 9% in the fatality rate between the two sets and, while this was not statistically significant, we chose to adjust our comparison of fatalities between control and treatment turbines accordingly.
During the treatment period, we estimated an average of 12.8 bats (95% CI: 9.5, 17.2) were killed per turbine at treatment turbines compared to 22.9 bats (95% CI: 18.0, 29.3) killed per turbine at control turbines. Bat fatalities per control turbine was estimated to be 1.8 times greater than per treatment turbine (95% CI: 1.22-2.64); in other words, 18-62% fewer bats killed per treatment turbines relative to control turbines during the treatment period. As stated above, however, we determined an approximate 9% inherent difference between treatment and control turbines and fatality per control turbine was estimated to be 1.09 times greater than per treatment turbine (95% CI: 0.74-1.61) prior to implementation of the treatment. Thus, the ratio of fatality per control turbine relative to treatment turbines after implementing the treatment was estimated to be 1.64 times greater than the pre-treatment period ratio (95% CI: 0.98, 2.76). In other words, between 2% more and 64% fewer bats were killed per treatment turbine relative to control turbines after accounting for inherent turbine differences prior to treatment implementation.
In 2010, prior to implementation of the deterrent treatment, we estimated 1.47 times as many hoary bats (95% CI = 0.39, 3.42) and 1.32 times as many silver-haired bats (95% CI = 0.47, 3.27) killed per control turbine than treatment turbine. Although we estimated nearly twice as many hoary bats ( x x = 1.88, 95% CI = 1.19, 2.82) and nearly 4 times as many silver-haired bats ( x x = 3.78, 95% CI = 1.12, 12.82; Tables 5 and 6) killed per control turbine than treatment turbine during the treatment period, these represented only about a 20% increase in fatality relative to the pre-treatment period. High variation among turbines, small numbers of carcasses found and frequent zero-counts of these and other species at each turbine prevented formal statistical tests of these ratios using the BACI design (Tables 6 and 7).

Discussion
Previous research has indicated difficulty when attempting to mask or ''jam'' bats' echolocation except under specific conditions [e.g., 17,33]. Indeed, bats can actually adjust their echolocation under jamming conditions [34,35]. Bats are, however, likely ''uncomfortable'' when broadband ultrasound is present because it forces them to shift their call frequencies to avoid overlap, which in turn will lead to suboptimal use of echolocation or they may not echolocate at all [14,34].
In contrast to previously tested acoustic ''repellers'' [36], the device we have developed and tested shows some promise for deterring bats from the surrounding airspace near wind turbines. This study represents the first field test of a deterrent device to reduce bat fatalities at wind turbines by comparing fatalities at treated and untreated turbines. Our findings generally corroborate our previous conclusions from unpublished laboratory and field experiments that a regime of presumably uncomfortable or disorienting ultrasound may deter bats from occupying such a treated airspace. While the treatment response we observed generally falls within the range of variation of fatalities among turbines we studied, nothing in the statistical evaluation of our data suggested that our random selection of the 10 treatment turbines somehow skewed mortality rates among the turbines we chose. We acknowledge that 3 of our treatment turbines had to be located on the Locust Ridge I portion of the facility where no control turbines were selected. While this could have influenced the results, we noted in 2009 that two of these three turbines had fewer mean fatalities relative to the overall mean for deterrent turbines, while in 2010, the mean fatalities of all three of these turbines generally were equal to or greater than the overall mean for treatment turbines. Fatalities at other turbines in both the control and treatment sets also varied from one year to the next and we do not believe data from the three turbines from LR I biased our findings.
In 2010, we examined potential inherent difference between the two sets of turbines and our findings suggested a marginal difference existed in fatalities between control and treatment turbines prior to implementation of the treatment. However, we caution that data from our pre-treatment period in 2010 was collected prior to migration of migratory tree roosting species and the ratio of migrant to non-migrant species was different between these two periods in our study. Thus, different levels of fatality, different species composition, and possibly different behaviors of the bats during the two phases may have influenced our findings regarding inherent differences between control and treatment turbines. Future field tests of deterrent devices should better account for potential differences in fatalities among different species when determining inherent variation among sample turbines.
The effectiveness of ultrasonic deterrents as a means to prevent bat fatalities at wind turbines is limited by the distance and area that ultrasound can be broadcast. Unfortunately, rapid attenuation of ultrasound in air, which is heavily influenced by humidity (Table S1), limits the effective range of broadcasts. Nightly humidity in this region of Pennsylvania averaged 86.5% in August  ). Ultrasound emission in the perpendicular plane of the rotor-swept area may be adequate to affect approaching bats, particularly those species influenced at the lower frequencies. However, it is clear that effective emissions in the parallel plane of the rotor-swept area will be difficult if not impossible to achieve based on sound attenuation in humid environments. The effective airspace would be different and larger in more arid environments, however (Table  S1). We also note that some devices were not operating all the time during our study, due to malfunctions. Although we were unable to account for this factor in our analysis, clearly the affected airspace was reduced when some devices were inactive, which further influenced our findings. We assume that when bats encounter a gradient of increasingly strong emissions as they approach the deterrent device, they will respond by flying opposite to that gradient to escape the effect of emissions. However, at present we know little about the general responses that various species of bats have upon entering a field of ultrasound emissions. It is therefore important to consider our assumptions when interpreting results of this study. Although our acoustic deterrent device could only generate a limited effective volume of uncomfortable airspace, bats could have detected the presence of such airspace from a greater range, possibly beyond the rotor swept area. Bats previously experiencing the discomfort of ultrasound broadcast may avoid approaching other treated towers, which they could detect as treated from beyond the zone of discomfort. In this way, ultrasound broadcast may effectively serve as acoustic beacons to direct bats away from wind turbines. Over time, bats may learn to avoid all turbines from their experience with those equipped with deterrents, similar to documented behavior of bats encountering other discomforting experiences such as mist nets. Just as bat capture success in mist nets declines on successive nights as bats apparently learn to detect the presence of nets and thereafter avoid them [37,38,39,40], we speculate that after experiencing a disagreeable encounter with ultrasound treated airspace bats may opt to subsequently avoid it. Other lines of evidence indicate that bats learn and remember spatial locations or stimuli associated with obstacles or threats. A study that modified experiments conducted by Griffin [14] challenged bats to maneuver through vertical wires, and they did so by tilting and scrunching their wings; the same bats continued these maneuvers at the locations of wires even after they were removed [41]. In practice, the actual decline of activity at any treated site will likely depend upon immigration of naïve bats into the area. We did not monitor bat activity with night vision or thermal imaging cameras [42] and, thus, were unable to assess activity patterns of bats simultaneous with fatality searches. It is also possible that insects preyed on by bats in this region were deterred from the turbines, which could represent the ultimate cause of avoiding treated turbines. Indeed, studies have demonstrated that ultrasound can repel insects [43] and influence their reproduction [44]. However, we did not assess insect abundance and suggest future studies should attempt to address causal factors of avoidance including effect on insect prey.
Conversely, bats may habituate to the presence of ultrasound emissions and acoustic deterrents may actually lose their effectiveness over time. However, in prior field tests of deterrents, bats did not appear to habituate or accommodate to the presence of ultrasound emitted from a previous prototype deterrent at least over short periods of 5-7 days (J. Szewczak, Humboldt State University and E. Arnett, Bat Conservation International, unpublished data). Habituation to deterrents should, however, continue to be investigated in future studies.
The effectiveness of acoustic deterrents will likely vary among different species of bats. Hoary bats, for example, employ the lowest frequency range of the species we studied (,20-25 kHz) and may be affected more so than other species that use higher frequencies and perhaps fly at further distances from the device. Hoary bats had significantly fewer fatalities at turbines with deterrents relative to those without them in both years, and silverhaired bats also had fewer fatalities at turbines with deterrents in 2010. In 2010, however, after accounting for inherent differences between turbine sets prior to treatment, hoary and silver-haired bats killed per control turbine relative to treatment turbines during the treatment period represented only a 20% increase in fatality over the pre-treatment period. Species-specific effectiveness warrants further investigation in a study with more power to detect differences among species. Future studies hopefully will also elucidate whether deterrents can eventually serve as a mitigation tool for minimizing or eliminating take of threatened or endangered species such as the Indiana bat (Myotis sodalis). The limited range of ultrasound broadcast from a wind turbine tower or nacelle might have only a moderate contribution toward reducing impacts of bats randomly flying through the rotor-swept area. However, for bats that may be drawn to and approach turbine towers as potential roosts or gathering sites [1,45], the combination of effective range and learned avoidance response to ultrasound broadcast may have longer term effects in reducing bat mortality at wind turbines. We also note that we only tested broadband ultrasound emission (20-100 kHz) and short pulses mimicking echoes of insects [16], for example, could prove to be more effective for some species and should be tested in future studies. Introducing ultrasound emissions into the environment could potentially yield negative environmental effects on other species of wildlife, but we do not feel this is of concern because the device we tested only had a limited effective range because of rapid attenuation of ultrasound with distance from its source. Within the effective range of the treated airspace, emissions could affect ultrasound-sensitive ) insects and disperse them, providing less reason for bats to occupy that airspace, assuming food sources attract bats to turbines [1,10]; and 2) passerines that may be deterred from the airspace, thus reducing strikes. If ultrasound could reach the ground from locations where deterrent devices are mounted, which was not the case in our study, ultrasound sensitive small mammals could be dispersed away from the base of the turbines, but this would likely be a positive impact that could reduce strikes of stooping raptors.

Conclusions
This study, and previous experiments with earlier prototypes, revealed that broadband ultrasound broadcasts may affect bat behavior directly by discouraging them from approaching the sound source, or indirectly by reducing the time bats spend foraging near a turbine if insects are repelled by ultrasound [e.g., 42, 43; also recognizing not all insects have ears to detect ultrasound] and ultimately reduce bat fatalities at wind turbines. However, variation among turbines yielded inconclusive evidence of a strong effect of deterrents on bat fatality and while the approach may hold some promise, further refinement and investigation is needed. We did experience technical issues in both years of the study, including water leakage that rendered some deterrents inoperable during portions of the study period which clearly influenced our findings. Thus, results from this study may reflect a more conservative estimate of potential fatality reduction achievable through application of the deterrent device we tested. Still, we caution that the response estimated in this study falls generally within the range of variation for bat fatalities among turbines in this and other studies in the region [2]. Additionally, deterrents resulted in lower reductions in bat fatality relative to curtailing turbine operations by increasing cut-in speeds (44-93%) [11,12]. We further caution that it would be premature and unwarranted to conclude or interpret from these initial results that this technology provides an operational deterrent device ready for broad-scale deployment at wind facilities. While we do not consider acoustic deterrents to be an acceptable mitigation strategy at this time, with further experimentation and modifications, this type of deterrent method may prove successful and broadly applicable for protecting bats from harmful encounters with wind turbine blades. Future research and development and field studies should attempt to improve the device and it's weatherproofing and emission performance, and optimizes the placement and number of devices on each turbine that would affect the greatest amount of airspace in the rotor-swept area to estimate potential maximum effectiveness of this tool to reduce bat fatalities. New studies also should test other emission types such as short ultrasonic pulses that mimic insects [16]. Finally, we did not attempt to develop comparative estimates of costs associated with our deterrent devices relative to lost revue of operational mitigation because current deterrent development costs are high and dynamic and operational costs to maintain them over a period of time have not been established. Future efforts should determine production and maintenance costs of newly developed deterrents that can be factored into a costbenefit analysis comparing different approaches for mitigating bat fatalities. Table S1 Calculated decibel level at different distances and frequencies at two different levels of relative humidity (10 and 40%) for acoustic deterrent devices used in this study. Calculations assume ambient temperature of 20uC and air pressure of 101.325 kPa (kilopascal). (DOCX) Table S2 The attenuation of sound in air due to viscous, thermal and rotational loss mechanisms is simply proportional to f 2 . However, losses due to vibrational relaxation of oxygen molecules are generally much greater than those due to the classical processes, and the attenuation of sound varies significantly with temperature, water-vapor content and frequency. A method for calculating the absorption at a given temperature, humidity, and pressure can be found in ISO 9613-1 (1993). The table and figure below gives values of attenuation in dB m 21 for a temperature of 20uC and an air pressure of 101.325 kPa. The uncertainty is estimated to be 610%. (DOCX) Table S3 Habitat visibility classes used during this study, following Pennsylvania Game Commission Protocol [21]. Data for Classes 3 and 4 were combined during our final analyses. (DOCX) this work. Finally, we wish to thank M. Jensen (Binary Acoustic Technology, Tucson, Arizona) for his initial work on the deterrent devices and Deaton Engineering (Georgetown, Texas) for their follow-up design work and manufacturing of the deterrents used in this study.