SciELO - Scientific Electronic Library Online

 
vol.42 issue4Pedigree traceability in whiteleg shrimp ( Litopenaeus vannamei ) using genetic markers: A comparison between microsatellites and SNPsRelative abundance of Antillean manatees in the Pantanos de Centla Biosphere Reserve in the coastal plain of Tabasco, Mexico author indexsubject indexsearch form
Home Pagealphabetic serial listing  

Services on Demand

Journal

Article

Indicators

Related links

  • Have no similar articlesSimilars in SciELO

Share


Ciencias marinas

Print version ISSN 0185-3880

Cienc. mar vol.42 n.4 Ensenada Dec. 2016

https://doi.org/10.7773/cm.v42i4.2628 

Articles

Spatial characterization of the fisheries in the Avilés Canyon System (Cantabrian Sea, Spain)

Antonio Punzón*  1 

Juan Carlos Arronte1 

Francisco Sánchez1 

Ana García-Alegre1 

1Instituto Español de Oceanografía [Spanish Institute of Oceanography], Promontorio de San Martín SN, 39080 Santander, PO Box 240, Spain.


Abstract:

Marine protected areas are one of the main tools in ecosystem-based fishery management. This approach requires analyzing the role of fisheries in the ecosystem by, for example, identifying fishing tactics and analyzing the spatial distribution with sufficient precision to allow us to determine whether fishing activities and the habitat or ecosystem to be protected can coexist spatially. In the Avilés Canyon System (Cantabrian Sea, Spain), we detected 6 different gears and 17 fishing tactics were identified by a non-hierarchical classification technique (CLARA). The spatial distribution of demersal fishing activity is strongly affected by the topography of the study area and the large presence of rocky ground. There is a clear spatial segregation between gear types, mainly due to the operating requirements of the gears and the distribution of target species. This minimizes conflicts between different components of the fleet. The fishing tactics targeting pelagic species were mainly located on the coast and shelf, while those targeting demersal species were used on the shelf and slope. Longline fishing tactics are those that had a more restricted spatial distribution.

Key words: fisheries; vessel monitoring systems; spatial distribution of fishing effort; Avilés Canyon System; fishing tactics

Resumen:

Las áreas marinas protegidas son una de las principales herramientas en la gestión pesquera basada en el ecosistema. Este enfoque requiere evaluar el papel de las pesquerías en el ecosistema mediante, por ejemplo, la identificación de tácticas pesqueras y el análisis de su distribución espacial, y esto permitirá identificar la coexistencia espacial de las actividades pesqueras y los hábitats o ecosistemas objeto de protección. En el Sistema de Cañones de Avilés (mar Cantábrico, España), hemos identificado el uso de 6 artes de pesca y 17 tácticas pesqueras, éstas últimas con una técnica de clasificación no jerárquica (CLARA). La distribución espacial de las actividades pesqueras demersales está fuertemente afectada por la topografía del área de estudio y la gran presencia de sustrato rocoso. Hay una clara segregación espacial entre tipos de artes de pesca, principalmente debido a los requerimientos operativos de las artes de pesca y la distribución de las especies objetivo. Esto minimiza los conflictos entre los distintos componentes de flota. Las tácticas pesqueras dirigidas a especies pelágicas se localizaron principalmente en la zona litoral y en la plataforma continental, mientras que aquellas dirigidas a especies demersales se distribuyeron por la plataforma continental y en el talud continental. Las tácticas pesqueras de palangre de fondo son las que tuvieron una distribución espacial más restringida.

Palabras clave: pesquerías; sistema de seguimiento de buques; distribución espacial del esfuerzo; Sistema de Cañones de Avilés; tácticas pesqueras

Introduction

Fishing activity has important consequences for both the exploited communities and the ecosystems where these activities are conducted. Fishing activity is part of the ecosystem, and its effects are numerous and widely studied (Hollingworth 2000, Cury et al. 2005). Among the most important effects are habitat destruction, mortality of non-target species, changes in population demographics, and changes in the function and structure of ecosystems (Jennings and Kaiser 1998, Pikitch et al. 2004, Sánchez and Olaso 2004).

This focus has changed the foundations of fisheries management and administration towards what has been called ecosystem-based fishery management (Sinclair et al. 2002, Link 2005). According to Pikitch et al. (2004), the main objective under this paradigm is to preserve the health of marine ecosystems and the fisheries that depend on them. Accordingly, the concepts of ecosystem approach to fisheries (EAF) and marine protected areas (MPAs, including marine reserves) have been developed as governance and management measures, respectively (Jones 2001, Pauly et al. 2002, Claudet and Pelletier 2004, García and Cochrane 2005). For the implementation of an EAF, it is necessary to obtain a better understanding of how the ecosystem functions (Gislason et al. 2000, Sinclair et al. 2002, Fréon et al. 2005), taking into consideration fisheries as part of the ecosystem and not only as a fishing activity to regulate. In this sense, knowing the behavior of fishing activities, the target species, and the space where they are working are essential aspects for the development of future management plans that minimize the socioeconomic effects on fishing activity and maximize the protection of habitats and target species.

In this context, the countries of the European Union have to create a network of MPAs to ensure the conservation and protection of habitats, resources, and breeding areas, connectivity between areas, etc. Continuous and intensive fishing occurs in many of the MPA candidate sites and any implementation of regulatory measures could bring about large social conflict (Baelde 2005). Therefore, it is essential to gain an understanding of the fishing activity and its relationship with the ecosystem in order to spatially characterize the possible effects on the ecosystems.

