Spatial management of the European hake Merluccius merluccius fishery in the Catalan Mediterranean : Simulation of management alternatives with the InVEST model

European hake (Merluccius merluccius) is an important commercial fisheries species that shows growth overfishing, with catches basically focused on juveniles. This study assesses the benefit of closing a coastal area (an essential habitat for European hake recruits) to fishing, in addition to other alternatives of spatial management, compared with traditional, non–spatial management scenarios on fishing grounds exploited by the bottom trawl fleets of Blanes and Palamós (Province of Girona, NE Spain). We use InVEST, a spatially explicit model of intermediate complexity that simulates the bioeconomic effects of management measures for decision making. The sensitivity analysis of the model results shows the high influence of some parameters, particularly the parameterization of the recruitment submodel and European hake’s fecundity coefficients. The results are also examined in the light of uncertainty on migration parameters: in the two cases analysed (considering migration patterns or not), the results of the indicators (catch and revenues, abundance, recruitment and spawning stock biomass) were qualitatively similar and all show that the application of a restricted fishing area in one particular fishing ground (Vol de Terra) is the best management alternative. Its bioeconomic effects are comparable to a reduction of fishing effort of up to 20%. With high levels of ontogenetic migration, fishing on a second fishing ground (Cul de Rec – El Pas) should be restricted to enhance the biomass of the European hake population.


INTRODUCTION
European hake [Merluccius merluccius (Linnaeus, 1758)] is one of the most important demersal species of the Mediterranean demersal ecosystem.This species of great commercial value shows growth overfishing, with catches basically focused on juveniles (Recasens et al. 1998, Lleonart et al. 2003).In the Mediterranean Sea, European hake catches have decreased from 52394 t in 1994 to 20370 t in 2014.In GSA 6, this species requires immediate attention as its stocks show a clear pattern of declining recruitment and a high exploitation rate focused on recruits (age group 0) and age group 1 (STECF 2015).
The protection of European hake spawning areas has been proposed as an effective measure for improving the size composition of catches (Caddy 1999).The implementation of new spatial management measures such as marine protected areas (MPAs) could be the best alternative for the recovery of the fishery and for aligning the exploitation of Mediterranean fisheries resources with the Common Fisheries Policy objectives in terms of fishing mortality levels compatible with the maximum sustainable yield (MSY) in 2020 (Goñi et al. 2004, Maynou 2014).Several studies of European hake populations in the Mediterranean show that the species recruits in stable nursery areas that are recurrent from year to year (Colloca et al. 2015, Fiorentino et al. 2003, Tserpes et al. 2008), and restricting fishing effort in these areas would make an important contribution to reducing fishing mortality for age 0 European hake, thus contributing to stock rebuilding.
Several quantitative models have been developed to support decision making in fisheries management.According to their level of complexity, the fisheries literature has investigated "simple" models such as VIT (Lleonart and Salat 1992) that reflect the traditional vision of fisheries management by considering a single species with fleets that interact technically only, and "complex" models such as ECOPATH (Polovina andOw 1983, Christensen andPauly 1992), ECOSIM (Walters et al. 1997) and ECOSPACE (Walters et al. 1999) that try to represent all trophic levels in an ecosystem in a balanced manner.Finally, the literature offers models of intermediate complexity (Plagányi et al. 2014) or minimally realistic models that represent only a subset of (fisheries-relevant) ecosystem processes, with a limited number of species most likely to have important interactions with one species or a group of species of interest.These include the InVEST model (Integrated Valuation of Ecosystem Services and Tradeoffs, Sharp et al. 2016) chosen in this study, and which includes economic aspects in addition to the biological ones.
The InVEST suite of models includes software tools used to map and value natural goods and services that sustain human economies (Sharp et al. 2016).The InVEST toolbox contains ecosystem models designed for freshwater, marine and coastal ecosystems, among which the InVEST fisheries model (Arkema et al. 2015) was selected for our analyses.InVEST focuses on the quantification of trade-offs associated with alternative management strategies (in the context of fisheries, see Smith 1994).
The main objective of this study is to evaluate the biological and economic benefits of closing a shallow area (an essential habitat of European hake recruits) to fishing and the consequences of various spatial management scenarios applied to the European hake fishery in the Catalan sea.The study emphasizes the importance of simulating the bioeconomic effects of such spatial management measures for decision making through the use of a management strategy evaluation (MSE) approach.MSE has been widely recognized as a valuable tool for testing the robustness of management procedures to uncertainties in the fishing system (Kell and Fromentin 2007).MSE can be applied by using simulations to compare the effectiveness of different combinations of management strategies for achieving the expected management objectives and to compare the strength of these strategies with regard to errors and uncertainty (Smith 1994).Here, we assess the trade-off of management strategies based on spatial restrictions compared with effort limitation strategies.

