Next Article in Journal
In Situ Gel with Silver Nanoparticles Prepared Using Agrimonia eupatoria L. Shows Antibacterial Activity
Next Article in Special Issue
European Rabbit Invasion in a Semi-Arid Ecosystem of Chile: How Relevant Is Its Role in Food Webs?
Previous Article in Journal
Hippocampal Metabolic Alterations in Amyotrophic Lateral Sclerosis: A Magnetic Resonance Spectroscopy Study
Previous Article in Special Issue
Abundance and Dynamics of Small Mammals in New Zealand: Sequential Invasions into an Island Ecosystem Like No Other
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Numerical Response of Owls to the Dampening of Small Mammal Population Cycles in Latvia

1
Department of Zoology and Animal Ecology, Faculty of Biology, University of Latvia, Jelgavas Iela 1, LV-1004 Riga, Latvia
2
Latvian Ornithological Society, Skolas Iela 3, LV-1010 Riga, Latvia
*
Author to whom correspondence should be addressed.
Life 2023, 13(2), 572; https://doi.org/10.3390/life13020572
Submission received: 14 December 2022 / Revised: 5 February 2023 / Accepted: 14 February 2023 / Published: 17 February 2023
(This article belongs to the Special Issue Abundance and Dynamics of Small Mammals and Their Predators)

Abstract

:
Strong numerical and functional responses of owls to voles in cyclic environments are well known. However, there is insufficient knowledge from the boreonemoral region in particular, with depleted populations of small mammals. In this study, we describe the dynamics of the small mammal population in Latvia from 1991 to 2016 and link them to owl population characteristics. We used food niche breadth, number of fledglings, and population trends to lay out the numerical response of six owl species to dampened small mammal population cycles. We found temporarily increasing food niche breadth in tawny and Ural owls. There were no other responses in the tawny owl, whereas the breeding performance of three forest specialist species—pygmy, Tengmalm’s, and Ural owls—corresponded to the vole crash years in Fennoscandia. Moreover, the populations of forest specialist owls decreased, and the change in the Ural owl population can be attributed to the depletion of small mammal populations. We found evidence of a carry-over effect in the eagle owl arising from a strong correlation of declining breeding performance with the small mammal abundance indices in the previous autumn. We conclude that dampening of the small mammal population cycles is an important covariate of the likely effects of habitat destruction that needs to be investigated further, with stronger responses in more specialized (to prey or habitat) species.

1. Introduction

Small mammals play an important role in various ecological processes. This role ranges from influence on a natural succession [1] through influence on plant and microorganism community composition and chemistry [2] to demographic processes of small mammal predators [3,4,5,6] and even population processes and the behavior of directly unrelated species [7,8,9]. Fluctuations in the availability of small mammals are known to affect the demography of several species of birds of prey, and the analysis of the breeding performance of the latter can reveal large-scale spatiotemporal patterns of population dynamics of the former [10]. In the boreal region, small mammal populations typically show strict cyclicity over long time periods and spanning hundreds of kilometers [10,11,12]. On a smaller spatial scale, the synchrony of cycles has also been documented in western Europe [13] and the boreonemoral region, specifically in the Baltic states [14,15,16]. Based on long-term rodent abundance studies in Fennoscandia [12], substantial evidence for predation being the main reason for small mammal population cyclicity has been gathered [17]. The predation hypothesis suggests that the rodent abundance gradient reflects the relative influence of destabilizing specialists and stabilizing generalists on vole dynamics, modulated by the presence of snow cover [17]. This hypothesis has also been supported by the characteristics of rodent dynamics in central and western Europe [13]. In the late 20th century, rodent cycles showed irregularities and dampening, particularly in the more northern latitudes [18,19,20,21,22]. Due to the large spatial extent, this dampening has been explained by climate forcing—a decrease in delayed density dependence caused by milder winter conditions [23,24,25]. Nonetheless, the generality of this hypothesis has been refuted due to the return of the vole cycle in southern Finland [26].
There is abundant evidence revealing the importance of small mammals, particularly voles, to owls. First, voles form an important part of the owl diet, ranging from more vole-specialized long-eared owl (Asio otus, hereafter ASIOTU) [27,28,29,30,31] and Tengmalm’s owl (Aegolius funereus, AEGFUN) [32] to more generalist species such as the Eurasian pygmy owl (Glaucidium passernium, GLAPAS) [33,34], the Ural owl (Strix uralensis, STRURA) [34,35], the tawny owl (Strix aluco, STRALU) [28,34,36], and the eagle owl (Bubo bubo, BUBBUB) [37]. Second, voles are one of the key elements ensuring higher breeding performance of owls (generally, [6,34]; ASIOTU [29]; AEGFUN [38,39]; GLAPAS [40]; STRURA [41]; STRALU [36]) and survival (ASIOTU [42]; AEGFUN [43]; GLAPAS [44]; STRURA [3,39]; STRALU [36,45]). Finally, vole density affects owl behavior [46,47], migratory decisions [32], and life history via the carry-over effect (events that occur in one season but influence individual success in the following season) [48]. Most of studies on the subject have been conducted in cyclic environments of Fennoscandia, and there have been very few studies conducted in the Baltic states (see [28,36] and references therein). We are not aware of other studies focusing on the boreonemoral region during prolonged periods of depleted small mammal population dynamics.
In this study, we described small mammal population densities in Latvia for the period from 1991 to 2016 and linked them to several population characteristics of six owl species, i.e., ASIOTU, AEGFUN, GLAPAS, STRURA, STRALU, and BUBBUB. First, we described the owl diet and its changes to establish the importance of small mammals and, in particular, voles for different owl species. Then, we divided the owl species based on food niche breadth and the overall proportion of voles in their diet. Next, we compared nationwide owl population trends and their slopes with cyclic vole populations and since dampening. Finally, we described owl breeding performance and related some demographic parameters to the observed population change among owls.

2. Materials and Methods

2.1. Location and Field Methods

This study was conducted in Latvia, northern Europe (Figure 1a). The country is located in the boreonemoral region [49], with a humid continental climate [50]. It lies within the continuous distribution of all six investigated owl species (ASIOTU, AEGFUN, GLAPAS, STRURA, STRALU, and BUBBUB) [51].

2.1.1. Small Mammal Monitoring

Monitoring of the relative abundance of small mammals was conducted with snap traps from 1991 to 2016. It consisted of two schemes: the first with 2 transects per site, with 1 in a forest and 1 in a meadow (1); and the second with 11 transects per site, with 1 in a meadow and 10 in different forest habitats (2).
The first scheme was officially run from 1991 to 2011. A total of 100 snap traps per transect (approx. 5 m between traps) were applied in autumn (August–September) for three to four days [52]. Volunteers partially repeated this monitoring in 2015 and 2016. This scheme was conducted at four sites, but not every site was monitored every year (Table 1).
The second scheme was conducted in autumns (August–September) during the period of 2012–2016 by volunteers. In this scheme, 20–25 snap trap transects (approx. 5 m between traps) were applied in four areas (Table 1), although not all the areas were monitored every year. Forest transects were stratified into 10 categories as follows (minimum rotation ages for dominant tree species in Latvia are provided in Table A1):
  • YP—young (clearcuts and stands <7 years old) stands on poor soils;
  • YF—young (clearcuts and stands <7 years old) stands on fertile soils;
  • MPU—medium-aged (between 8 years and 80% of rotation age) stands on poor soils without drainage;
  • MFU—medium-aged (between 8 years and 80% of rotation age) stands on fertile soils without drainage;
  • MPD—medium-aged (between 8 years and 80% of rotation age) stands on poor drained soils;
  • MFD—medium-aged (between 8 years and 80% of rotation age) stands on fertile drained soils;
  • OPU—older (≥80% of rotation age) stands on poor soils without drainage;
  • OFU—older (≥80% of rotation age) stands on fertile soils without drainage;
  • OPD—older (≥80% of rotation age) stands on poor drained soils;
  • OFD—older (≥80% of rotation age) stands on fertile drained soils.
The locations of all the small mammal monitoring sites are shown in Figure 1b.

2.1.2. Owl Diet during the Breeding Season

Owl diet analysis was based on prey remains and pellets found in nests or on the ground near the nest. Only the material from a single breeding occasion was used, based on annual nest-box and cavity inspections (GLAPAS, AEGFUN, STRALU, and STRURA) or based on the assumption that pellets cannot survive for many months in open nests or on the ground in the case of ASIOTU. The material was collected in autumn or winter from nest boxes and cavities and during chick ringing from ASIOTU nests. When collecting material, all the soft contents of nest boxes and cavities were removed. In the case of ASIOTU, all the useable material was collected. The distribution of the owl diet sampling sites in Latvia is shown in Figure 1c.
Analysis of the prey remains and identification of the minimal number of individuals were conducted as described by Vrezec et al. [35]. Insects were assumed to represent 1 g biomass, and amphibians and reptiles were assumed to weigh 16 g based on the average body mass of 100 measured individuals during chick ringing in 2016. In birds, the reference size group (i.e., woodpigeon, song thrush, chaffinch, and chiffchaff) weight from the general literature [53] was utilized. Region-specific weight of mammals from our trapping data or the literature [28,54] was utilized. We assumed young hare Lepus sp. to weigh 350 g.

2.1.3. Owl Population Change Monitoring

Monitoring of owl population change was conducted with traditional (territory mapping with playback broadcasting) methods [55,56,57,58,59] in permanent sample areas from 1991 to 2020, as well as with fully standardized point counts (with playback broadcasting) from national Breeding Birds of Prey Monitoring [60] from 2015 to 2022. The spatial distribution of monitoring sites is shown in Figure 1d.

2.1.4. Owl Breeding Performance

We used the number of fledglings per successful nest as a nest-level descriptor of the breeding performance strongly related to food availability. We used four data sources: nest-box inspections (1), information reported by ringers to the Latvian Ringing Centre (2), reports in the Breeding Birds of Prey Monitoring and the previous Monitoring for Owls (3), and citizen scientist reports on nature observation platform dabasdati.lv (4). The first two sources, as well as the third source (partially), covered information from the nests during ringing just before the young fledge. The citizen science and monitoring databases (partially) covered information on the number of young soon after fledging.
Most of the information before 2010 (apart from monitoring data) did not contain exact coordinates to be employed in spatial filtering and removal of duplicate records. Therefore, we combined nest-box inspection, monitoring, and ringing databases based on location attributes (indicated by nest name given by the ringer, which was most often also the person performing monitoring) to remove duplicates. We used citizen science reports only if there was no other information on the species in the particular year in the particular spatial reference. We used the national 1 km projected coordinate grid (epsg: 3059) if coordinates were known or the reported municipality otherwise as a spatial reference.

2.2. Data Analysis

We used R software [61] for data analysis, using the ‘tidyverse’ package [62] for data processing and visualizations and ‘sf’ [63] for spatial data. We treated results with p-values ≤ 0.05 as statistically significant but also reported insignificant results with full test descriptions.

2.2.1. Small Mammal Monitoring

We standardized the small mammal trapping data to the number of individuals per 100 trap days for further analysis and filtered only for autumn counts due to low representation of spring data. We used a graphical representation of the standardized counts per sampling area and habitat to compare variability between sites and habitats. We used generalized linear mixed-effects modelling (GLMM) to compare the differences in standardized densities between sampling areas, age classes, and soil fertility classes (including meadow habitats as a separate class in the latter two). We created two main effect models per comparison:
  • Random intercept per transect and the comparable variable in the fixed part;
  • Random intercept per transect and the comparable variable and year as a factor in the fixed part.
We utilized the Poisson family of distribution with a logarithmic link function and selected the best model based on the lowest value of sample-size-corrected Akaike information criterion value (AICc) [64]. We applied marginal means contrasting with Tukey’s p-value correction for a post hoc analysis of the comparable variable between groups. For mixed-effects modelling, we implemented R packages ‘lme4′ [65] ‘emmeans’ [66] for contrasting.
We found no differences in either peaks or depressions between the sample areas in graphical analysis or mixed-effects models. Therefore, we used information from all the areas to obtain the countrywide small mammal population change index with TRIM analysis implemented in R package ‘rtrim’ [67]. The baseline model in this tool can generally be written as:
ln(μij) = αi + βj,
where μij is an expected count, αi is a site parameter for site i, and βj is a time-point parameter for year j (for full explanation, see [68]). We created a model (model = 3) for the number of individuals pooled across species (small mammals) and separately for voles of genera Microtus and Clethrionomys. We evaluated serial correlation and overdispersion and selected the best model based on the lowest AIC value. We used the multiplicative slope of imputed values reported by TRIM to describe the overall population trend.
We utilized graphical evaluation and Pearson’s correlation analysis to compare the dynamics of both vole groups using yearly indices produced by TRIM analysis. We evaluated the presence of cyclicity with autocorrelation function analysis.

2.2.2. Owl Diet during the Breeding Season

We utilized only samples with at least 5 mammal prey individuals to avoid the influence of some very small samples. We used the small mammal data grouped according to genus as in [69] and other prey groups as described in Section 2.1.2. to calculate Levin’s niche breadth (FNB):
B = 1/∑pi2,
where pi is the fraction of a given prey in the total consumed biomass [70].
In Supplementary Materials, Table S1, we provide a description of the abundance and weight of the annual diet per owl species. For the description, we used the total number and cumulative biomass of the prey group, their arithmetical mean values, and proportions with Wilson’s 95% confidence intervals.
We utilized linear regression (LM) analysis to evaluate the overall temporal change in food niche breadth. In species with repeated samples from the same nest sites, we also created linear mixed-effects models (LMMs) with nest ID as a random intercept. In the case of singularity in the random part, we employed the result of LM [71,72], as there were too few replicates per nest to contribute to the model fit. Model coefficients in the fixed part did not differ in these models. We fitted simple regression models with niche breadth as a response and year as a regressor variable. Both types of models were parameterized for Gaussian residual distribution with the identity link function.
Then, we evaluated the effect of each previously described small mammal population index on niche breadth with the same LM and LMM approach. We fitted models with the independent variable of the year of breeding, admitting a possibly reduced effect as the small mammal yearly indices represent late summer to early autumn rather than late spring to early summer, when owl breeding occurs.
Additionally, we used the biomass proportion of bank voles and voles of genus Microtus for regression analysis, as well as estimated population indices and bank vole proportion, with the index of Microtus voles. In this set of analyses, we applied GLMMs with a random intercept for nest ID and generalized linear models (GLMs) if the identifier did not contribute to model performance due to the number of replicates being too low [71,72]. We used the binomial family of distribution with logistic link function to compare the proportion of prey in the diet with its annual population index value in nature the same year.
Due to the small number of samples, we only described the AEGFUN diet without further analysis.