The Avilés Canyon System (Fig. 1) is one of 10 areas chosen by the Ministry of Environment in Spain to be designated an MPA and integrated within the European Natura 2000 network. In the area of influence almost all types of fishing gear are used (trawl, longline, purse seine, and gill nets). This makes this area an excellent example of how a structurally and geographically complex zone (Gómez-Ballesteros et al. 2014) is exploited by mixed fisheries.

Figure 1 Study area (ACS, Avilés Canyon System). 

In the European Union, information on the spatial distribution of fishing activity can be obtained from 2 main sources: logbooks and vessel monitoring systems (VMS). The coupling of VMS and logbook data offers the possibility of performing a precise spatial analysis of fishing activity and its relationship with the affected ecosystems. The present study analyzes and characterizes for the first time the spatial distribution of Spanish fishing activity by fishing tactic (defined as a combination of gear, target species, and working area) in the southern part of the Bay of Biscay. This will allow us to analyze the fishing activity depending on the target species and relate each type of activity to the characteristics of the area. In addition, we can analyze the use of the area made by each gear and fishing tactic and the relationship with the other fishing activities. This information will be one of the main tools for the development of management plans for the Avilés Canyon MPA.

Materials and methods

Study area

The study area is located in the Cantabrian Sea (southern Bay of Biscay) and consequently falls within OSPAR Region IV. This whole area is characterized by a narrow platform andby deep-sea canyons. The Avilés Canyon System (Fig. 1) is a complex region of canyons and valleys comprising 3 main canyons of different morphostructural characteristics (Gómez-Ballesteros et al. 2014). The near continental shelf is generally narrow; the width varies from 12 km (where the head of the Avilés Canyon is incised) to 40 km. Maximum depths are approximately 200 m, but occasionally reach 300 m east of the canyon. The continental shelf displays rocky outcrops of different reliefs and morphologies, and a nearly horizontal surface corresponding to sedimentary areas with unconsolidated sediment cover or heavily eroded basement surfaces. The strong currents over the shelf are responsible for the minimal sedimentary cover in the area.

Data source

The bathymetric information presented here was obtained during two surveys carried out aboard the R/V Vizconde de Eza between 2010 and 2011 using a Simrad EM-300 multi-beam echosounder (Gómez-Ballesteros et al. 2014).

The spatial distribution of the activities of the industrial fleet was analyzed using logbook and VMS information. The logbook information corresponded to the period 2002-2010 for the entire Cantabrian Sea (southern Bay of Biscay), and it was used to identify the fishing gear and fishing tactics (Marchal et al. 2006). The VMS information covered the period 2006-2010, and only the data corresponding to the study area were used. VMS data were also used to carry out the spatial analysis of the fishing tactics previously identified in the logbooks. The VMS and logbook data were provided by the Spanish Ministry of Agriculture, Food and Environment (MAGRAMA, for its initials in Spanish).

A VMS is a tracking device that sends a signal (ping) every 2 h that contains the code of the emitting ship, date, time, position (latitude and longitude), speed in knots, course, and whether the vessel is carrying out fishing operations or not. Many filters and processing techniques are used to eliminate signals not related to fishing activity (Lee et al. 2010, Hintzen et al. 2011). In this study, the following processing technique was used: the time interval and the Euclidean distance between successive signals were obtained, and each of these values was associated with the first signal of each corresponding pair; when the time interval between signals was longer than 4 h, the beginning and end of each fishing expedition was determined; the average speed of the vessel was calculated using the interval between successive signals (pings); vessels with less than 10 signals in a year were eliminated; signals recorded within a distance of 3 miles or less from the closest fishing harbor were also eliminated. Each signal coinciding with a fishing trip registered in the logbooks (according to the vessel code and the date of capture) was associated with a fishing gear and a fishing tactic. Based on the distribution of the frequencies of average speeds, a working range for each fishing gear was defined, and all signals with associated velocities out of the working range were eliminated (Hintzen et al. 2011). For the study area, 51.6% of 46,075 VMS fishing days could be matched with their logbook records for the period 2006-2010 and both fishing gear and fishing tactic could be assigned to each matched day.

The frequency distribution of the average velocities was used to determine the average speed ranges at which we considered fishing activity to be carried out (Fig. S1). The identification of these ranges can be achieved by either locating changes in the tendency through the use of regression models (segmented regression) or using available information from the fleet and on-board observers. In the case of mobile fisheries (trawling and handline trawling), both methods are used, whereas in the case of static fisheries (purse seine, gillnet, longline, and handline fishing), it is better to use the information provided by on-board observers and the information obtained from the fleet according to the frequency distribution. The speed ranges by fishing gear can be seen in Table 1.

Table 1 Speed range (nm) by gear when it was considered that there was fishing. 

On wide continental shelves with sedimentary bottoms such as in the North Sea, the fishing areas are generally extensive. In these cases, after filtering and analyzing VMS data, the distribution of fishing efforts can clearly be observed. Therefore, areas containing points with no fishing activity are scarce (Hintzen et al. 2011). However, fishing areas located on narrow continental shelves with rocky mosaic bottoms or on continental slopes are generally narrow and short. In these cases, it is frequent to find areas where fishing activity is apparently being carried out when, in reality, there is no activity. To eliminate false fishing activity, it is necessary to apply thresholds to define when fishing activity is actually carried out.