Investigated fisheries
Our study area is located in the Catalan sea (GSA 6), including the area fished by the bottom trawl fleet from the ports of Blanes and Palamós (province of Girona, NE Spain, Fig. 1).In 2015, these two ports registered 27% of the total catches of European hake landed in Catalonia (1500 t) and are the most important ports in the autonomous community in terms of the operational fishing fleet: 40% of the operational fishing fleet in Girona (258 vessels) and more than 14% of the operational fleet in Catalonia (727 vessels).
In the study area, European hake is mainly exploited by bottom trawling, with 45 vessels (79% of the local fleet) with a size between 9 and 28 m length overall.It is also exploited, but to a lesser extent, by bottom longline, with 12 operational vessels (21% of the local fleet) and a ship size of between 6 and 11 m length overall.The number of fishing vessels has shown a decreasing trend in the last ten years (STECF 2015), as has the trawling effort (number of boats * days at sea), with a marked decline from 2009 to 2010 and from 2014 to 2015.The same trends have been reported for the Catalan fleet since 2005.The reasons for this decline are not well understood (Irazola et al. 1996), but probably originate from a combination of political measures aimed at limiting fishing effort, the removal of less efficient fishing units, and the general decrease in fisheries profitability in the Mediterranean owing to declining fish stocks.

