Next Article in Journal
Correction: Thess, A.D.; Lengsfeld, P. Side Effects of Wind Energy: Review of Three Topics—Status and Open Questions. Sustainability 2022, 14, 16186
Next Article in Special Issue
Scavenger Activity and Anti-Predator Behaviour in an Apennine Wolf Area
Previous Article in Journal
A Review of the Promotion of Sustainable Mobility of Workers by Industries
Previous Article in Special Issue
Red Deer in Lithuania: History, Status and Management
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Does Wolf Management in Latvia Decrease Livestock Depredation? An Analysis of Available Data

Latvian State Forest Research Institute Silava, Rigas Street 111, LV-2169 Salaspils, Latvia
*
Author to whom correspondence should be addressed.
Sustainability 2023, 15(11), 8509; https://0-doi-org.brum.beds.ac.uk/10.3390/su15118509
Submission received: 15 April 2023 / Revised: 21 May 2023 / Accepted: 22 May 2023 / Published: 24 May 2023
(This article belongs to the Special Issue Mammal Status: Diversity, Abundance and Dynamics Volume II)

Abstract

:
In Latvia, livestock depredation by wolves has increased during the last two decades. Most of the attacks occur in summer and autumn during wolf hunting season. Use of effective preventive measures in Latvia is low, and farmers primarily rely on wolf hunting as a depredation reduction measure. The total numbers of wolf attacks and number of affected sheep per year in regional forest management units were analyzed in relation to the estimated wolf density, extent of culling, and proportion of juveniles, as well as the sheep density and estimated number of wild prey animals. The response variables (number of attacks and affected sheep per year) were modelled using a negative binomial regression, testing the effects of every covariate separately and building models from the significant covariates. The depredation level was related to sheep density and estimated wolf population size. No reducing effect was found for culling, and an even greater depredation rate was expected when the proportion of culled wolves increased. In addition, no significant effect was associated with the other covariates. However, greater numbers of affected sheep were expected at higher red deer density, suggesting increased opportunistic livestock depredation when red deer locally outcompete roe deer, the preferred wolf prey in Latvia.

1. Introduction

For centuries humans have had diverse and complicated relationships with wolves [1], and livestock depredation is of the predominant sources of human–carnivore conflict [2,3]. Wolf attacks on livestock affect human attitudes towards them [2,3,4,5], leading to persecution and even complete eradication of this predator in many countries [1], even though more livestock are lost to diseases, harsh weather conditions, and other factors [2,6]. Due to the relatively recent recovery of wolf populations in many European countries [7] and increased depredation associated with prolonged livestock herding and breeding in the absence of wolves, derogation and more extensive application of lethal control is being reconsidered [8]. Mitigation of conflicts with wolves is important to ensure their conservation as an important part of the ecosystem, to maintain the habitual lifestyle and sources of income of local people, and to improve attitudes towards these carnivores [9,10,11].
Compared to other European countries, such as Sweden, Norway, Estonia, Poland, Slovenia, France, and Italy [2,12,13,14,15,16], livestock depredation by wolves in Latvia in the 21st century is rather low [17], with an average of 25 reported cases per year (varying from 9 to 79) [18], although it can cause significant damages to individual farmers. There are currently no subsidies for acquisition of preventive measures and no compensation paid for lost animals in Latvia [17]. As the wolf population is not endangered and wolves are hunted in Latvia, culling is seen as a management measure that can decrease the amount of livestock depredation. Lethal control of carnivore populations to reduce depredation and support the livestock industry is implemented in many countries [1,6,19]; however, the effectiveness of hunting is unclear, and has been questioned in some cases [9,10,20,21,22,23], as there are many factors (for example, lack of wild prey, number of wolves in the area, existence of wolves specializing in livestock depredation, social structure of the wolf population, stability of wolf packs, livestock density in the area, use of preventive measures, landscape characteristics) influencing the occurrence of depredation and the impact of hunting [9,13,14,19,20,24,25,26].
The aim of this study was to investigate the relationship between reported livestock depredation in Latvia and available data on the estimated density of wolves and their wild prey, as well as on culling from 2004 to 2022. Specifically, we searched for evidence of a negative impact of wolf hunting on the reported number of attacks and affected livestock to determine whether livestock depredation was locally minimized by wolf culling according to the implemented management approach. As wolf hunting can disrupt pack structure, and may cause juvenile individuals to resort to livestock depredation [19], we examined the relationship between livestock depredation and juvenile proportion, which was estimated according to the observed age structure among the culled individuals.

2. Materials and Methods

2.1. Wolves in Latvia

The study area, wolf management, and routine sampling of culled individuals in Latvia is described in more detail by Šuba et al. [27]. Wolves are distributed throughout the whole country [17,18]. A hunting season from 15 July until 31 March and a hunting quota were both introduced in Latvia in 2004. Before their implementation, wolves were hunted without any restrictions. The quota is set annually for the whole country, and the amount of livestock depredation is one of the considerations when deciding the size of the quota. The quota has been increased following its introduction, and for the last decade has been set to around 270–300 wolves per hunting season. In Latvia, the hunting pressure on wolves is rather high, with human-caused mortality estimated to be around 37% in the last two decades [27]. While the occurrence of illegal wolf hunting is plausible, its prevalence is unknown. The main prey species for wolves are roe deer (Capreolus capreolus), wild boar (Sus scrofa), and to a lesser extent red deer (Cervus elaphus) and Eurasian beaver (Castor fiber) [17,28].

2.2. Data Acquisition and Preparation

