The selection pressure from echolocating bats has driven the development of a diverse range of anti-bat strategies in insects. For instance, several studies have proposed that the wings of some moths absorb a large portion of the sound energy contained in a bat's ultrasonic cry; as a result, the bat receives a dampened echo, and the moth becomes invisible to the bat. To test the hypothesis that greater exposure to bat predation drives the development of higher ultrasound absorbance, we used a small reverberation chamber to measure the ultrasound absorbance of the wings of nocturnal (Bombycoidea: Saturniidae) and diurnal moths (Chalcosiinae: Zygaenoidea: Zygaenidae). The absorption factor of the nocturnal saturniids peaks significantly higher than the absorption factor of the diurnal chalcosiines. However, the wings of the chalcosiines absorb more ultrasound than the wings of some diurnal butterflies. Following a phylogenetic analysis on the character state of diurnality/ nocturnality in the Zygaenidae, we propose that diurnality in the Chalcosiinae is plesiomorphic (retained); hence, the absorbance of their wings is probably not a vestigial trait from an ancestral, nocturnal form but an adaptation to bat activity that overlaps their own. On a within-species level, females of the saturniids Argema mittrei and Samia cynthia ricini have significantly higher absorption factors than the males. In the female S. c. ricini, the higher absorption factor corresponds to a detection distance by bats that is at best 20-30% shorter than that of the male.

The prey-predator interaction between moths and bats has served for decades as an example of a coevolutionary arms race. Moths have developed a series of measures to avoid detection from the echolocation cries of bats (primary defences) and to promote survival once detected (secondary defences). Of all defence mechanisms, ultrasonic hearing is probably the one that has received the most attention (see reviews from Yack et al., 1999; Minet and Surlykke, 2003). Studies have reported at least three response mechanisms of eared moths exposed to echolocation cries or artificial ultrasound. The first is the evasive action of flying moths, which increases their chance of survival 40-50% over sympatric earless species (Acharya and Fenton, 1999; Roeder, 1998). The second is the adaptive silence of the non-aerial male lesser wax moth (Achroia grisella) that ceases its courtship song (Greenfield and Baker, 2003; Spangler, 1984). The third is the clicking sounds of tiger moths (Noctuoidea: Erebidae) that function as acoustic aposematism, Batesian or Müllerian mimicry, startle signals, or echolocation jamming (Conner, 2014; Conner and Corcoran, 2012; Miller and Surlykke, 2001).

Hearing clearly offers a survival advantage; however, the auditory system comes with energetic and behavioural costs and when it is no longer beneficial, such as in some moths that escape bat predation, it regresses (Miller and Surlykke, 2001). Several studies have measured directly the high energetic expenditure on the neural tissue of visual sensory systems (Niven and Laughlin, 2008), and although there are no equivalent direct data for auditory systems, the energy cost hypothesis is applicable throughout the nervous system (Niven and Laughlin, 2008). The potential behavioural costs for an eared moth emanate from the evasive response, which can result in the moth landing on a water surface (Guignion and Fullard, 2004) or missing a mating opportunity (Yager, 2012). These could be reasons why some moths under bat predation are deaf, yet rely on non-auditory anti-bat measures, such as their flight behaviour. Earless moths tend to fly less (Morrill and Fullard, 1992; Soutar and Fullard, 2004) or erratically (Lewis et al., 1993), and achieve acoustic crypsis by remaining close above vegetation in order to blend into the echoes of surrounding clutter (Fullard and Napoleone, 2001; Lewis et al., 1993; Rydell, 1998). Other earless species avoid bats more drastically: they emerge during seasons or times of the day when bat activity is low (Lewis et al., 1993; Soutar and Fullard, 2004; Yack, 1988).

Recent research revealed how some moths utilize the morphology of their wings towards two non-auditory defences: acoustic deflection and ultrasound absorption. The first is the acoustic equivalent of the visual deflection that several Lepidoptera species employ to deflect visually guided attacks from their body towards eyespots on their wings (Brakefield and Larsen, 1984; Lyytinen et al., 2003; Olofsson et al., 2010; Stevens, 2005; Vallin et al., 2011). Accordingly, luna moths (Actias luna) deflect bat attacks towards their long hindwing tails. Tailed moths have a ∼47% survival advantage over moths that have their hindwing tails ablated, an advantage similar to the one of eared over earless moths (Barber et al., 2015). In addition to diverting attacks, moth wings may absorb a large portion of the energy of an incidental echolocation cry. As a result, the intensity of the reflected echo and the detection distance of the moth decline. Several studies have proposed that the scales on moth wings may increase the absorbance of bat echolocation cries (Moss and Zagaeski, 1994; Roeder, 1963; Zeng et al., 2011). Zeng et al. (2011) designed a novel method for the measurement of absorbance in the ultrasonic spectrum and compared the absorption factors (measured absorption quantity) between the wings of nocturnal moths and diurnal butterflies. The latter absorbed significantly less ultrasound in the frequency range 40-60 kHz, presumably because they are not exposed to bat predation.

Accordingly, we compare the absorption factors between the wings of members of the Saturniidae (Bombycoidea) and of the Chalcosiinae (Zygaenoidea: Zygaenidae). Several studies have reported on the nocturnal activity of the saturniid study species Argema mittrei (Guérin-Méneville 1847) (Jolly et al., 1984), Automeris io (Fabricius 1775) (Fullard and Napoleone, 2001), and Samia cynthia ricini (Jones 1791) (Eguchi et al., 1982). Although there are no studies on the behaviour of the chalcosiines in this study [Amesia aliris analis (Jordan 1907), Campylotes burmanus (Hampson 1919), and Erasmia pulchera (Hope 1840)], the chalcosiines are primarily diurnal (Subchev, 2014; Yen et al., 2005). The purpose of this comparative study is to test the hypothesis that the between-families and between-species absorption factors of wings differ significantly due to different degrees of exposure to bat predation. Moreover, the saturniid study species are sexually dimorphic, and so the within-species absorption factors of males and females are also compared.

Measurement of moth wing absorption factors in a small reverberation chamber

Dried specimens and a small reverberation chamber (Fig. 1) were used to measure the random incidence absorption factors of the wings of each study species (separately for male and female saturniids) over the frequency range 20–100 kHz (Fig. 2). For each treatment n=6 sets were used, with each set comprising the maximum number of non-overlapping moth wings of same-sex specimens of a species that could fit on the reverberation chamber floor. The number of wings of each set depended on species size and availability of dried specimens.