The InVEST model
We used the InVEST model (Integrated Valuation of Ecosystem Services and Tradeoffs, https:// www.naturalcapitalproject.org/invest/), developed by the Natural Capital Project (Sharp et al. 2016), which can be used to quantify and map the values of ecosystem services to translate the understanding of the values of nature into management decisions.The software includes 17 models for analysing marine and terrestrial environments, including fish production, aquaculture, habitat quality, coastal protection, water quality and recreational fishing.The InVEST fisheries production model produces estimates of the volume and economic value of catches of the monospecific fisheries.It is a population model structured by age or stage, and it is presented as a generic model that can be adapted to most species and geographies.However, InVEST is not intended to provide an accurate prediction of these estimates, but rather is used as a tool to explore the consequences of decisions that could affect fish production.The particularity of the InVEST model is that it includes both biological and economic aspects.Furthermore, it is spatially explicit, using maps as a source of information and producing maps as outputs.Natural Capital provides software in Python for the execution of simulations.We used InVEST version 3.3.32and worked with the Windows 10 operating system.

Inputs
InVEST uses life-history information and survival parameters to estimate the volume and the value of harvests.In the case of a population structured by age, we estimated the parameters of the population that characterize the life history of the species (recruitment, ontogenic migrations, natural mortality, maturity, weight and fertility), the subregions of interest as defined by the fishing grounds within the study area, the attributes of each subregion (fraction of exploitation and larval dispersal), the behaviour of the fishery (vulnerability to fishing), the habitat dependencies (importance and availability of nursery habitat) and optionally, economic valuation (price per unit catch).

Outputs
The model is executed according to a number of time steps specified by the user and sufficient for the population to reach a state of equilibrium.The results of the model are estimates of the volume of the catch (in t) and economic value of the catches (in €) in the last step of execution and within the subregions established, in addition to estimates of abundance by age class and subregion in all years of simulation.The results of multiple executions of the model can be compared, each one representing different scenarios of habitat extension, environmental conditions and/ or fishing pressure.The results can be added from the folders to the ArcGIS document.

Sources of InVEST inputs
For the present study, the definition of the inputs related to the biological data, valuation data, commercialization of the catches and data on fishing activity were obtained from literature sources, fisheries production data of the Catalan Fisheries Directorate database and spatial distribution of effort from the Spanish Ministry of Agriculture, Marine Affairs and Environment by analysing vessel monitoring data (VMS) with the QGIS software.The VMS data provided the position of all the fishing boats of Blanes and Palamós for the years 2013-2015 at an average frequency of 60 minutes.Finally, data on the migratory behaviour (ontogenic changes in spatial distribution) of European hake were obtained by associating the existing bibliographic data (Doumenge 1966, Recasens et al. 1998) with the European hake distribution maps by age extrapolated in QGIS 2.10.

Areas of interest
The VMS data were used to establish the extent of fishing activity in the area of interest and to determine the fishing effort at the level of each fishing zone (de- fined in Fig. 2, with local names; see also Supplementary material Table S1).Fishing effort was calculated as number of hours of trawler activity in each cell of a 1×1 km grid covering the study area (Fig. 3), by filtering fishing vessel signals reporting speeds of between 2 and 4 knots (i.e.fishing tracks).Duplicated positions were deleted.The daily cleaned VMS data were then aggregated by year and averaged.All operations were carried out with QGIS 2.10.The fishing effort map (Fig. 3) indicates a heterogeneous distribution of fishing effort (average for 2013-2015) throughout the study area, which shows that the most exploited fishing grounds are particularly Gamba -St.Sebastià, Sot de la Malica, Barana, Vol de terra, Través, Abissínia, Cul de Rec -El Pas and L'Avió.

The European hake population model
We modelled the European hake population as 28 subpopulations interconnected via ontogenetic migration, simulating European hake movement from nursery shallow waters to greater depths with age (ages 2, 3 and 4).The age-structured population was projected from 2015 to 2025 in annual time steps.Hence, the parameters are based on annual rates and the model progresses in one-year increments.
The European hake population dynamics is given by: ∑ Harvest is assumed to occur at the beginning of the year, before the mortality from natural causes.Harvest (H x,t ) from each subregion in the final (equilibrium) time step is calculated on the basis of user-defined exploitation rate(s) and vulnerability.It is given by: where Ex x is exploitation, which is the proportion of the population vulnerable to harvest that is actually harvested, V a is vulnerability to harvest for age a, and w a is weight for age.
The number of spawners is the product of the number of individuals in each age class for the entire study region and the proportion that are mature by age: The biomass of spawners is the product of the number of individuals in each age class for the entire study region, the proportion that are mature at each age and their weight at a given age: where W a is weight by age.
We consider the hypothesis of fast growth for European hake, with a maximum age of approximately 15 years (García-Rodríguez andEsteban 2002, Mellon-Duval et al. 2010), but demographic data for the study area revealed only five age classes.Thus, five age classes of the species were considered, with a group 5+ that represents the maximum age class and includes all individuals of this age and individuals of more than 5 years.We pooled both sexes in the data, as fishing practices or sale of landings are independent of sex.Maturity at age were derived from the Data Collection Framework, in STECF (2015) The biological parameters used as inputs are presented in Supplementari material Tables S2-S5.
Survival from natural mortality is the proportion of individuals that continue to the next age.It is calculated from the instantaneous natural mortality rates (M a ) as where t is the duration of the temporary step over which the survival is calculated (one year).
Estimating that the largest age of European hake is 15 years, the vector of natural mortality by age was calculated from the formula of Caddy, using the PRO-BIOM Excel spreadsheet (Abella et al. 1997), and considering that natural mortality decreases with the older ages (Supplementary material Table S6).
The exploitation fraction (Ex) in InVEST refers to the proportion of the vulnerable and effectively captured population in each subregion with the effort applied (0 = 0% captured, 1 = 100% captured).Each subregion is treated independently.In this study, the calculation of this proportion was based on the fishing effort data available in the study area through the VMS data.By means of the spatial union of the two layers, "subregions" (Fig. 2) and "effort" (Fig. 3), with the QGIS software, the average effort in each subregion was calculated.The maximum exploitation fraction = 1 (100%) was attributed to the subregion Gamba -St.Sebastià which recorded the highest fishing effort.We calculated the exploitation fraction of the other subregions based on this fishing ground (Supplementary material Table S7).
Vulnerability to fishing may depend on size, specific behaviour of the species during a life stage, habitat use or regulations, and may change depending on fishing gear and fishing strategies.A value of 1.0 indicates that age is totally vulnerable to fishing, while values below 1 indicate vulnerability relative to fully vulnerable age.It is assumed that the vulnerability is the same in all subregions.On this basis, vulnerability is calculated by observing the values of fishing mortality applied to each age (Supplementary material Table S8).According to the stock assessment of 2015, age 2 represent the age totally vulnerable to fishing (STECF 2015).

Recruitment submodel
In the case of European hake, as for other Mediterranean fishery stocks, it is difficult to define the relationship between spawning stock biomass (SSB) and recruitment with any degree of certainty.In addition, the stock assessment figures in GSA 6 show considerable fluctuations in European hake recruitment in the last two decades, which may result from a predominance of environmental factors in the recruitment process, rather than a direct relationship between SSB and recruitment.We therefore examined the simulation results under two major hypotheses about the recruitment function assumed in the InVEST model.We first considered the hypothesis of a fixed recruitment function (H1), and then analysed the results assuming a linear recruitment (fecundity function: H2).The fecundity function assumes a constant reproductive rate for adults while the fixed recruitment function assumes that recruitment is constant and not dependent on the number of adults.The parameterization of each function is given by Fixed recruitment function: To determine the initial number of recruits, a virtual cohort analysis was carried out using VIT (Version 1.3).We incorporated the data of length frequency, maturity of the different age classes, the European hake catch for 2015 and the other required population parameters proportionally scaled to the study area (Supplementary material Table S9).
In the case of the linear recruitment function (fecundity), InVEST requires the age-specific fecundity values representing the number of descendants (offspring) per mature individual.These fecundity coefficients refer to the number of European hake eggs and larvae that survive and become recruits.Knowing the number of eggs per female of mature European hake (Recasens et al. 2008), we need to know the daily mortality rates of the eggs of the species to define the number of survivors corresponding to the recruits.
First, based on the observations of European hake individuals in MEDITS (2012MEDITS ( -2013) ) and in the commercial fishing data, we calculated the number of eggs per individual through the equation of Recasens et al. (2008): where y is fecundity (number of eggs, in thousands) and x is size in cm.
Length to age conversion (slicing) allowed us to estimate the average number of eggs for each age class, considering the MEDITS observations for ages 0, 1 and 2 and the commercial fishing observations for the higher ages (Supplementary material Fig. S1).Given the number of eggs (N 0 ), we applied an exponential model of decreasing abundance (average density) with respect to age to define the average number of recruits per European hake age class: where N a is abundance at age a, N 0 is estimated abundance at hatching, Z is instantaneous mortality coefficient (daily) and a is age in days.
The data presented in the literature indicate that mortality rates of eggs and larvae of demersal fish with pelagic eggs in general must be high and variable both between years and between populations.The range of mortality rates observed is 7% to 67% per day.Daan (1981) found a daily mortality rate of 22% for cod eggs in the North Sea.To be able to define this daily mortality rate of European hake eggs, we parameterized the model according to the catches observed in 2015 (of the order of 70 t), and considering an interval of 60 to 120 days of egg development.The best combination provided a daily mortality of 12% (0.1193) during a period of 90 days.Based on this, the result of the fecundity coefficients by age is presented in Supplementary material Table S5.
The combination of the sampling data of MED-ITS cruises with the commercial fishing series of the longline fleet (which exploits the older age classes) allowed us to draw a distribution of the different age classes and to exploit the available observations to make an interpolation of the spatial distribution of European hake throughout the study area.The interpolation of the available data provided the spatial distribution of each age class by depth ranges (Supplementary material Fig. S1).The interpolation method used is inverse distance weighting in QGIS.
As an input required by InVEST, the "larval dispersal" represents the proportion of the cumulative larvae pool that disperses into each subregion.For models with subregions (as defined by the fishing grounds, Fig. 2), the model assumes that adults from each subregion contribute to a common larval pool.Larvae are then distributed across subregions.To model the distribution of the larvae among the 28 subregions of interest, considering the lack of large-scale studies on the dispersal and behaviour of European hake larvae, we assumed the same spatial distribution of larvae and age 0 individuals throughout the study area.We used the QGIS software to interpolate the densities of age 0 over the entire study area [MEDITS cruises (2012[MEDITS cruises ( -2013) ) and commercial fishing data].Given that the dispersal across all subregions should add up to 1, a total density value is extracted in each subregion and then expressed as a relative density to the total density of age 0 individuals throughout the subregions of interest (Supplementary material Fig. S2).
The migration submodel of InVEST allows European hake movements between the subregions to be integrated and was applied here with high uncertainty because no study on European hake ontogenetic migration had been conducted in our study area.To limit the uncertainty of this submodel, the results were analysed in two parts, first without considering migration and second taking into account the ontogenetic migration of European hake, and the sensitivity of the outputs to this submodel are discussed.
To map the migration of European hake in this model, we consider the migration diagram proposed by Doumenge (1966), the existing literature and the European hake density maps by age resulting from the interpolation of the MEDITS abundance and commercial fishing data in the QGIS software.We conclude that migration occurs for three main ages: age 2, age 3 and age 4, with an increasing distribution towards deeper fishing grounds with age (Ontogenic migration).Combining all these data, we defined the probable direction of the migration between fishing grounds, and consequently the migration coefficients for the three ages (Supplementary material Table S10).
Valuation, V t,x reflects the earnings from the sale of the harvest.It is intended to give a rough idea of the current market value for an equilibrated population based on user-defined price parameters.It is simply where Price is the value in price per units (where units match those given by H x ) and FractionProcessed is the proportion of each harvest unit that remains to be sold after processing.To define the unit price of the catch, we referred to the European hake fishing statistics in the ports of Palamós and Blanes (Catalonia Fisheries Statistics, 2000-2015), giving an average price of 8.16 euros kg -1 .On the other hand, the "Fraction Kept After Processing" is 1, as the species is sold as a fresh product without processing.

Sensitivity analysis of the InVEST model
The various parameters used in this study were implemented with a certain level of uncertainty, so we included the effect of the uncertainty in the results for a better interpretation.The sensitivity of the results to the submodels and the most uncertain parameters in this study are analysed, namely: the European hake growth parameters (±10%) and by default natural mortality, daily mortality rate (Z) in the calculation of fecundity coefficients (±10%) and, finally, the submodels of recruitment and migration.For this analysis, we proposed new input matrices each time, executed the program again and collected the results of the new simulations.The results of sensitivity analysis are presented in the form of principal component analysis (PCA).

Simulation scenarios
The InVEST fisheries model facilitates comparison of fisheries production under different scenarios.In our study, we used InVEST to analyse the bioeconomic effect of a spatial management measure (MPA) and various cases of effort limitation in different fishing grounds of the study area, and then to compare the short-, medium-and long-term results of the management scenarios.The selected scenarios are consistent with the requirements of Council Regulation (EC) 1967/2006, regarding the inclusion of spatial aspects such as the establishment of protected fishing zones in order to protect nurseries and/or spawning areas for the sustainable exploitation of fishery resources in the Mediterranean Sea, and with the recommendations of the General Fisheries Commission for the Mediterranean that encourage the reduction of fishing effort in Mediterranean fisheries.

Scenario 1: Spatial closures
The first scenario consisted of the introduction of spatial closures in areas representing the essential habitats of European hake recruits in order to limit the fishing pressure at their level.Among the fishing areas with a high abundance of recruits, we selected the shallow fishing grounds located closer to the edge of the continental shelf (average depth ≤200 m) as the most likely to coincide with the habitat of European hake recruits.Then, we studied the effect of closing each of these six fishing zones: Vol de Terra (75 m), L'Avió (99 m), Planassa (111 m), Vol de fora (130 m), Bravada Dbf (141 m) and Cul de Rec -El Pas (182 m).The main scenario evaluated was the closure of the fishing ground Vol de Terra, which was proposed by the fishing sector during the participatory sessions organized by the MINOUW project (http://www.minouw-project.eu) in the period 2015-2016, as a step towards joint co-management of the demersal fishery in the area.The fishing sector and the fisheries administration perceived that closing this fishing ground to fishing would help protect European hake and other fisheries resources (e.g.red mullet) and sensitive habitats (Posidonia beds are located in the shallowest reaches of this fishing ground).

Scenario 2: Reduction of fishing effort
The second scenario consisted of the application of reductions in the fishing effort for the trawl fleet in our study area, taking into account the distribution map of the effort in the established subregions.

Sensitivity of the results to InVEST submodels
The PCA carried out to analyse the effect of the two InVEST submodels that we parameterized with most uncertainty, namely, the selected recruitment function (fixed or fecundity) and the migration matrix is shown in Figure 4A.The first two factorial axes of the PCA accounted for 99.51% of the total variability.It is noted that all the outputs of the model showed a positive correlation with the first component, which differentiates between the simulations with positive impact on the catch: the SSB and the abundance.The second component showed a positive correlation with the capture and the SSB in the short term and a negative correlation with the other outputs, but with a small contribution (1.18%).We conclude that the integration of the migration submodel when a constant recruitment is assumed has a limited effect on the results.The results were highly sensitive to the recruitment submodel when a linear SSB-R relation was assumed for European hake.In this case, the migration submodel also had an important effect on the results through a considerable increase in the total catch, abundance and SSB owing to the high sensitivity of the productivity to any change in the fecundity of the species.

Sensitivity to population parameters
In the PCA (Fig. 4B), the factorial axes F1 and F2 accounted for 99.54% of the total variation.The first axis had the highest contribution (+82%) and sets the parameters that positively affected all the outputs against those that had an opposite effect on all the output variables.The second axis contributed 17% to the total variation and sets the parameters with a positive correlation with the long-term capture and the abundance at the end of the simulation against the parameters that negatively affected the indicators in the short term.
This analysis shows that the parameters that had the greatest impact on the results were the fertility coefficients.In addition, it is noted that when one is working with a linear recruitment, varying the daily mortality rate (Z) or the growth parameters (G) of the species directly affects the fertility coefficients and generates quite different results.A low daily mortality (F+) (high fertility) or high growth parameters (G+F) (low natural mortality) lead to an increase in the assumed fertility and result in an increase in the total catch, the SSB and the abundance, especially in the long term.On the other hand, a slow growth (G-) or a higher daily mortality rate of European hake (F-) lead to a reduction of the productivity variables.
In the case of constant recruitment, the growth parameters affected the results positively when fast growth was assumed, and negatively when slow growth was assumed.However, this effect was limited compared with the sensitivity shown by the parameters under the issue of linear recruitment.

Impact of the management measures
The results of the present study were analysed in two parts (with or without the migration aspect) and under two main hypotheses (constant and linear recruitment).This allowed us to consider the effect of the two submodels of InVEST parameterized with more uncertainty in the representability of the results.Firstly, we explored the effect of the closure of fishing grounds (Supplementary material Tables S7, S8) and selected those with the most important effect on the results.Then we compared the effect of the MPAs selected (A1, A2, A3, A6) with the effect of effort limitation (-10% and -20%).
In the cases analysed, we conclude that it would be beneficial to close the Vol de terra fishing ground (A1) among the other MPA scenarios (Fig. 4C, E).The application of an MPA in Vol de terra seemed to have the best effect on the biological indicators representing abundance and SSB (Fig. 5A).In addition, the analysis of the results with linear recruitment showed the benefit of this MPA in the medium and long term and for all indicators, including total catch, which recovered over time (Fig. 5B).The recruitment function showed higher values for all the indicators because the time evolution of each indicator showed an increasing trend.Thus, all the results were presented relatively to the baseline scenario.
The beneficial effect of the Vol de terra MPA on the indicators was also maintained, when the ontogenetic migration of the European hake for ages 2, 3 and 4 was considered.In this case, we observed another fishing ground (Cul de Rec -El Pas) whose closure may also be beneficial for the fishery, with similar results to those of the Vol de terra MPA (Figs 6 and 7).In addition, considering European hake migration involved a change in the distribution of the catch and abundance in the different subregions, with an increase in catch and/ or abundance in some fishing grounds, and a decrease in others in response to the direction of displacement with migration (according to depth) (Table 1).
The closure of Vol de terra is approximately equal to a reduction in fishing effort of up to 20% in the entire area of interest (Fig. 4D, F), which would be equivalent to the interruption of fishing activity for one day per  week.When considering the migration of European hake, the results also revealed that the effect of this MPA is quite close to the effect of reducing fishing effort by 20% (Figs 6 and 7).Finally, InVEST produced distribution maps of catch and income in each fishing ground and for each management scenario.We present specifically the distribution of catch in a response to the closure of Vol de Terra as the incomes show the same distribution as catches (Fig. 8).

DISCUSSION
This work represents the first application of the In-VEST model to a specific case of demersal fishery in the Spanish Mediterranean.It is also the first attempt to pre-evaluate the effect of spatial management on an overexploited fishery of great interest for the fishing sector in Catalonia.
Throughout this study, we have evaluated the effect of the application of an MPA at the level of the Vol de Terra fishing ground, which seems to be an essential habitat for European hake recruits, together with other management scenarios in order to compare their efficiency and determine their effect on different biological indicators (abundance, SSB, recruitment) and economic indicators (catch or amount) of the fishery.

The InVEST model and the quality of inputs
The InVEST model has been treated with a level of uncertainty regarding the quality and availability of the information introduced (inputs).InVEST requires specific data on the biology and population dynamics of the species of interest and data on the fishing activity in the study area.This model allows us to integrate available information about the fishery and use it to explore various management scenarios (Sharp et al. 2016).A particularity of this model is that it needs a limited data volume in comparison with other spatially explicit models: ISIS-Fish (Ifremer 2007) and DIS-PLACE (Bastardie et al. 2014).
Some submodels and parameters were implemented with more uncertainty than others.Firstly, the relation between SSB and recruitment for European hake remains undetermined throughout GSA6 given the limitation of the historical series of these data.We therefore studied two hypotheses for the recruitment function (constant and linear), analysing the impact of each function on the results, which showed qualitatively similar results.Second, the fecundity coefficients  proposed in the case of linear recruitment are a source of considerable uncertainty and have a high influence on the results, as there are no scientific studies available on the survival and/or daily mortality of European hake eggs and larvae Furthermore, in relation to the migration submodel, the migration coefficients were considered with uncertainty due to the absence of mark-recapture studies or other specific studies that explore the migratory behaviour of European hake in the Catalan Mediterranean Sea.The integration of a migration matrix in this study has a more perceptible effect on the results, especially when we consider the hypothesis of linear recruitment, owing to the added uncertainty of the migration coefficients to that of the fertility coefficients.
Subsequently, the parameters of growth and, consequently, natural mortality (resulting from the change of these parameters) affect the results to a lesser extent than the previous parameters and do not change the general conclusions of the study.The uncertainty of the results for these parameters was considered in the representation of the results with confidence intervals, but not in the case of the linear recruitment function, given the high sensitivity shown.We note that European hake tagging studies (from de Pontual et al. 2003, 2012, Mellon-Duval et al. 2010) will be useful to improve knowledge about migration, natural mortality and growth of this species.
All the other parameters considered in this study were obtained from the most recent bibliographic sources and/or from the scientific data available to date (MEDITS), so their effect on the results is estimated as very limited.However, additional research to clarify the dynamics of European hake recruitment through the stations at the level of the study area will improve the information on the density distribution of the recruits and thus delimit the zones of nurseries of the species for better management.

InVEST results
On top of the problem of over-exploitation, the risks of collapse or depletion of the Mediterranean stocks seem limited (Caddy 1999) due to the high resilience of these stocks compared with Atlantic ones (Lleonart 2005).This is explained by the high exploitation rates of juveniles against the low mortality of adults with the occurrence of refuge zones for adults on the high seas.This high mortality of juveniles significantly affects the productivity and structure of the population (Hidalgo et al. 2009).The average capture size of these stocks is smaller than the size that would produce the maximum yield and suffers high variability due to instability of population dynamics and a decrease in the stocks' capacity to cope with environmental fluctuations (Anderson 2008).
Given the multi-specific peculiarity of Mediterranean fisheries, a combination of management tools based only on technical measures and effort control has not managed to guarantee the long-term sustainability of fisheries or the conservation of important habitats (Colloca et al. 2013, Tudela 2004).Therefore, the implementation of management measures to protect aggregation areas of juveniles during their first year of life has the potential to substantially improve the present fishing exploitation patterns and is a major prerequisite for the sustainability of trawl fisheries in the future (Colloca et al. 2015).
In addition to the results of our study, many studies have shown that the spatial closure of Mediterranean nurseries can provide important benefits to fisheries in terms of increased fishing resistance and improved yields (Apostolaki et al. 2002, Scarcella et al. 2014).It would also allow a better fulfilment of the landing obligation established by the new European Fisheries Policy and a reduction in the high fishing mortality of juveniles of demersal stocks in the Mediterranean (Colloca et al. 2013, Tsagarakis et al. 2014).Spatial restriction of fishing need not be permanent in order to be effective.Cinner et al. (2006) show that periodic closures can have a positive effect on fisheries resources in terms of both catches and biomass conservation.Other management schemes are also possible, such as the periodic rotation of area closures, whose positive effects are demonstrated in Hart and Rago (2006) or "pulse fishing", i.e. restricting fishing on given fishing grounds on alternate years (Da Rocha et al. 2012).
Outside the Spanish Mediterranean, the European hake fishery in the Gulf of Lion (GSA7) (Ifremer 2011) also shows over-exploitation by a mixed fleet, composed mainly of French (88%) and Spanish (12%) trawlers fishing small European hake individuals.The pre-evaluation of this fishery through a bioeconomic model (Macher et al. 2010) showed that interrupting the fishing activity of the trawls for one month seems to lead to an improvement in the catch and the SSB, and a decrease in the global and individual performance of the trawling fleet, without reaching the MSY.In addition, this study showed that a reduction of fishing effort by about 15% per year for four years allows the MSY to be reached with more beneficial results on catch and SSB.In this case, reducing the number of trawls (instead of fishing days) was the only scenario with a positive impact on the individual performance of all boats (trawls, longlines and others).
The analysis of another mixed European hake fishery in the Bay of Biscay using an ISIS-Fish simulation tool (Ifremer 2007) evaluated the effect of spatio-temporal management on the dynamics of this mixed fishery (Drouineau et al. 2006).The parameterization of the ISIS-Fish model requires the same inputs as InVEST, in addition to other additional parameters.The authors concluded that the application of MPAs in the European hake nursery areas improves not only the stock of European hake but also the total biomass, the catches, the SSB and the recruitment of the second species of interest for this fishery.Therefore, this study highlights the interest of MPAs implemented in areas of European hake nurseries and supports the conclusions of our work.In addition, it is suggested through this study that the application of an MPA in our area of interest may be useful for the other demersal species of this mixed fishery, a finding which could not be revealed with InVEST because of its monospecific peculiarity.
The understanding of spatial patterns in population dynamics is essential to protect critical habitats and ensure sustainable management of fisheries resources (Berkeley et al. 2004, Caddy 2000).In the specific case of European hake, such information on the essential habitats for the reproduction of the species in the Mediterranean will contribute to the good application of the spatial management of fisheries required by the EU and, in particular, to limiting fishing mortality of the recruits (Druon et al. 2015).
To identify the appropriate areas for closure against fishing, many authors have studied the spatial distribution of juvenile European hake at a regional scale and identified the main nursery areas where the highest concentrations of juveniles remain stable over the years (Colloca et al. 2015, Fiorentino et al. 2003, Tserpes et al. 2008).Other studies are based on the use of the ecological niche model for the modelling of potential European hake nurseries in the Mediterranean Sea (Druon et al. 2015).

CONCLUSION
The present study has shown that a spatially explicit tool such as InVEST can be very useful for preevaluating spatial management measures and choosing among the best possible alternatives for a fishery by integrating a limited number of non-complex data.The present work has allowed us to explore this model to guide decision making in the management of the European hake fishery, which is overexploited in the Mediterranean, and specifically at the ports of Blanes and Palamós.The results of the study coincide with other studies highlighting the benefits of spatial management and confirm the importance of protecting this fishery to sustain the stock of European hake, which is very sensitive to changes affecting juveniles.
The results of this work allow us to draw two main conclusions: First, the closure of the Vol de Terra fishing ground is the best spatial management alternative for the recovery of catches and biomass of European hake in the medium and long term, and would have a minor effect on fisheries performance in the medium term.Fishers are well aware of the importance of this area as a hake nursery, and our work helped confirm their preferences, as expressed during the participatory approach between scientists and stakeholders.According to our results, spatial closures in other fishing grounds would not attain a similar degree of stock protection.Second, the application of an MPA in this fishery would be equivalent to a reduction in fishing effort of 20% in the entire study area, but it would be easier to implement and would meet with less resistance from the sector than a fishing effort reduction measure.

Fig. 4 .
Fig. 4. -PCA results.A, sensitivity analysis of the indicators to the submodels Recruitment and Migration; B, sensitivity analysis of the indicators to the parameters (G+, fast growth, G-, slow growth, F+, low daily mortality, F-, high daily mortality); C, comparison between spatial management scenarios, constant recruitment; D, comparison between spatial and non-spatial management scenarios, constant recruitment; E, result of spatial management, linear recruitment; F, comparison between spatial and non-spatial management scenarios, linear recruitment.Legend: R, fixed recruitment; F, linear recruitment.Cc, short-term catch; Cm, medium-term catch; Cl, long-term catch; Nc and Nf, short-term and end-of-simulation abundance; Nl, long-term abundance; E0, baseline scenario; A1, MPA in Vol de Terra; A2, MPA in L'Avió; A3, MPA in Planassa; A4, MPA in Vol de fora; A5, MPA in Bravada Dbf; A6, MPA in Cul de Rec -El Pas; C1-10% and C1-20%, reduction of effort by 10% and 20% in the entire study area.

Fig. 8 .
Fig. 8. -Spatial distribution of catch under A1 (MPA in Vol de Terra) at the end of simulation (time t) according to two recruitment models.A, constant recruitment; B, linear recruitment.

Table 1 .
-Effect of MPA Vol de Terra on results, highlighting in bold face the main fishing grounds (Area) contributing with recruits to MPA Vol de Terra.

Table S7 .
-Exploitation fraction by subregion based on fishing effort.

Table S8 .
-Fishing mortality and vulnerability by age class.