Thresholds were applied to determine if the filtered VMS data corresponded to real fishing activity. Thus, the effort value was assigned to each corresponding point where fishing was detected (see "Estimation of relative spatial effort" section below). Subsequently, it was necessary to set a threshold for the effort value below which fishing activity was considered to be negligible or non-existent, for which there are numerous techniques (Jiménez-Valverde and Lobo 2007). The techniques that provided the best outcomes were based on applying quartile thresholds. Points at which the presence of fishing was detected were eliminated according to the distribution of the effort frequencies (Fig. 2). The use of other statistics based on generating pseudo-absences to distinguish points with real presence was not fruitful because efforts were located in specific places. Thus, when generating pseudo-absences and eliminating those that corresponded to real activity locations, pseudo-absences always occurred over zones where there was no activity; consequently, the threshold values were so low that we could not distinguish between real and false presences. Thresholds can be applied to fishing tactics or to fishing gear. If the relative importance of thresholds is small, it is better to apply them to fishing tactics; otherwise, it is more effective to apply them to fishing gear. In our case, comparing the distribution of effort with the distribution of pings (Fig. 2), when a 0.2 threshold was applied all those points that were outside the fishing area were eliminated without affecting the estimation of the effort of the main fishing areas. After inspecting the spatial distribution of the effort of all fishing gears the threshold was applied to all.

Figure 2 Threshold selection for the pair trawl fishery. 

Identification of fishing tactics

The terminology used in the present study was defined in ICES (2003) and used by Marchal et al. (2006): fleet, a group of vessels with similar characteristics, and fishing tactic, trips made targeting the same species/stock or a combination of these using similar gear in the same period of the year and in the same area. The identification of fishing techniques was performed by gear. To standardize the landings by trip for each of the gears, the species composition of each trip was expressed as a percentage of the total amount caught (Jiménez et al. 2004, Campos et al. 2007). Those species representing less than 1% of the total landings were incorporated into the matrix in a group called "other". Given the doubtful classification of deep-sea sharks conducted on board, this taxonomic group was separated into 2 categories: deep-sea sharks and other sharks, which included all remaining sharks.

The non-hierarchical cluster technique CLARA (Clustering Large Applications, Kaufman and Rousseeuw 1986) was used to identify and classify fishing tactic by fishing gear. The silhouette width (Rousseeuw 1987) was used to select the number of clusters, and to validate the consistency of the clustering we followed the criteria of Struyf et al. (1996), in which 4 levels are defined for the global silhouette coefficient and partial silhouette coefficient: 0.71-1.00, consistent pattern; 0.51-0.70, reasonable pattern; 0.26-0.50, weak pattern; <0.26, no pattern. A detailed explanation of the use of this technique in fisheries can be found in Castro et al. (2010).

The analyses were performed using R (R Development Core Team 2005). The correspondence and principal component analyses were performed using the "ade4" library (Dray and Dufour 2007), and the cluster analyses were performed with the "cluster" (Maechler et al. 2005) and "stats" (R Development Core Team 2005) libraries.

Estimation of relative spatial effort

The spatial effort was estimated using a grid spacing of 250 m. Thus, a 250-m search radius was chosen in each node such that adding up the time assigned to each ping emitted by the VMS corresponded to the presence detected within each radius. By taking a radius whose length was the same as the grid spacing, we achieved better geographic precision of effort, overestimating those locations where presence was more frequent.

Estimation of overlap between demersal fisheries

The overlap between demersal fishing gears was analyzed: bottom otter trawl (OTB), bottom pair trawl (PTB), set gillnet (GNS), and set longline (LLS). When analyzing the overlap between fishing activities, it is necessary to take into account that the VMS indicates the position of the boat, not the position of fishing gear. In order to avoid false positives or false negatives, values greater than 20% with respect to the maximum value of effort were taken into account and the resolution of layers was lowered to 1,000 m. To analyze the overlap, each pixel with activity of the layers analyzed was given the value 1, and the layers with the fishing techniques analyzed were added (PTB vs OTB, GNS vs LLS, and set gear vs mobile gear). The result was a final raster layer with the following values: NA, absence of activity; 1, activity of one of the 2 fishing techniques analyzed; and 2, overlap of the activity.

Results

Fishing tactics in the study area

All of the fishing gears were present in the study area: OTB, PTB, GNS, LLS, purse seine (PS), and handline (HLM). Longline targeting Dicentrarchus labrax (sea bass) and purse seine targeting Engraulis encrasicolus (anchovy) were not considered in the present study because they were poorly represented in the study area (3 and 6 fishing trips, respectively). The global silhouette coefficient by fishing gears showed values equal to or greater than 0.5; hence, all of the classifications could be considered robust (Table 2). A total of 17 fishing tactics were identified. For all fishing tactics the target species appear in 100% of the trips (frequency of occurrence, FO), except for mixed otter trawl (OTB MIX), purse seine targeting Trachurus spp. (PS_HOM), and handline targeting tuna (HLMTUN).

Table 2 Catch profile expressed as percent of total weight and frequency of occurrence (FO) by gear (GCS, global silhouette coefficient) and fishing tactic (PCS, partial silhouette coefficient). Only the species that comprised more than 1% of the catch are shown. 

Table 2 (Cont.) 

