Interspecific variation in the phenology of advertisement calling in a temperate Australian frog community

Abstract Spatial and temporal partitioning of resources underlies the coexistence of species with similar niches. In communities of frogs and toads, the phenology of advertisement calling provides insights into temporal partitioning of reproductive effort and its implications for community dynamics. This study assessed the phenology of advertisement calling in an anuran community from Melbourne, in southern Australia. We collated data from 1432 surveys of 253 sites and used logistic regression to quantify seasonality in the nightly probability of calling and the influence of meteorological variables on this probability for six species of frogs. We found limited overlap in the predicted seasonal peaks of calling among these species. Those shown to have overlapping calling peaks are unlikely to be in direct competition, due to differences in larval ecology (Crinia signifera and Litoria ewingii) or differences in calling behavior and acoustics (Limnodynastes dumerilii and Litoria raniformis). In contrast, closely related and ecologically similar species (Crinia signfera and Crinia parinsignifera;Litoria ewingii and Litoria verreauxii) appear to have staggered seasonal peaks of calling. In combination with interspecific variation in the meteorological correlates of calling, these results may be indicative of temporal partitioning of reproductive activity to facilitate coexistence, as has been reported for tropical and temperate anurans from other parts of the globe.


Introduction
Niche partitioning in ecological communities may entail differential responses to a large number of environmental axes (Hutchinson 1957). This partitioning is often considered in spatial terms, but the partitioning of resources through time can also facilitate the coexistence of ecologically similar species (Townsend et al. 2008). Evidence of temporal partitioning may be found at various scales, from differing activity times within the diel cycle (Kronfeld-Schor and Dayan 2003), to interspecific variation in reproductive effort between years (Silvertown 2004). However, differences in phenology are the most familiar forms of temporal partitioning within ecological communities.
Here, examples include sympatric mantids whose annual life cycles are staggered to allow coexistence (Hurd and Eisenberg 1989) and temporal segregation of larval development in odonate communities to reduce competitive and predatory interactions (Crowley and Johnson 1982).
In communities of amphibians, interspecific variation in reproductive phenology and its implications for coexistence has received some attention, particularly for species whose reproductive activity can be tracked by male advertisement calls (Wells 2007). Advertisement calling is exhibited by most species of frogs and toads. Within species, calls convey information about the location and fitness of males to prospective partners and competitors; among species, calls evolve distinct acoustics that act as an important premating isolation mechanism (Wells 2007). However, phenological differences also represent a key source of interspecific variation in advertisement calling in anuran communities, signaling considerable variation in the timing of reproduction.
In tropical environments, calling occurs primarily during the rainy season, but there may be considerable heterogeneity in the timing of calling linked to the wide variety of reproductive modes displayed. For example, among 31 species of frogs studied by Gottsberger and Gruber (2004) in French Guiana, species could be characterized as obligate early season callers (species that lay foam nests or display direct development), mid-or late season callers (species with embryonic development over water), sporadic and explosive callers (species that lay eggs directly into water), or generalists that call right through the rainy season (species with parental care of larvae). A comparable diversity of strategies has been reported for other tropical anuran communities, including species that call exclusively during the dry season (Aichinger 1987;Donnelly and Guyer 1994;Bertoluci 1998;Bertoluci and Rodrigues 2002). Temperate anurans may also vary substantially in calling seasonality. Among 13 species in eastern Texas, Saenz et al. (2006) documented calling seasons ranging from 2 to 12 months. In Ontario, de Solla et al. (2006) monitored a community of eight species and found little or no overlap in the calling seasons of several taxa.
Sympatric anurans can also vary markedly in their meteorological cues for calling. Even among aquaticbreeding amphibians from the same community, responses to rainfall can vary considerably, from explosive initiation of calling in response to rain, to dampening or even cessation of calling in response to rain (Oseen and Wassersug 2002;Gottsberger and Gruber 2004;Saenz et al. 2006). Positive and negative relationships between environmental temperature and calling activity are also known in anuran communities, with some species favoring warm conditions and others cooler conditions (Gottsberger and Gruber 2004;Steelman and Dorcas 2010). Similarly, while the fundamental effects of relative humidity and wind strength on water balance and thermoregulation in anurans (Tracy 1976) suggest that the former should be positively correlated with calling activity and the latter negatively, the relationship between these variables and calling may be idiosyncratic (Oseen and Wassersug 2002).
While it is clear that the phenology of calling can vary widely within anuran communities, it is also true that the available data on the subject are overwhelmingly derived from anuran communities in the Americas (Wells 2007). Few, if any, studies are available for large parts of Africa and Eurasia. The same is true for Australia. Of the 241 frogs so far described in Australia (Anstis 2013), only nine species (4%) have been the focus of dedicated phenological studies on calling and reproduction to our knowledge (MacNally 1984;Williamson and Bull 1992;Driscoll 1998;Brooke et al. 2000;Lemckert 2001;Hollis 2003;Wong et al. 2004;Hauselberger and Alford 2005;Canessa et al. 2012). Likewise, we are aware of only seven studies of the phenology of calling across an Australian frog community (Humphries 1979;Gillespie 2001;Morrison 2001;Lemckert and Mahony 2008;Lemckert and Grigg 2010;Dostine et al. 2013;Lemckert et al. 2013).
Here, we describe the phenology of advertisement calling in an anuran community from the temperate climes of southern Australia, focusing on interspecific variation in both the seasonality of calling and the meteorological cues for calling. We combined data collected over 11 years to develop models of the nightly probability of calling by the focal species. As has been documented for tropical and temperate anuran communities elsewhere, our data reveal considerable interspecific variation in the seasonality of calling and the weather conditions that favor calling, which may be indicative of temporal partitioning of reproductive activity in this community.