The data we analyzed in this study, referring to a period from 2004 until 2022, were obtained from the Latvian State Forest Service (SFS) and the Agricultural Data Center via the Central Statistical Bureau of Latvia. The SFS conducts a game census and investigates reports on livestock depredation by wild carnivores, while the Agricultural Data Center compiles data on livestock numbers per statistical region annually. Reports to the SFS about cases of damage are voluntary, and as there is no compensation paid for losses, not every farmer reports occurrences of depredation. The amount of unreported cases is unknown. From all reported cases, we analyzed only wolf attacks. Although wolves can be determined as a culprit in attacks by stray dogs or bears, such cases should not be substantial. Bear attacks are usually rather distinct from wolf attacks, and an ongoing study on DNA analyses showed that in Latvia, stray dogs rarely cause damage to livestock. For example, in 156 DNA samples taken from livestock depredation victims between 2018 and 2022, the presence of wolf DNA was confirmed in 88.5% cases [29].
In this study, we focused on wolf attacks on sheep, which comprise 90% of all reported livestock depredation cases. The SFS data corresponding to ten regional forestry units (Figure 1) were published on the SFS website [18,30] or made available upon request. Reports on livestock depredation included the date, location, number of killed, injured and lost sheep, the circumstances of the attack, and the preventive measures applied. In this study, victims were pooled into a single category (i.e., affected sheep). Data on applied livestock protection measures at sites where depredation had occurred (n = 506) were available for the years 2000–2020. The effectiveness of preventive measures was evaluated according to available recommendations [23,31,32,33].
Summary statistics on wolf attacks on sheep, wildlife abundance, and number of sheep in the SFS within the forestry units as well as the number of sheep in statistical regions, are provided in the Appendix A (Table A1). The estimated numbers of wolves, as well as red deer, roe dear, wild boar, and Eurasian beaver, provided by the SFS were used to account for wolf density and availability of wild prey. The number of culled wolves per forestry unit per year was used to account for hunting pressure. However, as this number was expected to correlate with estimated abundance, the culling intensity was expressed as the ratio between the number of culled individuals and the estimated number of wolves in the forestry unit.
Data on number of sheep in five statistical regions of Latvia (Figure 1) were used to account for regional variation in sheep density. A summary of sheep numbers from 2004 to 2022 is provided in the Appendix A (Table A2). As the borders of the statistical regions differ from the borders of the forestry units, data from individual or adjacent statistical regions corresponding to two to six neighboring forestry units within common or negligibly differing borders were pooled. The calculated number of sheep per km2 refers to all the forestry units within the array. Data on wolf age structure from 2004 until 2021 were obtained by the authors following examination of legally culled individuals and otherwise found carcasses (see Šuba et al. [27] for more details). Tooth samples for age assessment were prepared according to the methods described by Klevezal [34]. Age was determined via microscopic inspection and counting the increment lines in a cross-section of the extracted canine. For the purposes of this study, age was assigned to three age classes, namely, juveniles (i.e., individuals born in spring prior to opening of the hunting season), subadults (individuals aged one year), and adults (individuals aged two years and older). This allowed inclusion of adult individuals from which a tooth sample was unavailable and the precise age remained undetermined. In total, ages were recorded for 1902 individual wolves of known location. As the number of wolves with known age from individual forestry units per year was often insufficient for credible estimates, data from three to four neighboring forestry units were combined and the calculated proportion of juveniles among the sampled individuals was assigned to every forestry unit in the array.

2.3. Data Analysis

First, we obtained pooled information on reported livestock depredation cases and the number of affected sheep. In addition to the general description, the timing of attacks and the relationship with estimated wolf density within the country were analyzed. Afterwards, depredation cases within local forestry units were examined in relation to the available data.
The total number of reported attacks and affected sheep per year within a regional forestry unit were treated as response variables. Covariates are listed in Table A3 of the Appendix A and included nominal variables (i.e., taking value 0 or 1) which corresponded to forestry units and quantitative variables, namely, the mean number of sheep per 1 km2, the estimated number of wolves, the proportion of juvenile wolves, the proportion of culled individuals in the current and previous year, and the abundance of prey species in the current and previous year. Estimated prey density was expressed in thousands. Proportional variables (i.e., the proportions of culled individuals and of juveniles) were transformed using a logit function, as this provided a more feasible range. In one case, this variable assumed a value of −7 below 0.1%, and in two cases above 100% the value of 7 was assigned, as slightly more individuals per forestry unit were culled than were estimated to be present.
The relationship between the number of attacks or affected sheep and the covariates was investigated by means of a negative binomial regression. All the models were based on the same 58 records of response variables (total numbers per year) and associated values of the covariates. Initially, the effects of each covariate on the response variables were tested separately, comparing the model containing an intercept and a single covariate with a null model using the likelihood ratio test. Then, models which contained combinations of significant covariates were compared to each other. The statistical analyses were conducted using R software [35], and the MASS package was applied.

3. Results

During the study period, both the number of reported and verified wolf attacks on livestock and the number of affected sheep have fluctuated considerably, showing a slightly increasing trend (Figure 2).
The mean number of affected sheep per reported attack increased from 2.6 (in 2004–2009) to 5.5 (in 2017–2022) with a slope of 0.219 (SE = 0.077) per year (Figure 3). This increase was found to be statistically significant (linear regression analysis, F1,17 = 8.16, p = 0.011).
No livestock protection measures were used in 181 (35.8%) reported depredation cases. In 266 (52.6%) cases, the applied preventive measures were considered as inappropriate (electric fences with only one or two wire lines; electric, wood, or barbed-wire fences less than 1 m high; chained guard dogs). Only ten (2%) farms on which depredation occurred used more effective preventive measures (e.g., a shepherd or appropriate electric fencing at least 1.2 m high with five or six wire lines or mesh weave). In 49 (9.7%) cases, the reports contained no information about the use of preventive measures.
While wolf attacks were reported throughout the year, the majority occurred in summer and autumn (Figure 4). As wolf hunting season begins on 15 July, in most years it covered the period in which the majority of the attacks were reported.
Overall, the total number of wolf attacks and affected sheep increased with the estimated wolf density in the country (Figure 5), which was determined to be a significant factor according to negative binomial regression and likelihood ratio tests (for the number of the attacks, λLR = 5.911, df = 1, p = 0.015; for the number of affected sheep, λLR = 20.849, df = 1, p < 0.001). However, investigation at the level of SFS local forestry units revealed other relationships in which the estimated number of wolves no longer had such a significant effect.
By investigating the relationships between the number of wolf attacks on sheep and the available covariates via negative binomial regression, the local forestry unit (λLR = 18.17, df = 8, p = 0.02), mean number of sheep per km2LR = 7.724, df = 1, p = 0.005), and proportion of culled wolves in the current year (λLR = 6.74, df = 1, p = 0.009) had significant effects, while other covariates had no significant effect on the intercept (likelihood ratio tests, p > 0.05). The statistics of the negative binomial regression models containing combinations of these covariates and the estimated wolf density are provided in Appendix A (Table A4). However, the proportion of culled wolves had positive coefficient values, i.e., higher expected depredation rate, at a higher culling intensity. Other covariates, such as density of other wildlife species and proportion of juveniles, had an insignificant effect on the total number of depredation cases according to likelihood ratio tests (p > 0.05). The mean number of sheep per km2 had a significant effect on the cumulative number of affected sheep (λLR = 6.616, df = 1, p = 0.01). The most parsimonious models according to AIC values included the forestry unit and the estimated number of wolves and red deer among the factors (Table A5 of the Appendix A). Likelihood ratio tests revealed no significant effect of other covariates (p > 0.05).

4. Discussion