Fig. 1.

Small reverberation chamber for the measurement of absorption factors in the ultrasonic spectrum. (A) Experimental set-up. Funnels direct the ultrasound emitted by the transducers into the chamber through a 7 mm opening. The microphone can move freely up and down inside the chamber, allowing for spatial averaging of the reverberation times. (B) Overview of reverberation chamber showing the engraved diffusers. (C) Absorption factor of empty reverberation chamber (Eqn 1); the error bars indicate 95% confidence intervals (CIs) (n=27 microphone positions).

Fig. 1.

Small reverberation chamber for the measurement of absorption factors in the ultrasonic spectrum. (A) Experimental set-up. Funnels direct the ultrasound emitted by the transducers into the chamber through a 7 mm opening. The microphone can move freely up and down inside the chamber, allowing for spatial averaging of the reverberation times. (B) Overview of reverberation chamber showing the engraved diffusers. (C) Absorption factor of empty reverberation chamber (Eqn 1); the error bars indicate 95% confidence intervals (CIs) (n=27 microphone positions).

Fig. 2.

Absorption factors of moth wings measured with small reverberation chamber (Eqn 3). (A) Chalcosiines. (B-D) Female and male saturniids. The absorption factors of A. mittrei and S. c. ricini differ significantly according to sex for the frequencies indicated by circles. (E) Absorption factors grouped by preparation group. The absorption factors differ significantly from each other for frequencies within the range 20–40 kHz. Circles indicate significant difference between male and female saturniids, triangles between female saturniids and chalcosiines, and squares between male saturniids and chalcosiines. All pairwise comparisons were conducted with Tukey-Kramer test for each frequency separately. One, two, or three symbols correspond to 0.01<P<0.05, 0.001<P<0.01, and P<0.001, respectively. The error bars indicate the 95% CIs and n=6 sets for each sample.

Fig. 2.

Absorption factors of moth wings measured with small reverberation chamber (Eqn 3). (A) Chalcosiines. (B-D) Female and male saturniids. The absorption factors of A. mittrei and S. c. ricini differ significantly according to sex for the frequencies indicated by circles. (E) Absorption factors grouped by preparation group. The absorption factors differ significantly from each other for frequencies within the range 20–40 kHz. Circles indicate significant difference between male and female saturniids, triangles between female saturniids and chalcosiines, and squares between male saturniids and chalcosiines. All pairwise comparisons were conducted with Tukey-Kramer test for each frequency separately. One, two, or three symbols correspond to 0.01<P<0.05, 0.001<P<0.01, and P<0.001, respectively. The error bars indicate the 95% CIs and n=6 sets for each sample.

The absorption factors of all preparations follow similar patterns, with a peak occurring within the frequency range 20–25 kHz (Fig. 2A-D). However, there is a significant effect of preparation on absorption factor, meaning that the estimated marginal means across the frequency range 20–100 kHz differ significantly according to preparation (ANOVA, F8,45=23.78, P<0.001; Fig. 3A). Apart from preparation, frequency also has a significant effect on absorption factor, signifying that the absorption factor is frequency dependant (repeated measured ANOVA, F7,315=129.13, P<0.001; P-value computed with Greenhouse-Geisser correction). Furthermore, there is a significant preparation-frequency interaction indicating that different preparations have different absorption factors at different frequencies (repeated measures ANOVA, F56,315=3.87, P<0.001; P-value computed with Greenhouse-Geisser correction). Pairwise comparisons showed that the absorption factors of the chalcosiines do not differ significantly from each other (Tukey-Kramer test, A. a. analisC. burmanus, P=0.842; A. a. analisE. pulchera, P=0.823; C. burmanusE. pulchera, P=0.051). However, when it comes to within-species comparisons in saturniids, it appears that the absorption mechanisms of A. mittrei and S. c. ricini are sexually dimorphic; although the absorption factors between male and female A. io do not differ significantly (Tukey-Kramer test, P<0.05, P<0.001, and P=0.99 respectively). Specifically, the between-sexes absorption factors of A. mittrei and S. c. ricini differ significantly for frequencies within the range 20–40 kHz (Tukey-Kramer test for each frequency separately, see Fig. 2B and D for significant P-values), where females have higher absorption factors.

Fig. 3.

Estimated marginal means of absorption factors of moth wings across the frequency range 20–100 kHz. (A) Estimated marginal mean of each preparation. There is a significant effect of preparation on estimated marginal mean (ANOVA, F8,45=23.78, P<0.001). (B) Estimated marginal means grouped by preparation group. The estimated marginal means differ significantly according to preparation group (ANOVA, F2,51=19.08, P<0.001) and also differ significantly from each other (Tukey-Kramer test). The error bars show the 95% CIs [n=48 in A (6 sets×8 frequencies) and n=144 in B (3 preparations×6 sets×8 frequencies)].

Fig. 3.

Estimated marginal means of absorption factors of moth wings across the frequency range 20–100 kHz. (A) Estimated marginal mean of each preparation. There is a significant effect of preparation on estimated marginal mean (ANOVA, F8,45=23.78, P<0.001). (B) Estimated marginal means grouped by preparation group. The estimated marginal means differ significantly according to preparation group (ANOVA, F2,51=19.08, P<0.001) and also differ significantly from each other (Tukey-Kramer test). The error bars show the 95% CIs [n=48 in A (6 sets×8 frequencies) and n=144 in B (3 preparations×6 sets×8 frequencies)].

To compare by preparation group, the absorption factors of each group were calculated (Fig. 2E). The estimated marginal means of the absorption factors across the frequency range 20–100 kHz differ significantly according to preparation group (ANOVA, F2,51=19.08, P<0.001; Fig. 3B), and pairwise comparisons show that they significantly differ from each other (Tukey-Kramer test, Chalcosiinae–female Saturniidae, P<0.001; Chalcosiinae–male Saturniidae, P<0.05; female Saturniidae–male Saturniidae, P<0.001). The female saturniids have the highest estimated marginal mean, followed by the chalcosiines and the male saturniids. Apart from the significant group effect on absorption factor, there is a significant effect of frequency on absorption factor as well as a significant group-frequency interaction (repeated measures ANOVA, F7,357=123.13, P<0.001 and F14,357=7.46, P<0.001 respectively; P-values computed with Greenhouse-Geisser correction). Pairwise comparisons grouped by frequency show that the differences of each pair of preparation groups lie within frequency range 20–40 kHz (Tukey-Kramer test for each frequency separately, see Fig. 2E for significant P-values). The female saturniids generally have higher absorption factors than the other two preparation groups for this frequency range, and while the male saturniids have a significantly higher peak than the chalcosiines, their absorption factors decline more steeply, hence the abovementioned higher marginal mean of the chalcosiines.

