Next Article in Journal
Geographic-Scale Harvest Program to Promote Invasivorism of Bigheaded Carps
Previous Article in Journal
Use of Environmental DNA to Detect Grass Carp Spawning Events
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Catch per Unit Effort Dynamic of Yellowfin Tuna Related to Sea Surface Temperature and Chlorophyll in Southern Indonesia

1
Department of Fishery Resources Utilization, Faculty of Fisheries and Marine Sciences, IPB University (Bogor Agricultural University), Bogor 16680, Indonesia
2
Environmental and Conservation Sciences, College of SHEE, Murdoch University, Perth 6150, Australia
3
Wildlife Conservation Society–Indonesia Program, Bogor 16128, Indonesia
4
Professorship for Aquaculture and Sea-Ranching, University of Rostock, 18051 Rostock, Germany
5
Faculty of Fisheries and Marine Sciences, IPB University (Bogor Agricultural University), Bogor 16680, Indonesia
6
Masyarakat dan Perikanan Indonesia Foundation, Denpasar 80223, Indonesia
*
Author to whom correspondence should be addressed.
Submission received: 3 August 2020 / Revised: 24 August 2020 / Accepted: 25 August 2020 / Published: 31 August 2020

Abstract

:
Tuna fisheries are the most valuable fisheries in the world, with an estimated market value of at least US$42 billion in 2018. Indonesia plays an important role in the global tuna fisheries and has committed to improve its fisheries management; therefore, a pilot of long-term spatial-temporal data bases was developed in 2012, however none have utilized data to have better understanding for management improvement. In this study, the annual and seasonal variation of large (≥10 kg) Yellowfin Tuna (YFT) catch per unit effort (CPUE) have been investigated and the influence of sea surface temperature (SST) and chlorophyll-a on these variables examined. We used fish landing data from West Nusa Tenggara recorded every month between 2012 and 2017 and analyzed using generalized linear models and generalized additive models. We found a seasonal and annual pattern of tuna abundance affected by SST and chlorophyll-a (chl a) and related to upwelling and El Nino event. These results also suggest that a two-month closure to fishing in August and September in southern Lombok is worth considering by the Government to maximize conservation of stocks due to a high abundance of juveniles emerging during the upwelling months from June to August.

1. Introduction

Tuna is one of the globally most important fish, caught by fishing vessels of more than 85 countries and contributing 5.7% or 4.8 million t to the global seafood production in the year 2016 [1,2]. In 2016, the value of tuna at all landing sites in the world was estimated to be more than US$10 billion annually [2], with approximately an additional US$42.2 billion market value [3]. They provide food, key nutrients, and income to coastal countries and their local communities [4]. Indonesia’s tuna fisheries are amongst the most important, diverse, and complex fisheries in the world. More than 130,000 vessels from 3 m in length to more than 100 m are operated with a total catch of more than 1 million tons of tuna per year [2,5]. Tuna is one of the most valuable species in Indonesia, with an export value of US$677.9 million in 2017. Hence, tuna is an important generator of wealth across all fleets, making it the most important seafood category for the country. The species of tuna caught in Indonesian waters include yellowfin (Thunnus albacares), bigeye (Thunnus obesus), albacore (Thunnus alalunga), and skipjack (Katsuwonus pelamis).
Tuna fishing has grown significantly in the waters of the Indonesian archipelago so that in 2016, Indonesia was the top country globally for tuna landings [2,6]. In the last 10 years, Indonesia has become involved in several regional fisheries management organizations, such as Indian Ocean Tuna Commission, Commission for the Conservation of Southern Bluefin Tuna, Western and Central Pacific Fisheries Commission, and the Inter-American Tropical Tuna Commission [7]. The Government of Indonesia has shown its commitment to sustainable tuna fisheries by issuing a Tuna Fisheries Management Plan in 2015 [7]. One of the most important components of this plan is the collection of data to improve decision support systems for tuna fisheries management. A key issue for fisheries scientists is the determination of fishing effort from the fishing fleets to be evaluated in relation to fishing mortality and it is commonly used as a basic input in fish stock assessment models [8]. Knowing yield per unit area and catch per unit effort can be complementary indicators of the stock status [9]; the first refers to the total catch without any reference to the fishing pressure, while the second one links the catch to the fishing efforts. Numerous projects and organizations have been involved in supporting the Government of Indonesia to develop data collection systems for tuna, including the Indonesia Marine and Climate Support (IMACs) project of United States Agency for International Development (USAID) Indonesia that developed the original Ifish system in 2012, which is now administered by Provincial Data Management Committees [10]. Data from tuna fishing boats are collected by enumerators through standardized methods and then uploaded to an online database (Ifish data) that is accessible to fishers, academicians and researchers, industries, and government. The Ifish data are available for advanced studies, awareness, policy development, and decision-making [10].
The abundance and composition of tuna resources caught by fishers are closely related to human factors, such as accessibility, proxy of fishing activity (boats characteristic), and experience in the sea, and other natural factors, such as weather condition, SST, chl a concentration, and the state of fish resources. Structural changes of the fishery resources are highly related to the utilization intensity of those resources. The oceanic features such as current dynamics, fronts, eddies, and convergences [11], and SST have been used to investigate productive frontal zones [12], which can further be used to indicate potential tuna fishing grounds. Thermal (or color) gradients in satellite images that arise from the circulation of water masses often indicate areas of high productivity [13]. Meanwhile, chl a data can also be used as a valuable indicator of water mass boundaries and may identify upwelling that can influence tuna distribution within a region [14,15]. Chl a concentrations over 0.20 mg·m−3, indicate high phytoplankton abundance for tuna prey, such as small pelagic fishery resources [16,17].
Fisheries in West Nusa Tenggara, one of the 34 provinces in Indonesia (Figure 1), are characterized as a deep-sea handline tuna fishery. This province contributes significantly to the national production of tuna, with yellowfin tuna, T. albacares, the most abundant species landed, comprising ~70% of the total catch [18]. Despite the significance of tuna in this province, few studies related to tuna fisheries have been made in this region [19]. We utilize the ‘Ifish database’ of the Data Management Committee, West Nusa Tenggara from 2012 to 2017, to investigate how CPUE of large tuna varies among different temporal factors, such as year, month of the year, different fishing days, and boat lengths, as well as the use of fish aggregation devices. The variation in CPUE of large T. albacares (>10 kg) related to sea-surface temperature and chl a is also examined. Thousands of records of CPUE from the Ifish database are analyzed using General Additive Models (GAMs) and Generalized Linear Models (GLMs) to evaluate the importance of these difference factors. Finally, we use the selected most parsimonious model to understand the seasonal variation of tuna resources and relative abundance related to environmental factors, i.e., chl a and sea surface temperature, aiming to provide recommendations for sustainable utilization and management of the tuna fishery.

