Owls (Strigiformes) are top predators and thus crucial modulators of biodiversity (top-down regulation, Sergio et al. 2008). Due to their low densities and elusive behavior, owls are difficult to study, and therefore there is little information related to their ecology, particularly for Neotropical species (Rivera-Rivera et al. 2012). The habitat characteristics that determine the presence of Neotropical owls are not well known. Nonetheless, a few studies have suggested that certain forest structure elements, such as large decaying trees, snags, and understory vegetation can be important for owls’ breeding and foraging activities (Borges et al. 2004, Barros and Cintra 2009, Rivera-Rivera et al. 2012, Ibarra et al. 2014). Furthermore, landscape features, such as topography or distance to streams, can affect the abundance and detectability of owls’ prey (Carey et al. 1992, Rivera-Rivera et al. 2012).
Raptor populations of central Argentina have de- clined alarmingly in the last decades, mainly due to habitat loss (Sarasola et al. 2018). However, the res- ponse of nocturnal species (i.e., owls) to this habitat degradation is virtually unknown, probably because there is little background information related to their ecology in this region (but see Campioni et al. 2013). Here we report detections of two small owl species in a mountain forest of central Argentina in relation to forest structure and landscape characteristics. We focused our surveys on two relatively abundant and widespread species: The Tropical Screech-Owl (Me- gascops choliba; hereafter TRSO) and the Ferruginous Pygmy-Owl (Glaucidium brasilianum; hereafter FEPO) (Konig et al. 2008). We also report the densities of these two species in the area. This is the flrst study of owls in relation to habitat characteristics in the mountain forests of central Argentina. However, this preliminary information will be useful for designing future studies on the ecology of owls in central Argentina targeted to their conservation.
methods
Study area
The study was carried out in a 10 000 ha protec- ted area of the mountain forests of Córdoba, Argenti na: Sierras Chicas Corridor (hereafter SCC; 31°12’S, 65°22’W; Fig. 1). The average annual temperature of SCC is 18.9°C and annual rainfall is approximate- ly 670 mm (Gavier and Bucher 2004). We conducted surveys between 700 and 1000 masl. At this altitude, the vegetation is typical of the Chaco Mountain Forest district (Chaco Serrano) and it is dominated by slow- growing deciduous tree species (Cabido et al. 2018). The SCC has an extensive history of disturbances, many of which persist until today (e.g., selective-log- ging, fires, livestock grazing, exotic invasions), being this area highly degraded (Gavier and Bucher 2004).
Owl surveys
We conducted crepuscular-nocturnal surveys (19:00-2:00 h) between late-August and early-Octo- ber 2019, within the early part of the breeding season for both TRSO and FEPO (Carrera et al. 2008, Konig et al. 2008, Schaaf et al. 2019). Peaks of vocal activity of both species occur during these months, increasing detectability (Cerasoli and Penteriani 1992, Konig et al. 2008). On different nights, we surveyed five ran- domly selected transects of 2 km each along existing paths and recorded all spontaneous calls of both species; transects were separated by at least 250 m. Although low detectability could be a problem when studying cryptic species (MacKenzie et al. 2005), for the purposes of this baseline study we assumed per- fect detectability of all calls within ~125 m on either side of the transects, giving a surveyed area of ~250 ha (Fig. 1). To increase owl detections, we avoided ra- iny and windy nights and we used playbacks of both species. Due to the irregularity of the SCC terrain (si- nuous paths and numerous valleys and ravines), we used a minimum distance of 250 m between play- back-stations, shorter than that commonly used in owl studies (e.g., Campioni et al. 2013). We surveyed owls at 40 stations. We broadcasted calls (CD of calls: Narosky and Yzurieta 2010) during 1-min cycles per species using a portable amplifier (Brookstone floa- ting®, volume ~100 dB at 1 m), then we listened for 8 min until we completed 10 min at each station (Ibarra et al. 2014). Preliminary observations in the field did not indicate inhibition between our focal owl species, thus we randomized the order of species calls at each station (Borges et al. 2004, Ibarra et al. 2014). We de- termined the localization of each individual by com- pass triangulation, which also allowed us to reduce the risk of double-counting (Enríquez and Rangel-Sa- lazar 2001). We immediately stopped the playback as soon as an owl was detected, to avoid attracting owls toward playback stations, which could bias location estimations (Borges et al. 2004, Hausleitner 2006, Larson and Holt 2016). Once we identified points with owl presence, we selected at least twice as many "absence” points where owls were not detected, using the random points generator tool in QGIS 3.8.3 (QGIS Development Team 2019) and a Satellite image of the study area (Copernicus Sentinel-2B ESA, data 2019). All absence points were within 100 m of the transects and at least 250 m from any other point (presence or absence). We checked the absence points by revisi- ting them twice on different nights and repeating the playback procedure (Barros and Cintra 2009). Revi- sits were homogeneously distributed within the sur- veying period to reduce the effects of breeding stages on owl detections.
We estimated densities for each owl species fo- llowing Borges et al. (2004). We registered indivi- duals’ positions as inside or outside of a 50 m ima- ginary belt from the center point of the transects. We estimated density according to the following formula:
Density=(n1+n2)/(2 x r x l) x ln (n1+n2)/n2
where n1= number of individuals inside the 50 m belt, n2 = number of individuals outside of the 50-m belt, r = 50-m, and l = length of the transect (2000 m).
Habitat characteristics
We measured landscape characteristics around presence and absence points using QGIS 3.8.3 (QGIS Development Team 2019). Also, around these points, we staked out 30 x 30 m plots in which we measured forest structure characteristics in the field (Table 1; Barros and Cintra 2009). We selected habitat cha- racteristics considered to affect owls’ foraging and reproductive activities (Barros and Cintra 2009, Cam- pioni et al. 2013, Ibarra et al. 2014).
Statistical Analysis
To analyze the effect of habitat characteristics on the presence of each owl species we fitted GLMs (Zuur et al. 2009) with binomial error distribution and logit link function. For each model, we calculated Akaike’s Information Criterion corrected for small sample si- zes (AlCc; Burnham and Anderson 2002). To evaluate the strength of support for each model, we compared models based on AAICc and Akaike weight (w; Burn- ham and Anderson 2002). We assessed the multicolli- nearity of all explanatory variables with Pearson’s Co- rrelations (Graham 2003). The explanatory variables considered non-redundant (r < 0.6) and thus used to model owl presence were: proportion of small trees, number of snags, proportion of exotic trees, canopy cover, slope of terrain, and distance to stream (Table 2). All statistical procedures were performed using software R version 3.6 (R Core Team 2019).
results
The lowest AlCc model for TRSO included only the number of snags (Table 3). The probability of TRSO presence increased with the number of snags (b = 0.52, SE = 0.23, Z = 2.27, P = 0.02; Fig. 2A). The best model for FEPO included only the proportion of small trees (Table 3). The probability of FEPO presence in creased with the proportion of small trees (b = 14.32, SE = 7.66, Z = 1.87, P = 0.03; Fig. 2B).
were measured in the field in 30x30 m plots and landscape characteristics were measured remotely using QGIS. DBH: Diameter at breast height.
in a mountain forest of central Argentina. The best model for each species is indicated in bold letters. Predictor variables included in models appear in bold letters in Table 2, measurements descriptions are in Table 1.
discussion
We conducted the first estimation of TRSO and FEPO densities for the mountain forests of central Argentina. Two forest structure characteristics were associated with the presence of small owls in the mountain forest of central Argentina. Landscape cha racteristics had no effect on owls’ presence in this rea. TRSO was positively associated with the number of snags, and FEPO with the proportion of small trees.
The densities we estimated for TRSO (0.16 individuals/ha) and FEPO (0.054 individuals/ha) were within the range of those reported for other protected areas. Reported values for both species vary among studies, ranging from 0.07 to 2.4 individuals/ha in TRSO (Amaral 2007, Claudino 2013) and 0.06 to 0.23 individuals/ha in FEPO (Campioni et al. 2013, Menq and Anjos 2015). Apparently, densities of TRSO and FEPO in tropical and humid forests are higher than those reported for subtropical and drier forests (e.g., Borges et al. 2004, Campioni et al. 2013). Coinciden- tally, densities in the SCC were similar to lower values. Despite the extensive history of disturbance of the SCC, densities of TRSO and FEPO were within expec- ted values, suggesting that these species might be able to cope with disturbed environments, as proposed by Amaral (2007) and Menq and Anjos (2015). Our sur- veys did not consider imperfections in detectability,
which could lead to underestimation of density (MacKenzie et al. 2005). However, our results contribute with baseline information, which is useful to monitor and study these little-known owl populations.
Studies of Neotropical owl species indícate that forest structure can affect their presence, but expla- nations for these associations are still poorly unders- tood (Enríquez and Rangel-Salazar 2001, Borges et al. 2004, Barros and Cintra 2009, Rivera-Rivera et al. 2012, Menq and Anjos 2015). The association of TRSO with higher snag abundance in the SCC might be re- lated to the availability of nesting sites as this species sually nests in cavities of large living trees (diameter at breast height; DBH > 40 cm) and snags (Schaaf et al. 2019). In the mountains of central Argentina, the abundance of large living trees capable of forming na tural cavities by decay (approximately DBH > 50 cm; Cockle et al. 2011) is very low, suggesting that excava- ted cavities in snags might be the only suitable nes- ting site for TRSO (Gavier and Bucher 2004). Specific studies on the breeding behavior and ecology of TRSO in the SCC could contribute to testing this hypothe- sis, and thus confirm the importance of snags for this species in degraded forests.
The association of FEPO with a higher proportion of small trees could be related to this species’ diet, mainly based on birds that nest in the lower stratum of the forest (e.g., the Creamy-bellied Thrush Turdus amaurochalinus; Carrera et al. 2008, Batisteli et al. 2020). The association with the forest lower strata was also reported for other Glaucidium species (Rive ra-Rivera et al. 2012, Menq and Anjos 2015). Althou- gh further studies are needed to test this hypothesis, vertical forest structure could be important for stud- ying Glaucidium species ecology.
The scarce ecological information on Neotropical owls and the wide distribution of some species, such as TRSO and FEPO, hinder the detection of clear pat- terns that relate these birds with the characteristics of their habitats. Furthermore, detections of these pa- tterns can vary between breeding and non-breeding seasons, or even within each season (Cerasoli and Penteriani 1992). Our study in the mountain forests of central Argentina reveals habitat characteristics that could be important predictors of the occurrence of small owls that inhabit these forests in the early sta- ges of their breeding season. Specifically, snags and small trees should be considered as covariates when studying habitat use of TRSO and FEPO, respectively. Knowledge of how these predators use their habitats could be crucial to conserving their populations and their important ecological role.
acknowledgments
We thank the administrations of Reserva Municipal Natural Los Quebrachitos, Reserva Hídrica Natural Municipal Los Manantiales, and Reserva Provincial Dique La Quebrada for allowing us to conduct our study. We are deeply grateful to C. Schneider and C. Riachi for providing logistic and professional support during this study. We thank all the SCC neighbors for their hospitality. Suggestions and comments from K. L. Cockle and D. Vergara-Tabares greatly enhanced earlier versions of this manuscript.