Survival advantage of female S. c. ricini over males in terms of smaller detection distance by bats

The detection distance of a target by a bat depends on the transmission loss due to spherical spreading and atmospheric attenuation, and on the target strength (Møhl, 1988). The detection distances of male and female S. c. ricini were compared because they have wings of similar size and shape; consequently the comparison is feasible because their target strengths, and hence their detection distances, differ due to the different absorption factors of their wings alone. Detection distances of moths typically range between 1–10 m (Surlykke et al., 1999). Eqn 13 was used with three hypothetical detection distances within this range for male S. c. ricini in order to derive the respective detection distances of female S. c. ricini (Fig. 4). At 25 kHz, where the difference between the absorption factors of male and female S. c. ricini is maximized (Fig. 2D), the detection distance of female S. c. ricini is 20–30% smaller.

Fig. 4.

Percentage difference in the detection distance of female S. c. ricini compared to three hypothetical detection distances of male S. c. ricini [r(2)]. The difference between the detection distances of male and female S. c. ricini is maximized at 25 kHz, where the difference between the absorption factors of their wings is also maximized (Fig. 2D). The detection distances of female S. c. ricini were calculated with Eqn 13 for atmospheric attenuation at temperature 20°C and 70% relative humidity (Surlykke, 1988).

Fig. 4.

Percentage difference in the detection distance of female S. c. ricini compared to three hypothetical detection distances of male S. c. ricini [r(2)]. The difference between the detection distances of male and female S. c. ricini is maximized at 25 kHz, where the difference between the absorption factors of their wings is also maximized (Fig. 2D). The detection distances of female S. c. ricini were calculated with Eqn 13 for atmospheric attenuation at temperature 20°C and 70% relative humidity (Surlykke, 1988).

Studies have reported bat echolocation calls with dominant frequencies ranging from 11 kHz (Euderma maculatum; Fullard and Dawson, 1997) to 212 kHz (Cloeotis percivali; Fenton and Bell, 1981), though most insectivorous bats echolocate with dominant frequency 20–60 kHz (Fenton et al., 1998). The ultrasound absorbance of the wings of this study's moth species peaks at the lower end of this range (20–25 kHz; Fig. 2). Despite the similar patterns, there are significant differences not only between the nocturnal saturniids and the diurnal chalcosiines (Fig. 2E), but also between male and female A. mittrei (Fig. 2B) and S. c. ricini (Fig. 2D).

Zeng et al. (2011) proposed that the scales on the wings are responsible for the ultrasound absorption. The scales have interstitial spaces between them that create a network of interconnected pores similar to that found in porous sound absorbers (Zeng et al., 2011); when a sound wave propagates through this network, thermal and viscous effects cause the dissipation of its acoustic energy (Cox and D'Antonio, 2009a). In addition, the ultrastructure of the scales resembles a perforated panel backed by air (Zeng et al., 2011), which could act as a microperforated panel absorber (Cox and D'Antonio, 2009b).

The absorption of ultrasound is not the sole defensive function of the moth wing scales. For instance, a moth can release itself from a spider web by shedding some of its wing scales (Eisner et al., 1964). Furthermore, the microstructure of the scales is responsible for some of the colours found in moth wings (Brink and Lee, 1996; Ghiradella, 1991). In saturniids, the wing colours and patterns can play a defensive role by achieving crypsis, aposematism, or mimicry (Blest, 1957a,b), and in Callosamia securifera, C. promethea, and Eupackaria calleta, only the males are Batesian mimics of the unpalatable swallowtail Battus philenor and/or B. polydamas (Jeffords et al., 1979; Sternburg et al., 1977; Waldbauer and Sternburg, 1975). Such sexual dimorphism may have been driven by a more intense selection pressure on males by visual predators (Allen et al., 2011), because the males fly during the day in search of the sedentary, pheromone-releasing females (Morton, 2009).

Accordingly, we hypothesized that male saturniids will exhibit adaptive increase to bat predation compared to females, however, the female A. mittrei and S. c. ricini exhibit significantly higher absorption factors over certain frequencies (Fig. 2B and D). In the female S. c. ricini, the higher absorption factor translates to a detection distance by bats that is 20–30% shorter than the detection distance of the male, assuming a bat echolocates with a dominant frequency of 25 kHz (Fig. 4). Since the absorption mechanism serves no other function but to render the moth inconspicuous to echolocating bats, our results indicate that the female A. mittrei and S. c. ricini may be under more intense selection pressure by bats than males. Yet further studies on the behaviour of these two species are required to corroborate this hypothesis. These would not be the first examples of sexual dimorphism due to differential bat predation on the two sexes. The males of some species of moths and mantids are under heavier selection pressure from echolocating bats compared to the females, and as a result they have developed more sensitive ultrasonic hearing (moths: Cardone and Fullard, 1988; Rydell et al., 1997; mantids: Yager, 1990). However, it is possible that the absorption factors do not vary between the sexes as predicted because there are trade-offs among the functions of the wing scales. In that case, intense selection for e.g. the colour-producing properties of the wing scales (Brink and Lee, 1996; Ghiradella, 1991) of males would undermine their ultrasound absorption properties.