2. Results

2.1. Distribution and Variation of Nominal Catch per Unit Effort

In total, 2831 trips caught YFT using handlines from 2012 to 2017, with 596 of these trips (~21.1%) fishing around fish aggregating devices (FADs). A total of 41 different grids were fished during all trips, with FAD fishing taking place in 23 of these grids (Figure 1b). The average total catch of large YFT was 15 fish from fishing trips with an average duration of 11 days (Table 1). The lowest mean catch was 10 fish per trip in 2010 and the greatest mean was 25 per trip in 2012. The average catch using FAD was 13 tunas per trip, while without FAD was 15 individuals per trip. The SST and chl a concentration at all fishing grounds ranged from 25.7 to 32.0 °C and 0.07 to 1.80 mg·m−3, respectively. The lowest mean annual SST was 28.3 °C in 2017, while the highest mean SST was 30.5 °C in 2012. For chl a, the lowest and the highest were 0.17 mg·m−3 (2012) and 0.42 mg·m−3 (2017), respectively (Table 1, Figure 2 and Figure 3).
An initial comparison of the total catches of large YFT (nominal/unstandardized catches) using analysis of variance indicated that they differed significantly among years and months (p < 0.01, Table 2). The highest and lowest annual mean catches were recorded in 2012 (25 individuals) and 2017 (10 individuals), respectively, while the highest and lowest monthly mean catches were in February (7 individuals) and October (20 individuals), respectively. Average catch of YFT (nominal/unstandardized catches) at the areas (grids) with FADs was significantly higher (p < 0.001, Table 2) than the catch at the areas (grids) without FAD. We also detected some nonlinear patterns when we plotted the YFT catches against the duration of fishing trips (days at sea), boat length, SST, and chl a concentration (Figure 2 and Figure 3).

2.2. Distribution and Variation of Standardized Catch per Unit Effort

When all factors were included in a single GLM and GAM, the most parsimonious model to explain the CPUE of large YFT included year, month, days at sea, boat length, chl a, and SST. We excluded the FAD because it did not significantly affect the standardized YFT catch. Later, we found that GLM models with higher Akaike information criterion (AIC) values had residuals that were not homogeneous (Table 3, Appendix A). Meanwhile, GAM models showed lower AIC with no obvious pattern of residuals detected (Table 3, Appendix A), hence we selected the generalized additive model to predict our data. We found that the best fitted GAM was: GAM Negative Binomial (NB) Model: CPUE~Year + month + spline (boat length) + spline (days at sea) + spline (Chl a) + spline (SST).
The factors of month, days at sea, year, boat length, and chl a concentration significantly affected the catch of YFT as a parametric effect, while days at sea, SST, and Chl a were significant when included as non-parametric terms (Table 4, Figure 4).
The standardized CPUE according to the estimated marginal mean of YFT catch (using length of boat = 9 m, trip = 9 days at sea, SST = 29.25 °C, and chl a = 0.15 mg·m−3) showed a fluctuating pattern for both monthly and annual CPUE (Figure 5 and Figure 6). Monthly standardized CPUE ranged from 7.6 in February to 19.3 in October and the annual standardized CPUE ranged from 10.8 in 2017 to 21.2 in 2012. When estimating the standardized CPUE from chl a concentration and SST, higher CPUEs were estimated at lower chl a concentrations, while CPUEs were not linearly influenced by SST (Figure 4 and Figure 7).

3. Discussion