2.2.3. Owl Population Change Monitoring

Owl population monitoring started in 1990, but not all species had reliable information since the beginning of the period. Due to its preference for large forest massifs, STRURA was covered only since 1993. Due to a lack of knowledge on the monitoring of GLAPAS, its population change can be analyzed only since 2004. Due to the low population size, we did not have reliable data on the population change of BUBBUB.
We combined the data from both owl monitoring schemes if sites had all the planned census activities—for sample areas, sufficient coverage marked by an observer and for point count sites, four standardized visits to each point every year. To analyze population change, we drew on TRIM as described in Section 2.2.1. We exploited only the sites with information from at least two years and comparable effort according to the prerequisites [73].
We further calculated yearly indices and overall population change, as described by multiplicative slope with standard error [68,74], covering the whole available data period for the species. Then, we calculated two different slopes with relatively cyclic small mammal populations and since cycles had vanished. We used 2004 as the threshold for this division because it is:
  • The approximate time since when the small mammal populations did not recover to previous peaks;
  • The approximate midpoint of small mammal monitoring;
  • The approximate midpoint of STRURA monitoring;
  • The beginning of GLAPAS monitoring.
To calculate “before” and “after” trends, we selected the necessary parts of yearly indices and conducted linear regression on ln-transformed indices. We defined index 1991–2004 as “before” and index 2004–2016 as “after”. To obtain significance tests, we transformed the time to start with year 1 in each group to use in regression with interaction between time and period. We defined “before” as a reference level.

2.2.4. Owl Breeding Performance

We calculated the annual mean values of spatially cleaned results because in most cases, we were not able to match different breeding performance reports to an exact breeding location (nest or territory). To provide more generalizable information, we calculated bootstrapped 95% confidence intervals from 1000 bootstrap resamples and implemented them in visualization. We used the annual mean values to establish a temporal trend of overall change in breeding performance. We employed slightly different approaches for further processing because the amount of information varied among owl species.
For the two species (STRALU and STRURA) with most data available on a nearly annual basis, we compared the temporal trends before and after the dampening of small mammal cycles in 2004. We used Gaussian linear regression with the annual mean number of fledglings as a response and compared the trends between the periods as in population change analysis.
We used all the years available to correlate the annual mean breeding performance of STRALU and STRURA with small mammal population indices in the year of breeding and one year before to evaluate a possible carry-over effect. In the case of ASIOTU and BUBBUB, we used a reduced timeframe (from 2002 and 2001, respectively) to avoid possible artefacts due to irregular sampling and a small sample size. We used Spearman’s rank correlation analysis due to some outliers and a slightly curved scatter plot.
We harnessed the R package ‘Hmisc’ [75] for bootstrapping and base R for correlation and regression analysis.

3. Results

3.1. Small Mammal Monitoring

The number of small mammals per 100 trap days over time in different sample areas and habitats is illustrated in Figure 2a. As the figure depicts, the peaks and depressions matched well between areas over time, with only a slight stochasticity between habitats within the same areas. This was confirmed by GLMM analysis, showing no significant differences in marginal mean ratios of the number of small mammals per 100 trap days among sample areas when accounting for an individual transect in a particular year (Table A2). There were observable differences in the relative abundances of small mammals between habitats (Figure 2b). GLMM analysis revealed that meadow habitats had significantly lower abundances than any forest age group but with no differences between age groups (Table A3). Comparison of fertility groups revealed that meadows had significantly lower abundance and that forests on fertile soils had significantly higher abundance after accounting for multiple comparisons (Table A4). In every comparison, GLMM including the hierarchical random intercept of transect in year and only the variable of interest in the fixed part was the best-performing model (with the lowest AICc values).
We pooled all the results to conduct a population change analysis because there were no important differences between the areas. All three models suggest statistically significantly declining populations (total number of small mammals: S = 0.9671 ± 0.0083, p = 0.0007; Microtus voles: S = 0.9306 ± 0.0167, p = 0.0005; bank voles (Clethrionomys glareolus): S = 0.9706 ± 0.0128, p = 0.0325). Yearly abundance indices are shown in Figure 3a–c.
Some cyclicity in populations was apparent in the late 20th century and early 21st century, but it lost amplitude later (Figure 3a–c). Since 2008, the total number of small mammals (Figure 3a), as well as Microtus voles (Figure 3b), has remained relatively stable and at low density level. The same happened with bank voles (Figure 3c) earlier—around 2003. Graphical evaluation indicated a certain degree of matching between vole population indices in depression years, e.g., 1997, 2003, and, to some extent, 2014, with overall moderate (r = 0.5604) and statistically significant (t (24) = 3.3146, p = 0.0029) correlation. Autocorrelation function analysis showed no temporal cyclicity in any species group (Figure 3d).

3.2. Owl Breeding Season Diet

In total, 164 STRALU samples from 86 different locations covering 23 years, 56 STRURA samples from 38 different locations covering 15 years, 24 ASIOTU samples from 21 different locations covering 9 years, 7 GLAPAS samples from 7 different locations covering 7 years, and 2 AEGFUN samples from 2 different locations covering 2 years (Figure 4, Table S1) were analyzed. The description of annual food composition per species is presented in Supplementary Materials, Table S1.
The overall average food niche breadth (FNB) of STRALU was 5.125 (95% bootstrapped confidence interval (bCI), 4.867–5.423). FNB increased significantly (β = 0.0840 ± 0.0198; t (129.2401) = 4.249; p < 0.0001) from 1992 to 2016. There were notable differences between samples in any given year (Figure 5a), and inclusion of nest ID as a random effect provided some help in dispersion taming (LMM: R2conditional = 0.115, R2marginal = 0.101, ICC = 0.015), indicating some degree of territory-specific variability. Although the explained variances were low, we found a statistically significant negative effect of each of the small mammal population indices on FNB (Table 2). The proportion of voles (both groups) in owl prey correlated positively with the vole abundance indices, but STRALU showed preference for Microtus voles, as their abundance index had a significant negative correlation with the bank vole proportion in prey (Table 3). On average, Microtus voles accounted for 15.55% of biomass, whereas bank voles and voles in total accounted for 5.51% and 31.65% of biomass, respectively.
STRURA also showed large variability in the diet, as overall FNB was 4.485 (95% bCI, 4.201–4.758). We found a significant (β = 0.0499 ± 0.0227; t (54) = 2.194; p = 0.0325) increase in FNB from 1994 to 2016 (Figure 5). The overall variability was lower than in STRALU, but no nest-specific intercepts were found to improve the model, and LM could explain only about 6.5% of the variance (R2adj. = 0.06486). We found no correlation between FNB and small mammal population indices (Table 2). The species showed a strong preference for Microtus voles, the proportion of which in prey correlated positively with its abundance index, while higher abundance in nature led to a lower proportion of bank voles in prey (Table 3). We suggest preference as a reason for the negative correlation of the bank vole abundance index with its proportion in prey because both vole abundance indices were correlated (Table 3). On average, Microtus voles accounted for 15.10% of biomass, whereas bank voles and voles in total accounted for 6.07% and 31.78% of biomass, respectively.
The overall FNB of GLAPAS was 3.526 (95% bCI 2.355–4.756), showing a temporal increase (Figure 5). However, this increase was not found to be statistically significant (LM: β = 0.0420 ± 0.1136, t (5) = 0.369, p = 0.727; R2adj. = −0.1681), probably due to high intersample variability and small overall sample size. We found no correlation between FNB and small mammal population indices (Table 2). Results of prey proportion and relative abundance in nature were similar to those for STRURA, but the preference for Microtus voles was greater, serving as better explanator of the bank vole proportion in prey (Table 3). On average, Microtus voles accounted for 10.03% of biomass, whereas bank voles and voles in total accounted for 11.31% and 29.46% of biomass, respectively.
ASIOTU’s FNB was only 1.629 (95% bCI, 1.429–2.002), and LM did not suggest any temporal change (β = 0.0010 ± 0.0176, t (22) = 0.059, p = 0.953; R2adj. = −0.0453). We found no correlation between FNB and small mammal population indices (Table 2). The species showed a strong preference for Microtus voles over bank voles, and Microtus vole abundance in nature significantly correlated with its proportion in prey (Table 3). On average, Microtus voles accounted for 71.18%, whereas bank voles and voles in total accounted for 4.64% and 87.13% of biomass, respectively.
The two analyzed samples of AEGFUN had FNBs of 1.588 and 4.318, respectively (Table S1). On average, Microtus voles accounted for 62.89% of biomass, whereas bank voles and voles in total accounted for 13.81% and 76.69% of biomass, respectively.

3.3. Owl Population Change

Our results differed between species when comparing owl population changes before and after small mammal cycle depletion (Figure 6, Table 4).
STRALU, with an overall (1990–2021) stable population (S = 1.002 ± 0.005), exhibited no significant difference in population trends before and after depletion (Table 4). Although the population experienced a considerable depression during the period of 2010–2012, it has since recovered (Figure 6a).
We obtained similar results for ASIOTU, with an overall (1990–2021) stable population (S = 0.992 ± 0.010) and no significant difference between the periods (Table 4). However, a visual extension of the trend since mammal population depletion suggested a decline that might be obscured by fluctuating population (Figure 6c).
The results of STRURA were different; although overall (1993–2021), the population was classified as stable (S = 1.014 ± 0.012), there was a significant difference in trends (Table 4). This species had a strongly increasing population before 2004 and a steep decline since small mammal depletion (Figure 6b).
Population change information for GLAPAS was available only since the depletion, and its overall (2004–2021) negative population trend (S = 0.965 ± 0.017) was similar to that observed in 2004–2016, reflecting a significant decline (Figure 6e, Table 4).
The results for AEGFUN are interesting, as the overall (1990–2021) population had a steep decline (S = 0.934 ± 0.020) that fit with estimated yearly indices (Figure 6d). Nevertheless, the difference between slopes “before” and “after” depletion was significant (Table 4) and suggests a steeper decline during the pronounced small mammal dynamics than since the depletion of cycles. However, visually extending the trajectory of “after” revealed a pattern to that in the “before” period; thus, the difference could be an artefact due to an increased influence of some years.

3.4. Owl Breeding Performance

On average, STRALU had 2.32 (95% bCI 2.23–2.40; n = 934) fledglings per successful nest. Despite the appearance of reduced breeding performance since 2004 (Figure 7a), these differences were not statistically significant (Table 5). The overall trend of breeding performance was insignificant (β: −0.0041 ± 0.0088, t (26) = −0.467, p = 0.645; R2adj. = −0.0298; F (1;26) = 0.2179, p = 0.6445). We did not find a correlation with the small mammal abundance indices in the year of breeding or the previous autumn (Table 6).
On average, STRURA had 1.69 (95% bCI 1.58–1.80; n = 280) fledglings per successful nest. Despite the appearance of some differences in trends of breeding performance before and after small mammal cycle dampening (Figure 7b), they were not statistically significant (Table 5). The overall trend of breeding performance was insignificant (β: −0.0014 ± 0.0122, t (24) = 0.112, p = 0.912; R2adj. = −0.0411; F (1;26) = 0.0125, p = 0.912). We did not find correlations with the small mammal abundance indices in the year of breeding or the previous autumn (Table 6).
The average number of ASIOTU fledglings was 2.54 (95% bCI 2.38–2.72; n = 189) per successful nest. However, this population parameter declined over time (Figure 7c) by approximately one fledgling in 16 years (β: −0.0627 ± 0.0297, t (22) = −2.109, p = 0.0466; R2adj. = 0.1303; F (1;22) = 4.446, p = 0.0466). We did not find a correlation with the small mammal abundance in the year of breeding or the previous autumn (Table 6).
The average number of BUBBUB fledglings was 2.26 (95% bCI 2.04–2.46; n = 81) per successful nest. This population parameter declined over time (Figure 7c) by approximately one fledgling in 22 years (β: −0.0450 ± 0.0137, t (20) = −3.293, p = 0.0036; R2adj. = 0.3191; F (1;20) = 10.840, p = 0.00036). We found positive correlations with the abundance index of pooled small mammals in the breeding season and the previous autumn, with a stronger effect of the latter (Table 6). The effect of the Microtus vole abundance index in the previous autumn was also statistically significant and positive (Table 6).
We had too few reliable observations of GLAPAS and AEGFUN breeding performance for analysis; therefore, we provide only a description of the average values: r = 2.75 (95% bCI 2.00–3.50; n = 8) and 1.75 (95% bCI 0.50–2.75; n = 4), respectively.

4. Discussion

4.1. Small Mammal Monitoring