Three otter trawl fishing tactics were identified. The first was a mixed fishing tactic (OTB_MIX) targeting several demersal species, none of them accounting for more than 17% of the landings and almost 45% of the landings corresponding to 3 species (Lophius spp., monkfish; Lepidorhombus spp., megrims; and Merluccius merluccius, hake); the main species in FO were Lophius spp. The second tactic targeted Trachurus trachurus (horse mackerel) (OTB_HOM), which accounted for 80% of the landings, with some demersal species as bycatch. The third tactic (OTB_MAC) targeted Scomber spp. (mackerel), with Trachurus spp. as bycatch (FO = 43). Regarding pair trawling, only 2 fishing tactics were identified, the first targeting Micromesistius poutassou (blue whiting) (PTB_WHB), with Merluccius merluccius (hake) as the main non-targeted species (FO = 80), and the second targeting Scomber spp. (PTB_MAC) and blue whiting as main bycatch (FO = 43).

Four longline fishing tactics were identified. One tactic targeted Conger conger (conger) (LLS_COE) and another targeted Phycis blennoides (greater forkbeard) (LLS_GFB), the former capturing P. blennoides as bycatch and the latter capturing C. conger as bycatch. A third tactic targeted hake (LLSHKE), without any significant bycatch, and a fourth targeted Pollachius spp. (LLS_POL) with sparids as bycatch, which were unloaded during 58% of fishing trips (FO), thus accounting for only 7.81% of the landings by weight. In the case of all fishing tactics, bycatches were of little importance compared with the target species both in quantity and frequency.

Three gillnet fishing tactics were identified. The first tactic targeted hake (GNS_HKE) and the second targeted Lophius spp. (monkfish) (GNS MNZ), without any characteristic bycatch species. The third targeted species from the "other" group (GNS_MIX), with hake and monkfish as bycatch species.

In the area of influence of the Avilés Canyon System, purse seines and handlines are used to catch pelagic species. In the case of the purse seine net, 3 different fishing tactics were observed: the first targeted horse mackerel (PS_HOM), with Boops boops (FO = 17) and Sardina pilchardus (FO = 29) as main bycatch species; the second targeted mackerel (PS_MAC); and the third targeted Sardina pilchardus (sardine) (PS_PIL), with horse mackerel (FO = 19) and mackerel (FO = 18) as bycatch. Finally, 2 handline fishing tactics were identified, one targeting mackerel (HLM_MAC) and another targeting tuna with live bait (HLM_TUN).

Spatial distribution of fisheries by gear

Bottom otter trawling was carried out on sedimentary bottoms outside the 100-m isobath, the legal limit (Fig. 3a). In the western part of the study area, the fishing effort was located close, but below, 200 m depth, although fishing activity between 100 and 200 m and deeper than 500 m was observed. In the eastern part, the effort distribution was similar, but higher activity could be observed near 100 m. Bottom pair trawling was carried out in 3 areas: on the flanks of the Avilés Canyon and in the southern part of the rocky outcrop called El Agudo between depths of 200 and 500 m (Fig. 3b).

Figure 3 Spatial distribution by gear: (a) bottom otter trawl (OTB), (b) bottom pair trawl (PTB), (c) set gillnet (GNS), (d) set longline (LLS), (e) purse seine (PS), and (f) handline (HLM). Color scale indicates relative effort. 

The set gillnet fishery (Fig. 3c) was conducted throughout the study area except in the eastern part above 200 m. Two areas concentrate most of the effort: the eastern rocky area between 500 and 1,000 m and the region from west to southeast, around 100 m. It is necessary to take into account that we do not have environmental information for the latter zone (Fig. 1) and thus do not know if rocky areas are present there.

In the case of set longline fishing, the 3 main areas subjected to fishing efforts were located over rocky bottoms: the headwater region of the Avilés Canyon, El Agudo rocky outcrop, and the area in front of Cape Peñas (Fig. 3d). Additionally, fishing activity was carried out along the flanks of the 3 canyons, between 500 and 1,000 m.

Most of the purse seine fishing activity was located in the western area of Cape Peñas below 100 m depth (Fig. 3e). Regarding handlines, the main activity was located at the head of the canyon, although the fishing effort was scattered over the entire area (Fig. 3f).

Overlapping of demersal fishing activity

Overlapping between the 2 types of trawls occurred in the northwestern zone and slightly in the northeastern part of Avilés Canyon (Fig. 4a). The overlap between gillnets and longlines was limited to the head of the canyon and the El Agudo outcrop (Fig. 4b). When comparing the spatial distribution between the mobile and set types of gears, the overlap was very low and only occurred in the northeastern part of Avilés Canyon (Fig. 5).

Figure 4 Spatial distribution of the overlapping between demersal fisheries: (a) bottom otter trawl (OTB) and bottom pair trawl (PTB), and (b) set gillnet (GNS) and set longline (LLS). 

Figure 5 Spatial distribution of the overlapping between mobile gear (otter trawl and pair trawl) and set gear (gillnet and longline). 

Spatial distribution of fisheries by fishing tactic

The spatial distribution of each fishing tactic can be seen in Figures S2-S5. Regarding otter trawling, the main difference was that the tactic targeting horse mackerel (OTB_HOM) was used mostly in the western area and along 100 m in the eastern zone. This fishing tactic and the tactic targeting mackerel (OTB_MAC) were carried out below 200 m depth (Fig. S2a-b). The mixed fishing tactic (OTB_MIX) was observed in the whole area where this fishing gear was used (Fig. S2c). The OTBMAC tactic was carried out during the first 5 months of the year.