Strong seasonal and annual patterns were recorded in the catch data of large (>10 kg) Yellowfin Tuna (YFT) in Lombok Island waters from 2012 to 2017. Based on GAM analyses, we found significant changes in the CPUE of large YFT recorded both among months of the year and among years. These patterns of YFT catch were also influenced by chl a concentration and ocean-atmosphere phenomena, such as sea-surface temperatures (SST) and El Niño events. The catch rates of large YFT were greater during the transition period of SST and chl a concentration, i.e., May–July (from rainy to dry seasons) and October–November (from dry to rainy seasons). During May and July, the SST decreased, and the chl a concentration increased, while the converse of this pattern was seen from October to November, i.e., SST increased and chl a concentration decreased. Chl a concentrations are influenced by SST [20], which in turn influences the net primary productivity of plankton. The abundance of tuna increases in high-productivity zones (such as upwelling zones, around islands, or at current interfaces), which attracts tuna to areas with greater food supply [21]. Therefore, the relative densities of the fish usually follow the distribution of primary production, with some time delays or lags.
Despite the limitations of the CPUE method, which include catch data and information about population dynamics that can help interpret productivity of fisheries and identify contradictions [22], our study overcome the shortcomings by adding environmental factors when considering the variation in standardized CPUE of large YFT among years. The greatest mean catches were recorded in 2012 and 2015, years with cooler monthly mean SSTs and mostly greater mean monthly chl a concentrations than in other years. A coupled ocean-atmosphere mode, namely an El Niño event (cooler waters and higher primary production), took place in the tropical Pacific Ocean during 2012 and 2015–2016, which might have influenced Lombok waters through the Indonesian Throughflow that passes through the Lombok Strait. In 2012, the El Niño occurred in August and September and was categorized as a weak El Niño [23], while in 2015, the El Niño event started to develop from April to May 2015, reached its maximum between November and December 2015, and finished in April–May 2016 [24]. This event had an El Niño index of 3.4 and was classified as a strong event [25], one of the strongest events recorded during the last two decades [24]. During El Niño events, it is likely that the water movement from Makassar Strait has lower flow [26], thus it could have influenced the upwelling level in the coastal waters of Lombok and Sumbawa waters because cold water with higher productivity in deep water flows to the surface. During the El Niño events, the coastal upwelling level in terms of chl a concentration was more intense than in other years due to no disturbance from the Indonesia Throughflow waters from Lombok Strait [27]. The abundance of various economically important fish species such as tuna is influenced by environmental conditions and food availability [21,28]. To a certain degree, the water movement related to the El Niño event determines the distribution of water temperature [29]. The distribution and abundance of large YFT are closely associated with physical changes in the ocean, such as temperature changes, storms (length and magnitude or intensity), and the El Niño, as well as variations in the mixing depth layer and disturbances caused by the wind [21,30]. Our study emphasis on role of ENSO using longer time series data and investigating the spatial variability of key environmental variables SST and food distribution in the Indian Ocean in relation to tuna distribution and its relationship with climate on different time scales could further improve tuna habitat models, as also suggested by previous researchers [31].
According to the best selected GAM model, the standardized CPUE can be significantly predicted by specific factors, such as month, fishing trip (days at sea), year, boat size, and chl a. The SST and chl a concentrations are the main variables that contribute to the upwelling phenomenon that influences the abundance of fish in this area. The peak of chl a concentration occurred in July and August, which indicated annual upwelling that used to occur annually in this region from June (see Ningsih et al. [32]). This was followed by the high catch of the mature size of tuna in October, and slightly reduced catch in November–December during 2012–2017, about three months after the peak in chl a. The “time-lag” between the upwelling event and the season of tuna catch is due to effect of trophic conversion, where the tuna prey reaches the highest abundance 3–7 months after the peak of phytoplankton or pigment concentration [14,33,34]. The influence of chl a on the catch of tuna resulted in varying degrees of lag effect on trophic conversion. Ortega-Garcı´a and Lluch-Cota [34] found a 3- to 5-month time-lag between high chl a and high tuna abundance in the tropical Pacific, while Lehodey et al. [33] recorded a longer time-lag of 3 to 7 months. Meanwhile, the correlation between upwelling time may have a direct correlation with tuna juveniles only [35]. Although many studies concluded that SST has an influence on tuna abundance, we did not find a significant parametric effect of SST on the standardized CPUE of large YFT, only a non-linear effect. A plausible explanation for our finding is that SST contributed to a nonlinear effect on the abundance of YFT in Lombok Island waters. Dunstan et al. [20] stated that there is no consistent correlation and mechanism between SST and chl a variation, and warming water results in decreasing chl a concentration as well as leading to increased chl a concentration. Additionally, YFT can efficiently regulate their body temperature and allow themselves to make large horizontal and vertical migrations, hence reducing the effect of SST on the CPUE [14,36]. Lan et al. [14] found that subsurface water temperature was linearly affecting the abundance of YFT, rather than the SST itself.
This study revealed a significant correlation between environmental parameter fluctuation and regular seasonality and annual abundance of YFT, indicated by seasonal and annual patterns of variation in the standardized catch per unit of effort (CPUE) (Figure 5, Table 4). These dynamics of the large tuna influence the movements of tuna hand-line fishers, especially those who fish in grids without FADs. Currently, small-scale tuna management in West Nusa Tenggara has not yet adopted spatial management, e.g., using time/area closures on the tuna fishing grounds. These closures might be effective in conserving juvenile tuna and reducing bycatch. Our results suggest that the government could consider implementing a temporal closure in the waters of southern Lombok for two months after the chl a concentration reaches the highest level from August to September. At this time, the abundance of large YFT is low and the abundance of juvenile YFT is probably very high due to the upwelling phenomenon. A closure at this time would thus allow the juvenile tuna to grow without experiencing fishing mortality and could lead to greater recruitment of larger fish. A three-month closure on FAD fishing in the Atlantic appears to have been successful in reducing the catch of juvenile YFT and bigeye tuna [37]. The results of area and month closures for tuna fishing in the Banda Sea, of southern Maluku Province, also indicate that a combination of area and temporal closures would aid in the recovery of YFT stocks and could potentially be beneficial for fisheries [38]. However, the decision of location choices, sizes, shapes, and duration of closure (seasonal or for a number of years) should be carefully studied and characterized by environmental parameters.

4. Materials and Methods

The waters of Indonesia’s Exclusive Economic Zone (EEZ) are divided into 11 fisheries management areas (FMAs) that were formed in 1999 to optimize fisheries management in Indonesia [39]. The 11 FMAs fall into two major regional fishing areas recognised by the FAO: The Eastern Indian Ocean zone (FMAs 571, 572, and 573) and the Western and Central Pacific (FMAs 711, 712, 713, 714, 715, 716, 717, and 718). This study took place in Lombok Island, West Nusa Tenggara (WNT) Indonesia, part of FMA 713, which includes the relatively shallow waters of the Flores Sea and FMA 573 in the eastern Indian Ocean (Figure 1). The waters in this FMA are characterized as deep sea, with depths of up to 1000 m within 5 nautical miles of the coast in some locations. Upwelling events occur frequently in these waters due to its bathymetry, currents, and other oceanography features [32]. Small-scale fisheries predominate in this province and the most prominent landing sites in Lombok Island are Labuhan Lombok, Tanjung Luar, and various traditional landing sites, such as Bangko-bangko and Tanjung Ringgit (Figure 1a). Labuhan Lombok is the most important landing site for tuna on Lombok Island as this is where yellowfin tuna, skipjack tuna, and frigate tuna (Auxis thazard) are all landed at its coastal harbor (Figure 1).

4.1. Data Collection

4.1.1. Fisheries Data