Overall, the small mammal densities and trapping indices in our study (Figure 2) were similar to findings in the neighboring countries of Estonia [16,76] and Lithuania [77,78,79]. We found that meadow habitats had lower abundance of small mammals than forests (Table A3, Table A4). However, due to large within-class variation, no clear differences between age groups were found (Table A3). In Estonia [16] and Lithuania [78], increasing small mammal abundance has been recorded with increasing forest age in early meadow–forest succession. Additionally, studies from Finland [80] and Norway [81] reported that mature forests had the highest abundance of voles.
However, a study conducted in northern Sweden suggested that young stands have higher small mammal diversity and abundance if a large amount of felling remains is left [82]. Many authors have found that high vegetation complexity, habitat structural diversity, and abundance of coarse woody debris are important factors that can help to ensure high diversity and abundance of small mammal species in young stands and unmanaged habitats under natural succession [79,81,82,83,84,85].
The negative effects of intensive forestry have been found to be important at the landscape scale [83,84,85,86]. However, in a mosaic landscape, ecotones (with at least 100m buffer zone of habitat edges) have been found to contain the highest small mammal density and diversity [77,81,84].
Although the insufficient number of transects and occurrences in our study did not allow for statistical testing of forestry impacts, we found some declines in the small mammal numbers linked to forestry activities but unrelated to changes in other transects (Figure 2). Most of the small mammal monitoring transects were in intensively managed forests. However, two of the longer-term areas were in protected areas (Apasalas and Žūklis) and also showed dampening of the cycles, suggesting larger than local (or management) effects on the dampening of the cycles.
One of the most robust explanations of cyclicity was provided by Hanksi et al. [17] with further extensions for different systems (see [87] for overview). One of those extensions, modelling multispecies rodent assemblages, revealed transient dynamics that alternated between long time periods with cyclic and non-cyclic fluctuations [19]. These fluctuations were expected to cover relatively small spatial scales, yet the phenomenon of dampened cycles was more recently found to occur Europe-wide [25], suggesting broader environmental drivers, for example, climate change [22,23]. However, in some parts of Europe, the period of dampened vole population cycles has been shorter than in others, refuting the generality of the climate forcing hypothesis [26]. Our results also showed clearly dampened vole cycles in Latvia (Figure 3a–c). The fact that such a dampening of cycles has not been reported in neighboring counties [14,15,16] suggests some smaller-scale processes, such as those described by predator–prey models. Although Hanski et al.’s models were created for the Fennoscandian environment, their generality has also been shown in central and western Europe [13]. According to these models and previous studies (see [87] for an overview), generalist predators tend to stabilize rodent dynamics, and nomadic avian predators have a similar effect on rodents, although they also increase the regional synchrony, whereas specialist predators have been thought to maintain the fairly regular rodent cycles [19].

4.2. Numerical Response of Owls

We found a statistically significant relationship between the proportion of voles (Microtus and bank voles) in owl diet and their relative abundance indices in nature (Table 3). This means that although we drew on mammal abundance information from autumn, it was still able to represent their abundance in owl prey. It is known that small mammal densities increase during summer [52] and that spring counts represent winter survival and reproduction [88], but the relative value of the current year (spring or autumn) still represents part of the cycle in cyclic environments [41]. We found preference for Microtus voles in every analyzed owl species in terms of the proportion in owl diet; these voles also accounted for a higher biomass proportion than bank voles (Table 3). Generally, a higher proportion of Microtus voles than bank voles in owl diet can be related to different breeding biology of voles and dispersal between vole species groups and predator-escaping behavior (see [19] for an overview).

4.2.1. Long-Eared Owl

We found ASIOTU to have the narrowest FNB among the investigated species. The calculated values were slightly lower than in Lithuania [28], possibly due to the pooling of the results to genus level. This species is known to be a small-mammal specialist in Europe [27] with a high proportion of Microtus voles in their diet [28,89,90,91,92,93]. ASIOTU has been found to show strong functional responses of diet, breeding success, and dispersal to vole abundance [29,30,31,42]. It has even been suggested that species can adapt migratory behavior and breeding region selection during migration to account for vole abundance [94]. Moreover, this species may even exhibit repeated breeding attempts if the vole abundance is high [95,96].
Although the average breeding performance in Latvia was similar to the 2.94 ± 0.42 (μ ± SD; n = 1339) recorded in Finland [97] and 2.39 (n = 72) in the United Kingdom [95], we observed a significant decline in the number of fledglings per successful nest, i.e., more than one chick in three generations (5.7 years; [98]). The steepest decline occurred in the last two generations and matched the time of dampened populations of small mammals (Figure 7c). The declining breeding performance did not have an impact on the population change (2004–2016), but extension of the trend (2004–2021) showed a significant decline (β: −0.0530 ± 0.0191, t (16) = −2.772, p = 0.0136). We suspect that for a longer period of time, the ASIOTU population was supported by immigration of migrating individuals hatched elsewhere [29,90,99,100] and that a later decline implies a delayed response of returning individuals of Latvian origin.
Habitat composition and prey abundance have been found to be the most important factors shaping local ASIOTU populations [90,101,102,103] because the species shows no strong territorial defence and hunting grounds may largely overlap between neighboring pairs [104,105]. We found no correlation between the breeding performance and prey abundance indices in the year of breeding or the previous year (Table 6), likely due to selection of breeding territories with sufficient abundance of prey. This is supported by the knowledge of species benefitting from relatively small landscape elements, for example, flower strips [103].

4.2.2. Tengmalm’s Owl

Although we had only two samples of AEGFUN diet, its FNB suggested a high specialization, which was supported by a high proportion of voles in the diet. The observed proportions in Europe have shown a high importance of voles (overview in [32]), averaging 54.89% according to previous studies. This species has a strong functional response to vole abundance influencing habitat selection via hunting behavior [47,106,107,108], the timing of breeding and breeding performance [6,39], and survival [39,43]. Even with a certain degree of carry-over effect, the species has shown strong adaptability to fluctuating food conditions in terms of breeding performance [38].
The average long-term number of fledglings in Finland is 4.04 ± 0.62 (μ ± SD; n = 13,817) [38] and around 2 fledglings per successful nest in poor vole years [32,109]. The scarcely available data on the breeding performance in Latvia suggests that it is similar to that in vole depression years in Finland.
We found a steep decline of AEGFUN population throughout the studied period (Figure 6d), but it was slower with depleted population dynamics of small mammals (Table 4). We expected this to be an artefact of some better seasons or immigration from Fennoscandia and Russia [32] rather than an actual difference; therefore, we extended the period of analysis in the “after” group. Our results (β: −0.0618 ± 0.0074, t (16) = −8.302, p = <0.0001) showed a decline since 2004 closely matching the overall population decline and the slope of the period with pronounced population dynamics of small mammals. Some researchers have hypothesized potentially negative effects of increasing STRURA population on the AEGFUN population [32]. However, we did not find any AEGFUN as a prey of STRURA, although superpredation is known [34] and both species coexisted in the same study areas (authors’ personal observations). Furthermore, in central Europe, breeding in the proximity of STRURA has been found to protect AEGFUN against STRALU [110,111,112].
Population declines have also been reported in Finland, Sweden, and Estonia [113], suggesting larger-scale factors affecting the population. This species is a mature spruce and mixed forest specialist [47,101,108,114,115,116,117,118,119]. These are habitats with some of the highest densities of small mammals [79,80,81,82,83,84,85]. We consider the loss of species-specific habitats to be the most important factor in population decline, amplified by dampened dynamics of small mammal populations in Latvia. The forestry intensity, as measured by tree cover loss, is increasing in Latvia and, in particular, in priority sites for species conservation [120].

4.2.3. Eurasian Pygmy Owl

We found average level of specialization of GLAPAS and the preference for Microtus voles was strongest among the analyzed owl species (Table 3), although with a low proportion of voles in the diet. The vole proportion was similar to the breeding season diet in Finland [34] and in central Europe [121]. Masoero et al. found a strong numeric and functional response of GLAPAS to vole abundance in winter [33], suggesting not only age- and gender-specific preference for voles but also stronger migratory behavior during low-vole-density years in boreal Finland. During the years of higher vole population densities, breeding density and performance of GLAPAS also increase [6,122]. The dependency on voles has been found to be stronger in boreal than boreonemoral regions, with breeding both in low and peak vole years in the latter [40]. In the boreonemoral zone, the onset of breeding was later with no correlation with breeding performance, and the clutches were slightly smaller than in the boreal zone [40].
For the few documented records of successful breeding in Latvia, the values were markedly lower than 5.85 ± 0.55 (μ ± SD; n = 13,817) in Finland [97] and boreal Norway (6.9 ± 1.1) and somewhat lower than in boreonemoral Norway in vole crash years (3.7 ± 2.8) [40]. The difference relative to boreonemoral Norway indicates a possible cumulative effect of longer-term dampened population cycles of small mammals, which is supported by a declining GLAPAS population. The population in Estonia and Lithuania is increasing [113], but it is declining in Latvia (Figure 6e and Table 4) and in Finland [97,113]. It has to be noted that only Finland and Latvia were able to provide analytical assessment of the population (type: interval) in the last report on the Article 12 of the Birds Directive [113]; therefore, it cannot be ruled out, that the increase in the other Baltic states is more based on increased survey efforts and knowledge than a genuine change. Although irruptions linked with low rodent availability occur from time to time [123], it is unlikely that Finland and Latvia are source populations for neighboring countries with declining populations themselves, despite the increased distribution of the species [51]. This is supported by similar patterns of yearly indices in Latvia (Figure 6e) and Finland [97] but with a steeper decline in Latvia.
GLAPAS is known to be a structurally rich, mature spruce and mixed forest specialist species during the breeding season [122,124,125,126,127,128,129,130,131], and clearcuts and logging have been shown to affect habitat suitability [131], as well as population size [132]. Structurally rich mature forests are habitats with some of the highest densities of small mammals [79,80,81,82,83,84,85,133]. Latvia and Finland are the countries in Europe with the highest forestry activity, even in protected areas [134]. We suspect the loss of species-specific habitats to be the most important factor in the population decline, amplified by the dampened dynamics of small mammal populations in Latvia. The forestry intensity, as measured by tree cover loss, is increasing in Latvia and, particularly in priority sites for species conservation [120].

4.2.4. Ural Owl

One of the highest and temporarily increasing FNB values was found in STRURA. We found a low proportion of voles in the diet of this species. This proportion, when compared by count, was lower than in Finland [34,35,135,136], similar to that in Belarus [137,138,139], and higher than that in Slovenia [35]. When comparing the food niche as a whole, STRURA diet in Latvia was found to be similar to that in Finland during the low vole phase [35]. Although this species is known to be a generalist predator, a strong functional response to vole abundance has been proven in Fennoscandia, ranging from the timing of breeding and breeding performance [3,6,39,41,140] to winter survival [3,6,39,41,45] and even demonstrating a carry-over effect from the previous year [48] and a change in behavior [46,141,142].
Not only the food niche but also the breeding performance of STRURA in Latvia was similar to that in Finland in bad vole years. In Latvia, we observed, on average, 1.69 (95% bCI, 1.58–1.80; n = 280) fledglings per successful nest and no temporal trend. The corresponding overall value in Finland (1986–2016) is 2.59 (±0.43 SD, n = 18,901; [98]) and between 1.3 and 2 [3,41] in bad vole years, roughly matching our results. Given the strong numerical response to voles, we expected a declining trend in breeding performance, but we did not find it. We consider this an example of strong parental investment [142,143,144] as reflected by adjustments in hunting activity and possibly habitat selection [101], demonstrating the high plasticity of the species [35]. As Figure 2 shows, even with dampened small mammal cycles, there are habitats with high prey abundance, allowing prey to meet the demands of the young. The size of nest boxes in Latvia is similar to that in Finland [35] and cannot be suspected as a reason for lower breeding performance.
Increasing STRURA populations and expanding range, even increasing the niche of utilized habitats, was observed in many parts of Europe during the first decade of the 21st century [145,146,147,148]. This overlaps with the increase in Latvia and breeding occurrences in a mosaic landscape [101]. Given the extent of population increase, some unknown large-scale factors are most likely the explanation. Nevertheless, in Latvia, the period of steep decline in the species population overlapped with the dampening of the small mammal population dynamics. We consider the relative abundance of small mammals to be an important collider for overall habitat change, as species ecological niche analysis in Latvia suggests strong dependency on large forest massifs with dominance of mature forests and only some openings [102]. These are habitats found to be important for the same species elsewhere [146,149,150]. Although the range is still expanding in Latvia [151], the overall population size is declining [113]. The forestry intensity, as measured by tree cover loss, is increasing in Latvia, in particular in priority sites for species conservation [120]. We consider this as an argument for the conservation of mature forests important for the species and holding higher densities of its main prey, i.e., small mammals [81].

4.2.5. Tawny Owl

The highest FNB value was found in STRALU with a relatively low proportion of voles in the diet. The average FNB value was slightly lower than in Lithuania [28]. We observed a temporal increase in FNB, which was similar to the observation in Lithuania, with a declining proportion of Microtus voles [36]. This species is known to be a generalist [152]. Its food composition can considerably vary between breeding regions within the same year and between years in the same breeding territory [28,36,152,153,154,155,156,157,158,159,160]. Nevertheless, in the cyclic environment of Fennoscandia, a strong numerical response to vole abundance has been reported, including the timing of breeding [6,160], breeding performance [3,5,6,39,160,161], and winter survival [3,39].
Both the population change and breeding performance of STRALU were stable and showed no differences relative to pronounced and dampened vole cycles. Breeding performance was lower than 3.26 ± 0.41 (μ ± SD; n = 9668) in Finland, where the population was also stable [97], as well as in Lithuania, where an increasing trend of breeding performance (2002–2014) was observed and co-occurred with a decline in the number of breeding pairs [36]. We consider the relatively low breeding performance in Latvia to be related to the high population density (estimated to be around 16,604 in Latvia and below 4000 in Lithuania) [113]. The observed depression of the STRALU population from 2010 to 2012 in Latvia partially matched with Lithuania [36]. This was likely the consequence of two consecutive snow-rich winters, with multiple freeze—thaw events forming ice sheets in snow cover—factors reducing species survival [45,162,163]. This event did not affect breeding performance, and the population recovered quickly.
We think that the quick population recovery and overall stable breeding performance, even with increasing FNB values, was possibly due to the breeding habitat availability. Although this species is a well known generalist breeding from cities to large forest massifs in more southern latitudes [34], in the boreonemoral region, STRALU has been found to prefer forest edges over the interior [101,164]. With increasing forestry and forest fragmentation, more suitable landscapes for this species are created [101], probably overwhelming the negative effects of depleted small mammal populations.

4.2.6. Eagle Owl

