## PENENTUAN STATUS PEMANFAATAN DAN SKENARIO PENGELOLAAN
IKAN CAKALANG ( Katsuwonus pelamis ) YANG TERTANGKAP DI PERAIRAN BOLAANG-MONGONDOW SULAWESI UTARA
John  S. Kekenusa 1) , Victor N.R. Watung 2) , dan Djoni Hatidja 1)
1) PS Matematika FMIPA UNSRAT, Manado
2) PS Manajemen Sumberdaya Perikanan FPIK UNSRAT, Manado e-mail: johnskekenusa@yahoo.com; vnrwatung@yahoo.com; dhatidja@gmail.com

## ABSTRAK
Ikan  Cakalang  (Katsuwonus  pelamis  L.),  perlu  dikelola  dengan  baik  sebagai  sumber  daya  alam terbarukan,  tetapi  semakin  habis  atau  punah.  Salah  satu  pendekatan  dalam  pengelolaan sumberdaya  ikan  adalah  dengan  pemodelan.  Analisis  dilakukan,  bertujuan  untuk  mendapatkan model  terbaik  untuk  model  surplus  produksi  adalah  Maximum  Sustainable  Yield  (MSY),  dan untuk mendapatkan pemanfaatan cakalang dan tingkat usaha.  Data yang digunakan untuk model produksi  surplus  dikumpulkan  dari  data  hasil  tangkapan  cakalang  yang  disediakan  oleh  Dinas Kelautan  Dan  Perikanan  Kabupaten    Bolaang  Mongondow  Provinsi  Sulawesi  Utara.  Model terbaik Surplus Produksi dapat digunakan untuk menilai potensi hasil cakalang adalah model Fox. Upaya  optimal  per  tahun  secara  biologi  EMSY  adalah  3.876  perjalanan.  Hasil  optimal  secara biologi  CMSY  1,505.00  ton  per  tahun.  Tingkat  pemanfaatan  untuk  tahun  2005  adalah  124,48%, dengan tingkat usaha 110,96%. Hasil tangkapan untuk tahun 2005 sudah di atas nilai MSY, yang menunjukkan overfishing
Kata kunci: ikan cakalang, Model surplus produksi, Maksimum Sustainable Yield, Bolaang Monogndow
DETERMINATION OF THE STATUS OF UTILIZATION AND MANAGEMENT SCENARIOS SKIPJACK ( Katsuwonus pelamis L.) CAUGHT IN THE BOLAANG MONGONDOW BEACH, NORTH SULAWESI

## ABSTRACT
Skipjack (Katsuwonus pelamis L.), needs to be  managed well as a renewable natural resources, but  vunerable  to  depletion  or  extinction.  One  approach  in  fish  resource  management  is  by modeling. Analysis were carried out, aimed to get the best model for surplus production model, to assess the Maximum Sustainable Yield (MSY), and to get the skipjack utilization and effort level. Data  used  for  surplus  production  model  were  collected  from  skipjack  landings  data  provided  by Marine and Fisheries Service (Dinas Kelautan dan Perikanan) of  Bolaang-Mongondow Regency and North Sulawesi Province. The best Surplus Production Model that  can be used to assess  the skipjack potential yield is the Fox model. The optimum effort per year biologically E MSY is 3,876 trips.  The  optimum  yield  biologically  C MSY is  1,505.00  tons  per  year.  Utilization  level  for  2005 was 124.48 %, with effort level 110.96 %. The catch for 2005 had already above the MSY value, which shows an  overfishing.
Keywords : Skipjack, Surplus Production Model, Maximum Sustainable Yield, Bolaang Mongondow