In Latvia, numbers of reported livestock depredation cases and affected sheep were correlated with estimated wolf density, as most farms (88.4%) where wolf attacks occurred applied no or insufficient preventive measures against such attacks. In a survey on public attitudes towards large carnivores [15], most livestock farmers (73.4%) claimed that they do not use any preventive measures. Wolf hunting was deemed an effective means of reducing depredation by 84.1% of surveyed farmers, and 41.1% of farmers considered hunters to be responsible for the prevention and reduction of wolf depredation. Only 29% of surveyed farmers claimed personal responsibility for the prevention of depredation cases. Generally, prevention was introduced only after loss of livestock had been suffered due to wolf attacks.
Most attacks were reported during summer and autumn. Similar timing of wolf attacks on livestock has been observed in other countries [1,4,12,13,15,24,36,37,38,39]. In Latvia, unlike in neighboring Estonia [16] and Lithuania [40], wolf hunting season begins considerably earlier, on 15 July compared to 1 November and 15 October, respectively, coinciding with the majority of the observed attacks on livestock. Nevertheless, we found no indication that wolf hunting in the current or following year decreased the reported number of attacks or the number of affected sheep in SFS local forestry units. On the contrary, significantly more attacks were expected in the current year with a higher ratio between the number of culled wolves and the estimated number of wolves, as the coefficient was positive and significantly different from zero. As seen in previous studies, lethal predator control can be less effective than other preventive measures [14,20,23], and appropriate livestock protection can be more significant than a reduction in wolf numbers in decreasing the number of depredation cases [32,41].
While hunting can have a short-term positive effect on depredation reduction, it does not prevent attacks in the long term, as harvested animals are soon replaced by dispersing individuals [24]. In fact, in certain cases wolf hunting can increase the amount of depredation [19,20,25], as hunting impacts the demographic, territorial, and social structure of wolf populations, leading to potentially higher reproduction rates [42] and possible changes in animal behavior, including hunting habits [19,43,44]. As wolf hunting in Latvia begins when pups are very young and will continue to depend on adult animals for their survival for some time [45], the loss of parents or other adult pack members may make it more difficult for the remaining adults to provide for the pups [44], and as a result they may choose more vulnerable prey, e.g., livestock.
Theoretically, an increase in livestock depredation may be associated with disrupted pack structure and accidental removal of adults due to intensive hunting [19]. However, our analysis revealed no significant relationship between the number of attacks or number affected sheep and the observed proportion of juvenile wolves. In fact, juveniles are more likely to be removed from the population due to hunting than adults [46]. Additionally, having an abundant wild prey base decreases the possibility of depredation by juvenile wolves. The impact of culling on wolf pack structure needs to be evaluated in further studies involving existing kinship data, as individual circumstances in packs, such as the juvenile’s age when losing adult pack members or early dispersal from the natal pack, might be important factors leading juveniles to depredation.
Another significant factor influencing the number of reported depredation cases and affected sheep was the location in particular SFS forestry units. In addition to regional variation of the analyzed covariates (Appendix A, Table A1 and Table A2), local differences in the operation of the SFS and the activity of farmers in reporting cases may be relevant. However, qualitative or quantitative assessment of such characteristics is problematic.
No correlation was found between sheep depredation and the estimated numbers of most prey species. Although the numbers of roe deer and wild boar in Latvia have fluctuated [17,27,47], there is no reason to assume significant shortages at any point. However, the number of affected sheep was related to estimated numbers of red deer, and according to estimated coefficients, more affected sheep were expected at higher red deer densities. This may be associated with competition between the two deer species [48,49], as roe deer are more common in the wolf diet in Latvia [28], and may be affected by higher red deer density increasing opportunistic livestock depredation. In Europe, red deer are preferred in the wolf diet [50]; however, hunting them may require advanced hunting skills or greater pack size.
In Latvian society, various opinions exist concerning wolves [51,52]. Generally, livestock farmers and herders are the most negative in their attitudes towards wolves [53,54,55,56,57,58,59], as their income and lifestyle are affected by depredation. In addition, the wolf is sometimes seen as a symbol of the domination of the urban population over the lifestyle and needs of rural inhabitants. Therefore, negative attitudes towards this carnivore may originate from the symbolic meaning of the animal and general social and economic factors rather than from negative personal experiences with wolves [60,61]. Although it is often considered that attitudes should be improved in order to improve species conservation condition, in the case of livestock depredation practical measures that ensure successful coexistence may be of greater importance. The ability to accept the presence of wolves and coexist with them may be more important than engendering a positive attitude towards these carnivores. Acknowledgement of existing conflicts, hearing out of farmers and their problems, objective evaluation of the situation, and practical solutions for conflict mitigation might be more successful than attempts to improve knowledge and attitudes towards predators [62].

5. Conclusions

In Latvia, where the wolf population is not endangered and wolf culling is permitted, their lethal control is regarded as means of reducing livestock depredation. Nevertheless, as seen from this study, current wolf hunting practices in Latvia might not have the desired positive effect on depredation reduction. To minimize livestock depredation, improvements in applied livestock depredation measures must be considered rather than increasing culling quotas. Thus, the use of effective preventive measures and subsidies for their implementation are significant approaches in the context of sustainable coexistence with these carnivores.

Author Contributions

Conceptualization, G.B., G.D., A.O., J.O., D.P., A.S., J.Š. and A.Ž.; methodology, J.Š. and A.Ž.; software J.Š.; validation, J.Š. and A.Ž.; formal analysis, J.Š. and A.Ž.; investigation, G.B., G.D., A.O., J.O., D.P., A.S., J.Š. and A.Ž.; resources, G.B., G.D., A.O., J.O., D.P., A.S. and A.Ž.; data curation, G.B., G.D., A.O., J.O., D.P., A.S. and A.Ž.; writing—original draft preparation, J.Š. and A.Ž.; writing—review and editing, J.Š., A.Ž., G.B., D.P. and J.O.; visualization, J.Š. and G.D.; supervision, J.O.; project administration, J.O., J.Š. and A.Ž.; funding acquisition, J.O. and J.Š. All authors have read and agreed to the published version of the manuscript.

Funding

This study was conducted as part of the postdoctoral research project “Sub-population dynamics of grey wolf Canis lupus and Eurasian lynx Lynx lynx in Latvia and identification of depredation risk on livestock” (No. 1.1.1.2/VIAA/3/19/511), funded by the European Regional Development Fund (agreement No. 1.1.1.2./16/I/001). The sampling of culled wolves within investigation of variations in the large carnivore population in Latvia depending on hunting pressure was supported by the Game Management Development Fund of the Ministry of Agriculture of the Republic of Latvia. Part of the data was collected and analyzed with the support of the project “Strengthening of the capacity of doctoral studies at the University of Latvia within the framework of the new doctoral model” (No. 8.2.2.0/20/I/006).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Official data on game statistics in Latvia are available from the State Forest Service website.

Acknowledgments

We are grateful to the SFS personnel who prepared the requested data and to all who provided sample material for this study and contributed to the preparation and improvement of the manuscript.

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. Summary of Analyzed Data

Table A1. Summary statistics (minimum, maximum, and median values per year) on wolf attacks on sheep and applied wildlife data of Latvian State Forest Service local forestry units (Figure 1) from 2004 to 2022.
Table A1. Summary statistics (minimum, maximum, and median values per year) on wolf attacks on sheep and applied wildlife data of Latvian State Forest Service local forestry units (Figure 1) from 2004 to 2022.
DienvidkurzemeZiemeļkurzemeZemgaleRīga Regional
Reported wolf attacks on sheep1–9 (3)1–8 (3)51–4 (1)
Total number of affected sheep1–61 (20)1–50 (17.5)261–24 (12.5)
Estimated number of wolves95–337 (169)77–260 (158)11–129 (88)5–45 (14)
Number of culled wolves20–42 (28)17–55 (39)3–39 (14)1–11 (5)
Age structure175 juv
30 subad
130 ad
163 juv
36 subad
171 ad
85 juv
11 subad
58 ad
26 juv
4 subad
21 ad
Estimated numbers of other wildlife (thousands)red deer 4.9–14.3 (10.8)
roe deer 13.4–38 (22.4)
wild boars 2.4–12.7 (7.2)
beavers 4.9–13 (7.7)
red deer 6.6–13.6 (10.6)
roe deer 7.3–22.6 (9)
wild boars 1.2–11.4 (6.1)
beavers 3.0–8.4 (4.0)
red deer 1.9–13.9 (10.2)
roe deer 13.6–27.4 (22.9)
wild boars 2.3–11.1 (4.0)
beavers 3.4–12.3 (7.0)
red deer 1.3–4.3 (2.0)
roe deer 8.9–28.9 (17.4)
wild boars 0.8–5.8 (2.3)
beavers 4.1–10.6 (5.7)
SēlijaDienvidlatgaleAustrumlatgale
Reported wolf attacks on sheep1–6 (1)1–12 (2)1–11 (3.5)
Total number of affected sheep1–31 (11)7–68 (12)2–58 (16)
Estimated number of wolves63–190 (151)39–189 (124)68–180 (99)
Number of culled wolves5–37 (22)4–36 (15)6–53 (20)
Age structure101 juv
24 subad
82 ad
50 juv
10 subad
45 ad
57 juv
6 subad
56 ad
Estimated numbers of other wildlife (thousands)red deer 2.2–7.7 (5.2)
roe deer 11.1–28.3 (18.4)
wild boars 1.6–7.5 (3.7)
beavers 2–7.7 (6.6)
red deer 0.5–3.1 (1.4)
roe deer 10.1–25.4 (16.5)
wild boars 1.9–5.9 (3.0)
beavers 7.9–13.0 (11.6)
red deer 0.2–1.8 (0.7)
roe deer 9.2–15.1 (12.2)
wild boars 1.3–4.3 (2.7)
beavers 6.2–8.6 (7.3)
CentrālvidzemeZiemeļvidzemeZiemeļaustrumi
Reported wolf attacks on sheep1–28 (2)1–9 (3.5)1–9 (1)
Total number of affected sheep1–277 (14)3–76 (27)2–62 (5)
Estimated number of wolves11–205 (61)41–143 (79)33–169 (70)
Number of culled wolves4–47 (19)2–42 (28)3–45 (19)
Age structure71 juv
15 subad
68 ad
139 juv
17 subad
73 ad
103 juv
10 subad
65 ad
Estimated numbers of other wildlife (thousands)red deer 1–7.2 (4.5)
roe deer 8.9–21.5 (13.5)
wild boars 2.4–7.3 (4.2)
beavers 5.0–8.9 (6.4)
red deer 1.5–4.7 (3.7)
roe deer 12.6–41.2 (20.1)
wild boars 1.1–10 (4.4)
beavers 5.4–11.9 (6.5)
red deer 0.7–5.1 (2.7)
roe deer 6.1–23.3 (10.0)
wild boars 1.0–6.3 (2.6)
beavers 4.1–8.4 (4.9)
Table A2. Summary statistics on minimum, maximum, and median number of sheep in five Latvian statistical regions (Figure 1) from 2004 to 2022.
Table A2. Summary statistics on minimum, maximum, and median number of sheep in five Latvian statistical regions (Figure 1) from 2004 to 2022.
KurzemeZemgalePierīgaVidzemeLatgale
Non-urban area (km2)12,99510,678856215,75014,463
Number of sheep (thousands)4.4–19.0 (12.9)2.0–14.6 (10.7)3.8–17.8 (11.8)7.0–33.5 (23.1)20.0–28.5 (26.3)
Table A3. Investigated variables and coefficients corresponding to sheep depredation by wolves and available data.
Table A3. Investigated variables and coefficients corresponding to sheep depredation by wolves and available data.
Constant or VariableCoefficient
Cumulative number of wolf attacks on sheep per year1
Cumulative number of affected sheep in wolf attacks per year1
Interceptβint
Centrālvidzeme SFS forestry unitβforestry[CV]
Dienvidkurzeme SFS forestry unitβforestry[DK]
Dienvidlatgale SFS forestry unitβforestry[DL]
Rīga Regional SFS forestry unitβforestry[RR]
Sēlija SFS forestry unitβforestry[S]
Ziemeļaustrumi SFS forestry unitβforestry[ZA]
Ziemeļkurzeme SFS forestry unitβforestry[ZK]
Ziemeļvidzeme SFS forestry unitβforestry[ZV]
Mean number of sheep per 1 km2 in respective regionβsheep
Estimated number of red deer in current year in current year (in thousands)βwild[redd]
Estimated number of roe deer in current year in current year (in thousands)βwild[roed]
Estimated number of wild boars in current year in current year (in thousands)βwild[wildb]
Estimated number of beavers in current year in current year (in thousands)βwild[beav]
Estimated number of red deer in previous year in previous year (in thousands)βprevwild[redd]
Estimated number of roe deer in previous year in previous year (in thousands)βprevwild[roed]
Estimated number of wild boars in previous year in previous year (in thousands)βprevwild[wildb]
Estimated number of beavers in previous year in previous year (in thousands)βprevwild[beav]
Estimated number of wolves in current yearβwolf
Logit-transformed proportion of culled wolves in current yearβwcull
Logit-transformed proportion of culled wolves in previous yearβprevwcull
Logit-transformed proportion of juvenile wolves in current yearβwjuv
Table A4. Coefficients and statistics of investigated negative binomial regression models (overdispersion parameter, adjusted Akaike information criterion, difference, Akaike weight, and evidence ratio) describing cumulative number of wolf attacks on sheep per year in SFS local forestry units. Significant coefficients indicated by asterisks (*—p < 0.05, **—p < 0.01, ***—p < 0.001).
Table A4. Coefficients and statistics of investigated negative binomial regression models (overdispersion parameter, adjusted Akaike information criterion, difference, Akaike weight, and evidence ratio) describing cumulative number of wolf attacks on sheep per year in SFS local forestry units. Significant coefficients indicated by asterisks (*—p < 0.05, **—p < 0.01, ***—p < 0.001).
Coefficients (±SE)θ (±SE)AICcωER
βint = 0.742 (±0.366) *5.96 (±3.07)243.6500.5621
βsheep = 0.442 (±0.195) *
βwcull = 0.263 (±0.124) *
βint = 0.293 (±0.316)5.13 (±2.44)245.391.730.2362.4
βsheep = 0.55 (±0.196) **
βint = 1.501 (±0.163) ***4.71 (±2.09)246.372.710.1453.9
βwcull = 0.345 (±0.125) **
βint = 0.999 (±0.374)**12.6 (±11.0)250.026.370.02324.1
βforestry[CV] = −0.304 (±0.376)
βforestry[DK] = −1.086 (±0.368) **
βforestry[DL] = −1.279 (±0.391) **
βforestry[RR] = −1.234 (±0.669)
βforestry[S] = −0.858 (±0.335) *
βforestry[ZA] = −0.824 (±0.351) *
βforestry[ZK] = −1.173 (±0.377) **
βforestry[ZV] = −0.099 (±0.28)
βwolf = 0.007 (±0.003) *
βint = 1.127 (±0.101) ***3.81 (±1.52)250.967.310.01538.6
βint = 0.826 (±0.496)14.5 (±14.3)252.528.870.01246.4
βforestry[CV] = −0.545 (±0.367)
βforestry[DK] = −0.128 (±0.393)
βforestry[DL] = −0.867 (±0.39) *
βforestry[RR] = −1.436 (±0.638) *
βforestry[S] = −0.489 (±0.342)
βforestry[ZA] = −0.84 (±0.34) *
βforestry[ZK] = −0.359 (±0.352)
βforestry[ZV] = −0.14 (±0.276)
βsheep = 0.539 (±0.225) *
βwcull = 0.113 (±0.135)
βint = 1.735 (±0.221) ***8.58 (±5.69)252.478.820.00782.1
βforestry[CV] = −0.482 (±0.387)
βforestry[DK] = −0.736 (±0.344) *
βforestry[DL] = −1.042 (±0.389) **
βforestry[RR] = −1.735 (±0.649) **
βforestry[S] = −0.79 (±0.348) *
βforestry[ZA] = −1.042 (±0.355) **
βforestry[ZK] = −0.736 (±0.344) *
βforestry[ZV] = −0.253 (±0.289)
Table A5. Coefficients and statistics of negative binomial regression models (overdispersion parameter, adjusted Akaike information criterion, difference, Akaike weight, and evidence ratio) describing cumulative number of affected sheep in wolf attacks per year in SFS local forestry units. Significant coefficients indicated by asterisks (*—p < 0.05, **—p < 0.01, ***—p < 0.001).
Table A5. Coefficients and statistics of negative binomial regression models (overdispersion parameter, adjusted Akaike information criterion, difference, Akaike weight, and evidence ratio) describing cumulative number of affected sheep in wolf attacks per year in SFS local forestry units. Significant coefficients indicated by asterisks (*—p < 0.05, **—p < 0.01, ***—p < 0.001).
Coefficients (±SE)θ (±SE)AICcωER
βint = 1.486 (±0.41) ***1.505 (±0.287)474.9700.4231
βsheep = 0.833 (±0.226) ***
βwild[redd] = 0.075 (±0.031) *
βint = 1.564 (±0.42) ***1.462 (±0.277)476.681.710.182.4
βsheep = 0.674 (±0.221) **
βwolf = 0.005 (±0.002) *
βint = 1.414 (±0.427) ***1.512 (±0.289)476.982.010.1552.7
βsheep = 0.798 (±0.233) ***
βwild[redd] = 0.06 (±0.041)
βwolf = 0.002 (±0.003)
βint = 2.128 (±0.349) ***1.379 (±0.259)478.073.090.094.7
βsheep = 0.65 (±0.226) **
βint = 2.982 (±0.311) ***1.922 (±0.386)478.943.970.0587.3
βforestry[CV] = −1.048 (±0.586)
βforestry[DK] = −3.085 (±0.778) ***
βforestry[DL] = −0.521 (±0.438)
βforestry[RR] = −1.419 (±0.56) *
βforestry[S] = −1.366 (±0.504) **
βforestry[ZA] = −1.585 (±0.438) ***
βforestry[ZK] = −3.389 (±0.934) ***
βforestry[ZV] = −0.547 (±0.442)
βwild[redd] = 0.311 (±0.078) ***
βint = 2.291 (±0.463) ***2.026 (±0.411)479.054.070.0557.7
βforestry[CV] = −0.449 (±0.634)
βforestry[DK] = −2.657 (±0.78) ***
βforestry[DL] = −0.632 (±0.439)
βforestry[RR] = −0.826 (±0.629)
βforestry[S] = −0.99 (±0.507)
βforestry[ZA] = −1.21 (±0.465) **
βforestry[ZK] = −2.86 (±0.938) **
βforestry[ZV] = −0.111 (±0.473)
βwild[redd] = 0.225 (±0.087) **
βwolf = 0.007 (±0.004)
βint = 1.988 (±0.465) ***1.85 (±0.37)481.366.390.01724.4
βforestry[CV] = 0.712 (±0.497)
βforestry[DK] = −0.894 (±0.461)
βforestry[DL] = −0.602 (±0.458)
βforestry[RR] = −0.264 (±0.611)
βforestry[S] = −0.246 (±0.428)
βforestry[ZA] = −0.606 (±0.435)
βforestry[ZK] = −0.697 (±0.475)
βforestry[ZV] = 0.658 (±0.4)
βwolf = 0.011 (±0.003) ***
βint = 2.404 (±0.618) ***2.025 (±0.411)482.197.220.01136.9
βforestry[CV] = −0.489 (±0.646)
βforestry[DK] = −2.851 (±1.08) **
βforestry[DL] = −0.654 (±0.455)
βforestry[RR] = −0.875 (±0.636)
βforestry[S] = −1.082 (±0.599)
βforestry[ZA] = −1.254 (±0.486) **
βforestry[ZK] = −3.046 (±1.198) *
βforestry[ZV] = −0.15 (±0.508)
βsheep = −0.08 (±0.321)
βwild[redd] = 0.237 (±0.096) *
βwolf = 0.007 (±0.004)
βint = 3.114 (±0.121) ***1.241 (±0.228)482.537.560.0143.8