The largest European owl species, BUBBUB, is known to be a generalist predator, with the proportion of rodents in its diet ranging from 0 to 97.7%, with average of 49.7% among 182 studies (overview in [37]). We do not have reliable information on the diet of BUBBUB in Latvia, but during the ringing of the young, many bird feathers are found, as most of the known breeding sites are in close proximity to waterfowl lakes and landfills [101]. This species is resident with no known seasonal migrations in Europe [37,165], and breeding dispersal occurs mostly due to the loss of a mate [37]. We speculate that the BUBBUB population in Latvia depends highly on voles, at least in winter, when bodies of open water are typically frozen and most waterfowl and gulls in Latvia have moved to wintering areas [151]. This is supported by a study conducted in Finland evaluating the robustness of the alternative prey hypothesis for BUBBUB [166]. A correlation was found between vole abundance in nature and their proportion in the diet, and the proportion of alternative prey was found to be nearly independent of its abundance in the field [166]. Several other studies highlight the high proportion of voles in the BUBBUB diet [167,168,169].
The overall average number of fledglings per successful nest in Latvia was similar to that in Europe—around two (overview in [37]); however, we observed a declining trend, with a loss of more than one fledgling in two generations (generation length is 12.1 years [170]). We found that breeding performance was correlated with the abundance index of small mammals in field, and the correlation with the value from the previous autumn was stronger (Table 6). Other studies found that BUBBUB pairs with a diet based on high-value foods (rabbits and rodents) have comparatively larger broods and breed earlier [171,172], and higher reproductive productivity was associated with a higher proportion of the main prey (rats and rabbits in Spain) in the diet [173]. Our data do not lend themselves to examining such a relationship between breeding performance and diet. Nevertheless, we assume that the negative effect of a reduced abundance of small mammals highlights a carry-over effect, influencing adult fitness in spring and thereby reducing the breeding performance. This phenomenon is well known in STRURA [41,48,143,174,175] and has proven to be of increasing importance with the size of an owl species in Finland [6]. Ecological niche analysis in Latvia also suggests the importance of habitats with higher vole abundance [81] for BUBBUB [101].
Although we had annual data on a limited number of nests, they formed an important part of the whole population estimated at around 24 breeding territories, indicating a declining national population trend [113]. We consider our findings of a possible carry-over effect to be important in species conservation and to be linked to population decline via reduced breeding performance, as well as reduced winter survival, as both should be related via fitness, although this relationship needs to be studied more directly. Nevertheless, we consider the conservation of habitats important for breeding and winter feeding, together with nest site protection from ground predators, which is necessary to reduce the effects of dampened population dynamics of small mammals.

5. Conclusions

  • Small mammal relative abundance indices have shown depleted population cycles since approx. 2004. This has impacted the breeding performance, food niche breadth, and population trends of owl species to various degrees depending on the particular species;
  • The number of ASIOTU fledglings has declined since the depletion of small mammal populations. The population size of the species declined later and was significant for the period from 2004 to 2021. ASIOTU is the most specialized of the analyzed owl species in terms of the proportion of voles in the diet;
  • The breeding performance of the three forest specialist species AEGFUN, GLAPAS, and STRURA in Latvia was similar to vole depression years in the boreal and boreonemoral regions;
  • Populations of GLAPAS and AEGFUN declined in Latvia and showed no difference compared to periods with pronounced or depleted population dynamics of small mammals. In contrast, the population of STRURA has shown a significant decline since rodent depression. We suggest the depletion of the small mammal population dynamics to be an important negative contributing factor to more important effects of forestry, although the impact of forestry needs to be investigated further;
  • Neither the breeding performance nor population size of STRALU changed between the compared periods with pronounced and depleted population dynamics of small mammals. This suggests a strong plasticity of the species, as food niche breadth was temporarily increased;
  • We found evidence that suggests the dependency of BUBBUB on voles via a carry-over effect. The breeding performance of BUBBUB was significantly correlated with the abundance indices of small mammals in nature in the previous autumn.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/life13020572/s1, Table S1: Description of owl diet per year; Table S2: Population indices of small mammals; Table S3: Population indices of owls; Table S4: Description of owl breeding performance.

Author Contributions

Conceptualization, A.A. (Andris Avotins) and A.A. (Ainars Aunins); methodology, A.A. (Andris Avotins), A.A.S., and A.A. (Ainars Aunins); validation, A.A. (Ainars Aunins), A.A.S., and V.Ķ.; formal analysis, A.A. (Andris Avotins); data curation, A.A. (Andris Avotins); writing—original draft preparation, A.A. (Andris Avotins); writing—review and editing, A.A. (Ainars Aunins), V.Ķ., and A.A.S.; visualization, A.A. (Andris Avotins); supervision, A.A. (Ainars Aunins); funding acquisition, A.A.S., A.A. (Ainars Aunins), A.A. (Andris Avotins), and V.Ķ. All authors have read and agreed to the published version of the manuscript.

Funding

Small mammal monitoring in 1991–2011 was funded by the Teiči state reserve administration and the Nature Conservation Agency. Monitoring of owls in 1991–2007 was partially funded by the Teiči state reserve administration and the Nature Conservation Agency. Monitoring of owls in 2007–2009 was funded by the Latvian Environment, Geology, and Meteorology Center, and since 2014, it has been funded by the Nature Conservation Agency. Most of the owl prey remains analysis and small mammal monitoring conducted in 2016 was funded by the Administration of Latvian Environmental Protection Fund (1-08/158/2016).

Institutional Review Board Statement

Nest visits were conducted by professional ringers selected by the Latvian Ringing Centre and approved by the Nature Conservation Agency and do not require ethical review. The monitoring of the small mammals with snap-traps was carried out according to the regulation of the Cabinet of the Ministers from 8 January 2019 “Regulation on the protection of use of animals for scientific purposes” [Zinātniskiem mērķiem izmantojamo dzīvnieku aizsardzības noteikumi].

Informed Consent Statement

Not applicable.

Data Availability Statement

The aggregated (per species, per year) data presented in this study are available in the Supplementary Materials of this article. The raw data are not publicly available as we are not allowed to share raw monitoring information or nest locations.

Acknowledgments

We would like to thank the Latvian Ringing Centre and all the volunteer ringers for the availability of ringing data, particularly J. Lipsbergs, A. Kalvāns, G. Grandāns, O. Keišs, V. Boļšakovs, T. Ķeris, J. Jankavs, A. Freibergs, and P. Daknis. We also acknowledge the nature observation portal management team and volunteer citizen scientists reporting owl observations and breeding performance. We are grateful to people who provided us with materials of owl prey remains: I. Bauga, U. Bergmanis, P. Daknis, K. Funts, M. Funta, G. Grandāns, V. Ignatjevs, I. Jakovļevs, O. Keišs, J. Ķuze, U. Ļoļāns, I. Mārdega, M. Medne, J. Priednieks, E. Račinskis, R. Rekmanis, M. Rozenfelde, M. Strazds, D. Ūlands, and M. Zilgalvis. We are also grateful to those who helped with dissection: L. Dambeniece, V. Ignatjevs, M. Rozenfelde, D. Ūlands, and M. Zilgalvis. We thank A. Pupila, G. Done, G. Grandāns, and U. Bergmanis for sharing their historical small mammal trapping data and locations of transects.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

Appendix A

Table A1. Minimum legal rotation ages per site quality class in main tree species. Currently, there is no minimum rotation age in grey alder. We used 35 years, as it is the age of the youngest stand registered as “full grown” in FSR.
Table A1. Minimum legal rotation ages per site quality class in main tree species. Currently, there is no minimum rotation age in grey alder. We used 35 years, as it is the age of the youngest stand registered as “full grown” in FSR.
Dominant Tree SpeciesHighest QualityMedium QualityLowest Quality
Oaks101121121
Pines and larches101101121
Spruces, ashes, limes, elms, and maples818181
Birches717151
Black alders717171
Aspens414141
Table A2. Marginal mean comparison (tests were conducted on the log scale; p adjusted by Tukey’s method) of the number of small mammals per 100 trap days between sample areas. The model with a random intercept per transect and main effects of sample area and year was considered the best *.
Table A2. Marginal mean comparison (tests were conducted on the log scale; p adjusted by Tukey’s method) of the number of small mammals per 100 trap days between sample areas. The model with a random intercept per transect and main effects of sample area and year was considered the best *.
Sample Areas ContrastedRatio ± SEz-Ratiop-Value
Apsalas/Augstroze0.704 ± 0.459−0.5380.9983
Apsalas/Kaunata0.468 ± 0.307−1.1570.9099
Apsalas/Ķemeri1.816 ± 1.5870.6820.9936
Apsalas/Lubāna0.538 ± 0.35−0.9510.964
Apsalas/Pļaviņas0.566 ± 0.371−0.8680.9772
Apsalas/Žūklis1.329 ± 1.1760.3220.9999
Augstroze/Kaunata0.665 ± 0.236−1.1480.9132
Augstroze/Ķemeri2.578 ± 1.7521.3930.8057
Augstroze/Lubāna0.764 ± 0.265−0.7740.9874
Augstroze/Pļaviņas0.804 ± 0.286−0.6140.9964
Augstroze/Žūklis1.887 ± 1.3080.9170.9701
Kaunata/Ķemeri3.877 ± 2.651.9830.4257
Kaunata/Lubāna1.149 ± 0.4080.3930.9997
Kaunata/Pļaviņas1.209 ± 0.4390.5220.9985
Kaunata/Žūklis2.839 ± 1.9781.4980.7465
Ķemeri/Lubāna0.296 ± 0.201−1.790.5547
Ķemeri/Pļaviņas0.312 ± 0.213−1.7050.6123
Ķemeri/Žūklis0.732 ± 0.663−0.3440.9999
Lubāna/Pļaviņas1.052 ± 0.3730.1421
Lubāna/Žūklis2.47 ± 1.7111.3050.8498
Pļaviņas/Žūklis2.349 ± 1.6331.2280.8834
* AICc = 1603.663, R2conditional = 0.877, R2marginal = 0.177, ICC = 0.851.
Table A3. Marginal mean comparison (tests were conducted on the log scale; p adjusted by Tukey’s method) of the number of small mammals per 100 trap days between forest age classes. The model with a random intercept per transect and main effects of forest age class and year was considered the best *.
Table A3. Marginal mean comparison (tests were conducted on the log scale; p adjusted by Tukey’s method) of the number of small mammals per 100 trap days between forest age classes. The model with a random intercept per transect and main effects of forest age class and year was considered the best *.
Forest Age Groups ContrastedRatio ± SEz-Ratiop-Value
Young/Meadow7.078 ± 2.7724.997<0.0001
Young/Medium1.458 ± 0.4421.2440.599
Young/Old1.627 ± 0.4771.6610.3446
Meadow/Medium0.206 ± 0.073−4.4510.0001
Meadow/Old0.23 ± 0.079−4.2570.0001
Medium/Old1.116 ± 0.270.4520.9692
* AICc = 1579.355, R2conditional = 0.884, R2marginal = 0.409, ICC = 0.804.
Table A4. Marginal mean comparison (tests were conducted on the log scale; p adjusted by Tukey’s method) of the number of small mammals per 100 trap days between soil richness classes. The model with a random intercept per transect and main effects of soil fertility class and year was considered the best *.
Table A4. Marginal mean comparison (tests were conducted on the log scale; p adjusted by Tukey’s method) of the number of small mammals per 100 trap days between soil richness classes. The model with a random intercept per transect and main effects of soil fertility class and year was considered the best *.
Soil Fertility Classes ContrastedRatio ± SEz-Ratiop-Value
Meadow/Fertile0.153 ± 0.047−6.086<0.0001
Meadow/Poor0.309 ± 0.098−3.7170.0006
Fertile/Poor2.026 ± 0.3973.6040.0009
* AICc = 1572.443, R2conditional = 0.879, R2marginal = 0.529, ICC = 0.743.