## I.   PENDAHULUAN
Ikan  cakalang  merupakan  produk utama/unggulan  di  Sulawesi  Utara,  dan sebagai  komoditi  ekspor.    Total  produksi cakalang  di  Sulawesi  Utara  (termasuk perairan  Bolaang-Mongondow  untuk  tahun 2010  mencapai  60.190,3  ton,  dengan  nilai
milyaran  rupiah  (Anonim,  2011).    Penelitian tentang  ikan  cakalang  umumnya  membahas tentang  eksploitasi  untuk  meningkatkan produksi,  belum  banyak  yang  meneliti tentang  status  pemanfaatan  (termasuk  aspek kelestarian dan efisiensi) sumberdaya.
Penangkapan  cakalang  di  perairan Bolaang-Mongondow telah berlangsung
cukup  lama,  dengan  intensitas  yang  padat. Data  mengenai  tingkat  pemanfaatan  suatu sumberdaya ikan sangat penting, karena akan menentukan apakah pemanfaatan sumberdaya tersebut  kurang  optimal,  optimal,  atau berlebih. Pemanfaatan sumberdaya ikan yang berlebihan  akan  mengancam  kelestariannya . Dengan  mengetahui  tingkat  pemanfaatan sumberdaya  cakalang,  diharapkan  dapat dilakukan  pengelolaan  yang  terencana  dan lestari.
Model  yang  paling  sederhana  dalam dinamika populasi  ikan  ialah model produksi surplus, dengan memperlakukan ikan sebagai biomassa tunggal yang tak dapat dibagi, yang tunduk pada aturan-aturan sederhana kenaikan  dan  penurunan  biomassa.    Model ini, pada umumnya digunakan dalam penilaian stok ikan hanya dengan menggunakan data hasil tangkapan dan upaya tangkap yang umumnya tersedia.
Penelitian ini bertujuan untuk mendapatkan Model Produksi Surplus terbaik,  serta  mengetahui  berapa  besar  hasil tangkapan  maksimum  lestari  (MSY),  tingkat pemanfaatan,  dan  tingkat  pengusahaan  ikan cakalang di perairan Bolaang-Mongondow.

## Tujuan dan Manfaat Penelitian
Penelitian ini bertujuan untuk :
1. Mendapatkan  model  Model  Produksi Surplus  terbaik  untuk  diterapkan  pada perikanan  cakalang  di  perairan  Bolaang- Mongondow.
2. Mengetahui berapa besar hasil tangkapan maksimum lestari (MSY), tingkat pemanfaatan,  dan  tingkat  pengusahaan ikan cakalang.
3. Menetapkan alokasi upaya-tangkap (jumlah trip kapal) yang paling menguntungkan untuk daerah penelitian. Penelitian  ini  diharapkan  bermanfaat dalam:
(1).    Memberi  informasi    tentang  batas maksimum  jumlah  ikan  cakalang  yang  dapat ditangkap dan upaya-tangkap yang digunakan per  tahun,  agar  sumberdaya  ikan  cakalang tetap lestari maupun yang paling menguntungkan  secara  ekonomi;  dan  (2)
Memberi informasi tentang tingkat
pemanfaatan  dan  pengusahaan  perikanan cakalang  di  perairan  Bolaang-Mongondow,
apakah sudah terjadi tangkap-lebih ( overfishing ) atau tidak.

## II.   TINJAUAN PUSTAKA

## Model Produksi Surplus
Model  yang  paling  sederhana  dalam dinamika populasi  ikan  ialah model produksi surplus  yang  memperlakukan  populasi  ikan sebagai  biomassa  tunggal  yang  tidak  dapat dibagi,  yang  tunduk  pada  aturan-aturan sederhana  dari  kenaikan  dan  penurunannya. Model produksi ini tergantung pada 4 macam besaran, yaitu : biomassa populasi pada suatu waktu  tertentu  t  (B t ),  tangkapan  untuk  suatu waktu  tertentu  t  (C t ),  upaya  tangkap  pada waktu  tertentu  t  (E t ),  dan  laju  pertumbuhan alami  konstan  (r)  (Boer  dan    Aziz,    1995). Model  ini  pertama  kali  dikembangkan  oleh Schaefer,  yang  bentuk  awalnya  sama  dengan model pertumbuhan logistik.
Menurut Coppola dan Pascoe (1996), persamaan  surplus  produksi  terdiri  dari beberapa  konstanta  yang  dipengaruhi  oleh pertumbuhan alami, kemampuan alat tangkap,  dan  daya  dukung  lingkungan. Konstanta-konstanta  tersebut  diduga  dengan menggunakan model-model penduga parameter  biologi  dari  persamaan  surplus produksi, yaitu model : Equilibrium Schaefer, Disequilibrium Schaefer, Schnute, dan
Walter  –  Hilborn.    Berdasarkan  keempat model tersebut dipilih yang paling sesuai atau best fit dari pendugaan yang lain.
Menurut  Sparre  dan  Venema  (1999), rumus-rumus  model  produksi  surplus  hanya berlaku  apabila  parameter slope (b)  bernilai negatif,  yang  berarti  penambahan  upaya tangkap  akan  menyebabkan  penurunan  hasil tangkapan  per  upaya  tangkap.    Apabila parameter  b  nilainya  positip,  maka  tidak dapat  dilakukan  pendugaan  besarnya  stok maupun  upaya  optimum,  tetapi  hanya  dapat disimpulkan  bahwa  penambahan upaya tangkap masih memungkinkan untuk meningkatkan hasil tangkapan. Pendugaan upaya penangkapan optimum (E opt ) dan hasil tangkapan maksimum  lestari  (C MSY )  didekati  dengan Model  Produksi  Surplus.  Antara  hasil tangkapan  per  satuan  upaya  (CPUE t )  dan upaya tangkap ( effort ) dapat berupa hubungan linear maupun eksponensial (Gulland,  1983).    Model  Produksi  Surplus terdiri  dari  2  model  dasar  yaitu  Model Schaefer  (hubungan  linear)  dan  Model Gompertz  yang  dikembangkan  oleh  Fox
dengan bentuk hubungan eksponensial (Gulland, 1983).

## Model  Schaefer
Model  produksi  surplus  pertama  kali dikembangkan  oleh  Schaefer,  yang  bentuk awalnya  sama  dengan  model  pertumbuhan logistik. Model tersebut ialah sebagai berikut:
dt dB t =G( B t ) = r B t (1- K B t ) (1)
Persamaan  ini  belum  memperhitungkan pengaruh  penangkapan,  sehingga  Schaefer menuliskan kembali menjadi:
dt dB t =  r B t ( 1- K B t ) -  C t (2)
dengan  K  ialah  daya  dukung  lingkungan perairan,  dan  C t ialah  tangkapan  yang  dapat ditulis sebagai :
C t = q E t B t (3)
dengan  q    sebagai  koefisien  ketertangkapan ( catchability ),  dan  E t menunjukkan  upaya tangkap. Persamaan ini dapat ditulis menjadi:
t t E C =  q  B t =   CPUE (4)
Dari persamaan diferensial (2), tangkapan  optimum  dapat  dihitung  pada  saat
dt dB t =  0    atau  disebut  juga  penyelesaian pada  titik  keseimbangan  ( equilibrium ),  yang berbentuk  :
r B t ( 1 - K B t ) -  C t =   0 ,     atau
C t =  r B t ( 1 - K B t ) =  q E t B t (5)
Dari persamaan (3) dan (5) diperoleh nilai B t sebagai berikut :
B t =   K ( 1 - r qE t ) (6)
Sehingga persamaan (5) menjadi :
C t =  q  K E t ( 1 - r qE t ) = q  K E t - r K q 2 E t 2 (7)
Persamaan  (7)  disederhanakan  lagi  oleh Schaefer menjadi :
t t E C =  a   -  b  E t , atau C t =  a E t -  b E t 2 (8)
sedangkan  a  =  q  K  dan    b  = r K q 2 .
Hubungan  linear  ini  yang  digunakan  secara luas untuk menghitung C MSY melalui penentuan  turunan  pertama  C t terhadap  E t untuk  mencari  solusi  optimal,  baik  untuk tangkapan  maupun  upaya  tangkap.    Turunan pertama C t terhadap E t adalah:
t t dE dC =  a    -    2b  E t ,    sehingga  diperoleh dugaan  E opt (upaya  tangkap  optimum)  dan C MSY (tangkapan  maksimum  lestari)  masing- masing :
E opt = b a 2 = q r 2 (9)
dengan  memasukkan  nilai  E opt pada persamaan  (8),  akan  diperoleh  C MSY sebagai berikut :
C MSY =  a E t -  b E t 2 =  a ( b a 2 ) -  b ( b a 2 ) 2 = b a 4 2
dengan mensubstitusi a = qK dan  b = r K q 2 ,
akan diperoleh
C MSY = b a 4 2 = r K q K q / 4 2 2 2 = 4 rK (10)
Nilai-nilai  a  dan  b  diduga  melalui pendekatan  metode  kuadrat  terkecil  yang umum  digunakan  untuk  menduga  koefisien persamaan  regresi  sederhana.    Selanjutnya, dengan memasukkan nilai E opt pada persamaan (6) diperoleh   biomassa optimum (B MSY ) sebagai berikut :
B MSY =  K  - opt E r Kq =  K  - r Kq ( q r 2 ) =  K  - 2 K = 2 K (11)
Nilai-nilai  parameter  q,  K,  dan  r  dapat dihitung  dengan  menggunakan  algoritma Fox,  seperti  yang  diacu  dalam  Sularso (2005), sebagai berikut :
q t = ln ) /( 1 / 1 1 1 1 z b zU b zU t t                        (12)
dimana  z  =  -(a/b)/E*,    E*  =  (E t +E t+1 )/2,
t t t E C U  ,  dan  nilai  q  adalah  rata-rata
geometrik dari nilai q t .   Dari nilai a, b, dan q, selanjutnya dapat dihitung nilai K dan r.
Model   Fox Model Fox (1970) memiliki
beberapa  karakteristik  yang  berbeda  dari model Schaefer, yaitu pertumbuhan biomassa mengikuti  model  pertumbuhan  Gompertz. Penurunan CPUE terhadap upaya tangkap (E) mengikuti pola eksponensial negatif. C t =  E t . exp(a - b E t ) (13)
Upaya optimum diperoleh dengan menyamakan  turunan  pertama  C t terhadap  E t sama dengan nol.
E opt = b 1 (14)
Hasil  tangkapan  maksimum  lestari (C MSY )  didapat  dengan  memasukkan  nilai upaya  optimum  ke  dalam  persamaan  (13), dan diperoleh :
C MSY = b 1 e a-1 (15)

## Model  Schnute
Schnute  (1977),  mengemukakan  versi  lain dari  model  produksi  surplus  yang  bersifat dinamis serta deterministik.  Metode Schnute dianggap  sebagai  modifikasi  dari  model Schaefer  dalam  bentuk  diskret  (Roff,    1983, yang diacu dalam Tinungki 2005).
ln( t t U U 1  )   =  r - qK r ( 2 1   t t U U ) - q ( 2 1   t t E E )=  a  -  b  (
2 1   t t U U ) - c ( 2 1   t t E E ) (16)
dimana a = r,  b = qK r ,  dan  c = q,  adalah penduga parameter koefisien regresi berganda.
Model  Walter - Hilborn Walter  dan  Hilborn  (1976)  yang diacu dalam Tinungki (2005), mengembangkan  jenis  lain  dari  model produksi surplus, yang dikenal sebagai model regresi.    Model  Walter  –  Hilborn  ini,
menggunakan persamaan diferensial sederhana,dengan persamaan sebagai berikut:
t t U U 1  - 1  =   r  - Kq r U t - q E t = a - b U t -  c  E t (17)
dimana a = r,  b = Kq r ,  dan  c = q,  adalah
penduga parameter koefisien regresi berganda.

## Model  Clarke Yoshimoto Pooley (CYP)
Pendugaan parameter biologi untuk model  produksi  surplus  dapat  pula  dilakukan melalui teknik pendugaan yang dikemukakan oleh  Clarke,  Yoshimoto,  dan  Pooley  (Fauzi dan Anna 2005,  Tinungki 2005).  Parameter- parameter  yang  diduga  ialah  r,  K,  dan  q, dengan  model  yang  dinyatakan  sebagai berikut :
) ( 2 ) ln( 2 2 ) ln( 2 2 ) ln( 1 1                 t t t t E E r q U r r qK r r U (18)
di  mana  : r r a   2 2 ' , ) ln( ' qK a a  ,
r r b    2 2 , r q c   2
dengan  demikian  persamaan  (18)  dapat ditulis dalam bentuk :
) ( ) ln( ) ln( ' ) ln( 1 1       t t t t E E c U b qK a U ) ( ) ln( 1      t t t E E c U b a (19)

## III.   METODE  PENELITIAN

## Sumber Data
Data  yang  dikumpulkan  ialah  data sekunder  tentang  hasil  tangkapan  ikan cakalang  di  perairan  Bolaang-Mongondow dan  sekitarnya.    Data  produksi  dan  upaya tangkap  dikumpul  dari  Dinas  Kelautan  dan Perikanan  Kabupaten  Bolaang-Mongondow dan    Provinsi  Sulawesi  Utara  selama  tahun 1991 - 2010.
Data (variabel) yang digunakan untuk analisis  model  produksi  surplus  ialah  data hasil  tangkapan  (C t )  per  tahun  dan  upaya tangkap  (E t )  per  tahun,  serta  CPUE  ( ) t t E C .
Data (variabel) yang digunakan untuk analisis model  produksi  surplus  ialah  sebagai  berikut :
1. Hasil tangkapan  (C t )     :   berat ikan yang didaratkan (ton) pada tahun ke t 2. Upaya tangkap (E t ) :   jumlah kapal motor penangkap ikan yang mendaratkan
hasilnya  di  tempat  pendaratan  (trip)  pada tahun ke t
3. t t E C :    C t dibagi E t (ton /trip) pada tahun

## ke t

## Metode Analisis Data
Model  penduga  yang  dianalisis  dan dievaluasi  ialah  :  model  Schaefer,    model Fox,  model  Schnute,  model  Walter-Hilborn, dan  model  Clarke  Yoshimoto  Pooley  (CYP). Berdasarkan  hasil  evaluasi  secara  statistika, akan  diperoleh  suatu  model  yang  “terbaik” sebagai penduga.  Dari model terbaik tersebut
dapat  dihitung  nilai  C MSY ,  upaya  tangkap optimum,  tingkat  pemanfaatan,  dan  tingkat pengusahaan      sumberdaya  perikanan cakalang.

## IV.   HASIL DAN  PEMBAHASAN
Hasil  tangkapan  perikanan  cakalang di  perairan  Bolaang-Mongondow  dari  tahun ke  tahun  walaupun  berfluktuasi  tetapi cenderung  meningkat.    Data  hasil  tangkapan selang  tahun  1991-2010,  diringkas  dalam Tabel 1.
Tabel 1  Jumlah Hasil Tangkapan, upaya tangkap, dan CPUE ikan cakalang di perairan Bolaang-Mongondow  tahun 1991-2010
Tahun Tangkapan (ton) Upaya (trip) CPUE (ton/trip 1991 698,2 2.455 0,2844 1992 826,0 8.650 0,0955 1993 866,2 6.195 0,1398 1994 901,6 5.347 0,1686 1995 875,4 6.272 0,1396 1996 763,9 8.560 0,0892 1997 899,6 9.948 0,0904 1998 1.218,3 9.950 0,1224 1999 1.268,4 4.037 0,3142 2000 1250,3 3.027 0,4130 2001 1.304,3 1.102 1,1836 2002 1.461,3 3.366 0,4341 2003 2.137,1 4.368 0,4893 2004 1.995,7 6.464 0,3087 2005 1.873,4 4.301 0,4356 2006 1.203,7 1.020 1,1801 2007 1.246,9 7.544 0,1653 2008 1.369,0 6.410 0,2136 2009 1.379,0 6.400 0,2155 2010 1.631,0 6.452 0,2528 Rata-rata 1.258,465 5.593,40 0,3368 Sumber :  Diolah dari data Dinas Kelautan dan Perikanan  Kabupaten Bolaang- Mongondow Provinsi Sulawesi Utara Model  Schaefer
Dari hasil analisis diperoleh persamaan regresi
t t E C = 0,868 – 0,0000949
E t , dengan  nilai  koefisien  determinasi  (R 2 )  = 0,632  dan  tingkat  signifikansi  p  <  0,001. Dengan  demikian  model  produksi  penduga hasil  tangkapan  untuk  model  Schaefer  sesuai
persamaan  (8)  ialah  : C t =  0,886  E t – 0,0000949 E t 2 .

## Model  Fox
Dari hasil analisis regresi  diperoleh persamaan regresi :
Ln t t E C = 0,054 – 0,000258 E t , 58dengan R 2
= 0,808 (p < 0,001). Penduga hasil tangkapan untuk  model  Fox  sesuai  persamaan  (13) ialah: C t =  E t . e (0,054 – 0,000258 Et) .

## Model  Schnute
Untuk  metode  Schnute  diperoleh persamaan  regresi Ln(U t+1 /U t )  =  0,969  – 1,121  (U t+1 +U t )/2  -  0,00009201  (E t+1 +E t )/2 dengan  R 2 =  0,034,  dan  semua  koefisen regresi tidak ada yang signifikan  (p>0,05).
Model   Walter - Hilborn
Pada metode Walter–Hilborn diperoleh  persamaan  regresi (U t+1 /U t )-1  = 0,160 – 0,525 U t + 0,00004725 E t dengan R 2 =  0,135,  dengan  semua  variabel  tidak signifikan (p>0,05).
Model Clarke Yoshimoto Pooley (CYP)
Pada metode CYP diperoleh persamaan  regresi Ln(U t+1 )  =  0,265  -  0,357 ln(U t )  -  0,000192  (E t +E t+1 ) dengan  R 2 = 0,648,  dan  tidak  semua  koefisen  regresi signifikan.

## PEMBAHASAN
Hasil  perhitungan  dan  analisis  untuk kelima  model  produksdi  surplus,  disajikan pada Tabel 2.
Tabel 2.   Hasil perhitungan dan validasi model produksi surplus
Model Schaefer Model Fox Model Schnute Model Walter- Hilborn Model  CYP Kesesuaian Tanda Sesuai Sesuai Sesuai Tidak Sesuai Tidak Sesuai Nilai  R 2 0,632 0,808 0,034 0,135 0,648 Nilai Validasi 0,5706 0,2893 0,6462 2,6428 1,0292 Signifikansi Koefisien Regresi Individu Signifikan Signifikan Tidak Signifikan Tidak Signifikan Tidak Signifikan
Dari  hasil  perhitungan  pada  Tabel  2, terlihat bahwa yang paling sesuai ialah model Fox dengan nilai R 2 cukup besar (R 2 = 0,808) dan validasi (nilai residual) paling kecil.  Dari model Fox diperoleh nilai a = 0,054 dan nilai b = 0,000258, dengan persamaan (14) dan (15 dapat  dihitung  nilai  Upaya  optimum  (E opt ) dan  Tangkapan  maksimum  lestari  (C MSY ) sebagai berikut :
E opt =    9 96 , 3875 000258 , 0 1 1 b 3876
trip/tahun.
C MSY = b 1 e a-1 = 000258 , 0 1 e 0,054  -  1
= 1505,00  ton/tahun. Ini  berarti  bahwa  untuk  menjaga  kelestarian sumberdaya perikanan cakalang secara teknis dan  biologis,  dalam  setahun  jumlah  unit
penangkapan  tidak  boleh  melebihi  3876  trip. Untuk  menjaga  kelestarian  sumberdaya  ikan cakalang  di  perairan  Bolaang-Mongondow, maksimum ikan yang dapat ditangkap sebesar 1505,00 ton/tahun.
Selanjutnya  dari  nilai  E opt dan  C MSY dapat  dihitung  tingkat  upaya  penangkapan dan tingkat pemanfaatan ikan cakalang untuk tahun  tertentu  misalkan  tahun  2005,  sebagai berikut :
Tingkat  upaya tahun 2005=
opt E E 2005 x 100%
= 3876 4301 x   100%    =   110,96 %
Tingkat pemanfaatan tahun 2005
= MSY C C 2005 x100%
=
00 , 1505 4 , 1873 x   100%   =  124,48 %.
Dari  hasil  perhitungan,  ternyata  upaya tangkap  ikan  cakalang  di  perairan  Bolaang- Mongondow pada tahun 2005 sudah melebihi tingkat  upaya  maksimum  lestari,  demikian pula  tingkat  pemanfaatannya  sudah  melebihi 100  %.  Hal  ini  menunjukkan  bahwa  untuk tahun 2005 hasil tangkapan sudah melampaui batas  tangkapan  maksimum  lestari  atau terjadi overfishing (tangkap lebih).
Penelitian ini menjelaskan penggunaan beberapa kriteria statistika dalam memilih  model  produksi  surplus  terbaik.
Dengan menerapkan beberapa kriteria statistika  dalam  memilih  model  produksi surplus, akan diperoleh hasil yang lebih baik. Para peneliti di bidang perikanan mendapatkan  pedoman  dalam  menetapkan kriteria  pemilihan  model  produksi  surplus, sekaligus juga menghindari penerapan langsung satu model saja dalam menganalisis model produksi surplus di suatu perairan.
Karena adanya indikasi terjadi overfishing (tangkap  lebih),  disarankan  agar pemerintah  menerapkan  peraturan  mengenai pembatasan  ukuran  ikan  cakalang  yang  akan ditangkap.  Peraturan  tersebut  harus  diawasi pelaksanaannya  melalui  sistem monitoring, controlling,  dan  surveillance (MCS)  yang tegas.  Selain  itu,  untuk  mencegah  tingkat pemanfaatan  yang  mengganggu  kelestarian stok,  perlu  dilakukan  pengaturan  zonasi penangkapan  melalui  alokasi  penggunaan rumpon  yang  menjadi  daerah  penangkapan ikan cakalang.

## VI.  KESIMPULAN DAN SARAN
Kesimpulan 1.  Model    Produksi  Surplus  yang  dapat digunakan untuk menelaah hasil tangkapan  ikan  cakalang  di  perairan Bolaang-Mongondow  ialah  Model      Fox, dengan  persamaan    : C t =    E t .  e (0,054  – 0,000258 Et) .
2.    Hasil  tangkapan  maksimum  lestari  ikan cakalang  C MSY sebesar  1505,00  ton  per tahun,  diperoleh  pada  tingkat  upaya tangkap  E MSY 3876  trip.      Untuk  tahun 2005 besarnya tingkat pemanfaatan sebesar  124,48  %  dan  ini  tergolong overfishing (tangkap-lebih),
dengan tingkat  pengupayaan  sebesar  110,96  %
juga sudah melampaui batas nilai optimum.
Saran 1. Dalam menerapkan model produksi surplus  di  suatu  perairan,  tidak  hanya langsung menggunakan satu model tertentu saja, tetapi hendaknya menggunakan beberapa model yang dipilih  berdasarkan  kriteria  statistika. Kriteria  tersebut  menyangkut  antara  lain  : kesesuaian tanda, nilai koefisien determinasi  (R 2 ),  nilai  validasi,  dan signifikansi koefisien regresi.
2. Karena  ada  indikasi  terjadi overfishing (tangkap  lebih)  dan  besarnya  intensitas upaya  penangkapan  cakalang  di  perairan Bolaang-Mongondow, disarankan secepatnya  dilakukan  pengawasan  oleh pihak  berkompeten  untuk  menangani masalah ini. Terutama melakukan efisiensi terhadap upaya tangkap.
3. Perlu  adanya  perbaikan  sistem  penyajian dan  pengumpulan  data  perikanan  yang dilakukan  pemerintah,  agar  tersedia  data yang  sesuai  untuk  kebutuhan  analisis pengelolaan  sumberdaya.  Sebagai  contoh, untuk  data  upaya  tangkap  ( effort )  tidak hanya  dicatat  nilai  nominalnya  (jumlah kapal,  jumlah  alat,  dan  lain-lain),  akan tetapi  yang  lebih  penting  lagi  ialah  nilai efektifnya.    Misalnya  jumlah  hari-kapal ( boat-day )  yaitu  jumlah  kapal  dikali jumlah hari operasi.  Data seperti ini lebih dapat  menggambarkan  efektifitas  upaya tangkap,  dibanding  dengan  hanya  jumlah kapal  atau  jumlah  alat  tangkap  yang digunakan  yang  belum  tentu  semuanya efektif.

## DAFTAR  PUSTAKA
Anonim.  2011.  Buku  Tahunan  Statistik Perikanan  Tangkap  Sulawesi  Utara Tahun 2010. Manado.
Aziz, K. A.  1989.  Pendugaan Stok Populasi Ikan  Ikan  Tropis.    Bahan  Pengajaran (Tidak  Dipublikasi).    Pusat  Antar Universitas Ilmu Hayat, Institut Pertanian Bogor,  Bogor.  89 hal.
Aziz,  K.A.  dan  B.  Wahyudi.    1997.  Sistem
Pemantauan dan Basis Data Pengkajian Stok.  Makalah  dalam  Potensi  dan Penyebaran  Sumberdaya  Ikan  Laut  di Perairan Indonesia. Proyek Pengembangan dan Pemanfaatan Sumberdaya Perikanan Laut Direktorat Jenderal Perikanan, Departemen Pertanian,  Jakarta.  Halaman 13 – 32.
Boer,  M.,  dan  K.A.    Azis.    1995 . Prinsip- prinsip Dasar Pengelolaan Sumberdaya
Perikanan  Melalui  Pendekatan  Bio- Ekonomi. Jurnal  Ilmu-ilmu  Perairan dan Perikanan: 3(2):109-119.
Boer, M., K.A.  Azis, J. Widodo, A. Djamali, A.  Ghofar,  dan  R.  Kurnia.    2001 . Potensi,  Pemanfaatan  dan  Peluang
Pengembangan  Sumberdaya Ikan  Laut di  Perairan  Indonesia.    Pusat  Kajian Sumberdaya Pesisir dan Lautan, Institut Pertanian Bogor.
Coppola G., and S. Pascoe.  1996. A Surplus
Production  Model  with  a  non-linear Catch-Effort  Relationship .    (Research Paper  105)  Center  for  the  Economics and Managemant of Aquatic Resources University of Portsmouth.
FAO.    1994. Review  of  the  State  of  World Marine  Fisheries  Resources .    FAO Fisheries Technical Paper.  Rome.  136 p.
Fauzi,  A.  2002. A  Note  Surplus  Production Model .  IPB,  Darmaga Bogor.
Fauzi,  A.,  dan  S.  Anna.    2005.    Pemodelan
Sumberdaya  Perikanan  dan  Kelautan untuk Analisis Kebijakan. PT.
Gramedia Pustaka Utama,  Jakarta.
Fox,  W.W.    1970.    An  Exponential  Surplus Yield  Model  for  Optimazing  Exploited Fish  Population. Trans.  Am.  Fish  Soc . 99(1):80-88.
Gafa, B., I.G.S. Merta, H.R. Barus, dan E.M. Amin.    1993.    Penurunan  Hasil Tangkapan Ikan Tuna dan Cakalang  di Perairan  Sulawesi  Utara  dan  Faktor- faktor  yang  Mempengaruhinya. Jur. Pen. Perikanan Laut 72 : 11 – 19.
Gulland,  J.A.    1983.    Fishing  and    Stock  of
Fish  at  Iceland.    Mui.` Agric.  Fish
Food, Invest . (Ser.2) 23(4): 52 – 70.
Kekenusa,  J.S.    2006.  Pemodelan  Hasil
Tangkapan dan Evaluasi Model Produksi  Surplus  Ikan  Cakalang  yang Tertangkap  di  Perairan  Sekitar  Bitung Provinsi  Sulawesi  Utara.    Disertasi (Tidak  Dipublikasikan).    Program Pascasarjana  Universitas  Airlangga. Surabaya.  140 hal.
Kekenusa , J.S . 2007.  Analisis  Bio-ekonomi Ikan  Cakalang  yang  Tertangkap  di Perairan Sekitar Bitung Provinsi
Sulawesi  Utara. Pacific  Journal Vol.2 No.1 :71-76.
Kekenusa,  J.S.,  V.N.R.  Watung,  Dj.  Hatidja, dan  A.J.  Rindengan.  2008.  Penentuan Status  Pemanfaatan  dan  Skenario Pengelolaan
Ikan Cakalang
( Katsuwonus pelamis ) yang Tertangkap di  Perairan  Sulawesi  Utara.    Laporan Penelitian Hibah Bersaing.
Kekenusa,  J.S.,  V.N.R.  Watung,  dan  Dj. Hatidja. 2009. Penentuan Status Pemanfaatan dan Skenario Pengelolaan Ikan  Cakalang  ( Katsuwonus  pelamis ) yang  Tertangkap  di  Perairan  Sangihe- Talaud.    Laporan  Kegiatan  Penelitian Strategis Nasional Tahun 2009.
Kusnawan, I. 1999. Analisis Tingkat Pengusahaan  dan  Musim  Penangkapan Ikan  Cakalang  ( Katsuwonus  pelamis ) dengan  Pendekatan  Hasil  Tangkapan yang Didaratkan di Pelabuhan Perikanan  Pelabuhan  Ratu.  Skripsi (Tidak Dipublikasikan).  Program Studi Ilmu  Kelautan,  Fakultas  Perikanan  dan Ilmu Kelautan Institut Pertanian Bogor.
Bogor.  77 hal.
Menard,  F.,  A.  Fonteneau,  D.  Gaertner,  V. Nordstrom,  B.  Stequert,  and  E. Marchal.   2000.  Exploitation of Small Tunas  by  a  Purse-Seine  Fishery  With Fish  Aggregating  devices  and  Their Feeding Ecology in an Eastern Tropical Atlantic  Ecosystem. ICES  Journal  of
Marine Science 57:525-530. Meyer, W.J. 1987. Concepts of Mathematical  Modelling .    McGraw- Hill Inc. New York.  439p.
Monintja,  D.  R.,  dan  R.  Yusfiandayani.
1999.    Teknologi  Penangkapan  Ikan Cakalang  dan  Tuna.    Laboratorium Teknologi  Penangkapan  Ikan,  FPIK-
IPB.  Bogor.  27 hal.
Monintja,  D.  R.,  dan  R.  Yusfiandayani. 2001. Pemanfaatan Sumberdaya Pesisir Dalam Bidang Perikanan Tangkap. Prosiding Pelatihan Pengelolaan  Wilayah  Pesisir  Terpadu. IPB,  Bogor. Monintja,  D.  R.,  dan  Zulkarnain.    1995. Analisis Dampak Pengoperasian Rumpon  Tipe  Philipine  di  Perairan ZEE  terhadap  Perikanan  Cakalang  di Perairan  Teritorian  Selatan  Jawa  dan
Utara  Sulawesi.    Laporan  Penelitian. Fakultas  Perikanan  Institut  Pertanian Bogor.  Bogor.  70 hal.
Purwanto.  1988.  Bio-Ekonomi Penangkapan Ikan : Model Statik. Oseana. Vol. XIII No. 2 .  Departemen Pertanian,  Jakarta.
Schnute,   J. 1977.  Improved Estimates from the  Schaefer  Production  Models  : Theoretical  Considerations: J.  Fish. Res. Board Can ., 34:583-663.
Sparre  ,  P.  and  S.C.  Venema.    1999. Introduksi Pengkajian Stok Ikan Tropis.    Buku  1  Manual.  (Terjemahan J.  Widodo.  I.G.S.  Merta,  S.  Nurhakim, dan  M.  Badrudin).    Pusat  Penelitian dan  Pengembangan  Perikanan,  Badan Penelitian dan Pengembangan Pertanian (Kerjasama dengan Organisasi Pangan dan Pertanian
Perserikatan  Bangsa-bangsa).    Jakarta. 438 hal. Suharman,  I.    1981.    Suatu  Penelitian Tentang Penangkapan Skipjack dengan  Purse  Seine  di  Perairan  Utara Irian  Jaya.    Karya  Ilmiah,    Fakultas Perikanan  Institut  Pertanian  Bogor, Bogor. Sularso,  A.  2005.  Alternatif  Pengelolaan Perikanan  Udang  di  Laut  Arafura. Disertasi (Tidak Dipublikasikan). Sekolah  Pascasarjana  Institut  Pertanian Bogor.  Bogor.  130 hal. Tinungki,  G.  M.    2005.    Evaluasi  Model Produksi Surplus dalam Menduga Hasil Tangkapan  Maksimum  Lestari  untuk Menunjang Pengelolaan Perikanan Lemuru Di Selat Bali. Disertasi (Tidak Dipublikasikan).  Sekolah Pascasarjana Institut  Pertanian  Bogor.    Bogor.    207
hal.
Uktolseja,  J.C.B.    1997.    Laporan  Penelitian
Indeks  Kelimpahan  Ikan  Tuna  dan Cakalang  di  Sekitar  Rumpon  (Tidak Diterbitkan). Balai Penelitian Perikanan Laut,  Jakarta.  29 hal.
Uktolseja, J.C.B,  B. Gafa,  R. Purwasasmita, dan  B.  Iskandar.    1997.    Sumberdaya Ikan  Pelagis  Besar.    Makalah dalam Potensi  dan  Penyebaran  Sumberdaya Ikan  Laut  di  Perairan  Indonesia.
Proyek Pengembangan dan Pemanfaatan  Sumberdaya  Perikanan Laut  Direktorat  Jenderal  Perikanan
Laut, Departemen Pertanian,    Jakarta. Halaman 33 – 54.
Uktolseja,  J.C.B,  R.  Purwasasmita,  K. Susanto,  dan  A.B.  Sulistiadji.    1998.
Sumberdaya Ikan Pelagis Besar. Makalah
dalam
Potensi dan Penyebaran  Sumberdaya  Ikan  Laut  di Perairan  Indonesia.    Komisi  Nasional Pengkajian Stok Sumberdaya Ikan Laut LIPI,  Jakarta.  Halaman 40 – 88.
Wahyono,  U.    1999.    Potensi  Perikanan Indonesia.    Makalah dalam Prosiding Seminar  Validasi  Data  Inderaja  untuk Bidang  Perikanan  14  April  1999  di Jakarta.  BPPT Jakarta, hal. 18-32. Widodo,    J.    1987. Modified  Surplus Production Methods of Gulland (1961), and Schnute (1977) .  A Serial Seminars Published by Oceana XII(2):119-130. Zar,    J.H.    1984. Biostatistical    Analysis . Prentice-Hall,  New Jersey.