In the case of bottom pair trawling, we found 2 very different distribution patterns of effort for each fishing tactic. Thus, whereas the tactic targeting blue whiting (PTB_WHB) was spread across the entire area where pair trawling was carried out (Fig. S2d), the tactic targeting mackerel was centered on the western part of the continental margin of Avilés Canyon from 200 to 500 m depth (Fig. S2e). As observed for OTB_MAC, PTB_MAC was carried out during the first 5 months of the year.

Each of the 3 gillnet fishing tactics presented a different spatial distribution (Fig. S3). In the case of the mixed fishery (GNS_MIX), there did not appear to be any defined pattern, and this fishery operated throughout the entire area where gillnets were used. Gillnetting for hake (GNSHKE) was carried out in the western area of the continental shelf, between 100 and 200 m depth, and in the rocky area located in front of Cape Peñas. We did not have any information on the type of bottom in the western part of the platform and the eastern end at 100 m depth. In the case of the tactic targeting monkfish (GNS_MNZ), this activity was carried out at depths between 500 and 1,000 m and within the limits between the soft and rocky substrata surrounding the rocky area located in front of Cape Peñas.

The longline fishing tactic targeting forkbeard (LLS_GFB) was carried out primarily in the deep-sea areas located on the flanks of the three canyons, between 500 and 1,000 m depth (Fig. S4a). The tactic targeting conger (LLS_COE) was used at a shallower depth along the rocky area in front of Cape Peñas and at the 100 m depth contour line (Fig. S4b). The fishing effort of the tactic targeting hake (LLS_HKE) was concentrated in the headwaters of Avilés Canyon, in El Agudo, and in the outcrops located in front of Cape Peñas (Fig. S4c). The LLSPOL tactic was mainly used along the 100-m isobath (Fig. S4d).

The spatial distributions of the purse seine fishing techniques targeting horse mackerel and sardine were very similar because the fishing activity was carried out close to the shore and to the west of Cape Peñas (Fig. S5a-b). In the case of the fishery targeting mackerel, the fishing effort was concentrated in the middle of the continental shelf and close to the headwaters of the canyon (Fig. S5c). As for seasonality, all of the fishing tactics targeting mackerel were carried out during the first 4 months of the year.

In the case of handline fishing for mackerel, the effort distribution was very similar to purse seine targeting the same species (Fig. S5d). The tactic targeting tuna was carried out in the rocky area located in the waters of Avilés Canyon and in El Agudo (Fig. S5e).

Discussion

On narrow continental shelves subjected to strong fishing pressure, as in the case of the study area, fishing tactics are used in different areas to share the space available and avoid friction caused by the spatial overlapping of fishing activities. In this context, and considering the strong spatial segregations observed in this study, the filters used to eliminate signals not corresponding to fishing activity were effective. The validity of these filters was checked by analyzing the spatial distribution of bottom trawling because this fishing gear can be used only on sedimentary bottoms and is subject to specific legal regulations. Once filtering techniques were applied, trawling was not observed at depths above 100 m (closed zone to bottom trawling in Spanish territorial waters in the North Atlantic) or in harbor entrances. The fleet using bottom trawling cannot operate in rocky areas and, as we observed, bottom trawling fishing tactics were not used on rocky bottoms.

In the case of fishing gears already described in previous studies, the same structure was obtained for the fishing tactics (Castro et al. 2010, Punzón et al. 2011). Among those identified, mixed gillnet (GNSMIX) should not be considered a fishing tactic even though it presents a high partial silhouette coefficient value (Table 2). This high value can be attributed to the inclusion of a large number of under-represented species in the "other" group. In addition, it was necessary to take into account that the main bycatches were hake and monkfish and that the spatial distribution of GNS_MIX appeared to be the result of the GNS_HKE and GNS MNZ distributions, so it was probable that catches from the latter tactics were being included in the first.

Using statistical rectangles to define if a fishery should be considered a deep-water fishery, like LLS COE, may be problematic in narrow continental shelves (Punzón et al. 2011). As we have observed, the activity of this fishery takes place along the entire continental shelf, including the area close to the 500-m contour line, and for that reason it should be considered a deep-sea fishery. Therefore, the VMS is a powerful tool to characterize the fisheries in narrow continental shelves, especially for fisheries operating in areas of special interest where spatial management measures are to be implemented.

In the reviews conducted by Piñeiro et al. (2001) and Punzón et al. (2011), the authors identified 2 deep-water fisheries in the Northeast Atlantic region: GNSMNZ and OTBMIX. In the case of GNSMNZ, the fishing effort was located below 500 m depth, while in the case of OTBMIX, the fishing effort deeper than 500 m was carried out only in the western region and was not very intense. Although some fishing activity of PTBWHB and LLSHKE was observed close to the 500-m contour line, we could not determine whether these activities should be considered deep-water fisheries or not because it was not clear whether they were carried out below 500 m depth.

The spatial distribution of fisheries by type of fishing gear showed that they barely interacted in the study area, at least not in zones with the highest fishing effort. This small overlap implies that almost all the study area is subject to fishing effort; however, there is little competition between the fishing techniques for common resources. The almost non-existent overlap among the various types of activity will further serve to allocate fishing techniques with pings that were not used in the past due to the lack of logbook information (Russo et al. 2011).

