The octopus fishery on the Campeche Bank (Yucatán, Mexico) is considered the third largest in the world. In Yucatán, two fleets target this resource: an artisanal fleet and a semi-industrial fleet. The artisanal fleet only catches Octopus maya, while the semi-industrial fleet catches two species, O. maya and O. "vulgaris” Type II, because it operates at deeper waters (>30 m). Since there is no information on the abundance of O. "vulgaris” Type II, management is based only on O. maya. In order to generate information about the abundance of this species, four fishing research cruises were carried out in the northeastern area of the continental shelf off the Yucatán Peninsula. Four methods (a stratified random method, a swept area, geostatistics and a weighted swept area) were applied and compared to determine the instantaneous abundance and biomass of both species in the study area. The lowest potential biomass was calculated with the geostatistical method, with values between 18.5% and 36.7% lower than the other three methods. O. "vulgaris” Type II showed the lowest biomass (37.8±3.36 t) during May and July and the highest (189.56±11.6 t) in December. Our findings revealed that the total abundance of both species was similar in the study area, with a geographic overlap whose amplitude changed throughout the year according to the geographic position: O. maya dominated at approximately 88°W, while O. "vulgaris” Type II dominated towards the southeast at 87°W.
RESUMEN
La pesquería de pulpo en el banco de Campeche (Yucatán, México) se considera la tercera productora mundial de pulpo. En Yucatán, dos flotas capturan este recurso: una flota artesanal y una semi-industrial. La flota artesanal solo captura Octopus maya, mientras que la flota semi-industrial captura dos especies: O. maya y O. "vulgaris” Tipo II. Esto se debe a que la flota semi-industrial opera en aguas más profundas (>30 m). Dado que no existe información sobre la abundancia de O. "vulgaris” Tipo II, el manejo se basa únicamente en la evaluación de O. maya. Con objeto de generar información sobre la abundancia de esta especie, se realizaron cuatro campañas de prospección pesquera en la zona noreste de la plataforma continental de la Península de Yucatán. Se aplicaron cuatro métodos (método aleatorio estratificado, área barrida, geoestadística y área barrida ponderada) y sus resultados se compararon para determinar la abundancia y biomasa instantáneas de ambas especies en el área de estudio. La biomasa potencial más baja se obtuvo con el método geoestadístico, con valores de 18.5 a 36.7% más bajos que con los otros tres métodos. O. "vulgaris” Tipo II mostró una biomasa más baja (37.8±3.36 t) durante mayo y junio, mientras que la más alta (189.56±11.6 t) fue en diciembre. Nuestros hallazgos revelaron que la abundancia total para las dos especies fue similar; además ambas especies presentaron una superposición geográfica, cuya amplitud cambió con el período del año y la posición geográfica: O. maya fue dominante en la longitud de los 88°W, mientras que O. "vulgaris Tipo II dominó hacia el sudeste de la plataforma de Yucatán (aproximadamente 87°W).
KEYWORDSbiomassdistribution patternOctopus mayaOctopus "vulgaris” Type IICampeche BankYucatán PeninsulaMexicoPALABRAS CLAVEbiomasapatrón de distribuciónOctopus mayaOctopus "vulgaris” Tipo IIbanco de Campechepenínsula de YucatánMéxicoINTRODUCTION
The family Octopodidae contains more than 200 valid species. However, global octopus fisheries are based on only around 30 species, which are harvested in a range of fisheries from subsistence catches through to valuable, large-scale commercial fisheries (Norman et al. 2014). Since some of these catches are never properly identified, they are commonly reported as Octopus vulgaris. Due to the new findings in octopus taxonomy (Norman et al. 2014, Amor et al. 2017, Van Nieuwenhove et al. 2019), all landings reported under that name should be reviewed (e.g. Solís-Ramírez 1994, Norman et al. 2014, Emery et al. 2016). This confusion occurs in Mexico, one of the world’s main producers of octopus, where two species are reported to be caught, O. maya and O. "vulgaris”, the latter accounting for between 30% and 50% of the total landings (Fig. 1) (Velázquez-Abunader et al. 2013, Galindo-Cortes et al. 2014, Sauer et al. 2020). For many years, managers and researchers have assumed that O. maya and O. "vulgaris” occupy different areas on the continental shelf of Yucatán, separated by depth (Salas et al. 2008, Beléndez-Moreno et al. 2014, DOF 2016). For instance, areas of higher abundance of O. maya were described along the coast, decreasing towards deeper waters (30 m deep at approx. 88°W) and increasing towards the west of Campeche Bank (Solís-Ramírez and Chávez 1986, Gamboa-Álvarez et al. 2015, Avendaño et al. 2019).
Trend of capture of Octopus "vulgaris” Type II (dark-green line) (r=0.59) and Octopus maya (orange line) (r=0.58) during the period 1998 to 2018 (source: Federal delegation of the SAGARPA in Yucatán / Department of Fishery Administration and Velázquez-Abunader et al. 2013).
Despite its importance, there has been no stock assessment of O. "vulgaris” in the Yucatán Peninsula (Salas et al. 2008, Galindo-Cortes et al. 2014, Sauer et al. 2020), which is the largest octopus-producing region in Mexico. In addition, some authors have noted that O. "vulgaris” from Yucatán has been misidentified (Amor et al. 2017, 2019, Ritschard et al. 2019, Van Nieuwenhove et al. 2019). It will therefore be referred to as O. "vulgaris” Type II in this article.
The National Institute of Fisheries of Mexico (INAPESCA) has focused only on O. maya to determine management measures due to the difficulties of monitoring the fishing zones away from the coast where O. "vulgaris” Type II occurs (>30 m depth) (Gamboa-Álvarez et al. 2015, Ávila-Poveda et al. 2016). The assessment of O. maya included a total allowable catch (Salas et al. 2008, Velázquez-Abunader et al. 2013, DOF 2016). Other regulations for that fishery include a minimum landing size, a closed season for fishing and authorized fishing gear (DOF 2016), none of which consider the characteristics of O. "vulgaris” Type II.
Octopus production in Yucatán amounts to about 20000 t per year (Fig. 1) (Velázquez-Abunader et al. 2013), and it is considered that the O. maya catch has reached the maximum sustainable yield (Jurado-Molina 2010, Beléndez-Moreno et al. 2014, Galindo-Cortes et al. 2014). However, the production of O. "vulgaris” Type II has shown a tendency to increase over time (Fig. 1), and in recent years has accounted for almost 50% of the total catch (Velázquez-Abunader et al. 2013). The catch of O. "vulgaris” Type II has increased due to the expansion of fishing activities and the introduction of a new medium-scale fleet (vessels between 12 and 15 m in length), which operates in deeper zones (~30 m deep), where it overlaps with the large-scale fleet (>15 m in length) (Salas et al. 2019, Beléndez-Moreno et al. 2014). This new medium-scale fleet also overlaps with the artisanal fleet (boats of 8 to 12 m in length) that targets O. maya in shallow waters (DOF 2016, Salas et al. 2019).
Conducting fisheries management with limited knowledge of the biology of a species is risky (Pauly 1984, Leporati et al. 2015), and in some cases may result in overexploitation (Norman et al. 2014, Lima et al. 2017). In this regard, the Fishery Management Plan for the Octopus of the Gulf of Mexico and the Caribbean Sea (DOF 2016) underlines the need to enhance our scientific knowledge of O. "vulgaris” Type II because of the risk posed by the continuous growth in its landings.
There are several techniques for stock assessment, and their use depends on the information available and species characteristics (Hilborn and Walters 1992, Nevárez-Martínez et al. 2000, Hernández-Flores et al. 2015). Given that existing long time series and most of the information used for stock assessment of octopus from Yucatán are based on O. maya, the objective of this study was to assess the biomass status of O. "vulgaris” Type II and to determine its distribution on the northeast continental shelf of Yucatán at depths between 30 and 60 m. The data for performing this assessment were obtained from surveys conducted during scientific research cruises carried out independently from the fishery, so this information provides a preliminary analysis of the dynamics of a resource that has been catalogued with a status of potential exploitation. The analysis also explores the area of overlap between O. "vulgaris” Type II and O. maya.
MATERIALS AND METHODSStudy area
The study area is located in the northeastern zone of the Yucatán Peninsula, on the continental shelf known as Campeche Bank, between 30 and 60 m depth (Fig. 2). The zone is influenced by the Loop current, which enters the system through the Yucatán channel, circulating towards the Straits of Florida and generating seasonal upwellings from May to August that stratify the water column, forming colder zones at greater depth (Tester et al. 1991, Lee and Williams 1999, Paris et al. 2005). From October to February, it is influenced by cold and dry air masses coming from the north of the continent (Canada and United States of America), which, when merging with air masses of tropical origin, mainly from October to March, produce strong frontogenesis (gusts >100 km h) and temperature drops (Zavala-Hidalgo et al. 2002, Morey et al. 2006, Dubranna et al. 2011), locally known as nortes (Enríquez et al. 2010).
Study area and stations sampled in May-June, July and December 2016 and January 2017. Fishing area described for Octopus "vulgaris” Type II to the East of the Campeche Bank, Mexico (dark-green dashed line) and from Octopus maya (orange dashed line) on the east platform of Campeche Bank.
Field work
Four research cruises were conducted on a semi-industrial vessel, from the port of Progreso, Yucatán, before and after the fishing season (from 1 August to 15 December). Fishing operations targeted O. "vulgaris” Type II at 29 stations on average (±2) (Fig. 2). The cruises took place in five months: May-June and July 2016, and December 2016 and January 2017.
The method used to collect octopus samples was the traditional fishing technique known as "gareteo”. This technique consists of placing one bamboo pole (jimba) in the bow and another in the stern of a small wooden boat (<4 m long), known locally as alijo. These bamboo poles were approximately 4 m long (for the medium-scale fleet) (Salas et al. 2008, Sauer et al. 2020), and each had two monofilament lines with a specimen of Diplectrum sp. or Haemulon sp. tied to the tip as bait. Five alijos were deployed per cruise and each one had a GPS to trace the route and thus determine the swept area. The effective fishing time was standardized in one sampled by three hours per alijo (five) at each sampling station; thus, initial and final time of the fishing operations were recorded. At each fishing station or site, the total number of organisms caught (Nt), total catch weight (Wt) and total weight of each organism (wi) were recorded.
Determination of the area per sampling station
Thiessen polygons (Brassel and Reif 1979) were deployed to identify the spatial arrangement of the sample as a function of distance between the sampling points. This allowed us to determine the optimal (constant) area of representativeness of each sampling point to apply any of the three area-based methods: stratified random method, swept area and weighted swept area. This procedure allows us to assign a real significance to sampling point values and therefore tells us the area on which the octopus abundance calculation was made according to the captures at each sampling station (Avendaño et al. 2019). The Thiessen polygons and the area of each polygon were estimated using the ArcMap 9.2 software (Sawatzky et al. 2009).
Biomass estimates
Four statistical methods were applied: the stratified random method (Cochran 1978), the swept areas method of Rosenberg et al. (1990), the Geostatistics biomass model (Rivoirard et al. 2008) and the method of weighted swept areas (Avendaño et al. 2019). These methods increase the precision of biomass estimation, thus generating robust data concerning the availability of octopus in the deepest areas (30 to 60 m) in the eastern part of the platform of the Yucatán Peninsula.
The stratified random method (Cochran 1978) uses the frequency of the distribution of the total weight of the catch per station, which must be classified by strata. The formula for calculating the biomass by means of the stratified method was obtained from the average count (kg) in the i-th stratum i(ȳi):
(1)
Estimator of the variance of ȳi:
(2)
Estimator of the total population size expressed in kg (Nȳst):
(3)
Estimator of the variance of the total size of the population
(Nȳst):
(4)
For the estimation of 95% confidence intervals:
(5)
where Ni is the total number of units sampled in the i-th stratum, L is the number of strata, ni is the number of units sampled in the sample of the i-th stratum, ȳi is the average weight in the i-th stratum, and S2iis the variance of the counts in the i-th stratum.
The method of swept areas of Rosenberg et al. (1990) considered the catch in weight (biomass) obtained from the swept area, dividing the capture zone according to the observed abundances of the organisms:
(6)
with variance
(7)
where Yt is the total catch in the study area, At is the total area of the study zone, at is the total area swept by the five alijos, S2t is the variance of the total catch in the study area, mt is the number of fishing operations and
(BT) is the variance of the total biomass. In this case, ai represented the swept area covered by each alijo, so the total swept area at (expressed in km) for each fishing operation resulted from
(8)
ai was calculated using the following equation:
ai = Di* LJi,
(9)
where Di is the distance travelled by the i-thalijo, estimated using the tracking data recorded by the GPS, and LJi is the length of the bamboo rods or jimbas installed in the i-thalijo (LJi=8 m). Finally, the total abundance was calculated using the following equation:
(10)
where
is the average weight of the octopus obtained in biological samples.
The biomass calculated by the geostatistical method was determined by previous calculation of the catch per unit area (CPUA, expressed in number of octopus per km2), obtained by dividing the number of octopus captured by the area corresponding to each sampling station. The spatial correlation of the CPUA was determined by means of omnidirectional empirical variograms (isotropic), which measure the correlation between the semi-variance generated by all the differences between the data pairs separated by a previously established equal interval of distance (h) (Hernández-Flores et al. 2015). Later, a Kriging interpolation analysis was applied, giving the densities through interpolation nodes by the method of nearest neighbours (Cressie 1992), and generating a spatial structure that depends on the spatial arrangement of the population (Webster and Oliver 2007).
The empirical variograms were obtained using the following equation:
(11)
where γ(h) is the value of the variance for each distance h, N(h) is the number of pairs of observations separated by h, C(xi) is the CPUA observed at site xi and C(xi+h) is the CPUA observed at any other site separated by a distance h for each site xi. The obtained interpolations were divided into CPUA intervals, obtaining an average value for the i-th interval (
). The total abundance of each i-th interval (Ni) was obtained by multiplying the (
) by the total area covered per i-th interval; therefore, the total abundance (NT) was obtained using the equation
(12)
and the biomass was obtained using the equation
(13)
The method of weighted swept areas proposed by Avendaño et al. (2019) analyses the density of organisms as a specific datum of the area represented by each station. That is, the abundance at a station is representative of the area of that station. Later, the total biomass is calculated from the sum of the biomasses obtained in each sampled area using the following formula:
(14)
with a standard deviation of
(15)
where Yi is the total of the catch in the i-th stratum,
is the average of the catches, Ai is the total area in the i-th stratum, ai is the area swept in the i-th stratum and
(BT) is the standard deviation of the total biomass. The abundance is calculated using Equation 10.
Spatial distribution
The equation modified by Avendaño et al. (2019) from the method described by Guerra (1981) was used to determine the type of distribution, calculating the average probability of occurrence of octopus per fishing station, as well as, the type of distribution, and estimating the parameters of the negative binomial distribution p and k:
(16)
To determine whether the type of distribution of the octopus is homogeneous or aggregated, a random distribution was developed, assuming a negative binomial distribution (k), calculated as K1=
/S2–, subject to the following conditions: K/ >6 if is low, K>13 if is high or
if is moderate.
If none of these conditions can be verified, then it is considered that K1 is not suitable, and the parameter is estimated by means of the following equation:
(17)
For any case, p=/K . (18)
Once the parameters had been estimated, a comparison was made between the distribution function of the total of samples and the theoretical distribution assuming a negative binomial type distribution by means of a goodness-of-fit test to check whether there was in fact a grouping (Zar 1999).
Spatial distribution of abundances
Empirical distribution models were generated using the variograms estimated in the geostatistical analysis because there was no previous evidence about the distribution of O. "vulgaris” Type II on the continental shelf of Yucatán. The model that best described the observed data was selected using two criteria: 1) the model that most closely explained the spatial distribution of the observed data, and 2) the minimum value of the sum of squares of the residuals, as well as the spatial distribution of the variances.
Analysis of geographic overlap between species
The distribution of O. "vulgaris” Type II was compared with the distribution of O. maya obtained by Avendaño et al. (2019) during the same period in the same study area. This analysis should corroborate the fact that their distribution zones are separated by depth, as described elsewhere (DOF 2016). The geographic overlap (D) for both species was determined using the Schoener equation (Schoener 1968, Cornell 2011, Afkhami et al. 2014). The Schoener index computes the percentage of similarity of two species and estimates the range of overlap in the areas where both species occur, based on the abundance of each one. Therefore, the over-geographical space is defined with a range between 0 (no overlap) and 1 (full overlap), where pxi or pyi denotes the proportion assigned to the species x (or y) in cell i.
The total geographic superposition of the area (D) was calculated using the equation:
(19)
For both cases:
(20)
(21)
where xi and yi are the number of organisms of the species x or y respectively at the sampling station i. X and Y are the total number of observations per species x or y.
RESULTSPotential biomass analysis
The niche breadth per sampling station ranged between 60 and 940 km2, with an average of 242 km2, according to the Thiessen polygons. The lowest average biomass was estimated during the May-June cruise (37.8±3.36 t) while the highest was estimated in December (189.56±11.6 t) (Table 1). The coefficient of variation (CV expressed as a percentage) obtained for the biomass estimated by the four methods was lowest during the months of May-June and July (CV=13.3% and 13.9%, respectively), corresponding to the period before the fishing season, while biomass increased in December and January (CV=20.6% and 20.1%) after closure of the season.
Estimated values of the biomass per cruise±standard error (SE), density±SE and abundance±SE of the Octopus "vulgaris” Type II in the eastern sublittoral zone of Campeche Bank. Biomasses were standardized to a total area of 5000 km2. CV: coefficient of variation of the estimates by cruise of the four methods.
Method
Biomass (t)
± SE
Density (octopus km–2)
± SE
Abundance (N octopus)
± SE
May–June 2016
CV=13.3 %
CV=30.2 %
CV=30.2 %
Swept area
32.57
4.62
7.05
1.00
35256
5006
Stratified
44.66
1.99
3.82
0.17
19084
851
Geostatistic
36.40
2.21
8.09
0.49
40456
2453
Weighted
37.75
4.62
8.21
1.00
41030
5017
July 2016
CV=13.9 %
CV=13.4 %
CV=13.5 %
Swept area
125.99
11.79
27.39
2.56
136960
12817
Stratified
133.97
10.37
28.60
2.21
143002
11066
Geostatistic
95.94
10.40
20.86
2.26
104289
11308
Weighted
125.55
13.74
27.29
2.99
136466
14935
December 2017
CV=20.6 %
CV=20.5 %
CV=20.5 %
Swept area
232.89
18.77
57.52
4.64
287611
23178
Stratified
201.60
38.86
48.48
9.35
242418
46726
Geostatistic
139.34
12.91
34.40
3.19
172002
15933
Weighted
184.42
10.45
45.54
2.58
227679
12906
January 2017
CV=20.1 %
CV=20.4 %
CV=20.3 %
Swept area
58.64
6.88
13.94
1.63
69700
8172
Stratified
73.55
5.25
17.59
1.26
87946
6281
Geostatistic
46.42
9.78
11.04
2.33
55223
11637
Weighted
52.63
5.62
12.53
1.34
62652
6694
Among the four models, the lowest potential biomass was calculated using the geostatistical model, with values 18.5% to 36.7% lower than those obtained using the other methods. Similar patterns of variation were observed for density, but the differences from one month to another were not as evident as they were for biomass. Thus, the highest average density was recorded in December 2016 (46.5 octopuses per km–2), and the lowest average density was observed during the May-June 2016 cruise (6.8 octopuses per km–2). Finally, the weighted biomass method produced the second lowest values of biomass and abundance (Table 1), except during the May-June 2016 cruise.
Analysis of the distribution pattern
The value of the parameters k and p (k2=8, p=0.7) of the negative binomial distribution showed that individuals of O. "vulgaris” Type II tended to aggregate in patches. This species occupied the entire area, with density gradually increasing to a maximum in December and then decreasing to January. The variogram (Fig. 3) shows that as the distance between the sampling points increases, the O. "vulgaris” Type II CPUA (octopus km–2) increases towards shallow waters east of the Campeche Bank. This is due to the distribution of the species (Fig. 4).
Variogram analysis of the semi-variance of the abundances (CPUA; octopus km–2) of O. "vulgaris” Type II in the shallow sublittoral zone in the east of Campeche Bank in relation to the distance between sampling stations.
Abundance distribution (CPUA; octopus km–2) of O. "vulgaris” Type II (dark-green line) and O. maya (orange line) in the northeast of the Campeche Bank. The source of Octopus maya distribution was Avendaño et al. (2019).
Denser aggregations of O. "vulgaris” Type II were recorded towards the southeast of the Yucatán shelf. Cruises that took place before the fishing season (May-June and July 2016) showed the formation of high-density patches to the southeast of the study area (Fig. 4). At the end of the fishing season (December 2016 and January 2017), the octopuses tended to be distributed heterogeneously over the sampled area and, although patches continued to form, they were less extensive (Fig. 4).
Geographical overlay analysis
In order to identify the overlapping geographic areas of O. maya and O. "vulgaris” Type II, a distribution analysis of both species was conducted using data collected during the same period (Fig. 4). The distribution of O. maya was obtained from Avendaño et al. (2019), and the distribution of O. "vulgaris” Type II was determined in the current study. The geographical overlap according to the Schoener index (1968) was relatively low for the May-June and July samplings (D=0.30 and D=0.39, respectively), which corresponded to the period prior to the fishing season. These values increased after the end of the fishing season, with higher values of overlap in December 2016 and January 2017 (D=0.60 and D=0.62, respectively).
An overlap in the distribution of the two species was observed during all the research cruises. This was more evident towards the south of the study area, where medium to low overlap was observed (Di≈0.34) during the months prior to the opening of the fishing season; however, at those sites, O. maya was dominant in the west (approximately 88° W), while O. "vulgaris” Type II was dominant in the east (approximately 87° W). In the months following the fishing season, the heterogeneous distribution of smaller patches of both species caused the degree of overlap to increase throughout the study area (Di≈0.61).
DISCUSSION
A total allowable catch, which is a common measure in cephalopod fisheries, is calculated in the octopus fishery of the Mexican Atlantic coast (Jurado-Molina 2010, Norman et al. 2014). Although this measure is not respected, it is aimed at O. maya (Salas et al. 2008, Beléndez-Moreno et al. 2014, DOF 2016), while O. "vulgaris” Type II can be freely caught. This study analyses for the first time the distribution and relative abundance of O. "vulgaris” Type II on the continental shelf of the Yucatán Peninsula between 30 and 60 m depth. This analysis was carried out through systematic sampling, independent from fishery activities, covering the potential fishing area. Because any precursor of fisheries management measures must be based primarily on the identification of resource distribution (Pierce and Guerra 1994), this analysis provides an instantaneous estimation of relative abundance in that area in four different periods (Rosenberg et al. 1990, Nevárez-Martínez et al. 2000).
The results on the biomass calculated using the four methods were relatively consistent (CV<21%), and discrepancies found between the methods could be due to the differences in the representativeness of the area of influence assumed by each one (Pierce and Guerra 1994, Nevárez-Martínez et al. 2000, Hernández-Flores et al. 2015). For example, in three of the methods (the swept area method, the weighted swept area method and the geostatistical biomass method), the distance between sampling stations is a key element in the calculation. In contrast, in the stratified method, the number of intervals is the key element, so for this method it is necessary to apply a stratification or classification from the start of the analysis (Nevárez-Martínez et al. 2000, Avendaño et al. 2019). The precautionary approach recommends choosing the lowest value calculated from the models to determine the octopus abundance in the area of study for any potential management measure, which in this case was obtained using the geostatistical method.
As expected, our results showed that changes in the biomass and distribution of O. "vulgaris” Type II occurred before and after the fishing season in 2016. The increase in the relative biomass of O. "vulgaris” Type II from June to December (47%) and the change in its distribution pattern in those periods (Fig. 4) indicated that the stock reached its maximum abundance at the end of the fishing season (Table 1). This pattern differed from that observed for O. maya, which reached its highest value in October and gradually decreased to its lowest value in December (Table 2). This could be related to the fishing effort dynamics, as most fishermen target O. maya at the start of the season, but as its abundance decreases, they shift to O. "vulgaris” Type II.
Landings in tons (t) of the red octopus Octopus maya and the common octopus Octopus "vulgaris” Type II captured in the Yucatán Peninsula during the season of 2016. Source: National Institute of Fisheries and Aquaculture. Yucalpeten, Yucatán, Mexico.
Octopus maya
Octopus "vulgaris” Type II
August
4867
949
September
6421
1053
October
6923
704
November
5032
501
December
5134
1865
Total
28377
5072
The commercial landings of O. "vulgaris” during the 2016 season peaked in December (1865 t), which coincided with the highest relative biomass estimated in this study using data collected independently from the fishery. To explain the catch levels of O. "vulgaris” recorded by the commercial fleet, the instantaneous abundance calculated in this study (average 189.6 t) for December 2016 could be 12 or 15 times smaller than the value of the relative abundance in the entire fishing area, meaning that the biomass in that month would be in the order of magnitude of between 2300 and 2800 t.
The propensity for O. "vulgaris” Type II to form aggregations, as observed in this study, infers that the behaviour of the species in this area is similar to that observed in O. "vulgaris” Type III in South Africa (Oosthuizen and Smale 2003) and O. insularis in Brazil (Leite et al. 2009). The higher abundance of O. "vulgaris” Type II towards the southeast zone and the formation of zones of aggregation in the months of May–June and July, as well as the heterogeneous distribution for the months of December and January, could be related to changes in environmental conditions, as suggested for other coastal cephalopods (Pecl and Jackson 2008, Pierce et al. 2008, Gamboa-Álvarez et al. 2015). This hypothesis stems from evidence that the Loop current and the Yucatán current have a strong influence on some other resources in the area (i.e. Epinephelus morio and Octopus maya) (López-Rocha et al. 2013, Gamboa-Álvarez et al. 2015). It is also possible that O. "vulgaris” Type II paralarvae are recruited in the eastern part of the Yucatán shelf due to the oceanographic conditions in that area (i.e. upwelling and nutrients [Enríquez et al. 2010]).
In the current study the northeast of the Yucatán Peninsula showed a geographical overlap of both species, ranging from low to medium overlap; this area was previously considered exclusively as an O. "vulgaris” Type II area. It is possible that this speciesis in fact more abundant in deeper waters, since it has not been recorded in the landings of artisanal fishermen. However, O. maya can occupy deeper waters than previously thought (Avendaño et al. 2019). Our results show that there are no clear limits on the distribution of the two species in waters at depths greater than 30 m, with both species being able to occupy the same space. The dominance of one or other of the species will depend on intrinsic biological characteristics and the prevailing ecological conditions, as well as on competition for space (e.g. success in recruitment or faster growth) (Leite et al. 2009, Pierce et al. 2008, Rodhouse et al. 2014). This scenario is consistent with the fishers’ and researchers’ idea that O. maya and O. "vulgaris” Type II occupy different habitats on the continental shelf, according to depth (Salas et al. 2008, DOF 2016). Thus, most abundant fishing grounds for O. maya were located along the coast, decreasing towards deeper waters (30 m depth at approximately 88°W) and increasing towards the west of the peninsula (Solís-Ramírez and Chávez 1986, Gamboa-Álvarez et al. 2015, Avendaño et al. 2019).
In conclusion, our results showed the changes in the distribution and abundance patterns experienced by O. maya and O. "vulgaris” Type II over a short period of six months. The results suggest that O. "vulgaris” Type II occupied the entire sampled area and formed patches, especially during December and January. Moreover, our findings revealed a geographic overlap of both species, whose amplitude changed throughout the year according to the geographic position: O. maya dominated at approximately 88°W, while O. vulgaris Type II dominated towards the southeast at 87°W.
Although in general the total abundance of both species was similar in the entire area of study, Octopus "vulgaris” Type II showed higher abundance at depths of less than 40 m. Consequently, it would be necessary to calculate the maximum sustainable yield of this species to develop appropriate management measures. The two species must be managed jointly as they overlap in great part of the fishing grounds. We also propose that it would be advisable to identify and create protected areas, which would provide refuge for both species at depths greater than 40 m. These areas would protect O. maya and O. "vulgaris” Type II adults, resulting in improved reproductive success, as mentioned by Velázquez-Abunader et al. (2013) and Avendaño et al. (2019).
ACKNOWLEDGEMENTS
The authors thank the Federación de Sociedades Cooperativas de la Industria Pesquera del Centro y Oeste del Yucatán, SC de RL for logistical support and providing the boat to carry out the cruises. We also thank Miguel Ángel Cabrera, Jesús Miguel Soto Vázquez and Luis Enrique Ángeles-González for their support in the field work.
Financial support
The results of this study were obtained as part of the project "Distribución, reproducción, biomasa y patrones de movimiento del pulpo común Octopus vulgaris Cuvier, 1797 en la costa de Yucatán”, with financial support from the Consejo Nacional de Ciencia y Tecnología (grant number 237057).
REFERENCESAfkhamiM.E.McIntyreP.J.StraussS.Y.Mutualist-mediated effects on species' range limits across large geographic scales20141712651273https://doi.org/10.1111/ele.12332AmorM.D.NormanM.D.RouraA.Morphological assessment of the Octopus vulgarisspecies complex evaluated in light of molecular-based phylogenetic inferences201746275288https://doi.org/10.1111/zsc.12207AmorM.D.DoyleS.R.NormanM.D.Genome-wide sequencing uncovers cryptic diversity and mito-nuclear discordance in the Octopus vulgaris species complex2019573493https://doi.org/10.1101/573493AvendañoO.Velázquez-AbunaderI.Fernández-JardónC.Biomass and distribution of the red octopus (Octopus maya) in the north-east of the Campeche Bank20199913171323https://doi.org/10.1017/S0025315419000419Ávila-PovedaO.H.KouetaN.Benítez-VillalobosF.. Reproductive traits of Octopus maya (Cephalopoda: Octopoda) with implications for fisheries management2016362944https://doi.org/10.1080/13235818.2015.1072912Beléndez-MorenoL.F.J.Espino-BarrE.Galindo-CortesG.2014MexicoSecretaría de Agricultura, Ganadería. Desarrollo Rural, Pesca y Alimentación (SAGARPA)BrasselK.E.ReifD.A procedure to generate Thiessen polygons197911289303https://doi.org/10.1111/j.1538-4632.1979.tb00695.xCochranW.G.1978MexicoCompañía Editorial Continental513CornellH.Niche OverlapHastingsA.GrossL.J.2011California, USAUniversity of California Press489497CressieN.Statistics for spatial data19924613617https://doi.org/10.1111/j.1365-3121.1992.tb00605.xDiario Oficial de la Federación (DOF)2016Ciudad de MéxicoSAGARPAApril 13, 2016DubrannaJ.Pérez-BruniusP.LópezM.Circulation over the continental shelf of the western and southwestern Gulf of Mexico2011116C08009https://doi.org/10.1029/2011JC007007EmeryT.J.HartmannK.GardnerC.Management issues and options for small scale holobenthic octopus fisheries2016120180188https://doi.org/10.1016/j.ocecoaman.2015.12.004EnríquezC.Mariño-TapiaI.JHerrera-SilveiraJ.A.Dispersion in the Yucatan coastal zone: Implications for red tide events201030127137https://doi.org/10.1016/j.csr.2009.10.005Galindo-CortesG.Hernández-FloresÁ.Santos-ValenciaJ.Pulpo del Golfo de México Octopus maya y Octopus vulgarisBeléndez-MorenoL.F.J.Espino-BarrE.2014INAPESCA177207Gamboa-ÁlvarezM.Á.López-RochaJ.A.Poot-LópezG.R.Spatial analysis of the abundance and catchability of the red octopus Octopus maya (Voss and Solis-Ramírez, 1966) on the continental shelf of the Yucatán peninsula, México201534481492https://doi.org/10.2983/035.034.0232GuerraA.Spatial distribution pattern of Octopus vulgarisCuvier1981195133146https://doi.org/10.1111/j.1469-7998.1981.tb01897.xHernández-FloresA.CondalA.Poot-SalazarA.Geostatistical analysis and spatial modeling of population density for the sea cucumbers Isostichopus badionotus and Holothuria floridana on the Yucatan Peninsula, Mexico2015172114124https://doi.org/10.1016/j.fishres.2015.07.005HilbornR.WaltersC.1992New YorkChapman-Hall570Jurado-MolinaJ.A Bayesian framework with implementation error to improve the management of the red octopus (Octopus maya) fishery off the Yucatán Peninsula201036114https://doi.org/10.7773/cm.v36i1.1627LeeT.N.WilliamsE.Mean distribution and seasonal variability of coastal currents and temperature in the Florida Keys with implications on larval recruitment1999643556LeiteT.S.HaimoviciM.MatherJ.Habitat, distribution, and abundance of the commercial octopus (Octopus insularis) in a tropical oceanic island, Brazil: Information for management of an artisanal fishery inside a marine protected area2009988591https://doi.org/10.1016/j.fishres.2009.04.001LeporatiS.C.HartA.M.LarsenR.Octopus life history relative to age, in a multi-geared developmental fishery20151652841https://doi.org/10.1016/j.fishres.2014.12.017LimaF.D.Berbel-FilhoW.M.LeiteT.S.Occurrence of Octopus insularis Leite and Haimovici, 2008 in the Tropical Northwestern Atlantic and implications of species misidentification to octopus fisheries management201747723734https://doi.org/10.1007/s12526-017-0638-yLópez-RochaJ.A.Arreguín-SánchezF.Spatial dynamics of the red grouper Epinephelus morio (Pisces: Serranidae) on the Campeche Bank, Gulf of Mexico201377313322https://doi.org/10.3989/scimar.03565.13BMoreyS.L.Zavala-HidalgoJ.O'BrienJ.J.The seasonal variability of continental shelf circulation in the northern and western Gulf of Mexico from a high-resolution numerical modelSturgesW.Lugo-FernandezA.2006Washington, DC. USAGeophysical Monograph Series203218https://doi.org/10.1029/161GM16Nevárez-Martı́nezM.O.Hernández-HerreraA.Morales-BojórquezE.Biomass and distribution of the jumbo squid (Dosidicus gigas; d'Orbigny, 1835) in the Gulf of California, Mexico200049129140https://doi.org/10.1016/S0165-7836(00)00198-3NormanM.D.FinnJ.K.HochbergF.G.Family OctopodidaeJerebP.RoperC.F.E.3Octopods and vampire squidsFAO Species Catalogue for Fishery Puroises42014RomeFAO36215http://www.fao.org/3/a-i3489e.pdfOosthuizenA.SmaleM.J.Population biology of Octopus vulgaris on the temperate south-eastern coast of South Africa200383535541https://doi.org/10.1017/S0025315403007458hParisC.B.CowenR.K.ClaroR.Larval transport pathways from Cuban spawning aggregations (Snappers; Lutjanidae) based on biophysical modeling200529693106https://doi.org/10.3354/meps296093PaulyD.ICLARM Stud. Rev.19848325PeclG.T.JacksonG.D.The potential impacts of climate change on inshore squid: biology, ecology and fisheries200818373385https://doi.org/10.1007/s11160-007-9077-3PierceG.J.GuerraA.Stock assessment methods used for cephalopod fisheries199421255285https://doi.org/10.1016/0165-7836(94)90108-2PierceG.J.ValavanisV.D.GuerraA.A review of cephalopod-environment interactions in European Seas and other world areas20086124970https://doi.org/10.1007/s10750-008-9489-7RitschardE.A.Guerrero-KommritzJ.SanchezJ.A.First molecular approach to the octopus fauna from the southern Caribbean20197e7300https://doi.org/10.7717/peerj.7300RivoirardJ.SimmondsJ.FooteK.G.2008LondonBlackwell Science206RodhouseP.G.K.PierceG.J.NicholsO.C.Environmental effects on cephalopod population dynamics: implications for management of fisheriesVidalE.A.G.Adv. Mar. Biol.20146799233https://doi.org/10.1016/B978-0-12-800287-2.00002-0RosenbergA.KirkwoodG.CrombieJ.The assessment of stocks of annual squid species19908335350https://doi.org/10.1016/0165-7836(90)90003-ESalasC.CabreraM.PalomoL.2008Merida, MexicoCinvestav IPNSalasS.Torres-IrineoE.CoronadoE.Towards a métier-based assessment and management approach for mixed fisheries in Southeastern Mexico2019103148159https://doi.org/10.1016/j.marpol.2019.02.040SauerW.H.GleadallI.G.Downey-BreedtN.World Octopus Fisheries2020https://doi.org/10.1080/23308249.2019.1680603SawatzkyD.RainesG.Bonham-CarterG.2009http://arcscripts.esri.com/details.asp?dbid=15341SchoenerT.W.The Anolis lizards of Bimini: Resource partitioning in a complex fauna196849704726https://doi.org/10.2307/1935534Solís-RamírezM.J.1994MexicoSecretaria de Pesca, INPCENEDICUniversidad de ColimaCD MultimediaSolís-RamírezM.ChávezE.Evaluación y régimen óptimo de pesca del pulpo. de la península de Yucatán198613118TesterP.A.StumpfR.P.StumpfF.M.An expatriate red tide bloom: transport, distribution, and persistence19913610531061https://doi.org/10.4319/lo.1991.36.5.1053Van NieuwenhoveA.H.M.RatsimbazafyH.A.KochziusM.Cryptic diversity and limited connectivity in octopuses: Recommendations for fisheries management201914e0214748https://doi.org/10.1371/journal.pone.0214748Velázquez-AbunaderI.SalasS.CabreraM.A.Differential catchability by zone, fleet, and size: the case of the red octopus (Octopus maya) and common octopus (Octopus vulgaris) fishery in Yucatan, Mexico201332845854https://doi.org/10.2983/035.032.0328WebsterR.OliverM.A.2007Chichester. EnglandJWShttps://doi.org/10.1002/9780470517277ZarJ.H.1999Upper Saddle River, New Jersey. USAPrentice HallZavala-HidalgoJ.Parés-SierraA.OchoaJ.Seasonal variability of the temperature and heat fluxes in the Gulf of Mexico20021581104