We used landing data for tuna from the Data Management Committee (DMC) of West Nusa Tenggara province, which oversees and manages data collection and management for the tuna fishery in the province. The provincial fisheries office employs three trained enumerators to collect data from each trip in Labuhan Lombok during working days (Monday to Friday; at least 20 days a month) and one supervisor to oversee and ensure data quality. We extracted data 2012 to 2017 from the Ifish system for this region and only used data from fishing trips that used handlines (Table 5), as this is the dominant fishing gear used in the tuna fishery in Lombok Island [18]. Note that data collection in 2012 started in February and in 2017 covered the months from January to October. The analyses concentrated on large yellowfin tuna (≥10 kg) that have reached maturity [10,40].

4.1.2. Sea Surface Temperature and Chlorophyll Data

SST and chl a data were downloaded from the Environmental Research Division’s Data Access Program (ERDDAP; https://coastwatch.pfeg.noaa.gov/erddap/info/index.html). SST measurements were gathered by the Moderate Resolution Imaging Spectroradiometer (MODIS) that was carried aboard by the NASA’s Aqua spacecraft, while chl a measurements were also gathered using the MODIS and generated according to algorithm developed by Hu et al., 2012 [41]. We selected the following data from this source: grid code (total area of a grid = 1.2 million ha) (see Figure 1b), days at sea for each fishing trip, and the catch of large YFT. We used average SST and chl a for each fishing trip: for trips of ≤3 d = average over the days of the trip; trips of 4–10 days in length = 8–day mean value; and for trips of >10 d at sea, the 30 d mean value. When the daily SST or chl a data were not available, we used moving averages (see Johnston et al. [42]) to estimate the missing values.

4.2. Data Analyses

Prior to the statistical analyses, we plotted all data (catch versus all factors) and then used analysis of variance [43] to understand the patterns of variation in unstandardized catch and the factors contributing towards this variation. Then, we used Generalized Linear Models (GLM) and Generalized Additive Models (GAMs) to standardize CPUE and to understand the influence of month of the year (as a proxy for the season), SST, and chl a on the abundance of large YFT. The full models for both types of analyses were:
  • CPUE~Year + month + boat length + days at sea + Chl a + SST… fitted with a simple linear regression or Generalized Linear Model (GLM);
  • CPUE~Year + month + lo (boat length) + lo (days at sea) + lo (Chl a) + lo (SST)… fitted with a Generalized Additive Model (GAM) using a loess/local regression, lo [44];
  • CPUE~Year + month + s (boat length) + s (days at sea) + s (Chl a) + s(SST)… fitted with a Generalized Additive Model (GAM) using a smoothing function (s(x)) [44].
All models were run using R software [45] and the best model was selected using the dredge function through automated model selection with subsets of the full model and using the Akaike information criterion (AIC) values [46]. We then inspected the residuals from selected models using GLMs with different families and GAMs with different smoothing functions, and compared the fit of the models using the AIC values to determine the most parsimonious model [47]. Multi-collinearity was checked using the variance inflation factor [48]. We estimated the standardized CPUE for each month, year, range of SST, and range of chl a according to the mode of the set of data values.

Author Contributions

Conceptualization, B.W. and I.Y.; methodology, B.W., I.Y., U.M., J.P., W., P.S.T. and D.D.; software, I.Y. and U.M.; validation, B.W., N.L., U.M., S.K., P.I.W., J.P., W., P.S.T., D.D. and I.Y.; formal analysis; B.W., N.L., U.M., S.K., P.I.W., J.P. and I.Y.; writing—original draft preparation, B.W. and I.Y.; writing—review and editing, B.W., N.L., U.M., S.K., P.I.W., J.P., W., P.S.T., D.D. and I.Y.; visualization, I.Y and U.M.; supervision, B.W., N.L., and I.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

We are grateful to Lalu Wahyudi Adiguna and Nurjamil from Fisheries Co-Management Committee of West Nusa Tenggara Province and all staff of Masyarakat dan Perikanan Indonesia Foundation (MDPI) for supporting this research.

Conflicts of Interest

Deirdre Duggan, Wildan, and Putra Satria Timur are employee of Masyarakat dan Perikanan Indonesia Foundation (MDPI), who lead and organize the project on the Fisheries Improvement Program and Ifish data base for handline tuna fisheries in West Nusa Tenggara. Those data were used in this manuscript. However, we declare that MDPI affiliation do not influence the writing or interpretation of the results from the study reported in this manuscript. Other authors declare there is no conflict of interest related to this work.

Appendix A

Figure A1. Residual plot and AIC of all full models tested using different GLM distributions (ac) and different GAM smoothing lines (d,e).
Figure A1. Residual plot and AIC of all full models tested using different GLM distributions (ac) and different GAM smoothing lines (d,e).
Fishes 05 00028 g0a1

References

  1. FAO. Tuna a Global Perspective; FAO: Rome, Italy, 2018. [Google Scholar]
  2. FAO. The State of World Fisheries and Aquaculture (SOFIA)—Meeting the Sustainable Development Goals; FAO: Rome, Italy, 2018; ISBN 978-92-5-130562-1. [Google Scholar]
  3. Tidd, A.; Blanchard, J.L.; Kell, L.; Watson, R.A. Predicting Global Tuna Vulnerabilities with Spatial, Economic, Biological and Climatic Considerations. Sci. Rep. 2018, 8, 1–12. [Google Scholar] [CrossRef] [PubMed]
  4. McCluney, J.K.; Anderson, C.M.; Anderson, J.L. The Fishery Performance Indicators for Global Tuna Fisheries. Nat. Commun. 2019, 10, 1–9. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Sunoko, R.; Huang, H.W. Indonesia Tuna Fisheries Development and Future Strategy. Mar. Pol. 2014, 43, 174–183. [Google Scholar] [CrossRef]
  6. Coulter, A.; Cashion, T.; Cisneros-Montemayor, A.M.; Popov, S.; Tsui, G.; Le Manach, F.; Schiller, L.; Palomares, M.L.D.; Zeller, D.; Pauly, D. Using Harmonized Historical Catch Data to Infer the Expansion of Global Tuna Fisheries. Fish. Res. 2020, 221, 105379. [Google Scholar] [CrossRef]
  7. MMAF. Tuna Fisheries Management Plan; Ministry of Marine Affairs and Fisheries: Jakarta, Indonesia, 2015. [Google Scholar]
  8. Bordalo-Machado, P. Fishing Effort Analysis and Its Potential to Evaluate Stock Size. Rev. Fish. Sci. 2006, 14, 369–393. [Google Scholar] [CrossRef]
  9. Kantoussan, J.; Ecoutin, J.M.; Guy Fontenelle, G.; de Morais, L.T.; Lae, R. Catch Per Unit Effort and Yields As Indicators Of Exploited Fish Communities: Application To Two West African Reservoirs. Lake. Reserv. Res. Manag. 2014, 19, 86–97. [Google Scholar] [CrossRef]
  10. MDPI. Data Collection Protocol for Small Scale Handline Fisheries in Indonesia; MDPI: Denpasar, Indonesia, 2018. [Google Scholar]
  11. Polovina, J.J.; Howell, E.A. Ecosystem Indicators Derived from Satellite Remotely Sensed Oceanographic Data for the North Pacific. ICES J. Mar. Sci. 2005, 62, 319–327. [Google Scholar] [CrossRef]
  12. Zainuddin, M.; Saitoh, S.I.; Saitoh, K. Detection of Potential Fishing Ground for Albacore Tuna Using Synoptic Measurements of Ocean Color and Thermal Remote Sensing in the Northwestern North Pacific. Geophys. Res. Lett. 2004, 31, 1–4. [Google Scholar] [CrossRef]
  13. Saitoh, S.; Chassot, E.; Dwivedi, R.; Fonteneau, A.; Kiyofuji, H.; Kumari, B.; Kuno, M.; Matsumura, S.; Platt, T.; Raman; et al. Remote Sensing in Fisheries and Aquaculture; Forget, M.H., Stuart, V., Platt, T., Eds.; The International Ocean-Colour Coordinating Group: Dartmouth, NS, Canada, 2009; Volume 5, pp. 57–76. ISBN 1098-6030. [Google Scholar]
  14. Lan, K.W.; Lee, M.A.; Lu, H.J.; Shieh, W.J.; Lin, W.K.; Kao, S.C. Ocean Variations Associated with Fishing Conditions for Yellowfin Tuna (Thunnus Albacares) in the Equatorial Atlantic Ocean. ICES J. Mar. Sci. 2011, 68, 1063–1071. [Google Scholar] [CrossRef]
  15. Song, L.; Zhou, J.; Zhou, Y.; Nishida, T.; Jiang, W.; Wang, J. Environmental Preferences of Bigeye Tuna, Thunnus Obesus, in the Indian Ocean: An Application to a Longline Fishery. Environ. Biol. Fishes. 2009, 85, 153–171. [Google Scholar] [CrossRef]
  16. Polovina, J.J.; Howell, E.; Kobayashi, D.R.; Seki, M.P. The Transition Zone Chlorophyll Front, a Dynamic Global Feature Defining Migration and Forage Habitat for Marine Resources. Prog. Oceanogr. 2001, 49, 469–483. [Google Scholar] [CrossRef]
  17. Zainuddin, M.; Farhum, A.; Safruddin, S.; Selamat, M.B.; Sudirman, S.; Nurdin, N.; Syamsuddin, M.; Ridwan, M.; Saitoh, S.I. Detection of Pelagic Habitat Hotspots for Skipjack Tuna in the Gulf of Bone-Flores Sea, Southwestern Coral Triangle Tuna, Indonesia. PLoS ONE 2017, 12, 1–19. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Duggan, D.; Wildan; Riza Baroqi, A.; Satria Timur, P.; Juhrin; Nababan, N.; Kochen, M. Simple Economics in Small Scale Tuna Fisheries: Results of an Economic-Survey Conducted in Four MDPI Sites in Eastern Indonesia; MDPI: Denpasar, Indonesia, 2017. [Google Scholar]
  19. Setyadji, B.; Nugraha, B. Dynamics of Gears, Fleets, Catch and Fishing Season of Smallscale Tuna Fisheries in Labuhan Lombok, West Nusa Tenggara. Indones. Fish. Res. J. 2015, 21, 99–107. [Google Scholar] [CrossRef] [Green Version]
  20. Dunstan, P.K.; Foster, S.D.; King, E.; Risbey, J.; O’Kane, T.J.; Monselesan, D.; Hobday, A.J.; Hartog, J.R.; Thompson, P.A. Global Patterns of Change and Variation in Sea Surface Temperature and Chlorophyll, A. Sci. Rep. 2018, 8, 1–9. [Google Scholar] [CrossRef] [Green Version]
  21. Lehodey, P.; Bertignac, M.; Hampton, J.; Lewis, A.; Picaut, J. El Niño Southern Oscillation and tuna in the western Pacific. Nature 1997, 389, 715–718. [Google Scholar] [CrossRef]
  22. Maunder, M.N.; Sibert, J.R.; Fonteneau, A.; Hampton, J.; Kleiber, P.; Harley, S.J. Interpreting Catch Per Unit Effort Data to Assess the Status of Individual Stocks and Communities. ICES J. Mar. Sci. 2006, 63, 1373–1385. [Google Scholar] [CrossRef]
  23. BMKG. Data Data El Niño/La Niña; Meteorology and Geophysic Agency of Indonesia: Jakarta, Indonesia, 2019. [Google Scholar]
  24. Iskandar, I.; Utari, P.A.; Lestari, D.O.; Sari, Q.W.; Setiabudidaya, D.; Khakim, M.Y.N.; Yustian, I.; Dahlan, Z. Evolution of 2015/2016 El Niño and Its Impact on Indonesia. AIP 2017, 080001–080005. [Google Scholar] [CrossRef]
  25. Null, J. El Niño and La Niña Years and Intensities. Available online: https://ggweather.com/enso/oni.htm (accessed on 11 August 2019).
  26. Susanto, R.D.; Gordon, A.L.; Zheng, Q. Upwelling within the Indonesian Seas and its relation to Monsoon and ENSO. Geophys. Res. Lett. 2001, 28, 1599–1602. [Google Scholar] [CrossRef]
  27. Hendiarti, N.; Siegel, H.; Ohde, T. Investigation of Different Coastal Processes in Indonesian Waters Using SeaWiFS Data. Deep-Sea Res. Pt. II 2004, 51, 85–97. [Google Scholar] [CrossRef]
  28. Roger, C. Relationships among Yellowfin and Skipjack Tuna, Their Prey-fish and Plankton in the Tropical Western Indian Ocean. Fish. Oceanogr. 1994, 3, 133–141. [Google Scholar] [CrossRef]
  29. Suárez-Sánchez, J.; Ritter-Ortiz, W.; Gay-García, C.; Torres-Jácome, J. ENSO-Tuna Relations in the Eastern Pacific Ocean and Its Prediction as a Non-Linear Dynamic System. Atmosfera 2004, 17, 245–258. [Google Scholar]
  30. Blackburn, M. Conditions Related to Upwelling Which Determine Distribution of Tropical Tunas off Western Baja California. Fish. Bull. 1969, 68, 147–176. [Google Scholar]
  31. Kumar, P.S.; Pillai, G.N.; Manjush, U. El Nino Southern Oscillation (ENSO) Impact on Tuna Fisheries in Indian Ocean. SpringerPlus 2014, 3, 1–13. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Ningsih, N.S.; Rakhmaputeri, N.; Harto, A.B. Upwelling Variability along the Southern Coast of Bali and in Nusa Tenggara Waters. Ocean Sci. J. 2013, 48, 49–57. [Google Scholar] [CrossRef]
  33. Lehodey, P.; Andre, J.M.; Bertignac, M.; Hampton, J.; Stoens, A.; Menkes, C.; Memery, L.; Grima, N. Predicting Skipjack Tuna Forage Distributions in the Equatorial Pacific Using a Coupled Dynamical Bio-Geochemical Model. Fish. Oceanogr. 1998, 7, 317–325. [Google Scholar] [CrossRef]
  34. Ortega-García, S.; Lluch-Cota, S.E. Distribución de la abundancia de atún y su relación con la concentración de pigmentos fotosintéticos derivados de satélite en aguas al sur de México. Invest. Geogr. México 1996, 4, 85–93. [Google Scholar]
  35. Singh, A.A.; Sakuramoto, K.; Suzuki, N. Impact of Climatic Factors on Albacore Tuna Thunnus alalung in the South Pacific Ocean. Am. J. Clim. Chang. 2015, 4, 295–312. [Google Scholar] [CrossRef] [Green Version]
  36. Zagaglia, C.R.; Lorenzzetti, J.A.; Stech, J.L. Remote Sensing Data and Longline Catches of Yellowfin Tuna (Thunnus Albacares) in the Equatorial Atlantic. Remote. Sens. 2004, 93, 267–281. [Google Scholar] [CrossRef]
  37. Goujon, M.; Labaisse-Bodilis, C. Effets Des Plans de Protection Des Thonidés de l’Atlantique Depuis 1997 d’après Les Observations Faites Sur Les Thoniers Senneurs Gérés Par Les Armements Français. Collect. Vol. Sci. Pap. ICCAT 2001, 52, 575–589. [Google Scholar]
  38. Satrioajie, W.N.; Suyadi; Syahailatua, A.; Wouthuyzen, S. The Importance of the Banda Sea for Tuna Conservation Area: A Review of Studies on the Biology and the Ecology of Tuna. IOP Conf. Ser. Earth Environ. Sci. 2018, 184. [Google Scholar] [CrossRef]
  39. Ministry of Agriculture. Potensi Sumber Daya Ikan dan Jumlah Tangkapan yang Diperbolehkan (JTB); Ministry of Agriculture: Jakarta, Indonesia, 1999. [Google Scholar]
  40. Compean-Jimenez, G.A.; Dreyfus-Leon, M.J. Interaction between vessels fishing for yellowfin tuna (Thunnus albacares) in the northeastern and southeastern Pacific. In Status of Interaction of Pacific Tuna Fisheries in 1995, Proceeding of the Second FAO Expert Consultation on Interaction of Pacific Tuna Fisheries, Shimizu, Japan, 23–31 January 1995; Shomura, R.S., Majkowski, J., Harman, R.F., Eds.; FAO: Rome, Italy, 1996; pp. 339–349. [Google Scholar]
  41. Hu, C.; Lee, Z.; Franz, B. Chlorophyll a Algorithms For Oligotrophic Oceans: A Novel Approach Based On Three-Band Reflectance Difference. J. Geophys. Res. 2012, 117, 1–25. [Google Scholar] [CrossRef] [Green Version]
  42. Johnston, F.R.; Boyland, J.E.; Meadows, M.; Shale, E. Some Properties of a Simple Moving Average When Applied to Forecasting a Time Series. J. Oper. Res. Soc. 1999, 50, 1267–1271. [Google Scholar] [CrossRef]
  43. Chambers, J.M.; Freeny, A.E.; Heiberger, R.M. Statistical Models in S.; Routledge: New York, NY, USA, 2017; pp. 145–193. ISBN 9780203738535. [Google Scholar]
  44. Hastie, T. gam: Generalized Additive Models. Available online: https://cran.r-project.org/package=gam (accessed on 15 January 2020).
  45. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2018. [Google Scholar]
  46. Barton, K. MuMIn: Multi-Model Inference. Available online: https://cran.r-project.org/package=MuMIn (accessed on 15 January 2020).
  47. Zuur, A.; Ieno, E.N.; Walker, N.; Saveliev, A.A.; Smith, G.M. Mixed Effects Models and Extensions in Ecology With R.; Gail, A., Krickeberg, K., Samet, J., Tsiatis, A., Wong, W., Eds.; Springer Science & Business Media, LCC: New York, NY, USA, 2009; ISBN 978-0-387-87457-9. [Google Scholar]
  48. Dormann, C.F.; Elith, J.; Bacher, S.; Buchmann, C.; Carl, G.; Carré, G.; Marquéz, J.R.G.; Gruber, B.; Lafourcade, B.; Leitão, P.J.; et al. Collinearity: A Review of Methods to Deal with It and a Simulation Study Evaluating Their Performance. Ecography 2013, 36, 27–46. [Google Scholar] [CrossRef]
Figure 1. Map of (a) Indonesian Fisheries Management Areas (FMAs) showing the study region within West Nusa Tenggara and (b) fish landing and monitoring site, as well as the grids covering the fishing ground in the eastern Indian Ocean, Bali and Flores Seas, West Nusa Tenggara, Indonesia; combinations of letters (right) and numbers (below) are the grid codes for the fishing ground. N/A = not applicable as there are no records from these grids.
Figure 1. Map of (a) Indonesian Fisheries Management Areas (FMAs) showing the study region within West Nusa Tenggara and (b) fish landing and monitoring site, as well as the grids covering the fishing ground in the eastern Indian Ocean, Bali and Flores Seas, West Nusa Tenggara, Indonesia; combinations of letters (right) and numbers (below) are the grid codes for the fishing ground. N/A = not applicable as there are no records from these grids.
Fishes 05 00028 g001
Figure 2. Box plots of catch per trip, number of large (>10 kg) Yellowfin Tuna (YFT), Thunnus albacares, YFT individuals per trip by (a) year, (b) monthly average from all years, and (c) fish aggregating device. Dark horizontal line shows the median, boxes cover the 25th and 75th percentiles and the whiskers show the 5th and 95th percentiles.
Figure 2. Box plots of catch per trip, number of large (>10 kg) Yellowfin Tuna (YFT), Thunnus albacares, YFT individuals per trip by (a) year, (b) monthly average from all years, and (c) fish aggregating device. Dark horizontal line shows the median, boxes cover the 25th and 75th percentiles and the whiskers show the 5th and 95th percentiles.
Fishes 05 00028 g002
Figure 3. Plots of the total catches of large YFT versus (a) day at sea, (b) boat length, (c) chl a, and (d) sea surface temperature (SST).
Figure 3. Plots of the total catches of large YFT versus (a) day at sea, (b) boat length, (c) chl a, and (d) sea surface temperature (SST).
Fishes 05 00028 g003
Figure 4. Factor (year and month) and smoothing line (boat length, days at sea, chl a, and SST) of generalized additive model for number of large (>10 kg) Yellowfin Tuna, Thunnus albacares, individuals per trip by (a) year, (b) month, (c) boat length, (d) days at sea, (e) chl a, and (f) SST; degrees of freedom were determined using cross validation. Density of dots on X-axis reflects the number of points.
Figure 4. Factor (year and month) and smoothing line (boat length, days at sea, chl a, and SST) of generalized additive model for number of large (>10 kg) Yellowfin Tuna, Thunnus albacares, individuals per trip by (a) year, (b) month, (c) boat length, (d) days at sea, (e) chl a, and (f) SST; degrees of freedom were determined using cross validation. Density of dots on X-axis reflects the number of points.
Fishes 05 00028 g004
Figure 5. Seasonal (a) standardized catch per unit effort (CPUE) of large (>10 kg) Yellowfin Tuna Thunnus albacares, (b) seasonal SST, and (c) Chl a concentration from all fishing grounds plotted using boxplots.
Figure 5. Seasonal (a) standardized catch per unit effort (CPUE) of large (>10 kg) Yellowfin Tuna Thunnus albacares, (b) seasonal SST, and (c) Chl a concentration from all fishing grounds plotted using boxplots.
Fishes 05 00028 g005
Figure 6. Annual (a) standardized CPUE of large (>10 kg) Yellowfin Tuna, Thunnus albacares, (b) monthly SST, and (c) Chl a concentration from all fishing grounds plotted using boxplots.
Figure 6. Annual (a) standardized CPUE of large (>10 kg) Yellowfin Tuna, Thunnus albacares, (b) monthly SST, and (c) Chl a concentration from all fishing grounds plotted using boxplots.
Fishes 05 00028 g006
Figure 7. Standardized CPUE of large (>10 kg) Yellowfin Tuna, Thunnus albacares fitted using a General Additive Model by (a) variation of SST and (b) chl a concentration.
Figure 7. Standardized CPUE of large (>10 kg) Yellowfin Tuna, Thunnus albacares fitted using a General Additive Model by (a) variation of SST and (b) chl a concentration.
Fishes 05 00028 g007
Table 1. Summary of the variables used in the analysis of large yellowfin tuna (>10 kg), Thunnus albacares, caught from West Nusa Tenggara, Indonesia, showing the mean value (±1 SD) for Yellowfin Tuna catch, trip duration, length of boat, SST, and chlorophyll-a (chl a) from this study period (2012–2017).
Table 1. Summary of the variables used in the analysis of large yellowfin tuna (>10 kg), Thunnus albacares, caught from West Nusa Tenggara, Indonesia, showing the mean value (±1 SD) for Yellowfin Tuna catch, trip duration, length of boat, SST, and chlorophyll-a (chl a) from this study period (2012–2017).
YearCatch of YFT (Individual)Trip Duration (Days)Length of Boat (m)SST
(°C)
Chl a (mg m−3)
Mean ± SDMean ± SDMean ± SDMean ± SDMean ± SD
2012 (n = 131)25 ± 1810.45 ± 3.1117.30 ± 3.6330.48 ± 0.720.17 ± 0.03
2013 (n = 816)15 ± 1110.73 ± 3.6016.03 ± 2.9729.07 ± 1.050.24 ± 0.12
2014 (n = 796)14 ± 1211.84 ± 3.6914.34 ± 2.6528.48 ± 1.350.33 ± 0.25
2015 (n = 603)16 ± 1412.16 ± 3.3614.58 ± 2.8128.50 ± 1.420.33 ± 0.20
2016 (n = 400)14 ± 1211.61 ± 3.7313.62 ± 2.3329.55 ± 1.110.18 ± 0.06
2017 (n = 117)10 ± 911.65 ± 3.7713.23 ± 1.7528.29 ± 1.130.42 ± 0.31
2012–2017 (n = 2831)15 ± 1311.48 ± 3.6314.87 ± 2.9528.89 ± 1.330.28 ± 0.19
Table 2. Analysis of variance for nominal catch of large Yellowfin Tuna (>10 kg), Thunnus albacares, landed in West Nusa Tenggara, Indonesia.
Table 2. Analysis of variance for nominal catch of large Yellowfin Tuna (>10 kg), Thunnus albacares, landed in West Nusa Tenggara, Indonesia.
ParameterDegree of FreedomSum of SquareMean of SquareF Valuep-Value
Year518,740374824.57<2 × 10−16
Residuals2825430,887153
Month1138,558350524.04<2 × 10−16
Residuals2819411,069146
FAD117881787.911.290.000788
Residuals2829447,839158.3
Table 3. Description of the different models used to test the best fit model for large Yellowfin Tuna (≥10 kg), Thunnus albacares, and their over-dispersion and AIC values to select the “best” model (lo = loess/local regression; span = value to control loess smoothing (0.1–0.9); s() = smoothing spline; θ = df = degrees of freedom. Plots of residuals for each model are shown in Appendix A.
Table 3. Description of the different models used to test the best fit model for large Yellowfin Tuna (≥10 kg), Thunnus albacares, and their over-dispersion and AIC values to select the “best” model (lo = loess/local regression; span = value to control loess smoothing (0.1–0.9); s() = smoothing spline; θ = df = degrees of freedom. Plots of residuals for each model are shown in Appendix A.
NoFull ModelOver-DispersionAIC
Generalized Linear Models
1CPUE ~ Year + month + boat length + days at sea + Chl a + SST; family = GaussianNA22,263.00
2CPUE ~ Year + month + boat length + days at sea + Chl a + SST;
family = Poisson
9.837,384.00
3CPUE ~ Year + month + boat length + days at sea + chl a + SST;
family = negative-binomial
9.820,789.00
Generalized Additive Models
4CPUE~ Year + month + lo(boat length, span = 0.5) + lo(days at sea, span = 0.5) + lo(Chl a, span = 0.4) + lo(Chl a, span = 0.2);
family = negative binomial (θ = 1.75)
NA20,415.73
5CPUE ~ Year + month + s(boat length, df = 3) + s(days at sea, df = 5.5) + s(Chl a, df = 9.4) + s(SST, df = 9); family = negative binomial (θ = 1.8)NA20,403.25
Table 4. Summary for the best fitted Generalised Additive Model (GAM) for large Yellowfin Tuna (>10 kg), Thunnus albacares, landed in West Nusa Tenggara, Indonesia.
Table 4. Summary for the best fitted Generalised Additive Model (GAM) for large Yellowfin Tuna (>10 kg), Thunnus albacares, landed in West Nusa Tenggara, Indonesia.
ParameterDegree of FreedomSum of SquareMean of Square% MSF Valuep-Value
ANOVA for Parametric Effects
Month11309.9328.175235.6427.88>0.0001
Year565.6113.122716.6012.98521.63 × 10−12
Chl a15.325.32376.735.26790.022
Boat length16.856.84858.666.77680.009
SST10.080.08380.110.08290.773458
Days at sea124.4924.49230.9824.23549.02 × 10−07
Residuals27872816.511.01061.28
ANOVA for Nonparametric Effects
Month
Year
Chl a8.4 2.19230.02299
Boat length2 2.88340.05611
SST8.1 6.0727.65 × 10−08
Days at sea4.5 22.0429<2.2 × 10−16
Table 5. Summary of the variables extracted from the Ifish data base provided through the data Management Committee of West Nusa Tenggara province in eastern Indonesia. Data retrieved December 2017.
Table 5. Summary of the variables extracted from the Ifish data base provided through the data Management Committee of West Nusa Tenggara province in eastern Indonesia. Data retrieved December 2017.
VariableTypeDescription
YearCategoricalYear of trip (2012–2017)
MonthCategoricalMonth of trip (January–December)
Days at seaNumericTotal day of fishing trip
Boat lengthNumericLength of boat in meters
Fish aggregating device (FAD)CategoricalFishing on FAD or not
Fishing GroundCategoricalGrid code of fishing areas (Figure 1)
CatchNumericNumber of fish with weight ≥10 kg caught on each fishing trip

Share and Cite

MDPI and ACS Style

Wiryawan, B.; Loneragan, N.; Mardhiah, U.; Kleinertz, S.; Wahyuningrum, P.I.; Pingkan, J.; Wildan; Timur, P.S.; Duggan, D.; Yulianto, I. Catch per Unit Effort Dynamic of Yellowfin Tuna Related to Sea Surface Temperature and Chlorophyll in Southern Indonesia. Fishes 2020, 5, 28. https://0-doi-org.brum.beds.ac.uk/10.3390/fishes5030028

AMA Style

Wiryawan B, Loneragan N, Mardhiah U, Kleinertz S, Wahyuningrum PI, Pingkan J, Wildan, Timur PS, Duggan D, Yulianto I. Catch per Unit Effort Dynamic of Yellowfin Tuna Related to Sea Surface Temperature and Chlorophyll in Southern Indonesia. Fishes. 2020; 5(3):28. https://0-doi-org.brum.beds.ac.uk/10.3390/fishes5030028

Chicago/Turabian Style

Wiryawan, Budy, Neil Loneragan, Ulfah Mardhiah, Sonja Kleinertz, Prihatin Ika Wahyuningrum, Jessica Pingkan, Wildan, Putra Satria Timur, Deirdre Duggan, and Irfan Yulianto. 2020. "Catch per Unit Effort Dynamic of Yellowfin Tuna Related to Sea Surface Temperature and Chlorophyll in Southern Indonesia" Fishes 5, no. 3: 28. https://0-doi-org.brum.beds.ac.uk/10.3390/fishes5030028

Article Metrics

Back to TopTop