References

  1. McCaffery, R.; Jenkins, K.J.; Cendejas-Zarelli, S.; Happe, P.J.; Sager-Fradkin, K.A. Small Mammals and Ungulates Respond to and Interact with Revegetation Processes Following Dam Removal. Food Webs 2020, 25, e00159. [Google Scholar] [CrossRef]
  2. Moorhead, L.C.; Souza, L.; Habeck, C.W.; Lindroth, R.L.; Classen, A.T. Small Mammal Activity Alters Plant Community Composition and Microbial Activity in an Old-Field Ecosystem. Ecosphere 2017, 8, e01777. [Google Scholar] [CrossRef]
  3. Saurola, P.; Francis, C. Towards Integrated Population Monitoring Based on the Fieldwork of Volunteer Ringers: Productivity, Survival and Population Change of Tawny Owls Strix aluco and Ural Owls Strix uralensis in Finland. Bird Study 2018, 65, S63–S76. [Google Scholar] [CrossRef]
  4. Karell, P.; Ahola, K.; Karstinen, T.; Zolei, A.; Brommer, J.E. Population Dynamics in a Cyclic Environment: Consequences of Cyclic Food Abundance on Tawny Owl Reproduction and Survival. J. Anim. Ecol. 2009, 78, 1050–1062. [Google Scholar] [CrossRef]
  5. Solonen, T. Breeding of the Tawny Owl Strix aluco in Finland: Responses of a Southern Colonist to the Highly Variable Environment of the North. Ornis Fenn. 2005, 82, 97–106. [Google Scholar]
  6. Lehikoinen, A.; Ranta, E.; Pietiäinen, H.; Byholm, P.; Saurola, P.; Valkama, J.; Huitu, O.; Henttonen, H.; Korpimäki, E. The Impact of Climate and Cyclic Food Abundance on the Timing of Breeding and Brood Size in Four Boreal Owl Species. Oecologia 2011, 165, 349–355. [Google Scholar] [CrossRef]
  7. Blomqvist, S.; Holmgren, N.; Åkesson, S.; Hedenström, A.; Pettersson, J. Indirect Effects of Lemming Cycles on Sandpiper Dynamics: 50 Years of Counts from Southern Sweden. Oecologia 2002, 133, 146–158. [Google Scholar] [CrossRef]
  8. Summers, R.W. Breeding Production of Dark-Billed Brent Geese Branta bernicla bernicla in Relation to Lemming Cycles. Bird Study 1986, 33, 105–108. [Google Scholar] [CrossRef]
  9. Root-Bernstein, M.; Fierro, A.; Armesto, J.; Ebensperger, L.A. Avian Ecosystem Functions Are Influenced by Small Mammal Ecosystem Engineering. BMC Res. Notes 2013, 6, 549. [Google Scholar] [CrossRef]
  10. Sundell, J.; Huitu, O.; Henttonen, H.; Kaikusalo, A.; Korpimäki, E.; Pietiäinen, H.; Saurola, P.; Hanski, I. Large-Scale Spatial Dynamics of Vole Populations in Finland Revealed by the Breeding Success of Vole-Eating Avian Predators. J. Anim. Ecol. 2004, 73, 167–178. [Google Scholar] [CrossRef]
  11. Lindén, H. Latitudinal Gradients in Predator-Prey Interactions, Cyclicity and Synchronism in Voles and Small Game Populations in Finland. Oikos 1988, 52, 341–349. [Google Scholar] [CrossRef]
  12. Hansson, L.; Henttonen, H. Gradients in Density Variations of Small Rodents: The Importance of Latitude and Snow Cover. Oecologia 1985, 67, 394–402. [Google Scholar] [CrossRef]
  13. Lambin, X.; Bretagnolle, V.; Yoccoz, N.G. Vole Population Cycles in Northern and Southern Europe: Is There a Need for Different Explanations for Single Pattern? J. Anim. Ecol. 2006, 75, 340–349. [Google Scholar] [CrossRef]
  14. Balčiauskas, L.; Balčiauskienė, L. Long-Term Changes in a Small Mammal Community in a Temperate Zone Meadow Subject to Seasonal Floods and Habitat Transformation. Integr. Zool. 2022, 17, 443–455. [Google Scholar] [CrossRef]
  15. Balčiauskas, L.; Balčiauskienė, L. Small Mammal Diversity Changes in a Baltic Country, 1975–2021: A Review. Life 2022, 12, 1887. [Google Scholar] [CrossRef]
  16. Väli, Ü.; Tõnisalu, G. Community- and Species-Level Habitat Associations of Small Mammals in a Hemiboreal Forest-Farmland Landscape. Ann. Zool. Fenn. 2020, 58, 1–11. [Google Scholar] [CrossRef]
  17. Hanski, I.; Hansson, L.; Henttonen, H. Specialist Predator, Generalist Predator, and the Microtine Rodent Cycle. J. Anim. Ecol. 1991, 60, 353–367. [Google Scholar] [CrossRef]
  18. Henttonen, H.; Oksanen, T.; Jortikka, A.; Haukisalmi, V. How Much Do Weasels Shape Microtine Cycles in the Northern Fennoscandian Taiga? Oikos 1987, 50, 353. [Google Scholar] [CrossRef]
  19. Hanski, I.; Henttonen, H. Predation on Competing Rodent Species: A Simple Explanation of Complex Patterns. J. Anim. Ecol. 1996, 65, 220. [Google Scholar] [CrossRef]
  20. Hörnfeldt, B. Long-Term Decline in Numbers of Cyclic Voles in Boreal Sweden: Analysis and Presentation of Hypotheses. Oikos 2004, 107, 376–392. [Google Scholar] [CrossRef]
  21. Hörnfeldt, B.; Hipkiss, T.; Eklund, U. Fading out of Vole and Predator Cycles? Proc. R. Soc. B Biol. Sci. 2005, 272, 2045–2049. [Google Scholar] [CrossRef] [PubMed]
  22. Steen, H.; Ims, R.A.; Sonerud, G.A.; Ecology, S.; Dec, N. Spatial and Temporal Patterns of Small-Rodent Population Dynamics at a Regional Scale. Ecology 1996, 77, 2365–2372. [Google Scholar] [CrossRef]
  23. Kausrud, K.L.; Mysterud, A.; Steen, H.; Vik, J.O.; Østbye, E.; Cazelles, B.; Framstad, E.; Eikeset, A.M.; Mysterud, I.; Solhøy, T.; et al. Linking Climate Change to Lemming Cycles. Nature 2008, 456, 93–97. [Google Scholar] [CrossRef] [PubMed]
  24. Ims, R.A.; Henden, J.A.; Killengreen, S.T. Collapsing Population Cycles. Trends Ecol. Evol. 2008, 23, 79–86. [Google Scholar] [CrossRef]
  25. Cornulier, T.; Yoccoz, N.G.; Bretagnolle, V.; Brommer, J.E.; Butet, A.; Ecke, F.; Elston, D.A.; Framstad, E.; Henttonen, H.; Hörnfeldt, B.; et al. Europe-Wide Dampening of Population Cycles in Keystone Herbivores. Science 2013, 340, 63–66. [Google Scholar] [CrossRef]
  26. Brommer, J.E.; Pietiäinen, H.; Ahola, K.; Karell, P.; Karstinenz, T.; Kolunen, H. The Return of the Vole Cycle in Southern Finland Refutes the Generality of the Loss of Cycles through “Climatic Forcing”. Glob. Change Biol. 2010, 16, 577–586. [Google Scholar] [CrossRef]
  27. Birrer, S. Synthesis of 312 Studies on the Diet of the Long-Eared Owl Asio otus. Ardea 2009, 97, 615–624. [Google Scholar] [CrossRef]
  28. Balčiauskienė, L.; Jovaišas, A.; Naruševičius, V.; Petraška, A.; Skuja, S. Diet of Tawny Owl (Strix aluco) and Long-Eared Owl (Asio otus) in Lithuania as Found from Pellets. Acta Zool. Lith. 2006, 16, 37–45. [Google Scholar] [CrossRef]
  29. Village, A. The Diet and Breeding of Long-Eared Owls in Relation to Vole Numbers. Bird Study 1981, 28, 215–225. [Google Scholar] [CrossRef]
  30. Tome, D. Functional Response of the Long-Eared Owl (Asio otus) to Changing Prey Numbers: A 20-Year Study. Ornis Fenn. 2003, 80, 63–70. [Google Scholar]
  31. Tome, D. Changes in the Diet of Long-Eared Owl Asio otus: Seasonal Patterns of Dependence on Vole Abundance. Ardeola 2009, 56, 49–56. [Google Scholar]
  32. Korpimäki, E.; Hakkarainen, H. The Boreal Owl: Ecology, Behaviour, and Conservation of a Forest-Dwelling Predator; Cambridge University Press: Cambridge, UK, 2012; ISBN 9781107425323. [Google Scholar]
  33. Masoero, G.; Laaksonen, T.; Morosinotto, C.; Korpimäki, E. Age and Sex Differences in Numerical Responses, Dietary Shifts, and Total Responses of a Generalist Predator to Population Dynamics of Main Prey. Oecologia 2020, 192, 699–711. [Google Scholar] [CrossRef]
  34. Mikkola, H. Owls of Europe; A.D. & T. Poyser: Calton, UK, 1983. [Google Scholar]
  35. Vrezec, A.; Saurola, P.; Avotins, A.; Kocijančič, S.; Sulkava, S. A Comparative Study of Ural Owl Strix uralensis Breeding Season Diet within Its European Breeding Range, Derived from Nest Box Monitoring Schemes. Bird Study 2018, 65, S85–S95. [Google Scholar] [CrossRef]
  36. Grašytė, G.; Rumbutis, S.; Dagys, M.; Treinys, R. Breeding Performance, Apparent Survival, Nesting Location and Diet in a Local Population of the Tawny Owl Strix aluco in Central Lithuania Over the Long-Term. Acta Ornithol. 2016, 51, 163–174. [Google Scholar] [CrossRef]
  37. Penteriani, V.; del Mar Delgado, M. The Eagle Owl; T & A.D. Poyser: London, UK, 2019; ISBN 978-1-4729-0066-1. [Google Scholar]
  38. Korpimäki, E.; Hakkarainen, H. Fluctuating Food Supply Affects the Clutch Size of Tengmalm’s Owl Independent of Laying Date. Oecologia 1991, 85, 543–552. [Google Scholar] [CrossRef]
  39. Saurola, P.; Francis, C.M. Estimating Population Parameters of Owls from Nationally Coordinated Ringing Data in Finland. Anim. Biodivers. Conserv. 2004, 27, 403–415. [Google Scholar]
  40. Solheim, R. Breeding Biology of the Pygmy Owl Glaucidium passerinum in Two Biogeographical Zones in Southeastern Norway. Ann. Zool. Fenn. 1984, 21, 295–300. [Google Scholar]
  41. Brommer, J.E.; Pietiäinen, H.; Kolunen, H. Reproduction and Survival in a Variable Environment: Ural Owls (Strix uralensis) and the Three-Year Vole Cycle. Auk 2002, 119, 544–550. [Google Scholar] [CrossRef]
  42. Tome, D. Post-Fledging Survival and Dynamics of Dispersal in Long-Eared Owls Asio otus. Bird Study 2011, 58, 193–199. [Google Scholar] [CrossRef]
  43. Hakkarainen, H.; Korpimäki, E.; Koivunen, V.; Ydenberg, R. Survival of Male Tengmalm’s Owls under Temporally Varying Food Conditions. Oecologia 2002, 131, 83–88. [Google Scholar] [CrossRef]
  44. Masoero, G.; Laaksonen, T.; Morosinotto, C.; Korpimäki, E. Climate Change and Perishable Food Hoards of an Avian Predator: Is the Freezer Still Working? Glob. Change Biol. 2020, 26, 5414–5430. [Google Scholar] [CrossRef] [PubMed]
  45. Pavón-Jordán, D.; Karell, P.; Ahola, K.; Kolunen, H.; Pietiäinen, H.; Karstinen, T.; Brommer, J.E. Environmental Correlates of Annual Survival Differ between Two Ecologically Similar and Congeneric Owls. Ibis 2013, 155, 823–834. [Google Scholar] [CrossRef]
  46. Kontiainen, P.; Pietiäinen, H.; Huttunen, K.; Karell, P.; Kolunen, H.; Brommer, J.E. Aggressive Ural Owl Mothers Recruit More Offspring. Behav. Ecol. 2009, 20, 789–796. [Google Scholar] [CrossRef]
  47. Santangeli, A.; Hakkarainen, H.; Laaksonen, T.; Korpimäki, E. Home Range Size Is Determined by Habitat Composition but Feeding Rate by Food Availability in Male Tengmalm’s Owls. Anim. Behav. 2012, 83, 1115–1123. [Google Scholar] [CrossRef]
  48. Brommer, J.E.; Karell, P.; Pietiäinen, H. Supplementary Fed Ural Owls Increase Their Reproductive Output with a One Year Time Lag. Oecologia 2004, 139, 354–358. [Google Scholar] [CrossRef]
  49. Ahti, T.; Hämet-Ahti, L.; Jalas, J. Vegetation Zones and Their Sections in Northwestern Europe. Ann. Bot. Fenn. 1968, 5, 169–211. [Google Scholar] [CrossRef]
  50. Peel, M.C.; Finlayson, B.L.; McMahon, T.A. Updated World Map of the Köppen-Geiger Climate Classification. Hydrol. Earth Syst. Sci. 2007, 11, 1633–1644. [Google Scholar] [CrossRef]
  51. Keller, V.; Herrando, S.; Voříšek, P.; Franch, M.; Kipson, M.; Milanesi, P.; Martí, D.; Anton, M.; Klvaňová, A.; Kalyakin, M.V.; et al. European Breeding Bird Atlas 2: Distribution, Abundance and Change; European Bird Census Council & Lynx Edicions: Barcelona, Spain, 2020. [Google Scholar]
  52. Pupila, A.; Bergmanis, U. Species Diversity, Abundance and Dynamics of Small Mammals in the Eastern Latvia. Acta Univ. Latv. 2006, 710, 93–101. [Google Scholar]
  53. Demongin, L. Identification Guide to Birds in the Hand; Laurent Demongin: Beauregard-Vendon, France, 2016; ISBN 9782955501900. [Google Scholar]
  54. Tauriņš, E. Latvijas Zīdītājdzīvnieki; Zvaigzne: Rīga, Latvia, 1982. [Google Scholar]
  55. Avotiņš, A.; Ķemlers, A. Monitoring of Owls in Latvia. Ring 1993, 15, 104–114. [Google Scholar]
  56. Avotiņš, A.; Graubics, G.; Ķemlers, A.; Krēsliņš, V.; Ķuze, J.; Ļoļāns, U. Number and Breeding Densities of Owls in Latvia: Studies in Sample Plots. Vogelwelt 1999, 120, 333–337. [Google Scholar]
  57. Avotiņš, A. Owl Census in Sample Plots near Metriena and Laudona, Eastern Latvia. Bird Census News 1999, 12, 52–62. [Google Scholar]
  58. Avotiņš, A. Changes of Number and Structure in Population of Tawny Owl (Strix aluco) in Sample Plots at Eastern Latvia (1990–1994). Popul. Greifvogel Eulenarten 1996, 3, 377–386. [Google Scholar]
  59. Avotiņš, A. Tawny Owl’s Territory Occupancy in Eastern Latvia. Bird Census News 2004, 13, 167–173. [Google Scholar]
  60. Avotiņš, A.; Reihmanis, J. Plēsīgo Putnu Monitorings; Latvijas Ornitoloģijas Biedrība: Riga, Latvia, 2020. [Google Scholar]
  61. R Core Team. R: A Language and Environment for Statistical Computing; R Core Team: Vienna, Austria, 2022. [Google Scholar]
  62. Wickham, H.; Averick, M.; Bryan, J.; Chang, W.; McGowan, L.; François, R.; Grolemund, G.; Hayes, A.; Henry, L.; Hester, J.; et al. Welcome to the Tidyverse. J. Open Source Softw. 2019, 4, 1686. [Google Scholar] [CrossRef]
  63. Pebesma, E. Simple Features for R: Standardized Support for Spatial Vector Data. R J. 2018, 10, 439–446. [Google Scholar] [CrossRef]
  64. Burnham, K.P.; Anderson, D.R. Model Selection and Multimodel Inference: A Practical Information-Theoretic Approach, 2nd ed.; Springer: Berling/Heidelberg, Germany, 2002; ISBN 978-0-387-22456-5. [Google Scholar]
  65. Bates, D.; Mächler, M.; Bolker, B.M.; Walker, S.C. Fitting Linear Mixed-Effects Models Using Lme4. J. Stat. Softw. 2015, 67, 1–48. [Google Scholar] [CrossRef]
  66. Lenth, R.V. Emmeans: Estimated Marginal Means, Aka Least-Squares Means. Available online: https://CRAN.R-project.org/package=emmeans (accessed on 26 July 2022).
  67. Bogaart, P.; van der Loo, M.; Pannekoek, J. Rtrim: Trends and Indices for Monitoring Data. Available online: https://CRAN.R-project.org/package=rtrim (accessed on 26 July 2022).
  68. Pannekoek, J.; Bogaart, P.; van der Loo, M. Models and Statistical Methods in Rtrim; Statistics Netherlands: Heerlen, The Netherlands, 2018; pp. 1–34.
  69. Pucek, Z. Keys to Vertebrates of Poland. Mammals; PWN—Polish Scientific Publishers: Warsaw, Poland, 1981. [Google Scholar]
  70. Smith, E.P. Niche Breadth, Resource Availability, and Inference. Ecology 1982, 63, 1675–1681. [Google Scholar] [CrossRef]
  71. Oberpriller, J.; de Souza Leite, M.; Pichler, M. Fixed or Random? On the Reliability of Mixed-Effects Models for a Small Number of Levels in Grouping Variables. Ecol. Evol. 2022, 12, e9062. [Google Scholar] [CrossRef]
  72. Gomes, D.G.E. Should I Use Fixed Effects or Random Effects When I Have Fewer than Five Levels of a Grouping Factor in a Mixed-Effects Model? PeerJ 2022, 10, e12794. [Google Scholar] [CrossRef]
  73. Pannekoek, J.; Strien, A. Van TRIM 3 Manual (TRends & Indices for Monitoring Data); Statistics Netherlands: Heerlen, The Netherlands, 2005; pp. 1–49. [Google Scholar]
  74. Van Strien, A.J.; Pannekoek, J.; Gibbons, D.W. Indexing European Bird Population Trends Using Results of National Monitoring Schemes: A Trial of a New Method. Bird Study 2001, 48, 200–213. [Google Scholar] [CrossRef]
  75. Harrell, F.E., Jr. Hmisc: Harrell Miscellaneous. Available online: https://CRAN.R-project.org/package=Hmisc (accessed on 26 July 2022).
  76. Scott, D.M.; Joyce, C.B.; Burnside, N.G. The Influence of Habitat and Landscape on Small Mammals in Estonian Coastal Wetlands. Est. J. Ecol. 2008, 57, 279–295. [Google Scholar] [CrossRef]
  77. Balčiauskas, L. Small Mammal Communities in the Fragmented Landscape in Lithuania. Acta Zool. Litu. 2006, 16, 130–136. [Google Scholar] [CrossRef]
  78. Balčiauskas, L.; Čepukienė, A.; Balčiauskienė, L. Small Mammal Community Response to Early Meadow–Forest Succession. For. Ecosyst. 2017, 4, 11. [Google Scholar] [CrossRef]
  79. Mažeikytė, R. Small Mammals in the Mosaic Landscape of Eastern Lithuania: Species Composition, Distribution and Abundance. Acta Zool. Litu. 2002, 12, 381–391. [Google Scholar] [CrossRef]
  80. Savola, S.; Henttonen, H.; Lindén, H. Vole Population Dynamics during the Succession of a Commercial Forest in Northern Finland. Ann. Zool. Fenn. 2013, 50, 79–88. [Google Scholar] [CrossRef]
  81. Wegge, P.; Rolstad, J. Cyclic Small Rodents in Boreal Forests and the Effects of Even-Aged Forest Management: Patterns and Predictions from a Long-Term Study in Southeastern Norway. For. Ecol. Manag. 2018, 422, 79–86. [Google Scholar] [CrossRef]
  82. Ecke, F.; Löfgren, O.; Sörlin, D. Population Dynamics of Small Mammals in Relation to Forest Age and Structural Habitat Factors in Northern Sweden. J. Appl. Ecol. 2002, 39, 781–792. [Google Scholar] [CrossRef]
  83. Suchomel, J.; Šipoš, J.; Košulič, O. Management Intensity and Forest Successional Stages as Significant Determinants of Small Mammal Communities in a Lowland Floodplain Forest. Forests 2020, 11, 1320. [Google Scholar] [CrossRef]
  84. Panzacchi, M.; Linnell, J.D.C.; Melis, C.; Odden, M.; Odden, J.; Gorini, L.; Andersen, R. Effect of Land-Use on Small Mammal Abundance and Diversity in a Forest-Farmland Mosaic Landscape in South-Eastern Norway. For. Ecol. Manag. 2010, 259, 1536–1545. [Google Scholar] [CrossRef]
  85. Carey, A.B.; Harrington, C.A. Small Mammals in Young Forests: Implications for Management for Sustainability. For. Ecol. Manag. 2001, 154, 289–309. [Google Scholar] [CrossRef]
  86. Bowman, J.C.; Sleep, D.; Forbes, G.J.; Edwards, M. The Association of Small Mammals with Coarse Woody Debris at Log and Stand Scales. For. Ecol. Manag. 2000, 129, 119–124. [Google Scholar] [CrossRef]
  87. Henttonen, H.; Gilg, O.; Ims, R.A.; Korpimäki, E.; Yoccoz, N.G. Ilkka Hanski and Small Mammals: From Shrew Metapopulations to Vole and Lemming Cycles. Ann. Zool. Fenn. 2017, 54, 153–162. [Google Scholar] [CrossRef]
  88. Hanski, I.; Korpimäki, E. Microtine Rodent Dynamics in Northern Europe: Parameterized Models for the Predator-Prey Interaction. Ecology 1995, 76, 840–850. [Google Scholar] [CrossRef]
  89. Tome, D. Diet Composition of the Long-Eared Owl in Central Slovenia—Seasonal-Variation in Prey Use. J. Raptor Res. 1994, 28, 253–258. [Google Scholar]
  90. Sergio, F.; Marchesi, L.; Pedrini, P. Density, Diet and Productivity of Long-Eared Owls Asio otus in the Italian Alps: The Importance of Microtus Voles. Bird Study 2008, 55, 321–328. [Google Scholar] [CrossRef]
  91. Wijnandts, H. Ecological Energetics of the Long-Eared Owl (Asio otus). Ardea 1984, 72, 1–92. [Google Scholar] [CrossRef]
  92. Glue, D.E.; Hammond, G.J. Feeding Ecology of the Long-Eared Owl in Britain and Ireland. Br. Birds 1974, 67, 361–369. [Google Scholar]
  93. Romanowski, J.; Zmihorski, M. Effect of Season, Weather and Habitat on Diet Variation of a Feeding-Specialist: A Case Study of the Long-Eared Owl, Asio otus in Central Poland. Folia Zool. 2008, 57, 411–419. [Google Scholar]
  94. Heisler, L.M.; Somers, C.M.; Poulin, R.G. Owl Pellets: A More Effective Alternative to Conventional Trapping for Broad-Scale Studies of Small Mammal Communities. Methods Ecol. Evol. 2016, 7, 96–103. [Google Scholar] [CrossRef]
  95. Glue, D. Breeding Biology of Long-Eared Owls. Br. Birds 1977, 70, 318–331. [Google Scholar]
  96. Glue, D.; Nilsson, I.N. Long-Eared Owl Asio Otus. In The EBCC Atlas of European Breeding Birds: Their Distribution and Abundance; Hagemeijer, E.J.M., Blair, M.J., Eds.; T & A.D. Poyser: London, UK, 1997. [Google Scholar]
  97. Meller, K.; Björklund, H.; Saurola, P.; Valkama, J. Petolintuvuosi 2016, Pesimistulokset Ja Kannankehitykset. Linnut-vuosikirja 2017, 2016, 16–31. [Google Scholar]
  98. BirdLife International. Asio Otus. Available online: https://dx.doi.org/10.2305/IUCN.UK.2021-3.RLTS.T22689507A201150685.en (accessed on 9 December 2022).
  99. Korpimäki, E. Diet Composition, Prey Choice, and Breeding Success of Long-Eared Owls—Effects of Multiannual Fluctuations in Food Abundance. Can. J. Zool. Rev. Can. Zool. 1992, 70, 2373–2381. [Google Scholar] [CrossRef]
  100. Korpimäki, E. Population Dynamics of Fennoscandian Owls in Relation to Wintering Conditions and Between-Year Fluctuations of Food. In The Ecology and Conservation of European Owls; Joint Nature Conservation Committee: Peterborough, UK, 1992; pp. 1–10. [Google Scholar]
  101. Avotiņš Jun, A. Apodziņa Glaucidium Passerinum, Bikšainā Apoga Aegolius Funereus, Meža Pūces Strix Aluco, Urālpūces Strix Uralensis, Ausainās Pūces Asio Otus Un Ūpja Bubo Bubo Aizsardzības Plāns Plāns; LOB: Riga, Latvia, 2019. [Google Scholar]
  102. Martínez, J.A.; Zuberogoitia, I. Habitat Preferences for Long-Eared Owls Asio otus and Little Owls Athene noctua in Semi-Arid Environments at Three Spatial Scales. Bird Study 2004, 51, 163–169. [Google Scholar] [CrossRef]
  103. Aschwanden, J.; Birrer, S.; Jenni, L. Are Ecological Compensation Areas Attractive Hunting Sites for Common Kestrels (Falco tinnunculus) and Long-Eared Owls (Asio otus)? J. Ornithol. 2005, 146, 279–286. [Google Scholar] [CrossRef]
  104. Henrioux, F. Home Range and Habitat Use by the Long-Eared Owl in Northwestern Switzerland. J. Raptor Res. 2000, 34, 93–101. [Google Scholar]
  105. Galeotti, P.; Tavecchia, G.; Bonetti, A. Home-Range and Habitat Use of Long-Eared Owls in Open Farmalend (Po Plain, Northern Italy), in Relation to Prey Availability. J. Wildl. Res. 1997, 2, 137–145. [Google Scholar]
  106. Korpimäki, E. On the Ecology and Biology of Tengmalm’s Owl (Aegolius funereus) in Southern Ostrobothnia and Suomenselkä, Western Finland. Acta Univ. Oul. Ser. A Sci. Rer. Nat. 1981, 118, 1–84. [Google Scholar]
  107. Solheim, R. Breeding Frequency of Tengmalm’s Owl Aegolius funereus in Three Localities in 1974–1978. In Proceedings of the Third NOK, Ribe, Denmark, 3–8 August 1981; pp. 79–84. [Google Scholar]
  108. Hakkarainen, H.; Koivunen, V.; Korpimäki, E. Reproductive Success and Parental Effort of Tengmalm’s Owls: Effects of Spatial and Temporal Variation in Habitat Quality. Écoscience 1997, 4, 35–42. [Google Scholar] [CrossRef]
  109. Korpimäki, E. Clutch Size, Breeding Success and Brood Size Experiments in Tengmalm’ s Owl Aegolius funereus: A Test of Hypotheses. Ornis Scand. 1987, 18, 277–284. [Google Scholar] [CrossRef]
  110. Vrezec, A. Breeding Density and Altitudinal Distribution of the Ural, Tawny, and Boreal Owls in North Dinaric Alps (Central Slovenia). J. Raptor Res. 2003, 37, 55–62. [Google Scholar]
  111. Vrezec, A.; Tome, D. Habitat Selection and Patterns of Distribution in a Hierarchic Forest Owl Guild. Ornis Fenn. 2004, 81, 109–118. [Google Scholar]
  112. Kajtoch, Ł.; Matysek, M.; Figarski, T. Spatio-Temporal Patterns of Owl Territories in Fragmented Forests Are Affected by a Top Predator (Ural Owl). Ann. Zool. Fenn. 2016, 53, 165–174. [Google Scholar] [CrossRef]
  113. Eionet. Bird Population Status and Trends at the EU and Member State Levels 2013–2018. Available online: https://nature-art12.eionet.europa.eu/article12/ (accessed on 5 December 2022).
  114. Hakkarainen, H.; Korpimäki, E.; Laaksonen, T.; Nikula, A.; Suorsa, P. Survival of Male Tengmalm’s Owls Increases with Cover of Old Forest in Their Territory. Oecologia 2008, 155, 479–486. [Google Scholar] [CrossRef]
  115. Laaksonen, T.; Hakkarainen, H.; Korpimäki, E. Lifetime Reproduction of a Forest-Dwelling Owl Increases with Age and Area of Forests. Proc. R. Soc. B Biol. Sci. 2004, 271, 461–464. [Google Scholar] [CrossRef]
  116. Ravussin, P.-A.; Trolliet, D.; Willenegger, L.; Béguin, D.; Matalon, G. Choix Du Site de Nidification Chez La Chouette de Tengmalm Aegolius Funereus: Influence Des Nichoirs. Nos Oiseaux 2001, 5, 41–51. [Google Scholar]
  117. Sorbi, S. La Chouette de Tengmalm (Aegolius funereus) En Belgique. Synthèse et Mise à Jour Du Statut. Aves 1995, 32, 101–132. [Google Scholar]
  118. Shurulinkov, P.; Stoyanov, G. Some New Findings of Pigmy Owls Glaucidium passerinum and Tengmalm’s Owl Aegolius funerus in Western and Southern Bulgaria. Acrocephalus 2006, 27, 65–68. [Google Scholar]
  119. Brambilla, M.; Bassi, E.; Bergero, V.; Casale, F.; Chemollo, M.; Falco, R.; Longoni, V.; Saporetti, F.; Viganò, E. Modelling Distribution and Potential Overlap between Boreal Owl Aegolius funereus and Black Woodpecker Dryocopus martius Implications for Management and Monitoring Plans. Bird Conserv. Int. 2013, 23, 502–511. [Google Scholar] [CrossRef]
  120. Avotins, A. Modelling Owl and Woodpecker Habitat Suitability to Evaluate Forest Conservation in Latvia. In Proceedings of the 11th International Conference on Biodiversity Research, Daugavpils, Latvia, 20–22 October 2022; p. 17. [Google Scholar]
  121. Mikusek, R.; Kloubec, B.; Obuch, J. Diet of the Pygmy Owl (Glaucidium passerinum) in Eastern Central Europe. Buteo 2001, 12, 47–60. [Google Scholar]
  122. Morosinotto, C.; Villers, A.; Thomson, R.L.; Varjonen, R.; Korpimäki, E. Competitors and Predators Alter Settlement Patterns and Reproductive Success of an Intraguild Prey. Ecol. Monogr. 2017, 87, 4–20. [Google Scholar] [CrossRef]
  123. Lehikoinen, A.; Hokkanen, T.; Lokki, H. Young and Female-Biased Irruptions in Pygmy Owls Glaucidium passerinum in Southern Finland. J. Avian Biol. 2011, 42, 564–569. [Google Scholar] [CrossRef]
  124. Barbaro, L.; Blache, S.; Trochard, G.; Arlaud, C.; de Lacoste, N.; Kayser, Y. Hierarchical Habitat Selection by Eurasian Pygmy Owls Glaucidium passerinum in Old-Growth Forests of the Southern French Prealps. J. Ornithol. 2016, 157, 333–342. [Google Scholar] [CrossRef]
  125. Baroni, D.; Korpimäki, E.; Selonen, V.; Laaksonen, T. Tree Cavity Abundance and beyond: Nesting and Food Storing Sites of the Pygmy Owl in Managed Boreal Forests. For. Ecol. Manag. 2020, 460, 117818. [Google Scholar] [CrossRef]
  126. Baroni, D.; Masoero, G.; Korpimäki, E.; Morosinotto, C.; Laaksonen, T. Habitat Choice of a Secondary Cavity User Indicates Higher Avoidance of Disturbed Habitat during Breeding than during Food-Hoarding. For. Ecol. Manag. 2021, 483, 118925. [Google Scholar] [CrossRef]
  127. Pakkala, T.; Lindén, A.; Tiainen, J.; Tomppo, E.; Kouki, J. Indicators of Forest Biodiversity: Which Bird Species Predict High Breeding Bird Assemblage Diversity in Boreal Forests at Multiple Spatial Scales? Ann. Zool. Fenn. 2014, 51, 457–476. [Google Scholar] [CrossRef]
  128. Shurulinkov, P.; Ralev, A.; Daskalova, G.; Chakarov, N. Distribution, Numbers and Habitat of Pigmy Owl Glaucidium passerinum in Rhodopes Mts (S Bulgaria). Acrocephalus 2007, 28, 161–165. [Google Scholar]
  129. Sonerud, G.A. Risk of Nest Predation in Three Species of Hole Nesting Owls: Influence on Choice of Nesting Habitat and Incubation Behaviour. Ornis Scand. 1985, 16, 261. [Google Scholar] [CrossRef]
  130. Strom, H.; Sonerud, G.A. Home Range and Habitat Selection in the Pygmy Owl Glaucidium passerinum. Ornis Fenn. 2001, 78, 145–158. [Google Scholar]
  131. Avotins, A.; Kerus, V.; Aunins, A. National Scale Habitat Suitability Analysis to Evaluate and Improve Conservation Areas for a Mature Forest Specialist Species. Glob. Ecol. Conserv. 2022, 38, e02218. [Google Scholar] [CrossRef]
  132. Terraube, J.; Villers, A.; Poudré, L.; Varjonen, R.; Korpimäki, E. Increased Autumn Rainfall Disrupts Predator–Prey Interactions in Fragmented Boreal Forests. Glob. Change Biol. 2017, 23, 1361–1373. [Google Scholar] [CrossRef]
  133. Carey, A.B.; Johnson, M.L. Small Mammals in Managed, Naturally Young, and Old-Growth Forests. Ecol. Appl. 1995, 5, 336–352. [Google Scholar] [CrossRef]
  134. Tsiafouli, M.A.; Apostolopoulou, E.; Mazaris, A.D.; Kallimanis, A.S.; Drakou, E.G.; Pantis, J.D. Human Activities in Natura 2000 Sites: A Highly Diversified Conservation Network. Environ. Manag. 2013, 51, 1025–1033. [Google Scholar] [CrossRef]
  135. Jäderholm, K. Diets of the Tengalm’s Owl Aegolius funereus and the Ural Owl Strix uralensis in Central Finland. Ornis Fenn. 1987, 64, 149–153. [Google Scholar]
  136. Korpimäki, E.; Sulkava, S. Diet and Breeding Performance of Ural Owl Strix uralensis under Fluctuating Food Conditions. Ornis Fenn. 1987, 64, 57–66. [Google Scholar]
  137. Shamovich, D.I.; Shamovich, I.Y. Spatial and Temporal Differences in the Association of Forest Owls in Conditions of Different Landscape Types of Northern Belarus. In Owls of the Northern Eurasia; Volkov, S.V., Morozov, V.V., Sharikov, A.V., Eds.; MPGU: Moscow, Russia, 2005; pp. 121–136. [Google Scholar]
  138. Sidorovich, V.E.; Shamovich, D.I.; Solovey, I.A.; Lauzhel, G.O. Dietary Variations of the Ural Owl Strix uralensis in the Transitional Mixed Forest of Northern Belarus with Implications for the Distribution Differences. Ornis Fenn. 2003, 80, 145–158. [Google Scholar]
  139. Tishechkin, A.K. Comparative Food Niche Analysis of Strix Owls in Belarus. In Biology and Conservation of Owls of the Northern Hemisphere; USDA Forest Service: Washington, DC, USA, 1997; pp. 456–460. [Google Scholar]
  140. Pietiäinen, H. Breeding Season Quality, Age, and the Effect of Experience on the Reproductive Success of the Ural Owl (Strix uralensis). Auk 1988, 105, 316–324. [Google Scholar]
  141. Brommer, J.E.; Pietiäinen, H.; Kolunen, H. The Effect of Age at First Breeding on Ural Owl Lifetime Reproductive Success and Fitness under Cyclic Food Conditions. J. Anim. Ecol. 1998, 67, 248–258. [Google Scholar] [CrossRef]
  142. Karell, P.; Pietiäinen, H.; Siitari, H.; Pihlaja, T.; Kontiainen, P.; Brommer, J.E. Parental Allocation of Additional Food to Own Health and Offspring Growth in a Variable Environment. Can. J. Zool. 2008, 87, 8–19. [Google Scholar] [CrossRef]
  143. Pietiäinen, H.; Saurola, P.; Väisänen, R.A.; Ornis, S.; Scandinavian, S.; Dec, N. Parental Investment in Clutch Size and Egg Size in the Ural Owl Strix uralensis. Ornis Scand. 1986, 17, 309–325. [Google Scholar] [CrossRef]
  144. Brommer, J.E.; Kontiainen, P.; Pietiäinen, H. Selection on Plasticity of Seasonal Life-History Traits Using: Random Regression Mixed Model Analysis. Ecol. Evol. 2012, 2, 695–704. [Google Scholar] [CrossRef] [PubMed]
  145. Bylicka, M.; Kajtoch, Ł.; Figarski, T. Habitat and Landscape Characteristics Affecting the Occurrence of Ural Owls Strix uralensis in an Agroforestry Mosaic. Acta Ornithol. 2010, 45, 33–42. [Google Scholar] [CrossRef]
  146. Tutiš, V.; Radović, D.; Ćiković, D.; Barišić, S.; Kralj, J. Distribution, Density and Habitat Relationships of the Ural Owl Strix uralensis Macroura in Croatia. Ardea 2009, 97, 563–570. [Google Scholar] [CrossRef]
  147. Vrezec, A.; Mihelič, T. The Ural Owl, Strix Uralensis Macroura, in Slovenia: An Overview of Current Knowledge on Species Ecology. Riv. Ital. Ornitol. 2012, 82, 30–37. [Google Scholar] [CrossRef]
  148. Bashta, A.-T. Ural Owl Strix Uralensis Population Dynamics and Range Expansion in Western Ukraine. Ardea 2009, 97, 483–487. [Google Scholar] [CrossRef]
  149. Sonerud, G.A. Effect of Snow Cover on Seasonal Changes in Diet, Habitat and Regional Distribution of Raptors That Prey on Small Mammals in Boreal Zones of Fennoscandia. Holarct. Ecol. 1986, 9, 33–47. [Google Scholar] [CrossRef]
  150. Lõhmus, A. Do Ural Owls (Strix uralensis) Suffer from the Lack of Nest Sites in Managed Forests? Biol. Conserv. 2003, 110, 1–9. [Google Scholar] [CrossRef]
  151. Ķerus, V.; Dekants, A.; Auniņš, A.; Mārdega, I. Latvijas Ligzdojošo Putnu Atlanti 1980–2017; Latvijas Ornitoloģijas Biedrība: Riga, Latvia, 2021. [Google Scholar]
  152. Solonen, T.; Karhunen, J.; Kekkonen, J.A.; Kolunen, H.; Pietiäinen, H. Tawny Owl Prey Remains Indicate Differences in the Dynamics of Coastal and Inland Vole Populations in Southern Finland. Popul. Ecol. 2016, 58, 557–565. [Google Scholar] [CrossRef]
  153. Jędrzejewski, W.; Jędrzejewska, B.; Zub, K.; Andrzej, L.; Bystrowski, C. Resource Use by Tawny Owls Strix aluco in Relation to Rodent Fluctuations in Białowieża National Park, Poland. J. Avian Biol. 1994, 25, 308–318. [Google Scholar] [CrossRef]
  154. Goszczynski, J.; Jablonski, P.; Lesiński, G.; Romanowski, J. Variation in Diet of Tawny Owl Strix aluco L. along an Urbanization Gradient. Acta Ornithol. 1993, 27, 113–123. [Google Scholar]
  155. Lesiński, G.; Gryz, J.; Kowalski, M. Bat Predation by Tawny Owls Strix aluco in Differently Human Transformed Habitats. Ital. J. Zool. 2009, 76, 415–421. [Google Scholar] [CrossRef]
  156. Petty, S.J. Diet of Tawny Owls (Strix aluco) in Relation to Field Vole (Microtus agrestis) Abundance in a Conifer Forest in Northern England. J. Zool. 1999, 248, 451–465. [Google Scholar] [CrossRef]
  157. Romanowski, J.; Zmihorski, M. Seasonal and Habitat Variation in the Diet of the Tawny Owl (Strix aluco) in Central Poland during Unusually Warm Years. Biologia 2009, 64, 365–369. [Google Scholar] [CrossRef]
  158. Solonen, T. Are Vole-Eating Owls Affected by Mild Winters in Southern Finland? Ornis Fenn. 2004, 81, 65–74. [Google Scholar]
  159. Solonen, T.; Karhunen, J.; Kekkonen, J.; Kolunen, H.; Pietiäinen, H. Diet and Reproduction in Coastal and Inland Populations of the Tawny Owl Strix aluco in Southern Finland. J. Ornithol. 2017, 158, 541–548. [Google Scholar] [CrossRef]
  160. Solonen, T. Timing of Breeding in Rural and Urban Tawny Owls Strix aluco in Southern Finland: Effects of Vole Abundance and Winter Weather. J. Ornithol. 2014, 155, 27–36. [Google Scholar] [CrossRef]
  161. Solonen, T.; Ahola, K.; Karstinen, T. Clutch Size of a Vole-Eating Bird of Prey as an Indicator of Vole Abundance. Environ. Monit. Assess. 2015, 187, 588. [Google Scholar] [CrossRef]
  162. Karell, P.; Ahola, K.; Karstinen, T.; Valkama, J.; Brommer, J.E. Climate Change Drives Microevolution in a Wild Bird. Nat. Commun. 2011, 2, 208. [Google Scholar] [CrossRef]
  163. Solonen, T. Significance of Plumage Colour for Winter Survival in the Tawny Owl (Strix aluco): Revisiting the Camouflage Hypothesis. Ibis 2021, 163, 1437–1442. [Google Scholar] [CrossRef]
  164. Rumbutis, S.; Vaitkuvienė, D.; Grašytė, G.; Dagys, M.; Dementavičius, D.; Treinys, R. Adaptive Habitat Preferences in the Tawny Owl Strix aluco. Bird Study 2017, 64, 421–430. [Google Scholar] [CrossRef]
  165. Mitropolskiy, O.V.; Rustamov, A.K. Eagle Owl—Bubo Bubo Linnaeus, 1758. In The Birds of Middle Asia; Union of Protection of Birds of Kazakhstan: Almaty, Russia, 2007; pp. 423–431. [Google Scholar]
  166. Korpimäki, E.; Huhtala, K.; Sulkava, S. Does the Year-To-Year Variation in the Diet of Eagle and Ural Owls Support the Alternative Prey Hypothesis? Oikos 1990, 58, 47–54. [Google Scholar] [CrossRef]
  167. Randla, T. Eesti Röövlinnud [Estonian Birds of Prey]; Valgus: Tallinn, Estinia, 1976. [Google Scholar]
  168. Schweiger, A. Bubo Bubo Während Der Jungenaufzucht in Bayern. Ornitol. Anz. 2011, 50, 1–25. [Google Scholar]
  169. Frafjord, K. Population Dynamics of an Island Population of Water Voles Arvicola amphibius (Linnaeus, 1758) with One Major Predator, the Eagle Owl Bubo bubo (Linnaeus, 1758), in Northern Norway. Polar Biol. 2022, 45, 1–12. [Google Scholar] [CrossRef]
  170. BirdLife International. Species Factsheet: Bubo bubo. Available online: http://datazone.birdlife.org/species/factsheet/eurasian-eagle-owl-bubo-bubo (accessed on 8 December 2022).
  171. Penteriani, V.; Gallardo, M.; Roche, P. Landscape Structure and Food Supply Affect Eagle Owl (Bubo bubo) Density and Breeding Performance: A Case of Intra-Population Heterogeneity. J. Zool. 2002, 257, S0952836902000961. [Google Scholar] [CrossRef]
  172. Campioni, L.; Delgado, M.D.M.; Lourenço, R.; Bastianelli, G.; Fernández, N.; Penteriani, V. Individual and Spatio-Temporal Variations in the Home Range Behaviour of a Long-Lived, Territorial Species. Oecologia 2013, 172, 371–385. [Google Scholar] [CrossRef]
  173. Marchesi, L.; Sergio, F.; Pedrini, P. Costs and Benefits of Breeding in Human-Altered Landscapes for the Eagle Owl Bubo bubo. Ibis 2002, 144, E144–E177. [Google Scholar] [CrossRef]
  174. Karell, P.; Kontiainen, P.; Pietiäinen, H.; Siitari, H.; Brommer, J.E. Maternal Effects on Offspring Igs and Egg Size in Relation to Natural and Experimentally Improved Food Supply. Funct. Ecol. 2008, 22, 682–690. [Google Scholar] [CrossRef]
  175. Pietiäinen, H.; Kolunen, H. Female Body Condition and Breeding of the Ural Owl Strix uralensis. Funct. Ecol. 1993, 7, 726–735. [Google Scholar] [CrossRef]