Study sites and surveys
We conducted 1432 nocturnal surveys for frogs between 2001 and 2012 at 253 sites across the volcanic plains north and west of Melbourne, Victoria, Australia. Sites were distributed in the catchments of the Darebin, Merri, Moonee Ponds, and Kororoit Creeks. The landscape is undulating, rising to a maximum elevation of~300 m asl at the study sites. Upper catchments are dominated by grazing land, graduating to industrial and urban estates in the lower catchments. The climate is temperate, with cool winters (July mean daily temperature range = 5.4-13.1°C) and warm summers (February mean daily temperature range = 14.1-26.6°C) (BOM 2014). Rainfall averages 538 mm, with the highest monthly rainfall occurring in late winter and spring (August-November) in most years (BOM 2014). Sites included slow flowing pools along the main channel and tributaries of the streams listed above, as well as lentic wetlands such as farm dams, flooded quarries, swamps, and water treatment ponds. Site selection was haphazard, including both randomized protocols (mostly for selecting pools along streams) and nonrandomized protocols that aimed to maximize geographic coverage while accounting for both access and logistical constraints (for further details, see Canessa et al. 2012;Canessa and Parris 2013;Heard et al. 2015).
Surveys were conducted between September and April (inclusive) in each year, encompassing the Austral spring, summer, and autumn. Survey timing sought to optimize detection rates of Litoria raniformis in most cases; this endangered species having been studied intensively in this region since 2001 (Heard et al. 2012a(Heard et al. ,b, 2013(Heard et al. , 2015. However, aural surveys during the spring and early summer months optimize detection rates across the frog community in this region (Canessa and Parris 2013). Survey techniques were broadly consistent between sites and years. Surveys commenced at least 30 min after dark and began with a period of 5-10 min listening quietly for male advertisement calls. Frogs in the study area have distinctive calls that allow ready identification to species. Sites were then searched on foot with the aid of spotlights. Any additional calls heard during these searches were recorded. Surveys extended up to 265 min, with up to 10 repeat surveys at a site in any given year. Surveys were completed under permits 10001816, 10003005, and 10005649 issued by the Victorian Department of Environment, Land, Water and Planning.
To supplement the above survey dataset, we collated data from 263 nocturnal surveys for frogs completed at the same study sites in the same period by private consultants, other researchers, and local naturalists. All surveys were conducted during the Austral spring, summer, or autumn, and followed surveys protocols analogous to our own. Survey data were gathered from the relevant research report or thesis, or through personal communications. Full details are provided by Heard et al. (2015).

Meteorological variables
Six meteorological variables were considered as potential determinants of calling activity. Air temperature (dry bulb,°C) and relative humidity (%) were recorded during most surveys with the aid of a compact whirling hygrometer (Brannan P/L, Cumbria, UK). Surface water temperature was recorded using a digital probe thermometer (varying brands). Wind strength was recorded subjectively on an ordinal scale between 0 (still) and 3 (strong, regular gusts). Rainfall (mm) was derived from the Australian Bureau of Meteorology (BOM 2014), taking measure-ments from the Melbourne Airport weather station, located roughly in the center of the study area (-37.67 S, 144.83 E). Two measures of rainfall were used: rainfall on the day of the survey and cumulative rainfall over the preceding 7 days ("lagged" rainfall).

Modeling
We used Bayesian logistic regression to model the nightly probability of calling for those species with a sufficient number of detections. For each of these species, we began by filtering the nondetection data to exclude records from sites in years when no detections of calling were made. Hence, for each species, only detection and nondetection data from sites that were known to be occupied by reproductively active males in a given year were included. We also removed any surveys when the focal wetland was dry (which precludes calling among the focal species) and randomly filtered the remaining dataset for each species such that only one site was included for any given survey night. Repeat surveys of the same site in a given year were assumed to be independent. Hence, for each species, the detection (1) or nondetection (0) of calling on survey night i at a site known to be occupied by reproductively active males was treated as a Bernoulli variable with probability p i .
Covariate effects on p i were modeled using a logistic equation and linear link function: where a is the intercept, and b x is the regression coefficient for covariate X. The underlying seasonal variation in p i was modeled using a Fourier series approach following Burnham and Anderson (2002, p. 188). The approach treats phenomena as displaying a regular pattern that may be described using cosine functions of increasing complexity. We chose a simple two parameter version. The change in p i with days since the start of the Austral spring (September 1, days i ) was modeled as: With the exception of wind strength (which did not show any obvious seasonality), raw meteorological data for each survey were first adjusted to represent the anomaly from the seasonal mean. After Canessa et al. (2012), we fitted up to fourth degree polynomial regressions between each meteorological variable and survey date using maximum likelihood in R version 3.0.3 (R Core Team 2014). The best-fitting model for each meteorological variable was selected using Akaike's information criterion (Burnham and Anderson 2002), and the residuals from this model used as inputs for the logistic models.
Fourteen models of p i were fitted to the detection data for each species. Each model contained one of the two rainfall variables, plus one or more of temperature, relative humidity, and wind strength. Models were fitted twice for each species; once with temperature represented by air temperature and once with temperature represented by water temperature. The best-fitting temperature variable for each species was selected based on model selection statistics. Correlations between predictors were weak in the model set (Pearson's r correlations ranged from À0.13 to 0.25), except for a moderate negative correlation between air temperature and humidity (Pearson's r of À0.46). An effect of survey effort (survey duration multiplied by the number of observers) was included in all models to account for the fact that the probability of detecting males that are calling sporadically increases with survey duration, and multiple surveyors have a higher cumulative chance of hearing sporadic or weakly calling individuals. A random "year" effect was also included to account for annual variation in calling activity.
Models were fitted to the data using Markov Chain Monte Carlo (MCMC) sampling in OpenBUGS version 3.2.3 (Thomas et al. 2006), called from R. Uninformative priors were used for each parameter. With the exception of survey date, covariates were centered by subtracting the mean and dividing by two standard deviations. Parameter estimates and their 95% credible intervals (95% CI) were drawn from 40,000 MCMC samples after a burn-in of 60,000 samples. Models were compared using the deviance information criterion (Spiegelhalter et al. 2002). The best-fitting model for each species was used to predict their nightly probability of calling across 366 days, beginning and ending on September 1. In turn, these predictions were used to estimate the seasonal peak in the probability of calling for each species.
Crinia signifera and Lit. ewingii were recorded calling in all survey months, but most frequently in early spring (Fig. 1). For these species, the estimated seasonal peak of calling was in July (Fig. 3), with a steep decline in the nightly probability of calling from November to February, followed by a rapid recovery during autumn (Fig. 2). Limnodynastes tasmaniensis displayed a similar (although weaker) pattern, with observed calling in all survey months (Fig. 1), but an estimated peak in the nightly probability of calling in September and trough in March (Figs. 2, 3). Litoria verreauxii is primarily a spring breeder. For this species, the proportion of surveys in which calling was detected declined in a roughly linear fashion after September (Fig. 1); however, the estimated seasonal peak in the probability of calling was in October (Fig. 3). Thereafter, the nightly probability of calling declined sharply, reaching near zero in mid-autumn (Fig. 2). Limnodynastes dumerilii and Lit. raniformis displayed the closest pattern of calling activity. The nightly probability of calling was estimated to peak in early December for both species (Fig. 3), before declining to zero or slightly thereabouts during autumn and winter (Fig. 2). Crinia parinsignifera and Lit. peronii were only detected calling between November and early February.
Interspecific variation was also apparent in relationships between the nightly probability of calling and the meteorological variables considered (Tables 1, 2). Relatively high rainfall in the preceding 7 days increased the probability of calling by C. signifera, Lim. dumerilii, Lim. tasmaniensis, and Lit. ewingii (Table 2). However, there was a clear negative effect of lagged rainfall on the probability of calling by Lit. raniformis (mean estimate = À0.60, 95% CI = À1.10, À0.11; Table 2). Likewise, while the probability of calling by both Lim. dumerilii and Lit. raniformis increased as the temperature anomaly increased, the former responded more strongly to air temperature, while the later responded to water temperature (Tables 1, 2). Relative humidity was only related to the probability of calling for Lit. ewingii and Lit. verreauxii, having a positive effect in both cases (Lit. ewingii: mean estimate = 1.15, 95% CI = 0.39, 1.96; Lit. verreauxii: mean estimate = 1.1, 95% CI = 0.06, 2.22; Table 2).

Discussion
There are two obvious mechanisms through which temporal partitioning of reproductive effort may facilitate coexistence in anuran communities. The first is reducing acoustic competition. Intraspecifically, male anurans compete acoustically for mates in various ways, including extrinsic factors such as caller location (Littlejohn 1977;Parris 2002) and intrinsic factors such as call duration and acoustic frequency (Littlejohn 1977;Wagner 1992;Welch et al. 1998). Interspecifically, acoustic competition may arise purely through interference (Littlejohn 1977). For example, Littlejohn and Martin (1969)  Only surveys at sites that were known to be occupied in a given year are included for each species (calling was detected at least once at those sites in those years). Schwartz and Wells (1984) reported reductions in the attractiveness of calls of Hyla ebraccata in Panama when the primary notes overlapped with calls of the sympatric H. microcephala. In an extensive study of 39 species across three anuran communities in Brazil, Duellman and Pyles (1983) demonstrated considerable overlap in the acoustic properties of calls. They concluded that temporal partitioning of calling, both in terms of the seasonal timing and environmental conditions under which calling takes place, was an important mechanism by which acoustic interference is mitigated (and coexistence facilitated) in these communities. The second mechanism by which partitioning of reproductive phenology may facilitate coexistence in anuran communities is mitigation of larval competition. Competitive interactions between larval anurans are a key determinant of larval survival and metamorphic success (Alford 1999). Resource competition among larvae also has important carry over effects on the fitness and survival of postmetamorphic anurans (Chelgren et al. 2006). Hence, larval competition represents a key determinant of anuran population dynamics and community structure (Wells 2007).
The early work of Alford and Wilbur on the subject of larval competition in anuran communities produced evidence both for and against temporal partitioning as a mechanism for coexistence Wilbur and Alford 1985;Alford 1989a,b). Temporal partitioning was important for the outcomes of some interactions Wilbur and Alford 1985;Alford 1989a), but not others (Alford 1989b). Later work confirmed that the strength of larval competition between species varies intrinsically, and mediates the importance of temporal partitioning for larval fitness (Gascon 1992). It also showed that the temporal order in which interactions take place (so-called "priority effects") may be an important determinant of the strength of competitive interactions among anuran larvae (Lawler and Morin 1993).
We cannot demonstrate that the interspecific variation in calling phenology documented here is adaptive with regard to coexistence, but there are several interesting points to be made in this regard. In the focal community, two species pairs could be expected to benefit most from partitioning of reproductive phenology: Lit. ewingii & Lit. verreauxii and C. signifera & C. parinsignifera. These frogs are closely related, have analogous larval ecologies and similar advertisement calls (Loftus-Hills 1973; Watson et al. 1985;Anstis 2013). Litoria ewingii and Lit. verreauxii also occasionally hybridize (Watson et al. 1985;Smith et al. 2012). Our data suggest that Lit. ewingii and Lit. verreauxii display subtle differences in their seasonal calling peak, with Lit. ewingii being primarily a winter bree-der and Lit. verreauxii having a calling peak in early to mid-spring. However, there was considerable uncertainty surrounding our estimates of calling seasonality for Lit. verreauxii, and we note that an earlier, winter peak in calling activity has been recorded elsewhere for this species (Humphries 1979;Smith et al. 2003;Lemckert and Grigg 2010). Nevertheless, these studies of Lit. verreauxii generally took place in allopatry from Lit. ewingii. The lagged peak in calling for Lit. verreauxii suggested by our study could be indicative of temporal partitioning when in sympatry with Lit. ewingii. Our data are very limited for C. parinsignifera; however, they suggest that the peak calling period of C. parinsignifera occurs later than that of C. signifera, and this agrees with observations elsewhere (Humphries 1979;MacNally 1979). Indeed, MacNally (1979) reported that C. signifera truncates its spring calling period to minimize overlap with C. parinsignifera when the two species occur in sympatry.
It is also notable that the species pairs with the closest seasonal patterns of calling activity -C. signifera & Lit. ewingii and Lim. dumerilii & Lit. raniformisdisplay traits that should weaken competitive interactions. Tadpoles of C. signifera and Lim. dumerilii are bottom feeders, whereas those of Lit. ewingii and Lit. raniformis tend to feed higher in the water column (Anstis 2013; pers. obs.). As such, resource competition between larval C. signifera and Lit. ewingii and between larval Lim. dumerilii and Lit. raniformis may be limited. Likewise, while C. signifera and Lit. ewingii call over a similar frequency band (Parris et al. 2009), the former calls from concealed positions under vegetation or crevices, while Lit. ewingii generally calls from elevated positions in vegetation (pers. obs.). Elevated perches can significantly increase the effective distance of calls (Parris 2002); hence, acoustic interference of Lit. ewingii by C. signifera may be mitigated by microhabitat partitioning. Different caller positions may also mitigate acoustic interference between Lim. dumerilli and Lit. raniformis; Lim. dumerilli calls from hidden positions under banks and rocks, or in emergent vegetation (Anstis 2013;pers. obs.), whereas Lit. raniformis calls primarily from mats of aquatic vegetation (Heard et al. 2008;pers. obs.). These two species also favor slightly different meteorological conditions for calling, and there calls differ acoustically (dominant frequency of 680 Hz for Lim. dumerilii vs. 1200 Hz for Lit. raniformis; Loftus-Hills 1973).
Our data revealed both expected and unexpected relationships with regard to the meteorological correlates of calling in this community. The probability of calling increased with increasing rainfall over the preceding 7 days for C. signifera, Lim. dumerilii, Lim. tasmaniensis, and Lit. ewingii. None of these species displayed strong relationships with rainfall on the day of survey. Hence, the influence of rainfall on calling activity is most likely mediated by its effect on wetland water levels, as is often the case for aquatic-breeding anurans (Wells 2007). In contrast, rainfall curtailed calling activity by Lit. raniformis. Although somewhat counterintuitive, this relationship has been reported for ecologically analogous ranid frogs from North America. Saenz et al. (2006) demonstrated a negative relationship between calling activity and rainfall for both Lith. catesbeiana and Lith. clamitans in eastern Texas. Lit. raniformis, Lith. catesbeiana, and Lith. clamitans are all largely aquatic, summer breeding frogs that occupy permanent wetlands (Heard et al. 2008(Heard et al. , 2013Saenz et al. 2006;this study). They also all favor nights with warmer water temperatures for calling (Oseen and Wassersug 2002;Saenz et al. 2006;this study). While it may be supposed that water temperature is lower during rainy periods, and this could account for a negative effect of rainfall on calling activity by Lit. raniformis, lagged rainfall and water temperature were unrelated in our dataset (Pearson's r = À0.04). Saenz et al. (2006) hypothesized that acoustic interference from storms or rainfall itself may be the underlying mechanism. We can see two potential mechanisms for Lit. raniformis. Firstly, it may be that males of Lit. raniformis reduce calling activity following rainfall to reduce acoustic interference with other species in the community, with calling by four species being positively related to lagged rainfall (as above). Secondly, the negative relationship between calling activity and rainfall in Lit. raniformis may stem from the sensitivity of this species to rising water levels in lotic situations, where increased flows represent a threat to eggs and weakly swimming early-stage tadpoles. Of the frogs Table 1. The five top-ranked models of the nightly probability of calling for the six commonly detected species. All models include cosine effects of days since September 1 and a linear effect of survey effort (to account for the fact that the probability of detecting calls increases with effort). Dots show the meteorological variables included in each model. They are as follows: air or water temperature anomaly (Temp), humidity anomaly (Hum), wind strength (Wind), rainfall anomaly (mm) either for the day of survey (Rain) or for the preceding 7 days (Rain lag). Model selection statistics are the deviance information criterion (DIC), distance from the top model (DDIC), and model weight (w).

Species
Temp detected in our study area, Lit. raniformis was the most frequently observed in lotic habitats. In line with Lemckert and Grigg (2010), our data suggest that calling by Lim. dumerilii is stimulated by warmer than average air temperatures during the calling season. Lim. dumerilii is a terrestrial species that can range up to 1 km from wetlands during the nonbreeding season (Humphries 1979;Anstis 2013). We speculate that the tendency of male Lim. dumerilii to call on nights with comparatively warm air temperatures relates to the suitability of such nights for female migration to breeding sites. Finally, the positive relationship between the probability of calling and relative humidity demonstrated here for Lit. ewingii and Lit. verreauxii may be readily explained by their propensity to call from elevated perches or on land, rather than in the water (Anstis 2013; pers. obs.). Higher relative humidity facilitates hydroregulation in these circumstances (Tracy 1976), allowing prolonged calling.
Other ecological mechanisms could account for the variation in reproductive phenology observed during this study. For example, the focal species may vary in their susceptibility to terrestrial or aquatic predators, and time reproduction to minimize exposure to particular predators (e.g., Lips 2001). Invertebrates and fish are important larval predators for some of the frogs studied here; however, the available evidence suggests that they respond spatially rather than temporally to these predators Parris 2010, 2013). We encourage further research on the reproductive phenology of amphibians and the mechanisms underlying interspecific variation in this trait within communities. As well as insights into niche partitioning and community dynamics, research on reproductive phenology can inform assessments of the effects of anthropogenic disturbances on amphibian communities, including disruptions from climate change (Todd et al. 2011). In particular, we encourage research on the reproductive phenology of amphibians from Australia, Africa, and Eurasia, for which data are scarce.