Of the 2 trawling techniques used in the area, otter trawling was carried out in larger areas. This distinction was due to the fact that bottom otter trawls are smaller longitudinally and narrower at the mouth, so this fishing gear can be handled more easily. In other words, because of its better handling qualities and minor surface contact, this fishing gear can be used in narrower and rougher fishing grounds. Bottom otter trawls can also be used in the same fishing grounds where pair trawling is being carried out. However, in these areas, otter trawling activity was low. The main reason was that pair trawling targeted blue whiting, a species inhabiting the break shelf (Castro et al. 2010), but it ceased to be the target species at the beginning of the 21st century (Punzón et al. 2010). Since these fisheries did not compete for the same resource, they operated in different locations. Both trawling techniques were limited to specific areas, and they were carried out at a constant depth following the bathymetry to avoid affecting artisanal fisheries. The region located on the east side of the study area had a wider continental shelf and a smoother sea bottom. In this area, trawling patterns were changeable compared with the shelf break and rocky zones. This was possible because of the small number of artisanal fisheries operating in the western region of the continental platform, leaving this area to be used exclusively for trawling. However, the number of fisheries increased considerably in the shelf breaks, so the trawling patterns became more regular.

The main spatial overlap between fishing tactics was observed between gillnets targeting monkfish and longlines targeting hake in the rocky outcrop called El Agudo. The bulk of the fishing effort exerted in this area was related to gillnet fishing, whereas the activity related to longline fishing was of less importance. To determine the extent of the interaction between the different fishing gears in the area, we need further information about the characteristics of the local fishing fleet.

It is worth mentioning that the different fishing tactics were located in distinct areas depending on the sea bottom characteristics and, in the particular case of sandy bottoms, on the substrate type. Thus, the spatial distribution of the different tactics and fishing gears used was strongly influenced by this variable in the study area, affecting the way the fishing activity was interpreted. For instance, GNS_HKE was present in the rocky area where information about the sea bottom was available, whereas in the western part of the continental shelf, where we did not have information about the bottom type, its spatial distribution appeared to be related to a specific type of sea bottom. Thus, it was difficult to define the spatial interactions between fishing gears when this information was not available.

Hake, mackerel, and horse mackerel were targeted by more than one fishing tactic in the study area. In the case of hake, the spatial overlap between the two tactics (LLS_HKE and GNS_HKE) was not important. Therefore, hake was caught by longline gear in the head of Avilés Canyon, whereas gillnetting was used in the rocky area in front of Cape Peñas and in the western area between 100 and 200 m depth. Spatial overlap was only found at El Agudo, although the longline efforts were very low in this zone.

From the end of winter to early spring, the main fishing effort in the southern part of the Bay of Biscay is directed at mackerel during its reproductive migration (Uriarte and Lucio 2001, Villamor et al. 2011). Most of the fishing gears employed in the area are used to catch mackerel (Punzón et al. 2004). This pattern could have led to a decrease in the fishing efforts of other fishing tactics, but it was only observed for bottom trawling targeting mackerel from February to March. The fishing activity of all the tactics directed at mackerel was located on the continental shelf surrounding Cape Peñas. Fisheries targeting small and medium pelagic fish generally operated in offshore waters. Only purse seine fishing was carried out close to the coast, but the bulk of the fishing effort was located beyond 100 m. This pattern might be explained by the fact that the targeted species has to round Cape Peñas during its migration from east to west, resulting in higher abundances of this species on the mid and outer continental shelf.

The use of VMS information together with information from logbooks allowed us to study not only the spatio-temporal pattern of fishing activity but also fundamental aspects of fisheries management and the ecosystems affected by this activity. In this context, we can gain knowledge of the ecosystem structures with respect to commercial species, including information about non-accessible areas for research surveys carried out in the study area based on trawling hotspots (Sánchez and Serrano 2003), which were spatially limited to sedimentary bottoms and temporarily limited during autumn. Moreover, this research allowed us to study the behavior of commercial species because the VMS signals associated with a specific catch can be interpreted as the presence of the species. The use of this information in combination with environmental data obtained from other sources could lead to a better understanding of the environmental conditions in which these catches are carried out.

Finally, by knowing the distribution patterns of fishing activity and the biological and environmental variables determining these distributions, we could predict their responses to different management measures. Thus, we could simulate the effect of different fishing effort distributions in the area before establishing MPAs if fishing activity were modified or if mutually exclusive activities were carried out. This knowledge will allow us to estimate the social and economic impacts on fishing activity as well as on the mortality of the targeted species.

Acknowledgments

This study was partially funded by the TEPESCO project (2016, Instituto Español de Oceanografía) and by the European Commission LIFE+ "Nature and Biodiversity" program (INDEMARES project 07/NAT/E/000732); the Biodiversity Foundation of the Ministry of Environment was the institution responsible for coordinating the latter. The bathymetric information and ground characteristics were made available thanks to the invaluable work of María Druet, María Gómez-Ballesteros, and Araceli Muñoz-Recio.

References

Baelde P. 2005. Interactions between the implementation of marine protected areas and right-based fisheries management in Australia. Fish. Manage. Ecol. 12: 9-18. [ Links ]

Campos A, Fonseca P, Fonseca TJ, Parente J. 2007. Definition of fleet components in the Portuguese bottom trawl fishery. Fish. Res. 83: 185-191 [ Links ]

