Qinqin Lin, Yuying Zhang, Jiangfeng Zhu. Simulating the impacts of fishing on central and eastern tropical Pacific ecosystem using multispecies size-spectrum model[J]. Acta Oceanologica Sinica, 2022, 41(3): 34-43. doi: 10.1007/s13131-021-1902-3
Citation:
CUI Wei, YANG Jungang, MA Yi. A statistical analysis of mesoscale eddies in the Bay of Bengal from 22-year altimetry data[J]. Acta Oceanologica Sinica, 2016, 35(11): 16-27. doi: 10.1007/s13131-016-0945-3
Qinqin Lin, Yuying Zhang, Jiangfeng Zhu. Simulating the impacts of fishing on central and eastern tropical Pacific ecosystem using multispecies size-spectrum model[J]. Acta Oceanologica Sinica, 2022, 41(3): 34-43. doi: 10.1007/s13131-021-1902-3
Citation:
CUI Wei, YANG Jungang, MA Yi. A statistical analysis of mesoscale eddies in the Bay of Bengal from 22-year altimetry data[J]. Acta Oceanologica Sinica, 2016, 35(11): 16-27. doi: 10.1007/s13131-016-0945-3
Eddy properties in the Bay of Bengal are studied from 22 a archiving, validation and interpretation of satellite oceanographic (AVISO) data using a sea level anomaly (SLA)-based eddy identification. A geographical distribution and an eddy polarity, an eddy lifetime and propagation distances, eddy origins and terminations, eddy propagation directions and trajectories, eddy kinetic properties, the evolution of eddy properties, seasonal and interannual variabilities of eddy activities are analyzed in this area. Eddies exist principally in the western Bay of Bengal and most of them propagate westward. The polarity distribution of eddies shows cyclones prefer to occur in the northwest and south of the Bay of Bengal, while anticyclones mainly occur in the east of the bay. Five hundred and sixty-five cyclones and 389 anticyclones with the lifetime that exceeds 30 d are detected during the 22 a period, and there is a preference for the cyclones for all lifetime and propagation distances. The kinetic properties of all observed eddies show the average amplitude of the cyclones is larger than that of the anticyclones, whereas that is opposite for average radius, and their average velocities are basically the same. Moreover, the evolution of eddies properties reveals that the eddies with a long lifetime that exceeds 90 d have a significant double-stage feature of the former 50 d growth period and the dying period after 50 d. For the seasonal variability of the eddies, the cyclones occur more often in spring while the anticyclones occur more often in summer. The analysis of long-lived eddy seasonal distributions shows that there is the obvious seasonal variation of the eddy activities in the Bay of Bengal. The interannual variability of an eddy number shows an obvious negative correlation with the EKE variation.
Ecosystem modeling has been increasingly recognized as an essential part of fisheries management (Food and Agricultural Organization of the United Nations, 2001; Jacobsen et al., 2015; Bauer et al., 2019). There have been rapid development and implementation of ecosystem-based fisheries management (EBFM) over the past decades. Yet, the projection of a given fishing scenario varies with the choice of an ecosystem model (Collie et al., 2016). Also, marine ecosystems are often complicated because of the various life histories of species, as well as the direct and indirect interactions among species. One key characteristic of EBFM is to take ecosystem processes into account when developing fisheries management measures, including changes in fishing strategies and species interactions (Sissenwine and Murawski, 2004). Recently, size-spectrum models have become a promising tool for modeling fish communities (Jacobsen et al., 2015; Giacomini et al., 2016). It is a type of dynamic model based on individual-level biological processes, convenient to address fisheries management issues, such as the impacts of fishing on ecosystem structure and the effectiveness of management decisions (Polovina and Woodworth-Jefcoats, 2013; Jacobsen et al., 2014; Andersen et al., 2015). The size-spectrum model was first developed by Andersen and Beyer (2006), including three versions of increasing complexity: the community model, the trait-based model, and the multispecies model. They are a subset of physiologically structured models (Scott et al., 2014) and can be implemented to represent a realistic fish community. The application of the multispecies size-spectrum model in fishery management is that it provides a framework to evaluate various management measures and makes trade-offs. For example, Blanchard et al. (2014) evaluated the current and several alternative management strategies for exploited North Sea stocks, and concluded that reductions in fishing mortality would be needed to sustain biodiversity.
The Pacific Ocean is rich in fish resources and supports many productive pelagic fisheries, but it has experienced a long exploitation history. Tunas and billfishes are predominantly caught by longline, one of the most widespread types of fishing gear in the area (Sibert et al., 2006). Fishing activities not only directly affect the species’ abundance (Beverton and Holt, 1957; Jennings et al., 1999), but also lead to indirect impacts (Gislason and Rice, 1998; Law, 2000; Dulvy et al., 2004), such as the changes of species’ life history. Consequently, with the overexploitation of many ecologically and economically important fish populations, the Pacific Ocean ecosystems have been degraded (Pauly et al., 1998; Polovina et al., 2009). In addition to technological and operational developments, fishery management measures are arguably the most critical non-technical factors shaping the longline fishery (Shelley et al., 2014). To date, most fishery management tends to set limits for target species while neglecting the impacts of species interactions. Therefore, EBFM is increasingly applied in Tuna Regional Fisheries Management Organizations (tRFMOs) to develop strategies applicable to global tuna fisheries (Miyake et al., 2010).
For studies on coastal marine ecosystems, scientific survey data are usually available. However, previous ecosystem studies of the Pacific Ocean were primarily based on data from adjacent areas (Gerrodette et al., 2012; Zhu et al., 2012). Besides, much research on pelagic ecosystems remains at the theoretical stage, with few case studies of practical application in operational management. Zhang et al. (2016) developed a multispecies size-spectrum model for the Haizhou Bay, China, to study the ecosystem dynamics. They confirmed that multispecies size-spectrum models could be well applied to data-poor ecosystems. In this analysis, we used a multispecies size-spectrum framework (Hartvig et al., 2011) in assessing the dynamics of populations and ecosystem responding to fishing in the central and eastern tropical Pacific Ocean based on Chinese tuna longline observer data, the latest Western and Central Pacific Fisheries Commission (WCPFC) stock assessment reports (https://www.wcpfc.int/meetings/16th-regular-session-scientific-committee) and published data, such as FishBase (http://www.fishbase.org/). Twenty species were simulated in the size-spectrum model based on the simplified food web constructed in this area (Lin and Zhu, 2020), including common large-sized species with significant economic and social importance caught by longline. The main aim of this study is to (1) examine the direct and indirect fishing effects on community dynamics in the EBFM framework using common ecological indicators, and (2) project the biomass changes of nine main catch species in the longline fishery under five different fishing levels. This study informs appropriate management strategies regarding economic species in this area and contributes to the sustainable exploitation of fishery resources of WCPFC.
2.
Materials and methods
2.1
Model structure
According to the simplified food web constructed in the central and eastern tropical Pacific food web (Lin and Zhu, 2020) in terms of topological importance analysis and species’ economic value, we selected twenty species (Table 1) and projected their population dynamics under different fishing levels. These twenty species accounted for over 90% of the total catch of all species collected by China’s tuna longline observers in the central and eastern tropical Pacific from June to November, 2017 (Fig. 1). Trip information was listed in Table S1, and more details are available in China’s annual report provided to the WCPFC Commission (Dai et al., 2017). In addition to the relatively higher abundances in the ecosystem, these species play critical roles in maintaining food web structure and stability and reflecting population dynamics in a multispecies size-structured community (Lin and Zhu, 2020). Nine large-sized species are our research priority, not only because of the data availability but also the primary interests of longline fishery in the tropical Pacific Ocean, including Yellowfin tuna (Thunnus albacares), Bigeye tuna (Thunnus obesus), Swordfish (Xiphias gladius), Blue marlin (Makaira mazara), Albacore tuna (Thunnus alalunga), Silky shark (Carcharhinus falciformis), Blue shark (Prionace glauca), Skipjack (Katsuwonus pelamis ) and Striped marlin (Kajikia audax). If measurable, the size, sex, and other biological information of fish were recorded by on-board observers. Also, the stomach samples were sent to the laboratory of Shanghai Ocean University for further analysis. Information on stomach content samples, including length (mm), weight (g), and quantity of prey items, were measured. The detailed feeding relationships among the nine species in the ecosystem have been introduced in Feng et al. (2019) and are not repeated here.
Table
1.
Species-specific input parameters used in the size-spectrum model
Figure
1.
Location of Chinese longline fishery in the central and eastern tropical Pacific from June to November in 2017. Each dot represents a position where a longline set was deployed.
In the multispecies size-spectrum model, each species shows different change patterns driven by specific life history characteristics, including growth, predation, mortality, and reproduction (Andersen and Beyer, 2006; Hartvig et al., 2011; Scott et al., 2014). A governing principle in the size-spectrum model is that the available energy gets from predation by larger species on smaller ones and used for growth and reproduction can be linked to the asymptotic size of the predator (Ursin, 1973; Andersen et al., 2015).
The model framework builds on two generally accepted ecological assumptions. The first is that an essential trait of a species, i, is its body weight, w, and the model can be formulated with several general parameters as all the other parameters are relevant to the body weight. The population dynamics of a given species that are driven by growth and mortality can be formulated by the McKendrick-von Foerster Conservation Equation (M’Kendrick, 1925; von Foerster, 1959):
where gi(w) is the individual growth at weight w, and μi(w) is the total mortality of species i, t is time. The size spectrum Ni(w) of a particular size class in the model can be measured in the unit of numbers per weight per volume (Andersen and Pedersen, 2010). The second assumption is that the food preference (ϕ) depends on the species and the individual weight, described by the log-normal selection model (Ursin, 1973):
where w and wp are the bodyweights of the predator and prey, respectively. βi is the preferred predator-prey mass ratio, and σi is the width of the size selection function. Details of other equations (Hartvig et al., 2011; Andersen et al., 2016b) are listed in Table 2. A set of default parameters is listed in Table S2, and the interaction matrix (θ) estimated based on species’ vertical overlap is included in Table S3 (Woodworth-Jefcoats et al., 2019).
Note: $\kappa $ in Eq. (4) is the abbreviation of carrying capacity; p in Eq. (7) and Eq. (8) represents the exponent of standard metabolism and the default value in the model is 0.75; $\epsilon $ in Eq. (11) is the efficiency of offspring production; SS in Eq. (14) is the selectivity size; r0 in Eq. (15) is the coefficient of regeneration rate of resources, and the default values in the model are 0.1 and 4, respectively; $\text{λ} $ in the Eq. (16) is the exponent of resource spectrum following the equation 2−n+q, where n is the exponent for maximum food intake (2/3) and q is the exponent for volumetric search rate (0.8). Other default parameters are listed in Table S2.
The encounter rate of food Ei(w) (Eq. (3)) among individuals is based on the “Andersen-Ursin” encounter model (Andersen and Ursin, 1977). Both small-sized species and zooplankton which is part of the background resource can be considered as available food. The combination of food preference, volumetric search rate Vi(w) (Eq. (4)), and prey abundance determine the food encounter rate. The consumption of an individual after predation is related to the feeding level, fi(w) (Eq. (5)), ranging from 0 to 1, and the maximum consumption rate Imax, i(w) (Eq. (6)).
The energy obtained by effective assimilation of ingested food is used to meet the needs of standard metabolism, partly for somatic growth gi(w) (Eq. (7)) of juveniles and partly for reproduction and maturation $\psi $i(w) (Eq. (8)) of mature individuals. Species’ growth and reproduction will stop if the food supply fails to support the metabolic requirements. And the allocation proportion between somatic growth and reproduction is determined by the maturation function $\psi $i(w) (Eq. (9)) of individuals. The recruitment Ri (Eq. (10)) described by a Beverton-Holt type function (Beverton and Holt, 1957) is dependent on the amount of energy used for reproduction Rep (Eq. (11), calculated as the number of eggs); and the maximum recruitment Rmax, i (numbers of surviving recruits every year) as a key parameter essentially decides the total biomass of a population. The mortality of an individual has three categories: predation mortality μp, i(w) (Eq. (12)), background mortality μb, i(w) (Eq. (13)), and fishing mortality, Fi(w) (Eq. (14)), which ensures a mass balance in the size-spectrum model. Fishing was imposed with the default knife-edge selectivity in the size-spectrum framework. The catchability (Qi) was set to a constant value of 1 in the simulation.
The resource spectrum Nr(w) (Eq. (15)) acts as the food items of the smallest individuals, which prey on zooplankton or benthic production. The resource dynamics can be modeled by a semi-chemostat growth equation combined with a carrying capacity $\kappa(w) $ (Eq. (16)) defined as the control parameter for the resource productivity in the community.
The size-spectrum model can be implemented by the R package “mizer” (Scott et al., 2014).
2.2
Model parameterization and validation
Each of the twenty species is characterized by a range of parameters reflecting its species-specific physiological features and life history (Table 1). K is the growth coefficient; a and b are the parameters in weight-length relationships. Winf and Wmat represent the asymptotic body weight and mature weight of species, respectively. Species’ maximum size obtained from FishBase was converted to Winf following the standard exponential equation W =aLb. The above five growth parameters were publicly available from FishBase. Following previous studies, β, the preferred predator-prey mass ratio was assigned 100 for species fed mostly on fish, 400 for shark species, and 1000 for species with strong dependencies on zooplanktons (Blanchard et al., 2014; Reum et al., 2019). The initial Rmax was estimated from an empirical equation as 1011×$W_{\rm{inf}}^{-1.5} $ (Blanchard et al., 2014). The Rmax and default carrying capacity of background resources ($\kappa $) were then calibrated by minimizing the Root Mean Square Error (RMSE) between the log-transformed modelled and estimated spawning stock biomass (SSB) from the WCPFC single stock assessments of the nine large-sized species from 1993 to 2012 (20 years’ fishing period). RMSE is a frequently used measure of the differences between values predicted by a model and the values observed. According to Chai and Draxler (2014), the calculation of RMSE was expressed as:
The time-series fishing mortality (1952−2015) used for model validation were also derived from the most recent WCPFC stock assessment reports. Many of these species are non-target species in the longline fishery that are lacking fishing mortality data. The fishing mortality of bycatch species has increased about fourfold, owing to increases in effort in tuna fisheries in the Pacific Ocean (Shelley et al., 2014). Thus, we assumed that fish mortality of these species increases linearly (from 0 to 0.4) at this time period. The selectivity size (SS) parameters of each species are different. The size position susceptible to fishing mortality of each species was determined by the asymptotic size following the equation as Winf×0.05 (Scott et al., 2013). Cephalopods, Scomber, and Puffer are not sorted to species level, and the parameter estimates were averaged from several representative species. Besides, parameters of some species were obtained from relevant studies when biological data were not available in FishBase. Detailed sources were noted in Table 1.
The constructed size-spectrum model was first to run for 200 a before 1952 without fishing mortality for all species to get the virgin resources before the fishing efforts were introduced. The initial population abundances play a decisive role in determining the population dynamics in a community. We also run the model under historical fishing mortality from 1952 to 2015 to reach the current state of the community. In addition, after the calibration of Rmax and $\kappa $, we compared the simulated growth curves with empirical size-at-age data to ensure that the calibrated model could produce plausible stock dynamics. R-squared (R2) ranging from 0 to 1 was used to indicate the statistical measure of fit between two sets of data. The changing trend between predictions and observations can better demonstrate the goodness of fit than a specific value matched in a given time period (Jacobsen et al., 2017). Most species in the size-spectrum model began to be exploited in the 1950s. Some species started being fished in the 1970s, and the fishing mortality was set to 0 in the 1950s to 1970s, including Blue shark, Silky shark, Blue marlin, and Skipjack.
2.3
Scenario projections and analysis of model outputs
We used the calibrated model to predict the dynamic changes of nine species in the community from 2016 to 2050 under five fishing scenarios:
(1) unexploited Scenario: F= 0;
(2) recent F Scenario: F=Frecent, the average fishing mortality from 2013 to 2015;
(3) low F Scenario: F=0.5Frecent;
(4) medium F Scenario: F=2Frecent;
(5) high F Scenario: F=3Frecent.
The recent fishing mortality (Frecent) for nine species obtained from WCPFC stock assessment outputs was listed in the Table S4. These scenarios were chosen partly based on the fishing history of each species in the Pacific Ocean and can help us understand the fishery status of the target species relative to the unfished ecosystem. The total biomass, mean body weight (MeanW), and the slope of community size spectra (Slope) were used as ecological indicators to reflect the changes in species abundance and ecosystem structure.
3.
Results
3.1
Model Validation
Following previous studies, we calibrated the model by tuning two influential parameters Rmax and kappa, which are essential for setting species’ initial size spectra and determines the background resource in the system. The validation results were demonstrated in Fig. 2. The model produced strong correlations between the size-at-age data from our simulation and the WCPFC reports, with the R2 of nine species more than 0.9. Therefore, the model validation results indicated that the multispecies size-spectrum model could reflect the dynamics of a fish community in the tropical Pacific.
Figure
2.
Model validation results. Comparison of species’ size-at-age data modeled in this study with the outputs in the WCPFC stock assessment reports. The closer R2 is to 1, the better the fit is.
The multispecies size-spectrum model projects the dynamics of individual species. As shown in Fig. 3, all species showed different patterns driven by specific parameters, including biomass (Fig. 3a), predation mortality (Fig. 3b), and biomass spectra (Fig. 3c). It can be seen that predation mortality increased firstly and then varies with the increasing body size. There are two “bumps” of the predation mortality curve due to the trophic cascade effect. The enormous individuals in the community do not experience predation by other species. Species biomass fluctuated continuously. The economic species showed a slight downward trend under the historical fishing mortality, except that the biomass of Blue shark began to increase since 1990. The biomass of several species tended to reach a balance over time, such as Bigeye tuna and Blue marlin. The biomass density of large-sized species declined dramatically following historical fishing (Fig. 3c).
Figure
3.
The outputs of the multispecies size-spectrum model when historical fishing efforts from 1952 to 2015 were applied. Each line presents one of the species in the community. The biomass (a) change through time was demonstrated, and predation mortality (b) and biomass density (c) were the reflection of the last fishing year (2015).
Figure 4 demonstrated the projected biomass change of nine large-sized species caught by longline fishery from 2016 to 2050. The simulation results indicated that each species responded differently to fishing activities, and the decline in biomass was not proportional to the increase in fishing mortality. The biomass of nine species showed an upward trend over the projection period without fishing. The biomass of the main target tuna species in the community was still abundant by 2050 under the recent fishing level, except Albacore tuna. When the fishing mortality was increased to two or three times of Frecent, the biomass of these species decreased dramatically in a few years. We also found that the biomass of Silky shark remained at relatively low levels by 2050 when the fishing mortality was halved (0.5Frecent).
Figure
4.
Projected biomass of nine large-sized species under different fishing scenarios: F=0 (red), F=Frecent (yellow), F=0.5Frecent (green), F=2Frecent (blue), F=3Frecent (purple).
Figure 5 presents the changes of ecological indicators under five fishing scenarios from 2016 to 2050. The dashed lines (F=0) denote an unexploited community, where the MeanW and Slope are both stable. The direct effect of fishing on the community structure was the decrease of the abundance of large-sized individuals, which could be reflected in the declining MeanW and steepening of the size spectrum. Under the 2Frecent Scenario, species’ mean weight decreased about 1 500 g (21%) by 2050 compared with that under the Frecent Scenario. The community Slope indicated a more obvious downward trend with the increasing fishing mortality.
Figure
5.
Ecological indicators respond to different fishing levels over the projection time period (2016−2050), including mean body weight (MeanW in g) and slope of community size spectra (Slope). Dashed line, solid line, dot-dash line, dotted line, and longdash line represent fishing scenarios of F=0, F=Frecent, F=0.5Frecent, F=2Frecent, and F=3Frecent, respectively.
We developed a multispecies size-spectrum framework to analyze and predict the dynamics of species in the central and eastern tropical Pacific ecosystem. The methods provide a convenient way to explore the responses of the community to different fishing levels (Scott et al., 2013). The size-spectrum model has the advantage of maintaining ecological realism compared with traditional unstructured models (Andersen et al., 2016a; Zhang et al., 2018). Most parameters in the model can be obtained from a cross-species analysis (Hartvig et al., 2011) and the metabolic theory (Brown et al., 2004). The model is based on individual-level ecological processes. In Fig. 3, the biomass and mortality of each species responded differently to historical fishing mortality; the situation was mainly determined by size structure. The simulation results also confirmed the importance of species’ specific life history and their interactions with other species in sustaining ecosystem structure. Many species in our study are non-tuna bycatch species retained as byproducts in the commercial longline fishery. However, stock assessments have only been undertaken on very few stocks of non-tuna species, and the absence of information impacting the ability to make informed management decisions.
The tropical Pacific Ocean is large enough to ignore the migration effects for most species. Meanwhile, there are sufficient spatial overlaps among species to ensure regular species interactions, as shown in Table S3. Many previous studies on the Pacific Ocean were primarily based on data from adjacent areas due to the difficulties of obtaining samples in open-ocean (Worm et al., 2005; Gerrodette et al., 2012; Zhu et al., 2012). Thus, the acquisition of the Chinese tuna longline observer data can prompt the reliability of our study. However, our results may only reflect part of pelagic ecosystem dynamics because most of our species are top predators caught by pelagic longlines. In addition, we validated the SSBs and size-at-age data of nine species, but that may not guarantee a reliable projection of the whole ecosystem (Andersen et al., 2016b).
Another limitation of this study is that we only considered the effects of fishing on community structure. Climate-driven pressures such as changes in water temperature and ocean circulation may also threaten the survival of species, impact species distributions and change the community structure (Genner et al., 2010; Simpson et al., 2011). Secondly, we used the MeanW, Slope, and total biomass to reflect the responses of the community to different fishing levels. The study of Houle et al. (2012) revealed that length-based indicators weighted by biomass were more sensitive to fishing pressure. However, these indicators we analyzed may not reflect the population dynamics and ecosystem structure in all respects. The dynamic changes of ecosystem structure leading by fishing are nonlinear in the size-spectrum model. The response of community indicators is usually more sensitive at a smaller fishing effort level (Houle et al., 2012; Zhang et al., 2016). Therefore, as shown in Fig. 4, the species’ biomass showed a disproportionate changing pattern partly due to the different sensitivities to fishing mortality and indirect trophic effects in the community. This size-spectrum model would generate more complex results since a change in fishing effort generally does not result in a proportional change in relevant community indicators. Thus, more attempts to explore the sensitivity of different indicators in fishing are encouraged (Rochet and Trenkel, 2003; Houle et al., 2012).
The simulation results (Fig. 5) in this study are consistent with the conjecture that continuing fishing activities leads to a steepening of the spectrum (Rice and Gislason, 1996; Shin and Cury, 2004; Shin et al., 2005). A reason for this consistency is that fishing mainly targets large-sized individuals. Populations with small asymptotic sizes have larger abundances than those with large asymptotic sizes throughout the life history (Scott et al., 2014). We did not specify the species’ catchability. The default value is set to 1 in the simulation, which means that the fully-selected individuals of a species would be caught efficiently under fishing. Fishing reduces the abundance of small fishes, but compared with natural predation mortality, the fishing effects are relatively low (Andersen and Pedersen, 2010). Also, the reduction of prey resources affects the food availability of their predators; thus, large individuals would have a slower growth rate. In other words, the same fishing effort level had a limited impact on smaller-sized species compared to large-sized predators in our study. Combined with the projection results in Fig. 4, species will not become extinct under long-term fishing only by adopting the most practical management measures. That is to say, there must be sufficient large-sized individuals within a population. Therefore, the recent fishing efforts of Blue marlin and Silky shark need to be adjusted according to the evaluation results and actual situation. The study of Schindler and Hilborn (2015) also revealed that the size-spectrum model is more suitable for short-term predictions, which can be considered as a reference when deciding management strategies.
The fishing impacts triggered by the removal of top predators can also be described as a trophic cascade (Andersen and Pedersen, 2010). A trophic cascade is a sign of indirect impacts among species and can occur in various marine ecosystems when the species abundance at different sizes change due to fishing (Myers et al., 2007; Casini et al., 2008). The phenomenon can be explored by comparing the relative abundance changes of all individuals at different size ranges of F=Frecent, F=3Frecent, and unexploited community (Fig. 6). Compared with the recent fishing scenario, the abundance of large-sized species decreased sharply under the high fishing scenario, consistent with the results shown in Fig. 5. The constant F=0 prediction was used as the baseline, which is also a fundamental assumption proposed by Sheldon et al. (1972) about species size distribution in the marine ecosystem. The effects of fishing on relative abundance vary with the increasing individual size. The fishing pressure reduces the abundance of large-sized individuals. The abundance of smaller individuals will increase due to the release from predation by their predators, which is accordant with the results in the North Sea fish community (Daan et al., 2005). Similarly, this increases the predation mortality on their prey, thus reducing the abundances of the prey species. Subsequently, the fishery will begin to target medium-sized individuals, which will lead to changes in ecosystem structure.
Figure
6.
The relative abundance of different size individuals in the size-spectrum model when F=0 (dashed line), F=Frecent (solid line), and F=3Frecent (red line).
This study modelled a realistic fish community in the central and eastern tropical Pacific Ocean under different fishing scenarios. The response of nine ecologically and economically important species suggested that only the biomass of Bigeye tuna, Swordfish, Yellowfin tuna, and Skipjack could maintain a stable abundance under the recent fishing mortality. Current management strategies are insufficient to achieve an optimal balance between conservation objectives and exploitation of fishery resources, especially for sharks and billfishes. Ecosystem-based fisheries management is an effective way to minimize potential impacts on the ecosystem while realizing sustainable utilization of fishery resources. In an ecosystem context, it is needed to explore the community dynamics under simultaneous stressors such as fishing and environmental changes and consider various complementary indicators in the simulation analysis. In addition, biological research and rigorous stock assessments of non-tuna species captured in the longline fishery should be a high research priority. Thus, relevant studies can better inform and advise fisheries management.
Babu M T, Kumar P S, Rao D P. 1991. A subsurface cyclonic eddy in the Bay of Bengal. Journal of Marine Research, 49(3):403-410
Babu M T, Sarma Y V B, Murty V S N, et al. 2003. On the circulation in the Bay of Bengal during northern spring inter-monsoon(March-April 1987). Deep-Sea Research:Part Ⅱ. Topical Studies in Oceanography, 50(5):855-865
Chaigneau A, Gizolme A, Grados C. 2008. Mesoscale eddies off Peru in altimeter records:Identification algorithms and eddy spatiotemporal patterns. Progress in Oceanography, 79(2-4):106-119.
Chelton D B, Schlax M G. 2003. The accuracies of smoothed sea surface height fields constructed from tandem satellite altimeter datasets. Journal of Atmospheric and Oceanic Technology, 20(9):1276-1302
Chelton D B, Schlax M G, Samelson R M. 2011. Global observations of nonlinear mesoscale eddies. Progress in Oceanography, 91(2):167-216
Chen Gengxin, Wang Dongxiao, Hou Yijun. 2012. The features and interannual variability mechanism of mesoscale eddies in the Bay of Bengal. Continental Shelf Research, 47:178-185
Cheng Xuhua, Xie Shangping, McCreary J P, et al. 2013. Intraseasonal variability of sea surface height in the Bay of Bengal. Journal of Geophysical Research:Oceans, 118(2):816-830
Doglioli A M, Blanke B, Speich S, et al. 2007. Tracking coherent structures in a regional ocean model with wavelet analysis:application to Cape Basin eddies. Journal of Geophysical Research:Oceans, 112(C5):C05043, doi: 10.1029/2006JC003952
Dussurget R, Birol F, Morrow R, et al. 2011. Fine resolution altimetry data for a regional application in the Bay of Biscay. Marine Geodesy, 34(3-4):447-476
Eigenheer A, Quadfasel D. 2000. Seasonal variability of the Bay of Bengal circulation inferred from TOPEX/Poseidon altimetry. Journal of Geophysical Research:Oceans, 105(C2):3243-3252
Hacker P, Firing E, Hummon J, et al. 1998. Bay of Bengal currents during the northeast monsoon. Geophysical Research Letters, 25(15):2769-2772
Isern-Fontanet J, García-Ladona E, Font J. 2003. Identification of marine eddies from altimetric maps. Journal of Atmospheric and Oceanic Technology, 20(5):772-778
Isern-Fontanet J, García-Ladona E, Font J, et al. 2006. Non-Gaussian velocity probability density functions:an altimetric perspective of the Mediterranean Sea. Journal of Physical Oceanography, 36(11):2153-2164
Le Traon P Y, Dibarboure G. 2004. An illustration of the contribution of the TOPEX/Poseidon-Jason-1 tandem mission to mesoscale variability studies. Marine Geodesy, 27(1-2):3-13
Le Traon P Y, Faugère Y, Hernandez F, et al. 2003. Can we merge Geosat follow-On with TOPEX/Poseidon and ERS-2 for an improved description of the ocean circulation?. Journal of Atmospheric and Oceanic Technology, 20(6):889-895
Le Traon P Y, Nadal F, Ducet N. 1998. An improved mapping method of multisatellite altimeter data. Journal of Atmospheric and Oceanic Technology, 15(2):522-534
Lou Hao, Bracco A, Di Lorenzo E. 2011. The interannual variability of the surface eddy kinetic energy in the Labrador Sea. Progress in Oceanography, 91(3):295-311
McWilliams J C. 1990. The vortices of two-dimensional turbulence. Journal of Fluid Mechanics, 219:361-385
Nencioli F, Dong Changming, Dickey T, et al. 2010. A vector geometry-based eddy detection algorithm and its application to a high-resolution numerical model product and high-frequency radar surface velocities in the Southern California Bight. Journal of Atmospheric and Oceanic Technology, 27(3):564-579, doi: 10.1175/2009JTECHO725.1
Pascual A, Faugère Y, Larnicol G, et al. 2006. Improved description of the ocean mesoscale variability by combining four satellite altimeters. Geophysical Research Letters, 33(2):L02611, doi: 10.1029/2005GL024633
Penven P, Echevin V, Pasapera J, et al. 2005. Average circulation, seasonal cycle, and mesoscale dynamics of the Peru Current System:a modeling approach. Journal of Geophysical Research:Oceans, 110(C10):C10021, doi: 10.1029/2005JC002945
Qiu Yun, Li Li. 2006. Review of study on upper ocean circulation in the Bay of Bengal. Advances in Marine Science (in Chinese), 24(4):593-603
Sadarjoen I A, Post F H. 2000. Detection, quantification, and tracking of vortices using streamline geometry. Computers & Graphics, 24(3):333-341
Sangrà P, Pascual A, Rodríguez-Santana Á, et al. 2009. The Canary Eddy Corridor:a major pathway for long-lived eddies in the subtropical North Atlantic. Deep-Sea Research:Part I. Oceanographic Research Papers, 56(12):2100-2114, doi: 10.1016/j.dsr.2009.08.008
Sarma Y V B, Rao E P R, Saji P K, et al. 1999. Hydrography and circulation of the Bay of Bengal during withdrawal phase of the southwest Monsoon. Oceanologica Acta, 22(5):453-471
Schott F A, McCreary Jr J P. 2001. The monsoon circulation of the Indian Ocean. Progress in Oceanography, 51(1):1-123
Somayajulu Y K, Murty V S N, Sarma Y V B. 2003. Seasonal and interannual variability of surface circulation in the Bay of Bengal from TOPEX/Poseidon altimetry. Deep-Sea Research:Part Ⅱ. Topical Studies in Oceanography, 50(5):867-880
Souza J M A C, de Boyer Montégut C, Le Traon P Y. 2011. Comparison between three implementations of automatic identification algorithms for the quantification and characterization of mesoscale eddies in the South Atlantic Ocean. Ocean Science, 7(3):317-334, doi: 10.5194/os-7-317-2011
Sreenivas P, Gnanaseelan C, Prasad K V S R. 2012. Influence of El Niño and Indian Ocean dipole on sea level variability in the Bay of Bengal. Global and Planetary Change, 80-81:215-225
Vinayachandran P N, Masumoto Y, Mikawa T, et al. 1999. Intrusion of the southwest monsoon current into the Bay of Bengal. Journal of Geophysical Research:Oceans, 104(C5):11077-11085
Yu Lisen, O'Brien J J, Yang Jiayan. 1991. On the remote forcing of the circulation in the Bay of Bengal. Journal of Geophysical Research:Oceans, 96(C11):20449-20454
Pablo Couve, Nixon Bahamon, Cristian M. Canales, et al. Systematic Review of Multi-Species Models in Fisheries: Key Features and Current Trends. Fishes, 2024, 9(10): 372. doi:10.3390/fishes9100372
Qinqin Lin, Yuying Zhang, Jiangfeng Zhu. Simulating the impacts of fishing on central and eastern tropical Pacific ecosystem using multispecies size-spectrum model[J]. Acta Oceanologica Sinica, 2022, 41(3): 34-43. doi: 10.1007/s13131-021-1902-3
Qinqin Lin, Yuying Zhang, Jiangfeng Zhu. Simulating the impacts of fishing on central and eastern tropical Pacific ecosystem using multispecies size-spectrum model[J]. Acta Oceanologica Sinica, 2022, 41(3): 34-43. doi: 10.1007/s13131-021-1902-3
Note: $\kappa $ in Eq. (4) is the abbreviation of carrying capacity; p in Eq. (7) and Eq. (8) represents the exponent of standard metabolism and the default value in the model is 0.75; $\epsilon $ in Eq. (11) is the efficiency of offspring production; SS in Eq. (14) is the selectivity size; r0 in Eq. (15) is the coefficient of regeneration rate of resources, and the default values in the model are 0.1 and 4, respectively; $\text{λ} $ in the Eq. (16) is the exponent of resource spectrum following the equation 2−n+q, where n is the exponent for maximum food intake (2/3) and q is the exponent for volumetric search rate (0.8). Other default parameters are listed in Table S2.
DownLoad:
DownLoad:
DownLoad: