Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Cormorant Catch Concerns for Fishers: Estimating the Size-Selectivity of a Piscivorous Bird

Correction

2 Jan 2014: Troynikov V, Whitten A, Gorfine H, Pūtys Ž, Jakubavičiūtė E, et al. (2014) Correction: Cormorant Catch Concerns for Fishers: Estimating the Size-Selectivity of a Piscivorous Bird. PLOS ONE 9(1): 10.1371/annotation/b59e978e-722d-4835-9350-f2f92d6f3738. https://doi.org/10.1371/annotation/b59e978e-722d-4835-9350-f2f92d6f3738 View correction

Abstract

Conflict arises in fisheries worldwide when piscivorous birds target fish species of commercial value. This paper presents a method for estimating size selectivity functions for piscivores and uses it to compare predation selectivities of Great Cormorants (Phalacrocorax carbo sinensis L. 1758) with that of gill-net fishing on a European perch (Perca fluviatilis L. 1758) population in the Curonian Lagoon, Lithuania. Fishers often regard cormorants as an unwanted “satellite species”, but the degree of direct competition and overlap in size-specific selectivity between fishers and cormorants is unknown. This study showed negligible overlap in selectivity between Great Cormorants and legal-sized commercial nets. The selectivity estimation method has general application potential for use in conjunction with population dynamics models to assess fish population responses to size-selective fishing from a wide range of piscivorous predators.

Introduction

Ecological interactions between cormorants, fish, and fishers are problematic globally [1][7]. Among those studies describing the effects of increasing cormorant numbers on fish populations and fisheries in Europe [5], [8] and North America [9][11] many suggest cormorants have a negative impact, and that is certainly a common perception of fishers [7]. Other studies have demonstrated that cormorant impact can be negligible [12][14], and may not lead to declines of commercial and non-commercial fish species. Cormorants frequently target commercially and recreationally important fish [7], [15], [16] and increases in the size of piscivorous bird colonies can obviously place pressure on fish stocks. A key requirement to understanding the effects of cormorant populations on fish stocks is the size-specific selectivity of their predation, however estimating parameters to describe predation selectivity can be difficult.

Previous studies of cormorant predation have relied predominantly on qualitative and empirical assessments of predation selectivity and population level effects [17]. A number of studies identified length distributions for all prey fish (e.g. [18][21]), whereas others focused on the most important species in the diet (e.g. [7], [8], [22][24]). Several assessed those fish age groups most affected by cormorant predation (see [25], [26]) or considered seasonal differences in selectivity [27], [28], and one focussed on size and species selectivity effects on fish community structure [29].

Most studies on competition between cormorants and fisheries have concentrated on fish consumption rather than catch composition and size structure (e.g. [8], [13], [19], [20], [30]). Despite some research estimating potential or existing direct competition, to our knowledge, evaluation of the quantitative overlap of cormorant diet and fisheries has not been undertaken previously. Among dietary studies it has been found that cormorants exploit fish mostly smaller than fishery minimum sizes [31]. In a study comparing cormorant diets with anglers' catches in the Great Lakes USA, the authors concluded that each group captured different sizes of fish with only a minor overlap [32]. A similar study in Lithuania compared the size distribution of regurgitated fish with commercial gillnet catches and concluded that direct competition was low [7]. Examining the potential for direct competition with the commercial fishery in Lithuania's Curonian Lagoon, it was found that although legal-sized fish comprised more than one quarter of the cormorants' diet, the composition was dominated by species of low commercial value [24].

Conventional methods for directly estimating selectivity parameters from length-composition data restrict the estimation process to a single type of fishing gear: They are based on the assumption that data are collected by gear types with equal or comparable units of fishing effort (e.g. [33]). Integrated analysis methods allow selectivity parameters for multiple types of gear, and therefore, for multiple functional forms, to be estimated [34]. Such methods combine several sources of data into a single analysis, so that parameters to describe all or most components of a population and fishery can be estimated as part of a single likelihood. Conversely, direct methods of estimating population and fishery related parameters lead to a two-step approach: Data are analysed in their raw form to produce parameter estimates and those estimates are used in population or stock assessment models [35]. The latter approach, whilst being simpler and less computationally and data intense, also allows for selectivity parameters to be estimated as independent properties of fishing gear without statistical correlation with other parameters in a population model.

In this paper we describe methods to directly estimate the parameters of selectivity functions using length composition data; simultaneous estimation of selectivity parameters for both fishing nets and piscivorous predators is achieved without an assumption of equality in fishing power among nets or predators. The methods also allow the estimation of parameters for multiple functional forms of selectivity and provide a basis for incorporating predation selectivity into analyses of cohort dynamics or general size-age-structured models. Such models could in turn be used to determine the effects of piscivores on fishable stocks. This is essential for developing an understanding of actual, compared with perceived, levels of competition between cormorants and fishers.

The methods developed in this study were applied to data relating to commercial fishing and predation by Great Cormorants (Phalacrocorax carbo sinensis L. 1758) from a colony at Juodkrantė, on the western coast of the Curonian Lagoon in Lithuania. In an investigation of cormorant diets in this area, 76% of surveyed birds were observed to have consumed European perch [24]. European perch is the fourth most abundant fish species in the Curonian Lagoon [36] and one of the most important commercial species, with annual commercial landings for the combined Lithuanian and Kaliningrad (Russia) regions of the Lagoon of 100–140 tonnes live weight (Fisheries Dept., Ministry of Agriculture, Lithuania). Comparisons between the estimates of predation selectivity of cormorants and fishing selectivity of commercial gillnets offer useful insights into the degree of competition between commercial fishers and this increasingly abundant piscivorous bird.

Materials and Methods

Ethics statement

All animal work was conducted on public land and waterways, did not involve protected, threatened or endangered species, and complied with relevant national and international guidelines and legislation. Permits to perform fish surveys at the Curonian Lagoon, Permit No. 01 in 2009 and Permit No. 001 in 2010, were issued to the Nature Research Centre (NRC) by The Environmental Protection Agency under the Ministry of the Environment of the Republic of Lithuania. All live animals removed from sampling nets were humanely killed by spiking the head, with either a dissecting needle for smaller fish or a scalpel blade for larger specimens, to destroy brain function. Disposal of the carcasses of unpreserved specimens was in accordance with NRC Laboratory of Marine Ecology bio-security protocols. Institutional animal ethics considerations and approvals were met via the Research Council of Lithuania project application procedural requirements.

Sampling

Perch were sampled from the northern part of the Curonian Lagoon during the summer periods (July) of 2009 and 2010. The Lagoon is a 1584 km2 freshwater basin, located in the south-east Baltic Region and connected to the Baltic Sea through the 500 m-wide Klaipėda channel (Fig. 1). For a substantial part of the year the Lagoon is frozen, with a thick ice cover, and during autumn and spring some fish migrate outside the study area, making seasonal sampling logistically and biologically unrealistic. Sampling sites were chosen from among the deeper areas of the Lagoon, utilized as habitat by most perch populations, and to the north and south along the Nerija coast from the Juodkrantė Great Cormorant colony, as well as the southern part of the Lithuanian zone where cormorants are more active. The shallow banks and areas in the north-eastern part of the Lagoon were not sampled because they are less used by fisheries.

thumbnail
Figure 1. Study area, showing the Juodkrantė Great Cormorant (Phalacrocorax carbo L. 1758) colony (▪), and perch (Perca fluviatilis L. 1758) sampling sites (•).

https://doi.org/10.1371/journal.pone.0077518.g001

Samples were taken from 14 sites using a set of 11 spun nylon gillnets (mesh sizes: 14, 17, 21.5, 25, 30, 33, 38, 45, 50, 60 & 70 mm) [37], as described in Troynikov et al. (2011) [38]. A minimum size limit of 18 cm (LT) applies to the commercial fishery and the mesh size of gillnets used commercially in the Curonian Lagoon is limited to 40–45 mm (knot-to-knot).

Pellets of Great Cormorants at the Juodkrantė colony were collected for dietary analysis during the breeding season from the beginning of March to the end of July. In each instance 15–20 pellets were sampled at 10 day intervals in 2009 (n = 298) and 2010 (n = 276). Pellets were analysed according to the standard methodology described by Carss et al. (1997) [39]; for a detailed description see Pūtys and Zarankaitė (2010) [24]. Fish lengths and weights were reconstructed using allometric relationships, estimated from a reference collection. The length distribution data obtained from these pellets and the sampling nets are listed in Table 1.