Castro J, Punzón A, Pierce GJ, Marín M, Abad E. 2010. Identification of métiers of the northern Spanish coastal bottom pair trawl fleet by using the partitioning method CLARA. Fish. Res. 102: 184-190. [ Links ]

Claudet J, Pelletier D. 2004. Marine protected areas and artificial reefs: A review of the interactions between management and scientific studies. Aquat. Living Resour. 17: 129-138. [ Links ]

Cury PM, Mullon C, García SM, Shannon LJ. 2005. Viability theory for an ecosystem approach to fisheries. ICES J. Mar. Sci. 62: 577-584. [ Links ]

Dray S, Dufour AB. 2007. The ade4 package: Implementing the duality diagram for ecologists. J. Stat. Softw. 22: 1-20. [ Links ]

Fréon P, Drapeau L, David JH, Fernandez Moreno A, Leslie RW, Oosthuizen WH, Shannon LJ , van der Lingen CD. 2005. Spatialized ecosystem indicators in the southern Benguela. ICES J. Mar. Sci. 62: 459-468. [ Links ]

García SM, Cochrane KL. 2005. Ecosystem approach to fisheries: A review of implementation guidelines. ICES J. Mar. Sci. 62: 311-318. [ Links ]

Gislason H, Sinclair M, Sainsbury K, Boyle R. 2000. Symposium overview: Incorporating ecosystem objectives in fisheries management. ICES J. Mar. Sci. 57: 468-475. [ Links ]

Gómez-Ballesteros M, Druet M, Muñoz A, Arrese B, Rivera J, Sánchez F, Cristobo J, Parra S, García-Alegre A, González-Pola C, Gallastegui J, Acosta J. 2014. Geomorphology of the Avilés Canyon System, Cantabrian Sea (Bay of Biscay). Deep-Sea Res. (II): 106: 99-117. [ Links ]

Hintzen N, Bastardie F, Beare D, Piet G, Ulrich C, Deporte N, Egekvist J, Degel H. 2011. VMStools: Open-source software for the processing, analysis and visualization of fisheries logbook and VMS data. Fish. Res. 115: 31-43. [ Links ]

Hollingworth CE (ed.). 2000. Ecosystem Effects of Fishing. Proceedings of an ICES/SCOR Symposium held in Montpellier, France, March 1999. ICES J. Mar. Sci. 57: 465-791. [ Links ]

[ICES] International Council for the Exploration of the Sea. 2003. Study Group on the Development of Fishery-based Forecasts. ICES CM 2003/ACFM:08, Boulogne, France. [ Links ]

Jennings S, Kaiser M. 1998. The effects of fishing on marine ecosystems. Adv. Mar. Biol. 34: 201-352. [ Links ]

Jiménez MP, Sobrino I, Ramos F. 2004. Objective methods for defining mixed-species trawl fisheries in Spanish waters of the Gulf of Cádiz. Fish. Res. 67: 195-206. [ Links ]

Jiménez-Valverde A, Lobo J. 2007. Threshold criteria for conversion of probability of species presence to either-or presence-absence. Acta Oecol. 31: 361-369. [ Links ]

Jones PJ. 2001. Marine protected area strategies: Issues, divergences and the search for middle ground. Rev. Fish Biol. Fish. 11: 197-216. [ Links ]

Kaufman L, Rousseeuw PJ. 1986. Clustering large sets (with discussion). In: Gelsema ES, Kanal LN (eds.), Pattern Recognition in Practice II. Elsevier/North Holland, Amsterdam, pp. 405-416. [ Links ]

Lee J, South AB, Jennings S . 2010. Developing reliable, repeatable, and accessible methods to provide high-resolution estimates of fishing-effort distributions from vessel monitoring system (VMS) data. ICES J. Mar. Sci. 67: 1260-1271. [ Links ]

Link JS. 2005. Translating ecosystem indicators into decision criteria. ICES J. Mar. Science. 62: 569-576. [ Links ]

Maechler M, Rousseeuw PJ , Struyf A, Hubert M. 2005. Cluster: Cluster Analysis Basics and Extensions. R package version 1.10.2, http://cran.r-project.org/web/packages/cluster/. [ Links ]

Marchal P, Andersen B, Bromley D, Iriondo A, Mahévas S, Quirijns F, Rackham B, Santurtun M, Tien N, Ulrich C . 2006. Improving the definition of fishing effort for important European fleets by accounting for the skipper effect. Can. J. Fish. Aquat. Sci. 63: 510-533. [ Links ]

Pauly D, Christensen V, Guénette S, Pitcher TJ, Sumaila UR, Walters CJ, Watson R, Zeller D. 2002. Towards sustainability in world fisheries. Nature 418: 689-695. [ Links ]

Pikitch EK, Santora C, Babcock EA, Bakun A, Bonfil R, Conover DO, Dayton P, Doukakis P, Fluharty D, Heneman B, Houde ED, Link J, Livingston PA, Mangel M, McAllister MK, Pope J, Sainsbury KJ. 2004. Ecosystem-based fishery management. Science 305: 346-347. [ Links ]

Piñeiro C, Casas M, Bañón R. 2001. The deep-water fisheries exploited by Spanish fleets in the Northeast Atlantic: A review of the current status. Fish. Res. 51: 311-320. [ Links ]