References

  1. Fritts, S.H.; Stephenson, R.O.; Hayes, R.D.; Boitani, L. Wolves and Humans. In Wolves: Behavior, Ecology and Conservation; Mech, L.D., Boitani, L., Eds.; The University of Chicago Press: Chicago, IL, USA, 2003; pp. 289–316. [Google Scholar]
  2. Kaczensky, P. Large carnivore depredation on livestock in Europe. Ursus 1999, 11, 59–71. [Google Scholar]
  3. Boitani, L. Action Plan for the Conservation of Wolves in Europe (Canis lupus); Council of Europe Publishing: Strasbourg, France, 2000. [Google Scholar]
  4. Mishra, C. Livestock depredation by large carnivores in the Indian trans-Himalaya: Conflict perceptions and conservation prospects. Environ. Conserv. 1997, 24, 338–343. [Google Scholar] [CrossRef]
  5. Wilson, C.J. Could we live with reintroduced large carnivores in the UK? Mammal Rev. 2004, 34, 211–232. [Google Scholar] [CrossRef]
  6. Ginsberg, J.R. Setting priorities for carnivore conservation: What makes carnivores different. In Carnivore Conservation; Conservation Biology Series; Gittleman, J.L., Funk, S.M., Macdonald, D.W., Wayne, R.K., Eds.; Cambridge University Press: Cambridge, UK, 2001; pp. 498–523. [Google Scholar]
  7. Chapron, G.; Kaczensky, P.; Linnell, J.D.C.; von Arx, M.; Huber, D.; Andrén, H.; López-Bao, J.V.; Adamec, M.; Álvares, F.; Anders, O.; et al. Recovery of large carnivores in Europe’s modern human-dominated landscapes. Science 2014, 346, 1517–1519. [Google Scholar] [CrossRef] [PubMed]
  8. Meuret, M.; Moulin, C.-H.; Bonnet, O.; Garde, L.; Nozières-Petit, M.-O.; Lescureux, N. Missing shots: Has the possibility of shooting wolves been lacking for 20 years in France’s livestock protection measures? Rangel. J. 2020, 42, 401–413. [Google Scholar] [CrossRef]
  9. Fritts, S.H.; Paul, J.W.; Mech, L.D.; Scott, D.P. Trends and Management of Wolf–Livestock Conflicts in Minnesota; Resource Publication 181; United States Department of the Interior Fish and Wildlife Service: Washington, DC, USA, 1992. [Google Scholar]
  10. Wydeven, A.P.; Treves, A.; Brost, B.; Wiedenhoeft, J. Characteristics of wolf packs in Wisconsin: Identification of traits influencing depredation. In People and Predators: From Conflict to Coexistence; Fascione, N., Delach, A., Smith, M.E., Eds.; Island Press: Washington, DC, USA, 2004; pp. 28–50. [Google Scholar]
  11. Treves, A.; Martin, K.A.; Wydeven, A.P.; Wiedenhoeft, J. Forecasting Environmental Hazards and the Application of Risk Maps to Predator Attacks on Livestock. BioScience 2011, 61, 451–458. [Google Scholar] [CrossRef]
  12. Štrbenac, A. Wolf Management Plan for Croatia: Towards Understanding and Addressing Key Issues in Wolf Management Planning in Croatia; State Institute for Nature Protection: Zagreb, Croatia, 2005. [Google Scholar]
  13. Gula, R. Wolf Depredation on Domestic Animals in the Polish Carpathian Mountains. J. Wildl. Manag. 2008, 72, 283–289. [Google Scholar] [CrossRef]
  14. Krofel, M.; Cerne, R.; Jerina, K. Effectiveness of wolf (Canis lupus) culling as a measure to reduce livestock depredations. Zb. Gozdarstva Lesar. 2011, 95, 11–22. [Google Scholar]
  15. Widman, M.; Elofsson, K. Costs of Livestock Depredation by Large Carnivores in Sweden 2001 to 2013. Ecol. Econ. 2018, 143, 188–198. [Google Scholar] [CrossRef]
  16. Remm, J.; Hindrikson, M. Estonian Conservation and Management Plan of Large Carnivores 2022–2031; Environmental Board: Pärnu, Estonia, 2022. [Google Scholar]
  17. Ozoliņš, J.; Žunna, A.; Ornicāns, A.; Done, G.; Stepanova, A.; Pilāte, D.; Šuba, J.; Lūkins, M.; Howlett, S.J.; Bagrade, G. Action Plan for Grey Wolf Canis lupus Conservation and Management; LSFRI Silava: Salaspils, Latvia, 2017. [Google Scholar]
  18. Game Resources. Available online: https://www.vmd.gov.lv/lv/es-sfera-esoso-sugu-monitorings (accessed on 12 December 2022).
  19. Frank, L.G.; Woodroffe, R. Behaviour of carnivores in exploited and controlled populations. In Carnivore Conservation; Conservation Biology Series; Gittleman, J.L., Funk, S.M., Macdonald, D.W., Wayne, R.K., Eds.; Cambridge University Press: Cambridge, UK, 2001; pp. 419–442. [Google Scholar]
  20. Treves, A.; Krofel, M.; McManus, J. Predator control should not be a shot in the dark. Front. Ecol. Environ. 2016, 14, 380–388. [Google Scholar] [CrossRef]
  21. Berger, K.M. Carnivore-Livestock Conflicts: Effects of Subsidized Predator Control and Economic Correlates on the Sheep Industry. Conserv. Biol. 2006, 20, 751–761. [Google Scholar] [CrossRef] [PubMed]
  22. Harper, E.K.; Paul, W.J.; Mech, L.D.; Weisberg, S. Effectiveness of Lethal, Directed Wolf-Depredation Control in Minnesota. J. Wildl. Manag. 2008, 72, 778–784. [Google Scholar] [CrossRef]
  23. Bruns, A.; Waltert, M.; Khorozyan, I. The effectiveness of livestock protection measures against wolves (Canis lupus) and implications for their co-existence with humans. Glob. Ecol. Conserv. 2020, 21, e00868. [Google Scholar] [CrossRef]
  24. Musiani, M.; Mamo, C.; Boitani, L.; Callaghan, C.; Gates, C.C.; Mattei, L.; Visalberghi, E.; Breck, S.; Volpi, G. Wolf depredation trends and the use of fladry barriers to protect livestock in western North America. Conserv. Biol. 2003, 17, 1538–1547. [Google Scholar] [CrossRef]
  25. Eklund, A.; López-Bao, J.; Tourani, M.; Chapron, G.; Frank, J. Limited evidence on the effectiveness of interventions to reduce livestock predation by large carnivores. Sci. Rep. 2017, 7, 2097. [Google Scholar] [CrossRef]
  26. Kaartinen, S.; Luoto, M.; Kojola, I. Carnivore-livestock conflicts: Determinants of wolf (Canis lupus) depredation on sheep farms in Finland. Biodivers Conserv. 2009, 18, 3503–3517. [Google Scholar] [CrossRef]
  27. Šuba, J.; Žunna, A.; Bagrade, G.; Done, G.; Lūkins, M.; Ornicāns, A.; Pilāte, D.; Stepanova, A.; Ozoliņš, J. Closer to Carrying Capacity: Analysis of the Internal Demographic Structure Associated with the Management and Density Dependence of a Controlled Wolf Population in Latvia. Sustainability 2021, 13, 9783. [Google Scholar] [CrossRef]
  28. Žunna, A.; Ozoliņš, J.; Pupila, A. Food habits of the wolf Canis lupus in Latvia based on stomach analyses. Est. J. Ecol. 2009, 58, 141–152. [Google Scholar] [CrossRef]
  29. Bagrade, G.; Ruņģis, D.I. Latvian State Forest Research Institute Silava, Salaspils, Latvia. 2023; unpublished work. [Google Scholar]
  30. Population of Game Species. Available online: https://www.vmd.gov.lv/lv/medijamo-dzivnieku-populacijas (accessed on 1 March 2023).
  31. Linnell, J.D.C.; Odden, J.; Mertens, A. Mitigation methods for conflicts associated with carnivore depredation on livestock. In Carnivore Ecology and Conservation: A Handbook of Techniques; Boitani, L., Powell, R.A., Eds.; Oxford University Press: New York, NY, USA, 2012; pp. 314–332. [Google Scholar]
  32. Reinhardt, I.; Rauer, G.; Kluth, G.; Kaczensky, P.; Knauer, F.; Wotschikowsky, U. Livestock protection methods applicable for Germany—A Country newly recolonized by wolves. Hystrix 2012, 23, 62–72. [Google Scholar] [CrossRef]
  33. Salvatori, V.; Mertens, A.D. Damage prevention methods in Europe: Experiences from LIFE nature projects. Hystrix 2012, 23, 73–79. [Google Scholar] [CrossRef]
  34. Klevezal, G.A. Age-Related Structures in Zoological Studies of Mammals; Nauka: Moscow, Russia, 1988. (In Russian) [Google Scholar]
  35. R Core Team. R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing, Vienna, Austria, 2017. Available online: https://www.R-project.org/ (accessed on 2 December 2021).
  36. Ciucci, P.; Boitani, L. Wolf and Dog Depredation on Livestock in Central Italy. Wildl. Soc. B 1998, 26, 504–514. [Google Scholar]
  37. Blanco, J.C.; Cortes, Y. Wolf recolonization of agricultural areas in Spain. In Proceedings of the Beyond 2000: Realities of Global Wolf Restoration, Duluth, MN, USA, 23–26 February 2000. [Google Scholar]
  38. Adamič, M.; Kobler, A.; Korenjak, A.; Marinčič, A.; Zafran, J. The recovery of the wolf (Canis lupus) in Slovenia. Beitr. Jagd Wildforschung 2001, 26, 85–94. [Google Scholar]
  39. Balčiauskas, L.; Balčiauskienė, L.; Volodka, H. Preliminary assessment of damage caused by the wolf in Lithuania. Acta Zool. Lit. 2002, 12, 419–427. [Google Scholar] [CrossRef]
  40. Kaczensky, P.; Chapron, G.; von Arx, M.; Huber, D.; Andrén, H.; Linnell, J. (Eds.) Status, Management and Distribution of Large Carnivores—Bear, Lynx, Wolf and Wolverine—In Europe; Part 2—Species Country Reports; IUCN/SSC Large Carnivore Initiative for Europe; The Institute of Applied Ecology, Istituto di Ecologia Applicata: Rome, Italy, 2013. [Google Scholar]
  41. Kaczensky, P. Large carnivore-livestock conflicts in Europe; Munich Wildlife Society: Munich, Germany, 1996. [Google Scholar]
  42. Fryxell, J.M.; Sinclair, A.R.E.; Caughley, G. Wildlife Ecology, Conservation, and Management; John Wiley & Sons: New York, NY, USA, 2014. [Google Scholar]
  43. Adams, L.G.; Stephenson, R.O.; Dale, B.W.; Ahgook, R.T.; Demma, D.J. Population dynamics and harvest characteristics of wolves in the Central Brooks Range, Alaska. Wildl. Monogr. 2008, 170, 1–25. [Google Scholar] [CrossRef]
  44. Brainerd, S.M.; Andren, H.; Bangs, E.E.; Bradley, E.H.; Fontaine, J.A.; Hall, W.; Iliopoulos, Y.; Jimenez, M.D.; Jozwiak, E.A.; Liberg, O.; et al. The effects of breeder loss on wolves. J. Wildl. Manag. 2008, 72, 89–98. [Google Scholar] [CrossRef]
  45. Mech, L.D. The Wolf: The Ecology and Behaviour of an Endangered Species; University of Minnesota Press: Minneapolis, MN, USA, 1970. [Google Scholar]
  46. Ozoliņš, J.; Stepanova, A.; Žunna, A.; Bagrade, G.; Ornicāns, A. Wolf hunting in Latvia in the light of population continuity in the Baltics. In Beiträge zur Jagd- und Wildforschung, Band 36; Stubbe, M., Ed.; Gesellschaft für Wildtier- und Jagdforschung e.V.: Halle/Saale, Deutschland, 2011; pp. S. 93–104. [Google Scholar]
  47. Ozoliņš, J.; Žunna, A.; Howlett, S.J.; Bagrade, G.; Pilāte, D.; Ornicāns, A.; Pēterhofs, E. Population dynamics of large mammals in Latvia with an emphasis on prey-predator interactions. In Beiträge zur Jagd- und Wildforschung. Band 41; Stubbe, M., Ed.; Gesellschaft für Wildtier- und Jagdforschung e.V.: Halle/Saale, Deutschland, 2016; pp. S. 59–73. [Google Scholar]
  48. Latham, J.; Staines, J.B.W.; Gorman, M.L. Correlations of red (Cervus elaphus) and roe (Capreolus capreolus) deer densities in Scottish forests with environmental variables. J. Zool. 1997, 242, 681–704. [Google Scholar] [CrossRef]
  49. Richard, E.; Gaillard, J.M.; Saïd, S.; Harmann, J.-L.; Klein, F. High red deer density depresses body mass of roe deer fawns. Oecologia 2010, 163, 91–97. [Google Scholar] [CrossRef] [PubMed]
  50. Jędrzejewska, B.; Jędrzejewski, W. Predation in Vertebrate Communities: The Białowieża Primeval Forest as a Case Study; Springer: Berlin/Heidelberg, Germany, 1998; Volume 135. [Google Scholar]
  51. Andersone, Ž.; Ozolinš, J. Public perception of large carnivores in Latvia. Ursus 2004, 15, 181–187. [Google Scholar] [CrossRef]
  52. Žunna, A.; Bagrade, G.; Ozoliņš, J. Attitudes of the General Public and Hunters Towards Wolves in Latvia; Its Predictors and Changes Over Time. Proc. Latv. Acad. Sci. B Nat. Exact Appl. Sci. 2020, 74, 280–286. [Google Scholar] [CrossRef]
  53. Bath, A.J. Attitudes of Various Interest Groups in Wyoming Toward Wolf Reintroduction in Yellowstone National Park. Master’s Thesis, University of Wyoming, Laramie, WY, USA, 1987. [Google Scholar]
  54. Kellert, S.R. The Public and the Timber Wolf in Minnesota. Anthrozoös 1987, 1, 100–109. [Google Scholar] [CrossRef]
  55. Bath, A.J.; Buchanan, T. Attitudes of interest groups in Wyoming towards wolf restoration in Yellowstone National Park. Wildl. Soc. B 1989, 17, 519–525. [Google Scholar]
  56. Blanco, J.C.; Reig, S.; Cuesta, L. Distribution, status and conservation problems of the wolf Canis lupus in Spain. Biol. Conserv. 1992, 60, 73–80. [Google Scholar] [CrossRef]
  57. Williams, C.K.; Ericsson, G.; Heberlein, T.A. A quantitative summary of attitudes toward wolves and their reintroduction (1972–2000). Wildl. Soc. B 2002, 30, 575–584. [Google Scholar] [CrossRef]
  58. Røskaft, E.; Händel, B.; Bjerke, T.; Kaltenborn, B.P. Human attitudes towards large carnivores in Norway. Wildl. Biol. 2007, 13, 172–185. [Google Scholar] [CrossRef]
  59. Sponarski, C.C.; Semeniuk, C.; Glikman, J.A.; Bath, A.J.; Musiani, M. Heterogeneity among Rural Resident Attitudes Toward Wolves. Hum. Dimens. Wildl. 2013, 18, 239–248. [Google Scholar] [CrossRef]
  60. Ericsson, G.; Heberlein, T.A. Attitudes of hunters, locals and general public in Sweden now that wolves are back. Biol. Conserv. 2003, 111, 149–159. [Google Scholar] [CrossRef]
  61. Heberlein, T.A. Navigating Environmental Attitudes; Oxford University Press: New York, NY, USA, 2012. [Google Scholar]
  62. Sillero-Zubiri, C.; Laurenson, K. Interactions between carnivores and local communities: Conflict or co-existence. In Carnivore conservation; Conservation Biology Series; Gittleman, J.L., Funk, S.M., Macdonald, D.W., Wayne, R.K., Eds.; Cambridge University Press: Cambridge, UK, 2001; pp. 282–312. [Google Scholar]