thumbnail
Table 1. Observed numbers by length-class of European perch caught in gill nets of increasing mesh size (G (mm)), and by Great Cormorants (GC).

https://doi.org/10.1371/journal.pone.0077518.t001

Estimation procedure

This study extends the maximum likelihood estimation procedure for selectivity functions developed by Kirkwood and Walker (1986) [40]. The likelihood function is based on the assumption that for each gear type i and length-class j, the catches nij  =  fiμj Sij are independent observations from a Poisson distribution with mean μj Sij, where μj is the relative proportion in the population from length-class j, Sij is a selectivity function and fi is fishing power;(1)

Then, with an assumption that each gear has fishing power equal to 1, at the maximum of the log-likelihood function;(2)

and parameter μj  =  Σ nij/Σ Sij, i =  1,..., I (see [40]). The selectivity function S is proportional to the conditional probability P (“catch ”| j) of fish to be captured given size-class j. Thus, the fundamental relationship between the selectivity function Sij, the size distribution in the sampled fish population fj (prior distribution) and size distribution pij in a sample by gear type i (posterior distribution) can be formulated as the Bayesian theorem:(3)

Given that fishing power may only affect sample size, but not the size distribution in a sample, values pij are proportional to Pr (N =  nij) in model (1) and size-frequency data can be used in estimator (2). This consideration has an important practical implication; it allows for the use of normalised size-frequency data collected by different types of fishing gears (including Great Cormorants) that may have incomparable fishing efforts and different fishing powers. The only values that need to be estimated in this case are the parameters of our nominated selectivity curves. Following Kirkwood and Walker (1986) [40], gamma probability density functions (pdf), rescaled to modal values of one, were used to describe fishing mesh net selectivity. Accordingly, the functional form used to model selectivities as a function of length, l, is:(4)where g denotes fishing gear.

The length at maximum selectivity for net i is proportional to the mesh size mi, and variance θ2 is a constant over different nets (principle of geometrical similarity) so that:(5)

By adopting a specific functional form for the fishing gear selectivity function and solving equations (5) with respect to α and β, only two parameters, θ1 and θ2, need to be estimated. However, the Great Cormorant “gear” type cannot be rescaled to any mesh size and the selectivity parameters must have independent values. An investigation of the cormorant size-frequency catch data reveals a distribution with a heavy right-hand tail, indicative of a lognormal shaped selectivity function. Therefore a log-normal pdf, rescaled to a modal value of one, was used to describe cormorant predation selectivity. The functional form used to model cormorant selectivity as a function of length, l, is(6)

with mode at and variance;; and where c denotes cormorants.

Thus in distinction from Kirkwood and Walker (1986) [40], there are two extra parameters to estimate. The log-likelihood function of the whole data set is:(7)

Maximum likelihood estimates for parameters of interest, θ1, θ2, θ3, and θ4, were obtained using the Nelder-Mead simplex method [41]. All statistical analyses were performed using the R software environment (Version 2.13.0, R Development Core Team, 2011). The Nelder-Mead method was employed using the general-purpose optimization function optim, available in R standard libraries. All source code, together with example usage scripts and the data used in this paper are available with instructions for their use at https://www.github.com/awhitten/corselect.

Results

European perch catch numbers were converted to normalised frequency data before applying the maximum likelihood estimator (7). Initial values of the parameters of interest (a requirement for the optimization procedure) were approximated using the following heuristics: as starting points, (1) the mode of an estimated selectivity function can be set close to the value corresponding to the mode of the length-frequency distribution from any particular gear, and (2) the value of the variance can be set around two-thirds the width of that same frequency distribution. The maximum likelihood estimates of the four parameters of interest are shown in Table 2. Parameter estimates were used to determine a mode and variance for each of the fishing gears and cormorant predation selectivity functions according to the equations presented under Estimation Procedure. The variance for all fishing gear selectivities was 12.10. The mode and variance for cormorant predation were 8.00 cm and 10.8 respectively.

thumbnail
Table 2. Parameter and standard error estimates for theta one and two, used to describe gamma shaped selectivity curves for categories of fishing gear, and for theta three and theta four, used to describe a log-normal shaped selectivity curve for cormorant predation.

https://doi.org/10.1371/journal.pone.0077518.t002

Cormorant selectivity for the perch population in the Curonian Lagoon was similar to that of the 14 mm gillnet with around 70% area of overlap. However, apart from the differences in modes (7.1 cm for Great Cormorants compared with 9.9 for the 14 mm gillnet) the selectivity estimates differed in the tails of their distributions, with the cormorants having a much steeper left tail and a slightly steeper right tail (Fig. 2). Overlap with the 38 mm gillnet was <2% and there was negligible overlap with the 45 mm commercial sized nets showing cormorants rarely consumed fish in size classes targeted by commercial fishers. Indeed, the intersection point between cormorants and 38 mm selectivity curves occurs close to the minimum legal size limit of 18 cm total length.

thumbnail
Figure 2. Estimated relative selectivity as a function of length of European perch for eight gillnets with increasing mesh size (dashed lines, marked “G” with mesh size in mm) and Great Cormorants (solid line, marked “GC”).

https://doi.org/10.1371/journal.pone.0077518.g002

Discussion

Our results show there is limited direct interaction between Great Cormorant predation and commercial fishing in the Curonian Lagoon. Only a small proportion of perch (<2% of the area of intersection between the commercial gill net and cormorant selectivity functions) selected by commercial-sized gill nets were of a size that are regularly consumed by cormorants. Although no recreational catch data were available, the very small number of legal sized perch in the cormorants' diet suggests that direct interaction with this fishing sector should also be low.

Estimated selectivity curves indicate little to no direct competition between cormorants and fishers for commercial-sized perch, but they do reveal strong selection of sub-legal sized perch by cormorants. The estimated average annual consumption rate of perch by Great Cormorants in the Lithuanian part of the Curonian Lagoon is 118 t which is 2.5 times commercial landings and represents 18% of the fish biomass that cormorants from this area consume (Ž. Pūtys unpubl. data). Our findings indicate almost all this consumption is of juvenile perch, yet to reach the commercially targeted component of the stock.

Substantial consumption of juvenile perch by cormorants might be expected to affect recruitment to the fishery and therefore catches, but landings throughout the entire Lagoon increased by 22% during the period 2007–2011. Whilst increased catches do not necessarily reflect changes in abundance, others have shown substantial fish consumption by cormorants does not lead to stock declines in prey populations [7], [12]. There may be several mechanisms responsible for such effects: We outline one scenario based on the literature. If cormorant numbers have increased in proportion to prey abundance due to a reduction in other predators [42], [43], such as a depletion of pikeperch (Sander lucioperca L. 1758) by fishing [44], then reductions in fish population densities by cormorants may stimulate growth and fecundity rates in perch through compensatory mechanisms related to reduced intra-specific competition in semi-closed fishing grounds, thereby improving the rate of recruitment (e.g. [45]). Under this hypothesis, as cormorants reduce their prey to more stable densities, there will be a gradual attrition in their numbers commensurate to overall prey availability [46].

This study also provides a computational method for the simultaneous estimation of fish length selectivity parameters for multiple capture methods, no matter whether they are predators or types of fishing equipment. The approach has general ecological application potential, and could be used as a starting point when comparing the specificity and relative efficiency with which piscivorous birds or other animals predate upon fish. Although several methods exist for direct estimation of selectivity parameters (e.g. [33], [40]), such studies have previously been restricted to methods that deal only with single gear categories (and thus single functional forms). The method presented here allows multiple gear types and predators to be included within a single estimation procedure, avoids the need for assumed equality in fishing power, and enables direct estimation from length composition data alone. This study thus presents a key step toward modelling and understanding the impacts of competing predators on fish stocks, without the need to extend analyses to highly complex and data-hungry integrated analysis methods. Importantly, data collection costs can be reduced with this method as it permits the use of historical data that may not otherwise have been intended for gear selectivity estimation.

The methods developed in this study cannot, in isolation, test the assertion by fishers that increased numbers of Great Cormorants in the colony at Juodkrantė are negatively affecting commercially important fish stocks in the Curonian Lagoon. However, by providing important inputs to age- or size-structured population models, these methods enable a two-step approach to the examination of population level effects of cormorant predation on European perch stocks. Age- or size-structured population models, including simulation models, could enable better estimates of total mortality on young fish and thus improve predictions of future recruitment to the fishery. Indeed, models that explicitly consider both gear selectivity and predator induced natural mortality can significantly improve stock assessment models, especially those that address ecosystem concerns [47]. Thus future research and stock assessment projects would benefit from incorporating the methods presented in this study. Such analyses would likely improve understanding of the full effects of size-selective perch predation by Great Cormorants.

Perceived versus actual effects of cormorant populations on fishable stocks can cause unnecessary concern for commercial fishers. Understanding the wide range of ecological processes at play is critical to future management of their fishery. In addition to continuing to monitor the abundance and age-structure of fish populations, concurrent surveys of cormorant numbers and their diets would permit estimates of annual perch consumption and support population modelling. Together these efforts would help determine the effect of cormorant predation on the dynamics of the fishable stock, and whether it might be important to reduce that predation in order to protect future recruitment to the fishery.

The selectivity estimation procedure in this study provides a useful link in determining the relative impacts of fishing and natural predation on fisheries, not only for the Curonian Lagoon, but for fisheries generally. Further development of this method is recommended to allow for the estimation of different selectivity curves, and for methods to determine the best selectivity functions relating to different fishing gears or predators of interest.

Acknowledgments

The assistance of R. Bukontaitė, Laboratory of Marine Ecology of the Nature Research Centre, in drafting Figure 1 is gratefully acknowledged. This manuscript was also improved following the helpful and insightful comments provided by reviewers.

Author Contributions

Conceived and designed the experiments: LL ZP HG. Performed the experiments: LL ZP HG EJ JD. Analyzed the data: VT AW HG. Contributed reagents/materials/analysis tools: VT AW ZP. Wrote the paper: VT HG AW. Designed the software used in analysis: AW.

References

  1. 1. Wires LR, Carss DN, Cuthbert FJ, Hatch JJ (2003) Transcontinental connections in relation to cormorant – fisheries conflicts: perceptions and realities of a “bete noire” (black beast) on both sides of the Atlantic. Vogelwelt 124: 389–400.
  2. 2. Grémillet D, Enstipp M, Boudiffa M, Liu H (2006) Do cormorants injure fish without eating them? An underwater video study. Mar Biol 148: 1081–1087
  3. 3. Kloskowski J (2011) Human-wildlife conflicts at pond fisheries in eastern Poland: perceptions and management of wildlife damage. Eur J Wildl Res 57: 295–304
  4. 4. Krohn WB, Allen R, Moring JR, Hutchinson AE (1995) Double-Crested Cormorants in New England: Population and management histories. Waterbirds 18: 99–109.
  5. 5. Steffens W (2010) Great Cormorant – substantial danger to fish populations and fishery in Europe. Bulg J Agric Sci 16: 322–331.
  6. 6. Yorio P, Quintana F, Dell'Arciprete P, Gonzalez-Zevallos D (2010) Spatial overlap between foraging seabirds and trawl fisheries: implications for the effectiveness of a marine protected area at Golfo San Jorge, Argentina. Bird Conserv Int 20: 320–334
  7. 7. Žydelis RR, Kontautas AA (2008) Piscivorous birds as top predators and fishery competitors in the lagoon ecosystem. Hydrobiologia 611: 45–54
  8. 8. Stewart DC, Middlemas SJ, Gardiner WR, Mackay S, Armstrong JD (2005) Diet and prey selection of cormorants (Phalacrocorax carbo) at Loch Leven, a major stocked trout fishery. J Zool 267: 191–201
  9. 9. Fielder DG (2008) Examination of factors contributing to the decline of the yellow perch population and fishery in Les Cheneaux Islands, Lake Huron, with emphasis on the role of Double-crested Cormorants. J Great Lakes Res 34: 506–523
  10. 10. Fielder D (2010) Response of yellow perch in Les Cheneaux Islands, Lake Huron to declining numbers of Double-crested Cormorants stemming from control activities. J Great Lakes Res 36: 207–214
  11. 11. Dorr BS, Moerke A, Bur M, Bassett C, Aderman T, et al. (2010) Evaluation of harassment of migrating Double-crested Cormorants to limit depredation on selected sport fisheries in Michigan. J Great Lakes Res 36: 215–223
  12. 12. Dalton CM, Ellis D, Post DM (2009) The impact of Double-crested Cormorant (Phalacrocorax auritus) predation on anadromous alewife (Alosa pseudoharengus) in south-central Connecticut, USA. Can J Fish Aquat Sci 66: 177–186
  13. 13. Engström H (2001) Long-time effects of cormorant predation on fish communities and fishery in a freshwater lake. Ecography 24: 127–138
  14. 14. Linn IJ, Campbell KI (1992) Interactions between White-breasted Cormorants Phalacrocorax carbo (Aves: Phalacrocoracidae) and the fisheries of Lake Malaŵi. J App Ecol 29: 619–634. Available: http://www.jstor.org/stable/2404470. Accessed 2013 Mar 13.
  15. 15. Švagždys A (2006) Conditions and opportunities of recreational fishing in the lower Nemunas river and the Curonian Lagoon. Žuvininkystė Lietuvoje (Fishery and Aquaculture in Lithuania) 6: 195–206.
  16. 16. Rudstam LG, Vandevalk AJ, Adams CM, Coleman JTH, Forney JL, et al. (2004) Cormorant predation and the population dynamics of walleye and yellow perch in Oneida Lake. Ecol Appl 14: 149–163
  17. 17. Irwin BJ, Treska TJ, Rudstam LG, Sullivan PJ, Jackson JR, et al. (2008) Estimating walleye (Sander vitreus) density, gear catchability, and mortality using three fishery independent data sets for Oneida Lake, New York. Can J Fish Aquat Sci 65: 1366–1378
  18. 18. Dirksen S, Boudewijn TJ, Noordhuis R, Marteijn ECL (1995) Cormorants Phalacrocorax carbo sinensis in shallow eutrophic freshwater lakes: prey choice and fish consumption in the non-breeding period and effects of large-scale fish removal. Ardea 83: 167–184. Available: http://www.avibirds.com/pdf/a/aalscholver22.pdf. Accessed 2013 Aug 1.
  19. 19. Britton JR, Cowx IG, Harvey JP, Davies JM, Holden T, et al.. (2003) Key factor analysis to assess cormorant depredation on inland fisheries in the UK. In: Cowx IG, editor. Interactions between fish and birds: Implications for management. 14–27.
  20. 20. Eschbaum R, Veber T, Vetemaa M, Saat T (2003) Do cormorants and fishermen compete for fish resources in the Väinameri (Eastern Baltic) Area? In: Cowx IG, editor. Interactions between fish and birds: Implications for management. 72–83.
  21. 21. Lilliendahl K, Solmundsson J (2006) Feeding ecology of sympatric European shags Phalacrocorax aristotelis and Great Cormorants P. carbo in Iceland. Mar Biol 149: 979–990
  22. 22. Suter W (1995) The effect of predation by wintering cormorants Phalacrocorax carbo on grayling Thymallus thymallus and trout (Salmonidae) populations: two case studies from Swiss rivers. J Appl Ecol 32: 29–46. Available: http://www.jstor.org/stable/2404413. Accessed 2013 July 13.
  23. 23. Coutin PC, Reside J (2003) Fish predation by Great Cormorants, Phalacrocorax carbo carboides, in the Gippsland Lakes, south-eastern Australia. In: Cowx IG, editor. Interactions between fish and birds: Implications for management. 196–210.
  24. 24. Pūtys Ž, Zarankaitė J (2010) Diet of the Great Cormorant (Phalacrocorax carbo sinensis) at the Juodkrantė colony, Lithuania. Acta Zool Lituan 20: 179–189.
  25. 25. Wires LR, Cuthbert FJ, Trexel DR, Joshi AR (2001) Status of the Double-crested Cormorant (Phalacrocorax auritus) in North America. Final Report to USFWS. 376 p.
  26. 26. Leopold MF, Van Damme CJG, van der Veer HW (1998) Diet of cormorants and the impact of cormorant predation on juvenile flatfish in the Dutch Wadden Sea. J Sea Res 40: 93–107
  27. 27. Lehikoinen A (2005) Prey-switching and diet of the Great Cormorant during the breeding season in the Gulf of Finland. Waterbirds 28: 515–515.
  28. 28. Čech M, Čech P, Kubečka J, Prchalova M, Draštik V (2008) Size selectivity in summer and winter diets of Great Cormorant (Phalacrocorax carbo): does it reflect season-dependent difference in foraging efficiency? Waterbirds 31: 438–447.
  29. 29. Engström H, Jonsson L (2003) Great Cormorant Phalacrocorax carbo sinensis diet in relation to fish community structure in a freshwater lake. Vogelwelt 124: 187–196.
  30. 30. Takahashi T, Kameda K, Kawamura M, Nakajima T (2006) Food habits of Great Cormorant Phalacrocorax carbo hanedae at Lake Biwa, Japan, with special reference to ayu Plecoglossus altivelis altivelis. Fisheries Sci 72: 477–484
  31. 31. Andersen SM, Teilmann J, Harders PB, Hansen EH, Hjollund D (2007) Diet of harbour seals and Great Cormorants in Limfjord, Denmark: interspecific competition and interaction with fishery. ICES J Marine Sc 64: 1235–1245
  32. 32. Diana JS, Maruca S, Low B (2006) Do increasing cormorant populations threaten sportfishes in the Great Lakes? A case study in Lake Huron. J Great Lakes Res 32: 306–320
  33. 33. Millar RB, Holst R (1997) Estimation of gillnet and hook selectivity using log-linear models. ICES J Mar Sci 54: 471–477
  34. 34. Methot RD Jr, Wetzel CR (2013) Stock synthesis: A biological and statistical framework for fish stock assessment and fishery management. Fish Res 142: 86–99
  35. 35. Maunder MN, Punt AE (2013) A review of integrated analysis in fisheries stock assessment. Fish Res 142: 61–74
  36. 36. Ådjers K, Appelberg M, Eschbaum R, Lappalainen A, Minde A, et al. (2006) Trends in coastal fish stocks of the Baltic Sea. Boreal Environ Res 11: 13–25.
  37. 37. Thoresson G (1993) Guidelines for coastal monitoring-fishery biology. Kustrapport 1993(1): 1–35.
  38. 38. Troynikov VS, Gorfine HK, Ložys L, Pūtys Z, Jakubavičiūtė E, et al. (2011) Parameterization of European perch Perca fluviatilis length-at-age data using stochastic Gompertz growth models. J Fish Biol 79: 1940–1949
  39. 39. Carss DN, Bevan RM, Bonetti A, Cherubini G, Davies J, et al.. (1997) Techniques for assessing cormorant diet and food intake: towards a consensus view. Supplementi di ricerche biologia selvaggina XXVI: 197–230. Proceedings 4th European Conference on Cormorants (1–3 November 1995, Bologna, Italy).
  40. 40. Kirkwood GP, Walker TI (1986) Gill net selectivities for gummy shark, Mustelus antarcticus Günther, taken in south-eastern Australian waters. Aust J Mar Freshwater Res 37: 689–697
  41. 41. Nelder JA, Mead R (1965) A simplex method for function minimization. Comput J 7: 308–313.
  42. 42. De Nie H (1995) Changes in the inland fish populations in Europe in relations to the increase of the Cormorant Phalacrocorax carbo sinensis. Ardea 83: 115–122.
  43. 43. Van Eerden MR, Koffijberg K, Platteeuw M (1995) Riding on the crest of the wave: possibilities and limitations for a thriving population of migratory cormorants Phalacrocorax carbo in man-dominated wetlands. Ardea 83: 1–9.
  44. 44. Mustamäki N, Bergström U, Ådjers K (2013) Pikeperch (Sander lucioperca (L.)) in decline: high mortality of three populations in the Northern Baltic Sea. Ambio (online; unassigned) doi: 10.1007/s13280-013-0429-z.
  45. 45. Rose KA, Cowan JH, Winemiller KO, Myers RA, Hillborn R (2001) Compensatory density dependence in fish populations: importance, controversy, understanding and prognosis. Fish Fish 2: 293–327
  46. 46. Van Eerden MR, Van Rijn S, Nordhuis R (2003) How zebra mussels Dreissena polymorpha, smelt Osmerus eperlanus and commercial fisheries interact with Great Cormorants Phalacrocorax carbo sinensis in Lake Ijsselmeer: the food web cascade revisited. Vogelwelt 124: 277–290.
  47. 47. Hollowed AB, Ianelli JN, Livingston PA (2000) Including predation mortality in stock assessments: a case study for Gulf of Alaska walleye pollock. ICES J Marine Sc 57: 279–293