The second significant finding of our study regards the differences between the absorption factors of the nocturnal saturniids and of the diurnal chalcosiines. Comparative studies on moth ears that function as bat detectors have reported on the degenerate state of audition in diurnal moths compared to nocturnal ones. Such findings suggest that diurnal activity could be an apomorphic (derived) trait that has allowed some moths to escape bat predation (Fullard et al., 1997; Muma and Fullard, 2004). Accordingly, we hypothesised that if the ultrasound absorption mechanism is the result of bat predation pressure, then the diurnal chalcosiine study species should have significantly lower absorption factors compared to the nocturnal saturniids. Indeed, the absorption factor of the chalcosiines peaks at approximately 0.4, whereas the absorption factors of the male and female saturniids peak around 0.6 (Fig. 2E). However, the peak of 0.4 may be significant compared to the peaks of 0.1–0.2 reported for butterfly wings and nocturnal moth wings without scales (Zeng et al., 2011). An explanation is that diurnality is an apomorphic trait in chalcosiines; hence, their absorption mechanism is a case of evolution in reverse (Porter and Crandall, 2003). This means that the absorption mechanism is a vestigial trait that the chalcosiines have retained from their ancestral, nocturnal form. The peak absorption factor that is higher than in butterflies but lower than in nocturnal moths may signify that the absorption mechanism has degenerated, but not enough evolutionary time has passed for it to disappear completely. It is possible however that the difference in peak absorbance between the diurnal chalcosiines and the nocturnal saturniids is not due to degeneration of the absorption mechanism, but due to an adaptive increase of the latter to the ongoing selection pressure from bats. Consequently, the chalcosiines could be a window to the absorption mechanism in its primitive form. In cases such as this one, where a vestigial trait persists, there is probably no intense selection pressure against the trait, which may eventually disappear by genetic drift. A relevant example is the fully functioning or degenerate hearing of some moths that are completely isolated from bats, spatially or temporally (review: Miller and Surlykke, 2001).

The hypothesis that diurnality in chalcosiines is apomorphic was tested with phylogenetic analysis, by mapping the character state of nocturnality/ diurnality on the phylogenetic tree of the Zygaenidae (Niehuis et al., 2006). Of the four zygaenid subfamilies, the Zygaeninae, the Procridinae, and the Chalcosiinae are predominately diurnal, at least when it comes to mating, but there are not enough data to make a conclusion regarding the Callizygaeninae (Subchev, 2014). Since three out of the four zygaenid subfamilies are diurnal, it is highly probable that diurnality is a pleisiomorphic (ancestral), not apomorphic, trait in the Chalcosiinae. Therefore, the fact that the peak absorbance of chalcosiine wings is higher than that of butterfly wings (Zeng et al., 2011) requires an alternative explanation.

Fullard et al. (2000) suggested that the terms ‘diurnal’ and ‘nocturnal’ can be misleading, especially when describing diel flight activity at a higher taxonomic level (e.g. subfamily). Diurnality and nocturnality are not mutually exclusive, and most so-called diurnal moths exhibit mixed diurnal/nocturnal flight activities (Fullard and Dawson, 1999; Fullard and Napoleone, 2001; Fullard et al., 2000), though some are exclusively diurnal (Fullard et al., 2000; Muma and Fullard, 2004). The fact that the Chalcosiinae are well adapted to diurnal activity – their mimetic and aposematic wing patterns deter predators that hunt using visual cues (Endō and Kishida, 1999; Owada and Tinh, 2002) – does not preclude them from potential nocturnal activity, and thus exposure to bat predation. Ultimately, considering the phylogenetic analysis, it is more likely for the absorption mechanism of the chalcosiines to be an adaptive measure to counter occasional bat encounters than to be an apomorphic, vestigial trait. On the other hand, butterflies tend to be exclusively diurnal (Fullard and Napoleone, 2001; Fullard et al., 2000), a fact that is reflected by the very low absorption factors of their wings (Zeng et al., 2011).

Dried specimens

Male and female saturniid specimens (A. io from USA, A. mittrei from Madagascar, and S. c. ricini from China) and male chalcosiine specimens (A. a. analis and C. burmanus from Thailand, and E. pulchera from Burma) were obtained from the Lepidoptera Breeders Association (Bourne, UK).

Measurement of moth wings absorption factor with the reverberation chamber method

The time period for which it takes the sound pressure level to drop 60 dB after a sound stops is termed reverberation time and depends on the total absorption inside a reverberation chamber. Consequently, by comparing the reverberation time before and after the introduction of an absorbent material (e.g. moth wings) inside the chamber, the random incidence absorption coefficient of the introduced material can be derived. The standard for the measurement of the random incidence absorption coefficient in a reverberation chamber essentially concerns frequencies below 20 kHz (ISO, 2003). For this frequency spectrum, the measurement procedure requires large chambers and 10-12 m2 of absorbent material (Cox and D'Antonio, 2009c); however, such a quantity of moth wings is practically infeasible. Besides, the shorter wavelengths of the ultrasonic spectrum render feasible the measurement of the absorption coefficient with a smaller chamber and smaller quantities of absorbent material. Still, since this study did not follow the standard faithfully, the measured absorption quantity has been called absorption factor. The same term was employed by the only other documented study of sound absorption in the ultrasonic spectrum (Zeng et al., 2011).

A small reverberation chamber for ultrasonic applications (Fig. 1) was fabricated for the measurement of the random incidence absorption factors of the moth wings. The objective of the reverberation chamber is to create a spatially uniform acoustic field where the sound pressure level is ideally the same for every point within the chamber. In order to achieve the diffuse field, an acoustic wave should be equally probable to propagate towards any direction. For this reason, Cox and D'Antonio modified primitive root diffusers (Cox and D'Antonio, 2000), which promote diffusion of the incident wave towards all directions and reduce the number of standing waves, and these were engraved on the walls and ceiling of the chamber (Fig. 1A,B). In addition, a reverberation chamber must have hard, non-absorbing surfaces that reflect most of the acoustic energy back into the chamber. Our reverberation chamber was fabricated with Polymethyl 2-methacrylate (Perspex) that secured low ultrasound absorbance (Fig. 1C).

Even with these measures, the reverberation time is spatially dependent within the reverberation chamber. To reduce the effect of non-diffuseness, two ultrasonic transducers equipped with funnels were utilised to direct the ultrasound into the chamber through 7 mm openings, scattering multi-directionally the acoustic wave that entered the chamber. In addition, the recording device, an ultrasonic microphone, could move freely up and down inside the chamber, allowing for spatial averaging of the measured reverberation times (Fig. 1A).