Figure 1. Study site locations: (a) Latvia in Europe; (b) small mammal monitoring areas; (c) owl diet sampling sites; (d) owl population change monitoring areas.
Figure 1. Study site locations: (a) Latvia in Europe; (b) small mammal monitoring areas; (c) owl diet sampling sites; (d) owl population change monitoring areas.
Life 13 00572 g001
Figure 2. Observed small mammal (pooled across species) population densities per 100 trap days: (a) variation over time (1991–2016) in different sample areas and habitats; (b) differences in the observed number of small mammals (per 100 trap days) between habitats over time (2012–2016).
Figure 2. Observed small mammal (pooled across species) population densities per 100 trap days: (a) variation over time (1991–2016) in different sample areas and habitats; (b) differences in the observed number of small mammals (per 100 trap days) between habitats over time (2012–2016).
Life 13 00572 g002
Figure 3. Population indices with standard errors: (a) pooled small mammals; (b) voles of genus Microtus; (c) bank voles; and (d) autocorrelation function analysis of the yearly small mammal indices; dotted lines indicate the threshold of significance.
Figure 3. Population indices with standard errors: (a) pooled small mammals; (b) voles of genus Microtus; (c) bank voles; and (d) autocorrelation function analysis of the yearly small mammal indices; dotted lines indicate the threshold of significance.
Life 13 00572 g003
Figure 4. Number of owl diet samples (with at least five mammal individuals) over time.
Figure 4. Number of owl diet samples (with at least five mammal individuals) over time.
Life 13 00572 g004
Figure 5. Temporal change in Levin’s food niche breadth in four owl species: (a) STRALU; (b) STRURA; (c) ASIOTU; (d) GLAPAS. Grey points represent individual prey samples. The black line is the linear regression trend, and the grey ribbon is the 95% confidence interval. Numbers above the X-axis represent the number of samples. The Y-axis range differs depending on the facet.
Figure 5. Temporal change in Levin’s food niche breadth in four owl species: (a) STRALU; (b) STRURA; (c) ASIOTU; (d) GLAPAS. Grey points represent individual prey samples. The black line is the linear regression trend, and the grey ribbon is the 95% confidence interval. Numbers above the X-axis represent the number of samples. The Y-axis range differs depending on the facet.
Life 13 00572 g005
Figure 6. Population changes in owl species: (a) STRALU; (b) STRURA; (c) ASIOTU; (d) AEGFUN; and (e) GLAPAS. Black dots with error bars are yearly indices with standard errors, and trend lines and ribbons with different colors represent different population trends: black—overall trend, blue—1991–2004, and red—2004–2016. The Y-axis range differs depending on the facet.
Figure 6. Population changes in owl species: (a) STRALU; (b) STRURA; (c) ASIOTU; (d) AEGFUN; and (e) GLAPAS. Black dots with error bars are yearly indices with standard errors, and trend lines and ribbons with different colors represent different population trends: black—overall trend, blue—1991–2004, and red—2004–2016. The Y-axis range differs depending on the facet.
Life 13 00572 g006
Figure 7. Annual number of fledglings of four owl species: (a) STRALU; (b) STRURA; (c) ASIOTU; (d) BUBBUB. Grey points represent individual nest performance. Red crosses are bootstrapped 95% confidence intervals around the annual mean value. The black line is the linear regression through mean values, and the grey ribbon is the 95% confidence interval. Numbers above the X-axis represent the annual number of samples. The Y-axis range differs depending on the facet.
Figure 7. Annual number of fledglings of four owl species: (a) STRALU; (b) STRURA; (c) ASIOTU; (d) BUBBUB. Grey points represent individual nest performance. Red crosses are bootstrapped 95% confidence intervals around the annual mean value. The black line is the linear regression through mean values, and the grey ribbon is the 95% confidence interval. Numbers above the X-axis represent the annual number of samples. The Y-axis range differs depending on the facet.
Life 13 00572 g007
Table 1. Descriptions of small mammal monitoring efforts.
Table 1. Descriptions of small mammal monitoring efforts.
Sample AreaPeriodDescription
Apsalas1991–2011; 2015–20162 habitats: meadow and forest (OFD); 100 traps per transect
Lisiņa2001–20052 habitats: meadow and forest (OFU); 100 traps per transect
Žūklis2003–2011; 2015–20162 habitats: meadow and forest (OFD); 100 traps per transect
Ķemeri2006–2010; 2015–20162 habitats: meadow and forest (OFU); 100 traps per transect
Kaunata2012–201611 habitats: 1 meadow and 10 forest classes; 20–25 traps per transect
Lubāna2012–201611 habitats: 1 meadow and 10 forest classes; 20–25 traps per transect
Pļaviņas2012; 201611 habitats: 1 meadow and 10 forest classes; 20–25 traps per transect
Augstroze2012–201611 habitats: 1 meadow and 10 forest classes; 20–25 traps per transect
Table 2. Description of the influence of the prey abundance index on the FNB of owl species.
Table 2. Description of the influence of the prey abundance index on the FNB of owl species.
Owl
Species
Prey (Index)β ± SETest
Statistic
df *p-ValueAICcR2adj./R2marg. **R2cond.ICC
STRALUSmall mammals−0.6127 ± 0.1850−3.312144.7780.0012662.3080.0630.0680.005
Microtus voles−0.3886 ± 0.1551−2.506156.3030.0132667.0010.0370.0490.013
Bank voles−0.2268 ± 0.0812−2.795161.4490.0058666.9130.0460.0520.007
STRURASmall mammals−0.0039 ± 0.25140.015540.9880177.536−0.019
Microtus voles−0.1304 ± 0.1820−0.716540.4770177.007−0.009
Bank voles−0.0050 ± 0.0837−0. 06540.9520177.532−0.019
GLAPASSmall mammals−0.4948 ± 0.7738−0.63950.550739.803−0.109
Microtus voles−0.5290 ± 0.5286−1.00150.362839.0750.0003
Bank voles−0.4110 ± 0.2885−1.42550.213637.9670.147
ASIOTUSmall mammals0.0171 ± 0.15310.112220.912056.975−0.045
Microtus voles−0.0672 ± 0.1497−0.449220.658056.769−0.036
Bank voles0.0236 ± 0.07520.314220.757056.881−0.041
* Satterthwaite’s degrees of freedom in LMM. ** R2adjusted reported in the case of LM; R2marginal in the case of LMM.
Table 3. Description of prey weight proportion in owl diet relative to the abundance index in nature (first two rows per owl species) or their cross correlation.
Table 3. Description of prey weight proportion in owl diet relative to the abundance index in nature (first two rows per owl species) or their cross correlation.
Owl SpeciesPrey (Index) *β ± SETest Statisticp-ValueAICcR2MF/R2marg. **R2cond.ICC
STRALUBank voles0.0248 ± 0.00574.359<0.00015693.8220.00050.1340.133
Microtus voles0.1302 ± 0.006121.520<0.000114934.8230.0040.1910.189
Bank~Microtus−0.0138 ± 0.0116−1.1890.23405711.112<0.00010.1340.134
STRURABank voles−0.1116 ± 0.0072−15.570<0.00014693.1870.0581
Microtus voles0.1275 ± 0.010012.720<0.00016132.0770.0268
Bank~Microtus−0.1658 ± 0.0175−9.497<0.00014865.0010.0208
GLAPASBank voles−0.4168 ± 0.0717−5.812<0.0001261.0480.178
Microtus voles0.2136 ± 0.05563.8390.0001220.5910.073
Bank~Microtus−0.9097 ± 0.0811−11.220<0.0001145.4460.613
ASIOTUBank voles−0.5294 ± 0.0497−10.640<0.0001761.3180.215
Microtus voles0.0611 ± 0.01384.419<0.0001599.7000.044
Bank~Microtus−1.2718 ± 0.1208−10.530<0.0001705.7940.276
* Prey type Bank~Microtus represents the proportion of bank voles in the diet depending on Microtus sp. vole abundance in nature. ** R2McFadden reported in the case of GLM; R2marginal reported in the case of GLMM.
Table 4. Description of owl population change trends with small mammal cycles (“before”) and since they vanished (“after”). Model coefficients are in log-odds scale.
Table 4. Description of owl population change trends with small mammal cycles (“before”) and since they vanished (“after”). Model coefficients are in log-odds scale.
Owl SpeciesParameterβ ± SETest Statisticp-Valuedf *R2adj. **
STRALUIntercept0.1178 ± 0.03962.9780.006226−0.1096
Time0.0016 ± 0.00430.3680.7159
Before reference
After<0.0001 ± 0.0059<0.00011
Time:After<0.0001 ± 0.0006<0.00011
STRURAIntercept0.1389 ± 0.13611.0200.319210.5837
Time0.0319 ± 0.01791.7830.089
Before reference
After1.1350 ± 0.19255.898<0.0001
Time:After−0.1354 ± 0.0253−5.361<0.0001
GLAPAS ***Intercept0.0069 ± 0.074700.0930.9279110.4368
Time−0.0339 ± 0.0106−3.2100.0083
ASIOTUIntercept−0.8450 ± 0.2688−3.1440.004623−0.0302
Time−0.0143 ± 0.0328−0.4360.6666
Before reference
After0.3087 ± 0.38010.8120.4250
Time:After−0.0080 ± 0.0464−0.1720.8651
AEGFUNIntercept−0.5246 ± 0.0331−15.855<0.0001230.9740
Time−0.0607 ± 0.0040−15.033<0.0001
Before reference
After−0.7508 ± 0.0468−16.043<0.0001
Time:After0.0183 ± 0.00573.2110.0039
* df is the same in each variable. ** Values represent the whole model. *** GLAPAS had data only since 2004.
Table 5. Description of owl breeding performance trends with small mammal cycles (“before”) and since their depletion (“after”).
Table 5. Description of owl breeding performance trends with small mammal cycles (“before”) and since their depletion (“after”).
Owl SpeciesParameterβ ± SETest Statisticp-Valuedf *R2adj. **
STRALUIntercept2.6221 ± 0.159416.452<0.0001230.0765
Time0.0174 ± 0.01940.8960.380
Before reference
After−0.3597 ± 0.2254−1.5960.124
Time:After0.0216 ± 0.02750.7840.441
STRURAIntercept1.8291 ± 0.170910.705<0.000121−0.0812
Time0.0106 ± 0.02240.4730.641
Before reference
After−0.2022 ± 0.2416−0.8370.412
Time:After0.0293 ± 0.03170.9230.366
* df is the same in each variable. ** Values represent the whole model.
Table 6. Spearman’s correlation analysis results of the annual mean number of fledglings and small mammal abundance indices in the year of breeding and the previous autumn (annotated as −1).
Table 6. Spearman’s correlation analysis results of the annual mean number of fledglings and small mammal abundance indices in the year of breeding and the previous autumn (annotated as −1).
Owl SpeciesPrey (Index)rsp-ValueNumber of YearsS
STRALUSmall mammals−0.11520.5737263262
Microtus voles−0.13910.4962263332
Bank voles0.03380.8700262826
Small mammals−10.13000.5341252262
Microtus voles−10.13770.5100252242
Bank voles−10.19460.3496252094
STRURASmall mammals−0.01790.9340242341.1
Microtus voles0.05450.8005242174.8
Bank voles−0.08630.6886242498.4
Small mammals−1−0.17880.4145232385.8
Microtus voles−1−0.21040.3351232449.9
Bank voles−10.02180.9214231979.9
ASIOTUSmall mammals0.08050.775515514.92
Microtus voles0.12520.656615489.87
Bank voles0.07870.780415515.92
Small mammals−10.30620.287114315.69
Microtus voles−10.05290.857614430.95
Bank voles−1−0.26430.361214575.26
BUBBUBSmall mammals0.53290.040815261.59
Microtus voles0.38170.160415346.28
Bank voles0.34020.214615369.49
Small mammals−10.64380.013014162.09
Microtus voles−10.55270.040414203.50
Bank voles−10.20200.488614363.09
−1 Abundance index in the previous autumn.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Avotins, A.; Avotins, A., Sr.; Ķerus, V.; Aunins, A. Numerical Response of Owls to the Dampening of Small Mammal Population Cycles in Latvia. Life 2023, 13, 572. https://doi.org/10.3390/life13020572

AMA Style

Avotins A, Avotins A Sr., Ķerus V, Aunins A. Numerical Response of Owls to the Dampening of Small Mammal Population Cycles in Latvia. Life. 2023; 13(2):572. https://doi.org/10.3390/life13020572

Chicago/Turabian Style

Avotins, Andris, Andris Avotins, Sr., Viesturs Ķerus, and Ainars Aunins. 2023. "Numerical Response of Owls to the Dampening of Small Mammal Population Cycles in Latvia" Life 13, no. 2: 572. https://doi.org/10.3390/life13020572

APA Style

Avotins, A., Avotins, A., Sr., Ķerus, V., & Aunins, A. (2023). Numerical Response of Owls to the Dampening of Small Mammal Population Cycles in Latvia. Life, 13(2), 572. https://doi.org/10.3390/life13020572

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop