Fisheries Management under Ecological Interdependence : The Case of European Hake and Blue Whiting Fishery

In mixed fisheries where there is an increased ecological interdependence between two or more target species, the possibilities of fishing these species must be determined jointly, since catches of one species will impact the natural growth not only of that species but of the others, as well. The objective of this paper is to develop a predator-prey model for two major species caught by the EU fishing fleet on European Union fishing grounds. The predator and prey’s population dynamics follow the Lotka-Volterra equation and are assumed to be logic-based, and a lineal interaction between the predator and prey populations is assumed. Optimal single owner multispecies fishery management is analysed, and the applied model is solved, obtaining the equilibrium value of biomasses, catches, and net benefits of the mixed fishery. The results show that the MSY of the predator is higher than that estimated by ICES, while MSY for prey is lower. The sensitivity analysis of the results reveals that the biomass levels of both species decrease as the discount rate increases, while catches increase.


Introduction
The complex series of interactions among a fishery's different species has led an increasing number of scientists to recommend multispecies approaches to fishery management [1] [2].In practice, fisheries are exploited as a multispecies resource [3].Biological interactions among fish stocks can dominate the structuring of marine ecosystems and so can strongly affect the population dynamics economic model by considering the implications of dependencies, both ecological and technological, that arise in multispecies fisheries [12].
Along these lines, attempts have been made to account for ecological and technological interdependencies that result from the extractive activities of fishing fleets, which differ not only in terms of their technological characteristics but also in terms of their targeted species and fish population components [8] [12] [13] [14].Technological interdependence occurs when fleets of different characteristics-for example, as regards fishing power and/or types of skills usedtarget different components (e.g., juveniles and adults) of the same stock, or different species of a mixed stock [15], and therefore have distinct effects on species populations [16].Another example of technological interdependence is the fleet that catches species coexisting within the same space regardless of whether or not they are interdependent on an ecological level.Ecological interdependence occurs when there is a certain relationship between two exploited stocks [17] [18] or between two different population segments of the same stock.We define two types of such ecological relationships: competition and predation.Competition is defined as the negative effects on one organism when another organism consumes (or interferers with access to) a needed resource that is in limited supply [19].This interaction leads to the reduced growth, survival, and/or reproduction of competitors; especially those that are phylogenetically subordinate [20] [21].Predation is the consumption of one organism (the prey) by another organism (the predator); the latter may be herbivorous, carnivorous, or omnivorous [21].
Fisheries that host a greater number of targeted species face a complicated management problem because more biological parameters must be incorporated in their estimates [12] [22] [23] [24] [25].Practically speaking, sustainable performance is hard to achieve at the ecological level among all the species in a multispecies fishery [26].Furthermore, the fact that the management of a multispecies fishery must also incorporate economic variables [11]  Natural Resources [28], only increases the degree of complexity.Yet it is important to continue developing this approach and urging decision makers to consider it, since catching a greater or lesser number of tonnes of a given species can have serious repercussions for the biomass of other species with which there is a trophic interaction (altering the tonnage would also affect, at least some, within-fishery species that were not "targets").Such catch differences can ultimately affect the sustainability of the ecosystem and of the enterprises that depend on it.
In this paper, we examine the case of a multispecies fishery with two target stocks that exhibit significant trophic interaction (ecological interdependence) as well as predation: the stock of blue whiting (Micromesistius poutassou) and the stock of European hake (Merluccius merluccius).From [11], we introduce costs dependent on the level of biomass into the net benefits function and estimate the optimum levels for biomass and catches of both species.In Section 2 we describe this fishery, where hake is the predator and blue whiting is the prey.In Section 3 we analyse a management model.It is applied to this mixed fishery in Section 4, and primary results are summarized in Section 5.

Hake and Blue Whiting Fisheries
In European waters, hake (Merluccius merluccius) and blue whiting (Micromesistius poutassou) are common commercial species.Both are primarily distri- buted across the continental shelf, where they spawn during the winter months.As we will see, these two species have a clear predator-prey relationship, as blue whiting is the hake's primary prey.
European hake is widely distributed throughout the North-eastern Atlantic Ocean.It is a demersal benthopelagic species that can reach depths that vary from 30 to 1000 meters, although they generally dwell at depths ranging between 70 to 370 meters.They form schools that stick close to the coastline in the summer and keep a greater distance during the winter.Egg-laying season lasts from January to May in the Bay of Biscay, and from May to July in the Celtic Sea.The European hake, a predator at the top of the Northeast Atlantic demersal trophic pyramid, primarily feeds on other species of fish like anchovies (Engraulisencrasicholus), sardines (Sardina pilchardus), blue whiting (Micromesistiuspoutassou), horse mackerel (Trachurustrachurus) and mackerel (Scomber scombrus).
Hake is caught in mixed fisheries-along with blue whiting, megrim, monkfish, and Norway lobsters-by a multi-rigged fleet, using the following fishing methods: bottom trawling with doors, bottom pair trawling, bottom-set longline, and small fixed gillnets.Figure 1 shows that this species has recovered over the last several years after nearly reaching biosafety level limits in the 1990s and after European Union (EU) implemented management plans for this species that were based on the progressive reduction of catches and fishing effort [29] [30].Total landings have been on an upward trend since 1998, in harmony with the abundance of biomass.
Blue whiting is found throughout the North Atlantic from the Southern Barents Sea and Eastern Norwegian Sea up to Cape Bojador, on the African coast.It is a demersal species of the gadus family.Its habitat is the ocean, and as a benthopelagic species, it lives mainly along the continental shelf and slope, distributed vertically at depths of between 150 and 1000 meters, though it is most commonly found at 300 -400 meters below the surface.After laying its eggs, it migrates north in the summer (Faroe Islands, eastern Iceland, and Norway) and returns to the area where it laid eggs in January and February.The eggs are pelagic and hatch between February and June.The main area where eggs are laid is to the west of the British Isles, and their growth is quite rapid.The blue whiting's diet-clearly seasonal-is composed mainly of crustaceans, namely copepods, krill, the larvae of decapods and a decapod known as white glass shrimp (Pasiphaeasivado).Bottom pair trawling is the primary method used to catch blue whiting (it is one of this method's target species) and individual bottom trawling, which picks up the by-catch.The blue whiting's biomass has been significantly reduced over the past several years (see Figure 2).This species plays a major role in fishing quota exchanges between European countries to obtain higher-value species.One such example is the exchanges made between Norway and the EU for higher cod quotas.
The ecological interdependence of the hake and blue whiting populations in European waters has been highlighted by various authors.Reference [31] as well as [32], study the diet of the European hake in the Cantabrian Sea by analysing its seasonal and bathymetric variations and according to the size, which is based on the stomach contents of the sample population.The blue whiting is noted as the hake's primary prey.The regression analysis of the predator-prey size relationship confirms a significant correlation between hake size and blue whiting size [9] [33].These authors also observed that hake is far more dependent upon blue whiting in the Cantabrian Sea than in the northern part of the Bay of Biscay [2] confirmed these results by showing that the blue whiting is the species upon which the hake most commonly preys in the Cantabrian Sea, though its importance as prey is moderate in the northern part of the Bay of Biscay and the Celtic Sea.Reference [10] also studied the diet of the blue whiting, the hake, horse mackerel and the mackerel in Portuguese waters by analysing the stomach contents of samples from these species taken from along the Portuguese coast.They concluded that the blue whiting is the hake's primary prey, followed by crustaceans (Processa spp.and mysidacea).In conclusion, in the Eastern Atlantic Ocean, the blue whiting is one of the hake's main prey, although this is somewhat less true in northern waters.

Predator-Prey Multispecies Model
The model employed is based on the logic equations used by [11], which, in turn, are based on Lotka-Volterra's predator-prey model [34] [35].By reordering it we get the following expressions: ( ) ( ) where X is the predator's biomass, Y is the prey's biomass; r m and r L are the intrinsic population growth rates for both species, respectively; X and Y the environmental charged capacities for both species; h m and h L the respective catches, and α and β the fish stock interaction coefficients.
Unlike [11], and given that trawl fleet generally depend on stock density [36], we can assume unit costs that depend on the corresponding biomass.Thus, the function of fishery's net benefits at moment t, is defined as: where P m and P L are the unit prices of hake and blue whiting, and C m (X) and C L (Y) the respective catch unit costs.The objective function selected for the control problem of fisher under centralised management (single owner) is: where ρ stands for the social discount rate.The optimal feasible controls, h m (t) = m h * and h L (t) = * L h , will be the ones that maximise the objective function while satisfying the conditions of the problem: ( ) do not yet know if the dynamics of both stocks will actually fit the quadratic form and in order to simplify calculations, we opted for the generic functions F(X) and G(Y) in problem ( 5) by incorporating the term of interaction between the two species.The Hamiltonian function H. corresponding to the management problem is based on the following expression:

H X t Y t h t h t t t t P C X h t P C Y h t t F X h t XY t G Y h t YX
where λ 1 and λ 2 represent the respective shadow prices of the predator and prey.
The first order conditions associated with (6) are based on: Following [37] [38] [39] and [40], and adapting it to the multispecies case, we obtain the following equilibrium equations: give us the steady state optimal biomass levels for both species, X u = X* and Y u = Y*.

Application to Hake and Blue Whiting Fisheries
In order to apply the previous multispecies model to the case study, we must first determine the functional form of the dynamics of each fish stock.To do so, econometric regressions were made using the ordinary least squares (OLS) method based on hake and blue whiting biomass and catch data for the 1988-2014 period, which are published by ICES [41] [42].In addition to the standard quadratic form, the following exponential and potential expressions, respectively, were considered for both species: e X XY X h The results of the econometric estimates are shown in Table 1.Overall, these expressions provide significant parameters for the levels that are normally accepted in statistical inference hypothesis testing (see Table 1).To compare these expressions, we used the adjusted coefficient of determination (R 2 adjusted), which covers the econometric model's explicative capacity (the maximum value of all of them), and the Akaike criteria (AIC).
The results indicate that the quadratic form is most appropriate for the growth functions of hake and blue whiting, with a value of R 2 adjusted from 0.992661 and 0.984233, respectively (where this value is below 0.87 in all the other cases), and higher AIC.In addition, only the quadratic forms have a positive predator-prey interaction coefficient for hake while at the same time negative for blue whiting, which would be typical of a trophic relationship.Therefore the following quadratic expressions were selected to create growth functions for hake and blue whiting: 1.49868 0.0000000821632 0.00000312827 As for the cost functions, we obtained average catch cost data for both species for the entire European fleet involved in this fishery, based on Scientific Technical and Economic Committee for Fisheries (STECF) annual economic reports on European fleets during 2011-2014.The expressions that best fit the data are as follows: By including the parameters of functions ( 15)- (18) in expressions ( 10) and (11) of the previous section and performing the calculations, we get the following expressions: And levels corresponding to the maximum sustainable yield (MSY) are given by the following expressions: to fish prices, by changing the price of a single species or for both stocks, and under a discount rate of 5%, are shown in Table 3.In general, the greatest increase in net benefits is obtained when the prices of both species are raised simultaneously.

Conclusions
A better understanding of trophic relationships among fish stocks allows to develop multispecies fishery management models, which is the category under which most European fisheries fall.This paper studies a fishery with two of the European fleet's target species that display significant ecological predatory interdependence.The trophic interaction between the predator (hake) and prey (blue whiting) is added to the management problem to jointly determine fishing quotas for both species.This interaction may be considered as a biological externality.By adding that externality to the fishery management problem, the  results show that both species are being over-fished, although hake is more so than blue whiting.This could place the ecological equilibrium that exists between the two species at risk, as well as alter the ecosystem overall.Furthermore, as a consequence of ecological interdependence between the two species, the MSY level of the predator (hake) is higher in this model than the estimated level by ICES (about 200 thousand tonnes; [41]), while the MSY for the prey is slightly lower (approx.3800 thousand tonnes; [42]).The TACs set by the EU in 2015 were approximately 105 and 1260 thousand tonnes.The results obtained from the predator-prey model show that catches should fall below the TAC levels set for each species.It is important to keep these interactions in mind when establishing the catch levels, given that the catch level for one species will have repercussions on the growth of the other and, therefore, on the sustainability of the marine environment and fishing activity in the long run.

Table 1 .
Results from econometric estimations for stock dynamic.
DOI: 10.4236/nr.2017.88036576 Natural Resources As for prices, we did not obtain any data on the European fleet.Based on the annual price data for first sale in the Spanish market for both species during the 2001-2014 period, we have estimated the average catch unit price in constant 2014 monetary units (P m = 4619.23€/tonne;PL = 399.13€/tonne).Given that these prices may differ from those obtained from other European countries, various price scenarios have been proposed to obtain the equilibrium values.The results for biomass, catches and net benefits using different discount rate values are shown in Table2.The optimal level of biomass of the predator (hake) is below the MSY level, while the optimal biomass of prey (blue whiting) is exceeding the MSY biomass.Likewise, as one can see, as the discount rate increases, the biomass level of both species decreases, while catches increase.The net benefits derived from fishing both species would increase by increasing the discount rate, although the net benefits from blue whiting fishing increase in greater proportion.The results of the sensitivity analysis with regard

Table 2 .
Optimum levels under different discount rates.