The absorption factor of the empty chamber was determined by (Cox and D'Antonio, 2009c):
(1)
where α0 is the average absorption factor of the empty chamber, V (m3) is the volume of the chamber, f (Hz) is the sound frequency, T0 (s) is the reverberation time before the introduction of moth wings inside the chamber, m1 (m−1) is the air volume absorption coefficient, c (m s−1) is the speed of sound in the air, and S (m2) is the surface area of the chamber.
The absorption factors of the empty chamber and of the moth wings were measured over nine 1/3 octave bands that have centre frequencies from 20 to 100 kHz. The parameter m1 was calculated by:
(2)
where fl and fh (Hz) are the low and high frequency limits of the 1/3 octave band, μ (Pa s) is the dynamic viscosity of air, ρ (kg m−3) is the density of air, γ is the ratio of specific heats of air, Pr is the Prandlt number of air, Cp (J kg−1 K−1) is the heat capacity of air at constant pressure, and k (W m−1 K−1) is the thermal conductivity of air (Table S1).
The absorption factors of the moth wings were determined by (Cox and D'Antonio, 2009c):
(3)
where αw is the absorption factor of the moths wings introduced inside the chamber, Sw (m2) is their surface area, and T1 (s) is the reverberation time after the introduction of the moth wings.

For the measurement of a reverberation time with respect to frequency f, a function generator (Agilent 33220A) generated a sinusoidal signal of frequency f that was amplified with an amplifier (Ultrasound Advice S55) before transmission inside the chamber by two ultrasonic transducers (Ultrasound Advice S56); the transmission stopped when the response inside the chamber reached the steady state. The signal was recorded with a microphone (Bruel & Kjaer 4138) using a sampling frequency of 625 kHz and was amplified with a pre-amplifier (Nexus 2690). The reverberation times were derived from the steady state decay curves. To determine α0(f) (Fig. 1C), T0(f) was measured from n=27 microphone positions for each 1/3 octave band with centre frequency f. To determine αw(f) (Fig. 2), n=6 sets per treatment were used and T1(f) of each set was averaged from four microphone positions. The term set refers to each set-up of non-overlapping moth wings of same-sex specimens of a species that cover the maximum possible surface area of the reverberation chamber floor. The number of wings of each set depended on species size and availability of dried specimens.

Estimation of moth wings surface area

Sw was estimated with an image processing method (Fig. S1). The wings of each set were placed on a white background of known surface area and an image was obtained. The image was smoothened with edge-preserving bilateral filter (Tomasi and Manduchi, 1998) and was subsequently converted to binary via k-means++ clustering (Arthur and Vassilvitskii, 2007). In the binary image, one class of pixels represents the wings and the other represents the background. Sw was estimated with
(4)
where Nw is the number of pixels representing moth wings, N the total number of pixels, and Sb the surface area of the white background.

Derivation of reverberation time from steady state decay curve

First, the recorded signal of main frequency f was filtered with a bandpass digital filter designed with a Kaiser window. The filter passed frequencies between fl and fh in order to retain only frequencies within the 1/3 octave band, frequencies outside this range were attenuated 80 dB (Fig. S2A).

The decay curve of the filtered signal has many fluctuations that render difficult the identification of the point where the response drops 60 dB after the offset of sound; therefore, the curve has to be smoothed. The first step of the smoothing process is to obtain the filtered signal's envelope by using the Hilbert transform:
(5)
where t (s) is the time, s (V) is the signal, s^ (V) is its Hilbert transform, and E is its envelope. The envelope was further smoothed with a 625-samples long moving average filter that corresponds to a time interval of 1 ms (Fig. S2B).
The steady-state reverberation time was derived from the smoothed curve by Schroeder's integration (Schroeder, 1965):
(6)
where L [dB] is Schroeder's curve, Es (V) is the smoothed envelope, τoff (s) is the time point where the sound stops, and τN (s) is the time point where the signal merges with the noise level. The parameter τN was estimated with the method proposed by Lundeby et al. (1995).
L is a maximally flat curve, meaning that its gradient is constantly negative. As a result, the curve crosses the −60 dB point once. This property makes Schroeder's integration method useful for the estimation of the reverberation time. However, in most cases L(τN)> –60[dB]. This means that the signal merges with the noise level before reaching the −60 dB point. Therefore, the reverberation time has to be estimated using linear regression. L is approximated with the linear model, (t)=A·t+B, and the reverberation time can then be estimated with
(7)
where T60 is the reverberation time (Fig. S2C).

Statistical analysis

To compare the absorption factors among the nine preparations (3 chalcosiines and 3×2 male and female saturniids; Fig. 2A-D and Fig. 3A), a repeated measures model, specifically a subject-by-treatment model (Doncaster and Davey, 2007), was fitted. In the model, the absorption factors are the responses, preparation is the between-subjects factor, which is used as the predictor variable, and frequency is the within-subject factor. ANOVA was used to test if the absorption factors differ significantly according to preparation (Fig. 3A), and repeated measures ANOVA to test if there is a significant effect of frequency on absorption factor as well as significant preparation-frequency interaction. The P-values of the repeated measures ANOVA were computed using a Greenhouse-Geisser correction (ε=0.79) because Mauchly's test for sphericity indicated that that the sphericity, hence the compound symmetry assumption, does not hold (χ2=144.80, P<0.001). To do pairwise comparisons between preparations, post hoc Tukey-Kramer tests were used (Fig. 2A-D and Fig. 3A).

Accordingly, a repeated measures model was used for the comparisons among the three preparation groups (chalcosiines, male and female saturniids; Fig. 2E and Fig. 3B), albeit with the preparation group as the between-subjects factor. The same tests were carried out as with the above model. Again, the sphericity did not hold (χ2=142.61, P<0.001), and a Greenhouse-Geisser correction (ε=0.81) was used for the calculation of the repeated measures ANOVA P-values. All statistical analysis was conducted in MATLAB (Mathworks, UK), and a P-value of <0.05 was considered statistically significant.

Survival advantage of female S. c. ricini over males in terms of smaller detection distance by bats

The relationship between the detection distances of male and female S. c. ricini was derived using the sonar equation (Møhl, 1988):
(8)
where D is the auditory detection threshold of the bat, H is the source level, K is the one-way transmission loss, M is the target strength, and e is the noise. The equation is in dB form and all quantities are measured in dB.
K is given by:
(9)
where r (m) is the detection distance of the target (e.g. moth), and λ (dB m−1) is the atmospheric attenuation factor. The first term of the equation accounts for loss due to spherical spreading and the second one for loss due to atmospheric attenuation.
M is given by:
(10)
where Ir and Ii (W m−2) are the returned and incident sound intensities respectively.
The ratio of incident to returned sound intensity and the absorption factor of the moth wings are related as follows (Cox and D'Antonio, 2009c):
(11)
Since male and female S. c. ricini have wings of similar size and shape, their ratios of incident to returned sound intensity depend on the absorption factors of their wings alone. Specifically, they are related as follows:
(12)
where the superscripts (1) and (2) denote quantities for female and male S. c. ricini, respectively.
Expressing Eqn 8 for female S. c. ricini and subtracting the respective equation for male S. c. ricini, then solving the resulting equation with respect to r(1), gives:
(13)
where W() is the Lambert W function (Corless et al., 1993, 1996), which in this study returns exactly one real solution. Eqn 13 was used to compare the detection distances of female and male S. c. ricini over bat sonar frequency range 20–100 kHz (Fig. 4).

We thank colleagues in the laboratories of the Centre of Ultrasonic Engineering at the University of Strathclyde for their advice and assistance during this work.

Author contributions

A.N. conceived the study, which was designed by A.N., F.G. and J.F.C.W. Experiments were carried out by A.N. and F.G.; data analysis was performed by A.N. A.N., F.G. and J.F.C.W. interpreted the findings. A.N. drafted the article. All authors read and approved the final version of the manuscript.

Funding

The research leading to these results has received funding from the European Research Council under the European Union's Seventh Framework Programme (FP/2007-2013)/ERC grant Agreement no. 615030; and was also supported by the Engineering and Physical Sciences Research Council (EPSRC) grant EP/L026511/1.

Data availability

All data created during this research are openly available from the University of Strathclyde Pure/KnowledgeBase at http://dx.doi.org/10.15129/a1acdd88-f8aa-4ac4-8d6c-c3de5b672773.

Acharya
,
L.
and
Fenton
,
M. B.
(
1999
).
Bat attacks and moth defensive behaviour around street lights
.
Can. J. Zool.
77
,
27
-
33
.
Allen
,
C. E.
,
Zwaan
,
B. J.
and
Brakefield
,
P. M.
(
2011
).
Evolution of sexual dimorphism in the Lepidoptera
.
Annu. Rev. Entomol.
56
,
445
-
464
.
Arthur
,
D.
and
Vassilvitskii
,
S.
(
2007
).
k-means++: the advantages of careful seeding
. In
Proceedings of the Eighteenth Annual ACM-SIAM Symposium on Discrete Algorithms
, pp.
1027
-
1035
.
Philadelphia, PA
,
USA
:
Society for Industrial and Applied Mathematics
.
Barber
,
J. R.
,
Leavell
,
B. C.
,
Keener
,
A. L.
,
Breinholt
,
J. W.
,
Chadwell
,
B. A.
,
McClure
,
C. J. W.
,
Hill
,
G. M.
and
Kawahara
,
A. Y.
(
2015
).
Moth tails divert bat attack: evolution of acoustic deflection
.
Proc. Natl. Acad. Sci. USA
112
,
2812
-
2816
.
Blest
,
A. D.
(
1957a
).
The function of eyespot patterns in the Lepidoptera
.
Behaviour
11
,
209
-
258
.
Blest
,
A. D.
(
1957b
).
The evolution of protective displays in the saturnioidea and sphingidae (Lepidoptera)
.
Behaviour
11
,
257
-
309
.
Brakefield
,
P. M.
and
Larsen
,
T. B.
(
1984
).
The evolutionary significance of wet and dry season forms in some tropical butterflies
.
Biol. J. Linn. Soc.
22
,
1
-
12
.
Brink
,
D. J.
and
Lee
,
M. E.
(
1996
).
Ellipsometry of diffractive insect reflectors
.
Appl. Opt.
35
,
1950
-
1955
.
Cardone
,
B.
and
Fullard
,
J. H.
(
1988
).
Auditory characteristics and sexual dimorphism in the gypsy moth
.
Physiol. Entomol.
13
,
9
-
14
.
Conner
,
W. E.
(
2014
).
Adaptive sounds and silences: acoustic anti-predator strategies in insects
. In
Insect Hearing and Acoustic Communication
(ed.
B.
Hedwig
), pp.
65
-
80
.
New York
:
Springer
.
Conner
,
W. E.
and
Corcoran
,
A. J.
(
2012
).
Sound strategies: the 65-million-year-old battle between bats and insects
.
Annu. Rev. Entomol.
57
,
21
-
39
.
Corless
,
R. M.
,
Gonnet
,
G. H.
,
Hare
,
D. E. G.
and
Jeffrey
,
D. J.
(
1993
).
Lambert's W function in Maple
.
Maple Tech. Newsl.
9
,
12
-
22
.
Corless
,
R. M.
,
Gonnet
,
G. H.
,
Hare
,
D. E. G.
,
Jeffrey
,
D. J.
and
Knuth
,
D. E.
(
1996
).
On the LambertW function
.
Adv. Comput. Math.
5
,
329
-
359
.
Cox
,
T. J.
and
D'Antonio
,
P.
(
2000
).
Acoustic phase gratings for reduced specular reflection
.
Appl. Acoust.
60
,
167
-
186
.
Cox
,
T. J.
and
D'Antonio
,
P.
(
2009a
).
Porous absorption
. In
Acoustic Absorbers and Diffusers: Theory, Design and Application
, pp.
175
-
195
.
Abingdon
:
Taylor & Francis
.
Cox
,
T. J.
and
D'Antonio
,
P.
(
2009b
).
Resonant absorbers
. In
Acoustic Absorbers and Diffusers: Theory, Design and Application
, pp.
196
-
229
.
Abingdon
:
Taylor & Francis
.
Cox
,
T. J.
and
D'Antonio
,
P.
(
2009c
).
Measurement of absorber properties
. In
Acoustic Absorbers and Diffusers: Theory, Design and Application
, pp.
70
-
109
.
Abingdon
:
Taylor & Francis
.
Doncaster
,
C. P.
and
Davey
,
A. J. H.
(
2007
).
Repeated-measures designs
. In
Analysis of Variance and Covariance: How to Choose and Construct Models for the Life Sciences
, pp.
179
-
227
.
Cambridge
:
Cambridge University Press
.
Eguchi
,
E.
,
Watanabe
,
K.
,
Hariyama
,
T.
and
Yamamoto
,
K.
(
1982
).
A comparison of electrophysiologically determined spectral responses in 35 species of Lepidoptera
.
J. Insect Physiol.
28
,
675
-
682
.
Eisner
,
T.
,
Alsop
,
R.
and
Ettershank
,
G.
(
1964
).
Adhesiveness of spider silk
.
Science
146
,
1058
-
1061
.
Endō
,
T.
and
Kishida
,
Y.
(
1999
).
Day-flying Moths: Chalcosiinae, Epicopeia
.
Tokyo
:
Endless Science Information
.
Fenton
,
M. B.
and
Bell
,
G. P.
(
1981
).
Recognition of species of insectivorous bats by their echolocation calls
.
J. Mammal.
62
,
233
-
243
.
Fenton
,
M. B.
,
Portfors
,
C. V.
,
Rautenbach
,
I. L.
and
Waterman
,
J. M.
(
1998
).
Compromises: sound frequencies used in echolocation by aerial-feeding bats
.
Can. J. Zool.
76
,
1174
-
1182
.
Fullard
,
J. H.
and
Dawson
,
J. W.
(
1997
).
The echolocation calls of the spotted bat Euderma maculatum are relatively inaudible to moths
.
J. Exp. Biol.
200
,
129
-
137
.
Fullard
,
J. H.
and
Dawson
,
J. W.
(
1999
).
Why do diurnal moths have ears
.
Naturwissenschaften
86
,
276
-
279
.
Fullard
,
J. H.
and
Napoleone
,
N.
(
2001
).
Diel flight periodicity and the evolution of auditory defences in the Macrolepidoptera
.
Anim. Behav.
62
,
349
-
368
.
Fullard
,
J. H.
,
Dawson
,
J. W.
,
Otero
,
L. D.
and
Surlykke
,
A.
(
1997
).
Bat-deafness in day-flying moths (Lepidoptera, Notodontidae, Dioptinae)
.
J. Comp. Physiol. A
181
,
477
-
483
.
Fullard
,
J. H.
,
Otero
,
L. D.
,
Orellana
,
A.
and
Surlykke
,
A.
(
2000
).
Auditory sensitivity and diel flight activity in Neotropical Lepidoptera
.
Ann. Entomol. Soc. Am.
93
,
956
-
965
.
Ghiradella
,
H.
(
1991
).
Light and color on the wing: structural colors in butterflies and moths
.
Appl. Opt.
30
,
3492
-
3500
.
Greenfield
,
M. D.
and
Baker
,
M.
(
2003
).
Bat avoidance in non-aerial insects: the silence response of signaling males in an acoustic moth
.
Ethology
109
,
427
-
442
.
Guignion
,
C.
and
Fullard
,
J. H.
(
2004
).
A potential cost of responding to bats for moths flying over water
.
Can. J. Zool.
82
,
529
-
532
.
ISO
(
2003
).
ISO 354 Acoustics-Measurement of Sound Absorption in a Reverberation room
.
Geneva
:
ISO
.
Jeffords
,
M. R.
,
Sternburg
,
J. G.
and
Waldbauer
,
G. P.
(
1979
).
Batesian mimicry: field demonstration of the survival value of pipevine swallowtail and monarch color patterns
.
Evolution
33
,
275
-
286
.
Jolly
,
A.
,
Oberlé
,
P.
and
Albignac
,
R.
(
1984
).
Key environments: Madagascar
.
Oxford
:
Pergamon Press
.
Lewis
,
F. P.
,
Fullard
,
J. H.
and
Morrill
,
S. B.
(
1993
).
Auditory influences on the flight behaviour of moths in a Nearctic site. II. Flight times, heights, and erraticism
.
Can. J. Zool.
71
,
1562
-
1568
.
Lundeby
,
A.
,
Vigran
,
T. E.
,
Bietz
,
H.
and
Vorländer
,
M.
(
1995
).
Uncertainties of measurements in room acoustics
.
Acta Acust. united with Acust.
81
,
344
-
355
.
Lyytinen
,
A.
,
Brakefield
,
P. M.
and
Mappes
,
J.
(
2003
).
Significance of butterfly eyespots as an anti-predator device in ground-based and aerial attacks
.
Oikos
100
,
373
-
379
.
Miller
,
L. A.
and
Surlykke
,
A.
(
2001
).
How some insects detect and avoid being eaten by bats: tactics and countertactics of prey and predator
.
Bioscience
51
,
570
-
581
.
Minet
,
J.
and
Surlykke
,
A.
(
2003
).
Auditory and sound producing organs
. In
Handbook of Zoology/Handbuch der Zoologie. Lepidoptera, Moths and Butterflies, Morphology and Physiology
(ed.
N. P.
Kristensen
), pp.
289
-
323
.
Berlin
:
Walter de Gruyter
.
Møhl
,
B.
(
1988
).
Target detection by echolocating bats
. In
Animal Sonar: Processes and Performance
(ed.
P. E.
Nachtigall
and
P. W. B.
Moore
), pp.
435
-
450
.
New York
:
Springer US
.
Morrill
,
S. B.
and
Fullard
,
J. H.
(
1992
).
Auditory influences on the flight behaviour of moths in a Nearctic site. I. Flight tendency
.
Can. J. Zool.
70
,
1097
-
1101
.
Morton
,
E. S.
(
2009
).
The function of multiple mating by female promethea moths, Callosamia promethea (Drury) (Lepidoptera: Saturniidae)
.
Am. Midl. Nat.
162
,
7
-
18
.
Moss
,
C. F.
and
Zagaeski
,
M.
(
1994
).
Acoustic information available to bats using frequency-modulated sounds for the perception of insect prey
.
J. Acoust. Soc. Am.
95
,
2745
-
2756
.
Muma
,
K. E.
and
Fullard
,
J. H.
(
2004
).
Persistence and regression of hearing in the exclusively diurnal moths, Trichodezia albovittata (Geometridae) and Lycomorpha pholus (Arctiidae)
.
Ecol. Entomol.
29
,
718
-
726
.
Niehuis
,
O.
,
Yen
,
S.-H.
,
Naumann
,
C. M.
and
Misof
,
B.
(
2006
).
Higher phylogeny of zygaenid moths (Insecta: Lepidoptera) inferred from nuclear and mitochondrial sequence data and the evolution of larval cuticular cavities for chemical defence
.
Mol. Phylogenet. Evol.
39
,
812
-
829
.
Niven
,
J. E.
and
Laughlin
,
S. B.
(
2008
).
Energy limitation as a selective pressure on the evolution of sensory systems
.
J. Exp. Biol.
211
,
1792
-
1804
.
Olofsson
,
M.
,
Vallin
,
A.
,
Jakobsson
,
S.
and
Wiklund
,
C.
(
2010
).
Marginal eyespots on butterfly wings deflect bird attacks under low light intensities with UV wavelengths
.
PLoS ONE
5
,
e10798
.
Owada
,
M.
and
Tinh
,
T. H.
(
2002
).
Biodiversity and mimicry complex of diurnal insects in Vietnam
. In
Proceedings of the third and fourth Symposia on Collection Building and Natural History Studies in Asia and the Pacific Rim
(ed.
T.
Kubodera
,
H.
Masanobu
and
R.
Miyawaki
), pp.
179
-
188
.
Tokyo
:
National Science Museum
.
Porter
,
M. L.
and
Crandall
,
K. A.
(
2003
).
Lost along the way: the significance of evolution in reverse
.
Trends Ecol. Evol.
18
,
541
-
547
.
Roeder
,
K. D.
(
1963
).
Echoes of ultrasonic pulses from flying moths
.
Biol. Bull.
124
,
200
-
210
.
Roeder
,
K. D.
(
1998
).
Nerve Cells and Insect Behavior. Revised
.
Cambridge, MA
:
Harvard University Press
.
Rydell
,
J.
(
1998
).
Bat defence in lekking ghost swifts (Hepialus humuli), a moth without ultrasonic hearing
.
Proc. R. Soc. London B Biol. Sci.
265
,
1373
-
1376
.
Rydell
,
J.
,
Skals
,
N.
,
Surlykke
,
A.
and
Svensson
,
M.
(
1997
).
Hearing and bat defence in geometrid winter moths
.
Proc. R. Soc. Lond. B Biol. Sci.
264
,
83
-
88
.
Schroeder
,
M. R.
(
1965
).
New method of measuring reverberation time
.
J. Acoust. Soc. Am.
37
,
409
-
412
.
Soutar
,
A. R.
and
Fullard
,
J. H.
(
2004
).
Nocturnal anti-predator adaptations in eared and earless Nearctic Lepidoptera
.
Behav. Ecol.
15
,
1016
-
1022
.
Spangler
,
H. G.
(
1984
).
Silence as a defense against predatory bats in two species of calling insects
.
Southwest. Nat.
29
,
481
-
488
.
Sternburg
,
J. G.
,
Waldbauer
,
G. P.
and
Jeffords
,
M. R.
(
1977
).
Batesian mimicry: selective advantage of color pattern
.
Science
195
,
681
-
683
.
Stevens
,
M.
(
2005
).
The role of eyespots as anti-predator mechanisms, principally demonstrated in the Lepidoptera
.
Biol. Rev.
80
,
573
-
588
.
Subchev
,
M.
(
2014
).
Sex pheromone communication in the family Zygaenidae (Insecta: Lepidoptera): a review
.
Acta Zool. Bulg.
66
,
147
-
157
.
Surlykke
,
A.
(
1988
).
Interaction between echolocating bats and their prey
. In
Animal Sonar: Processes and Performance
(ed.
P. E.
Nachtigall
and
P. W.
Moore
), pp.
551
-
566
.
New York
:
Springer
.
Surlykke
,
A.
,
Filskov
,
M.
,
Fullard
,
J. H.
and
Forrest
,
E.
(
1999
).
Auditory relationships to size in noctuid moths: bigger is better
.
Naturwissenschaften
86
,
238
-
241
.
Tomasi
,
C.
and
Manduchi
,
R.
(
1998
).
Bilateral filtering for gray and color images
. In
Sixth International Conference on Computer Vision
, pp.
839
-
846
.
New Delhi
:
Narosa Publishing House
.
Vallin
,
A.
,
Dimitrova
,
M.
,
Kodandaramaiah
,
U.
and
Merilaita
,
S.
(
2011
).
Deflective effect and the effect of prey detectability on anti-predator function of eyespots
.
Behav. Ecol. Sociobiol.
65
,
1629
-
1636
.
Waldbauer
,
G. P.
and
Sternburg
,
J. G.
(
1975
).
Saturniid moths as mimics: an alternative interpretation of attempts to demonstrate mimetic advantage in nature
.
Evolution
29
,
650
-
658
.
Yack
,
J. E.
(
1988
).
Seasonal partitioning of atympanate moths in relation to bat activity
.
Can. J. Zool.
66
,
753
-
755
.
Yack
,
J. E.
,
Scudder
,
G. G. E.
and
Fullard
,
J. H.
(
1999
).
Evolution of the metathoracic tympanal ear and its mesothoracic homologue in the Macrolepidoptera (Insecta)
.
Zoomorphology
119
,
93
-
103
.
Yager
,
D. D.
(
1990
).
Sexual dimorphism of auditory function and structure in praying mantises (Mantodea; Dictyoptera)
.
J. Zool. Lond.
221
,
517
-
537
.
Yager
,
D. D.
(
2012
).
Predator detection and evasion by flying insects
.
Curr. Opin. Neurobiol.
22
,
201
-
207
.
Yen
,
S.-H.
,
Robinson
,
G. S.
and
Quicke
,
D. L. J.
(
2005
).
The phylogenetic relationships of Chalcosiinae (Lepidoptera, Zygaenoidea, Zygaenidae)
.
Zool. J. Linn. Soc.
143
,
161
-
341
.
Zeng
,
J.
,
Xiang
,
N.
,
Jiang
,
L.
,
Jones
,
G.
,
Zheng
,
Y.
,
Liu
,
B.
and
Zhang
,
S.
(
2011
).
Moth wing scales slightly increase the absorbance of bat echolocation calls
.
PLoS ONE
6
,
e27190
.

Competing interests

The authors declare no competing or financial interests.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/3.0), which permits unrestricted use, distribution and reproduction in any medium provided that the original work is properly attributed.

Supplementary information