Punzón A , Hernández C, Abad E , Castro J , Pérez N, Trujillo V. 2010. Spanish otter trawl fisheries in the Cantabrian Sea. ICES J. Mar. Sci. 67: 1-13. [ Links ]

Punzón A , Serrano A, Castro J , Abad E , Gil J, Pereda P. 2011. Deep-water fishing tactics of the Spanish fleet in the Northeast Atlantic. Seasonal and spatial distribution. Sci. Mar. 75: 465-476. [ Links ]

Punzón A , Villamor B, Preciado I. 2004. Analysis of the handline fishery targeting mackerel (Scomber scombrus, L.) in the north of Spain (ICES Division VIIIbc). Fish. Res. 69: 189-204. [ Links ]

R Development Core Team. 2005. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. ISBN 3-900051-07-0, URL http://www.R-project.org. [ Links ]

Rousseeuw PJ . 1987. Silhouettes: A graphical aid to the interpretation and validation of cluster analysis. J. Comput. Appl. Math. 20: 53-65. [ Links ]

Russo T, Parisi A, Cataudella S. 2011. New insights in interpolating fishing tracks from VMS data for different métiers. Fish. Res. 108: 184-194. [ Links ]

Sánchez F, Olaso I. 2004. Effects of fisheries on the Cantabrian Sea shelf ecosystem. Ecol. Model. 172: 151-174. [ Links ]

Sánchez F , Serrano A . 2003. Variability of groundfish communities of the Cantabrian Sea during the 1990s. ICES Mar. Sci. 219: 249-260. [ Links ]

Sinclair M , Arnason R, Csirke J, Karnicki Z, Sigurjonsson J, Rune Skjoldal H, Valdimarsson G. 2002. Responsible fisheries in the marine ecosystem. Fish. Res. 58: 255-265. [ Links ]

Struyf A , Mia H, Rousseeuw PJ . 1996. Clustering in an object-oriented environment. J. Stat. Softw. 1(4): 1-30. [ Links ]

Uriarte A, Lucio P. 2001. Migration of adult mackerel along the Atlantic European shelf edge from a tagging experiment in the south of the Bay of Biscay in 1994. Fish. Res. 50: 129-139. [ Links ]

Villamor B , González-Pola C , Lavín A, Valdés L, Lago de Lanzón A, Franco C, Cabanas JM, Bernal M, Hernández C , Iglesias M, Carrera P, Porteiro C. 2011. Environmental control of Northeast Atlantic mackerel (Scomber scombrus) recruitment in the southern Bay of Biscay: Case study of failure in the year 2000. Fish. Oceanogr. 20: 397-414. [ Links ]

Supplementary material

Figure S1 Distribution of the mean speed by gear: bottom otter trawl (OTB), bottom pair trawl (PTB), set gillnet (GNS), set longline (LLS), purse seine (PS), and handline (HLM). 

Figure S2 Spatial distribution of fisheries by fishing tactic for bottom otter trawl (OTB) and bottom pair trawl (PTB): (a) OTB_HOM, (b) OTB_MAC, (c) OTB_MIX, (d) PTB_WHB, and (e) PTB_MAC. Color scale indicates relative effort. 

Figure S3 Spatial distribution of fisheries by fishing tactic for gillnet (GNS): (a) GNS_MIX, (b) GNS_HKE, and (c) GNS_MNZ. Color scale indicates relative effort. Color scale indicates relative effort. 

Figure S4 Spatial distribution of fisheries by fishing tactic for set longline (LLS): (a) LLS_GFB, (b) LLS_COE, (c) LLS_HKE, and (d) LLS_POL. Color scale indicates relative effort. 

Figure S5 Spatial distribution of fisheries by fishing tactic for purse seine (PS) and handline (HLM): (a) PS_HOM, (b) PS_PIL, (c)PS_MAC, (d) HLM_MAC, and (e) HLM_TUN. Color scale indicates relative effort. 

Supplementary material

Figura S1 Distribución de la velocidad media por arte de pesca: arrastre de fondo con puertas (OTB), arrastre de fondo a la pareja (PTB), enmalle (GNS), palangre de fondo fijo (LLS), cerco (PS) y línea de mano (HLM). 

Figura 2 Distribución espacial de la actividad pesquera por táctica pesquera para el arrastre de fondo con puertas (OTB) y el arrastre de fondo a la pareja (PTB): (a) OTB_HOM, (b) OTB_MAC, (c) OTB_MIX, (d) PTB_WHB y (e) PTB_MAC. La escala muestra el esfuerzo relativo. 

Figura S3 Distribución espacial de las pesquerías por táctica pesquera para el enmalle (GNS): (a) GNS_MIX, (b) GNS_HKE y (c) GNS_MNZ. La escala muestra el esfuerzo relativo. 

Figura S4 Distribución espacial de las pesquerías por táctica pesquera para el palangre de fondo fijo (LLS): (a) LLS_GFB, (b) LLS_COE, (c) LLS_HKE y (d) LLS_POL. La escala muestra el esfuerzo relativo. 

Figura S5 Distribución espacial de las pesquerías por táctica pesquera para el cerco (PS) y la línea de mano (HLM): (a) PS_HOM, (b) PS_PIL, (c) PS_MAC, (d) HLM_MAC, (e) HLM_TUN. La escala muestra el esfuerzo relativo. 

Received: December 2015; Accepted: October 2016

* Corresponding author: antonio.punzon@st.ieo.es

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License