Figure 1. Borders of local forestry units of the Latvian State Forest Service (red) and statistical regions (black).
Figure 1. Borders of local forestry units of the Latvian State Forest Service (red) and statistical regions (black).
Sustainability 15 08509 g001
Figure 2. Number of reported wolf attacks and affected livestock (a) and number of affected (i.e., killed, injured, or lost) sheep per attack (b), showing minimum, maximum, median, inter-quartile range, and number of cases in Latvia from 2004 to 2022 (data from the Latvian State Forest Service).
Figure 2. Number of reported wolf attacks and affected livestock (a) and number of affected (i.e., killed, injured, or lost) sheep per attack (b), showing minimum, maximum, median, inter-quartile range, and number of cases in Latvia from 2004 to 2022 (data from the Latvian State Forest Service).
Sustainability 15 08509 g002
Figure 3. Increase in mean number of affected sheep per reported wolf attack in Latvia from 2004 to 2022 (solid and dashed lines indicate the linear trend and 95% confidence intervals, respectively; data from the Latvian State Forest Service).
Figure 3. Increase in mean number of affected sheep per reported wolf attack in Latvia from 2004 to 2022 (solid and dashed lines indicate the linear trend and 95% confidence intervals, respectively; data from the Latvian State Forest Service).
Sustainability 15 08509 g003
Figure 4. Timing of reported wolf attacks on livestock throughout the year (earliest, latest, median, inter-quartile range, and number of cases). Gray shading indicates period of closed hunting season for wolves.
Figure 4. Timing of reported wolf attacks on livestock throughout the year (earliest, latest, median, inter-quartile range, and number of cases). Gray shading indicates period of closed hunting season for wolves.
Sustainability 15 08509 g004
Figure 5. Observed and predicted number of reported attacks by wolves (a) and number of affected sheep (b) per year according to estimated wolf density (the solid and dashed lines indicate the expected number and 95% confidence intervals according to a negative binomial regression, respectively; data from the Latvian State Forest Service).
Figure 5. Observed and predicted number of reported attacks by wolves (a) and number of affected sheep (b) per year according to estimated wolf density (the solid and dashed lines indicate the expected number and 95% confidence intervals according to a negative binomial regression, respectively; data from the Latvian State Forest Service).
Sustainability 15 08509 g005
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

Šuba, J.; Žunna, A.; Bagrade, G.; Done, G.; Ornicāns, A.; Pilāte, D.; Stepanova, A.; Ozoliņš, J. Does Wolf Management in Latvia Decrease Livestock Depredation? An Analysis of Available Data. Sustainability 2023, 15, 8509. https://0-doi-org.brum.beds.ac.uk/10.3390/su15118509

AMA Style

Šuba J, Žunna A, Bagrade G, Done G, Ornicāns A, Pilāte D, Stepanova A, Ozoliņš J. Does Wolf Management in Latvia Decrease Livestock Depredation? An Analysis of Available Data. Sustainability. 2023; 15(11):8509. https://0-doi-org.brum.beds.ac.uk/10.3390/su15118509

Chicago/Turabian Style

Šuba, Jurģis, Agrita Žunna, Guna Bagrade, Gundega Done, Aivars Ornicāns, Digna Pilāte, Alda Stepanova, and Jānis Ozoliņš. 2023. "Does Wolf Management in Latvia Decrease Livestock Depredation? An Analysis of Available Data" Sustainability 15, no. 11: 8509. https://0-doi-org.brum.beds.ac.uk/10.3390/su15118509

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