Biến động nguồn lợi và một số đặc điểm sinh học cá phèn khoai (Upeneus japonicus Houttuyn, 1782) ở vùng đánh cá chung vịnh Bắc Bộ Việt Nam - Trung Quốc (giai đoạn 2013-2015)

Cá phèn khoai (Upeneus japonicus) là một trong những đối tượng hải sản khai thác chính bằng nghề lưới kéo đáy và nghề rê đáy ở vùng đánh cá chung vịnh Bắc Bộ. Cá phèn khoai là loài đồng sinh trưởng, hệ số đồng hóa (b) dao động trong khoảng từ 2,07 đối với cá đực đến 2,10 với cá giống cái và hệ số dị hóa (a) dao động trong khoảng 0,013×104 - 0,026×104.

Tạp chí Khoa học Cơng nghệ Biển; Tập 17, Số 2; 2017: 186-197 DOI: 10.15625/1859-3097/17/2/8052 http://www.vjs.ac.vn/index.php/jmst BIẾN ĐỘNG NGUỒN LỢI VÀ MỘT SỐ ĐẶC ĐIỂM SINH HỌC CÁ PHÈN KHOAI (Upeneus japonicus Houttuyn, 1782) Ở VÙNG ĐÁNH CÁ CHUNG VỊNH BẮC BỘ VIỆT NAM - TRUNG QUỐC (GIAI ĐOẠN 2013 - 2015) Mai Công Nhuận*, Nguyễn Khắc Bát Viện Nghiên cứu Hải sản, Bộ Nông nghiệp Phát triển Nơng thơn * E-mail: mcnhuan@rimf.org.vn Ngày nhận bài: 4-4-2016 TĨM TẮT: Cá phèn khoai (Upeneus japonicus) đối tượng hải sản khai thác nghề lưới kéo đáy nghề rê đáy vùng đánh cá chung vịnh Bắc Bộ Cá phèn khoai loài đồng sinh trưởng, hệ số đồng hóa (b) dao động khoảng từ 2,07 cá đực đến 2,10 với cá giống hệ số dị hóa (a) dao động khoảng 0,013×104 - 0,026×104 Chiều dài thành thục Lm50 cá 126 mm cá đực 116 mm Trong khoảng thời gian nghiên cứu từ năm 2013 đến năm 2015 suất khai thác trung bình dao động từ 0,12 kg/h đến 0,40 kg/h Mật độ phân bố trung bình dao động từ 4,22 kg/km2 đến 12,40 kg/km2 Trữ lượng nguồn lợi tức thời lồi cá ước tính dao động từ 1.134 đến 3.360 tấn, trung bình khoảng 2.271 khả khai thác cho phép khoảng 1.135 Từ khóa: Năng suất, mật độ, trữ lượng, kéo đáy, rê đáy MỞ ĐẦU Cá phèn khoai (Upeneus japonicus) thuộc họ cá phèn (Mullidae) loài cá phân bố khu vực Tây Thái Bình Dương: Từ Đài Loan, Hồng Kơng, Nam Triều Tiên đến vùng biển Nam Trung Quốc phía nam vùng biển Nhật Bản Ở Việt Nam, cá phèn khoai phân bố ven biển từ vịnh Bắc Bộ vùng biển Đông-Tây Nam Bộ chúng phân bố tập trung chủ yếu vùng biển vịnh Bắc Bộ [1] Chúng phân bố tương đối rộng dải độ sâu < 200 m chủ yếu độ sâu < 100 m Thức ăn chúng chủ yếu loại động vật cỡ nhỏ cá giáp xác Chiều dài lớn bắt gặp loài 20 cm (TL), phổ biến bắt gặp nhóm chiều dài từ - 14 cm (TL) [1] Đây lồi cá đáy có giá trị kinh tế cao đối tượng khai thác chủ yếu nghề lưới kéo đáy vịnh 186 Bắc Bộ Hiện nay, với gia tăng áp lực khai thác, trữ lượng nguồn lợi cá phèn khoai thể xu hướng suy giảm Vì vậy, nghiên cứu biến động quần thể lồi từ góp phần cung cấp sở khoa học cho việc phục hồi phát triển nguồn lợi loài cá cần thiết Bài viết sử dụng số liệu điều tra nguồn lợi thuộc dự án Điều tra liên hợp Việt - Trung với 12 chuyến điều tra, nói số liệu lồi cá phong phú có độ tin cậy cao PHƯƠNG PHÁP NGHIÊN CỨU Đối tượng nghiên cứu Đối tượng nghiên cứu loài cá phèn khoai (Upeneus japonicus) thuộc họ cá phèn (Mullidae) Vị trí phân loại loài cá phèn khoai hệ thống phân loại: Biến động nguồn lợi số đặc điểm… Lớp: Actionopteri Bộ: Perciformes Họ: Mullidae Giống: Upeneus Loài: Upeneus japonicus (Houttuyn,1782) Tên tiếng Anh: Bensasi goatfish Tên tiếng Việt: Cá phèn khoai Tên đồng danh: Mullus bensasi (Temmink & Schlegen,1843) Mullus japonicus (Houttuyn, 1782) Upeneoides bensasi (Temmink & Schlegen,1843) Upeneoides tokisensis (Doderlein, 1883) Upeneus bensasi (Temmink & Schlegen,1843) Cá phèn khoai lồi cá sống đáy, thể thn nhỏ, kích thước tối đa khoảng 15 cm (SL) Thân có màu hồng nhạt pha lẫn đen sậm Thùy vây có cá vạch trắng, vàng xen lẫn thùy vây khơng có Tia vây lưng thứ có tia cứng, tia vây ngực có 14 tia, số lượng mang 7+18 Số vảy đường bên từ 29 - 30 [2] Thời gian địa điểm nghiên cứu Trong thời gian nghiên cứu vùng biển đánh cá chung vịnh Bắc Bộ năm tiến hành chuyến điều tra vào tháng 1, 4, tháng 10 Hệ thống trạm nghiên cứu thiết kế cố định với thống nhóm chuyên gia nguồn lợi thủy sản nước Việt Nam Trung Quốc, gồm 35 trạm thực bên phía Việt Nam 30 trạm thực phía Trung Quốc Phía Việt Nam thực chuyến điều tra vào tháng tháng 10 tàu BV9262 TS, Trung Quốc thực chuyến điều tra vào tháng tàu Bắc Ngư Tàu ngư cụ sử dụng Số liệu thu thập sử dụng viết thực tàu nghiên cứu Tàu BV 9262 TS (phía Việt Nam) có cơng suất máy chính: 600 HP, ngư cụ sử dụng lưới kéo đơn, kích thước mắt lưới đụt 2a = 20 mm, chiều dài giềng phao: 26,48 m, chiều dài giềng chì 55,5 m Tàu Bắc Ngư (phía Trung Quốc) có cơng suất máy chính: 600 HP, ngư cụ sử dụng lưới kéo đơn, kích thước mắt lưới đụt 2a = 20 mm, chiều dài giềng phao: 37,7 m Bảng Số lượng mẻ lưới thời gian thực chuyến điều tra Chuyến điều tra Hình Sơ đồ trạm vị chuyến điều tra Số mẻ lưới 2013 2014 2015 Tổng Tháng 30 30 30 90 Tháng 35 35 35 105 Tháng 30 30 30 90 Tháng 10 35 35 35 105 Tổng 130 130 130 390 187 Mai Công Nhuận, Nguyễn Khắc Bát Thu thập số liệu Trữ lượng khả khai thác Sau kết thúc mẻ lưới, tiến hành phân tích số tiêu sinh học (Đo chiều dài cá thể, cân khối lượng, xác định độ chín muồi tuyến sinh dục) cá phèn khoai trường đồng thời thu mẫu mang phân tích phịng thí nghiệm Trữ lượng vùng biển nghiên cứu tính tổng trữ lượng dải độ sâu theo phương pháp Gulland (1971) [5] Khi cơng thức tính trữ lượng phương sai cho dải độ sâu là: Đo chiều dài toàn thân (FL, Folk Length) cân khối lượng cá thể (BW, Body Weight) theo hướng dẫn FAO (1995) Đơn vị đo chiều dài mm khối lượng cá thể gram Xác định giới tính độ chín muồi tuyến sinh dục theo thang bậc Nikolsky (1963) Cá có độ chín muồi tuyến sinh dục từ giai đoạn IV trở lên coi thành thục; từ giai đoạn III trở xuống coi chưa thành thục [3] STT Nội dung Đơn vị Số lượng Đo tần suất chiều dài Đo tương quan chiều dài khối lượng Phân tích độ chín muồi tuyến sinh dục Đánh giá nguồn lợi cá thể 1.850 cá thể 1.060 cá thể 996 mẻ lưới 390 Năng suất khai thác trung bình tính theo cơng thức (1) mật độ phân bố trung bình tính theo cơng thức (2) theo Pennington (1983): n n Ci  v w t i 1 i i i n Trong đó: C sản lượng đánh bắt mẻ lưới (kg); v vận tốc kéo lưới trung bình mẻ; w độ mở ngang lưới; t thời gian kéo lưới mẻ; n số trạm CPUA tính theo phương pháp Pennington (1983) [4] 188 MSY = 0,5*M*Bj (3) Phân tích tương quan chiều dài khối lượng (4) Trong đó: W: Khối lượng cá thể; L: Chiều dài cá thể; a: Hệ số sinh trưởng; b: Hệ số quan hệ Ước tính Lm50 Chiều dài Lm50 chiều dài có 50% số cá thể quần thể chín muồi tuyến sinh dục tham gia vào quần đàn sinh sản lần đầu Lm50 tính phương pháp hồi quy phi tuyến tính lặp theo cơng thức Jenning nnk., (2001) [6] P n CPUA( kg / km )  Khả khai thác tối ưu (MSY) tính theo cơng thức: (1) i 1 (2) Trong đó: Bj: Trữ lượng lồi; Aj: Diện tích vùng biển; CPUA j mật độ trung bình lồi; X1 hệ số lưới (X1=0,5) W = a × Lb Năng suất khai thác CPUE (kg/h) mật độ CPUA (tấn/km2) CPUE (kg / h)  X1 Tương quan chiều dài khối lượng cá phèn khoai ước tính phương pháp hồi quy lặp (Iterative Non-linear Regression), theo phương trình tuyến tính sau: Xử lý số liệu  Ci Aj Trong đó: MSY khả khai thác tối ưu; M: hệ số chết tự nhiên; Bj: trữ lượng loài Bảng Số lượng mẫu vật thu thập phân tích sử dụng nghiên cứu B j  CPUA j 1  exp  r   Lt  Lm 50    (5) Trong đó: Lm50: Chiều dài 50% cá thể quần đàn sinh sản lần đầu; P: Tỷ lệ cá thể thành thục chiều dài Lt; Lt: Chiều dài cá thể Ước tính trữ lượng quần đàn cá bố mẹ (SSB) Trữ lượng quần đàn cá bố mẹ (Spawning Stock Biomas - SSB) tính dựa tổng trữ lượng tức thời, tương quan chiều dài - khối Biến động nguồn lợi số đặc điểm… lượng, phân bố chiều dài tỷ lệ thành thục nhóm chiều dài Về bản, tổng trữ lượng quần đàn cá bố mẹ loài tính cơng thức: SSB   Bi * Lmi (6) i KẾT QUẢ VÀ THẢO LUẬN Năng suất khai thác Năng suất khai thác cá phèn khoai vùng đánh cá chung vịnh Bắc Bộ từ năm 2013 đến 2015 có biến động tương đối lớn qua năm trung bình đạt 0,26 kg/h Năng suất khai thác trung bình cao năm 2014 đạt 0,40 kg/h, năm 2013 đạt 0,26 kg/h thấp năm 2015 giảm xuống 0,12 kg/h (bảng 3) Bảng Năng suất khai thác (kg/h) cá phèn khoai vùng đánh cá chung vịnh Bắc Bộ qua tháng điều tra, giai đoạn 2013 - 2015 Chuyến điều tra Năm Trung bình Tháng Tháng Tháng Tháng 10 2013 - 0,26 0,04 0,50 0,26 2014 0,55 0,44 0,18 0,42 0,40 2015 0,01 0,09 0,04 0,36 0,12 Trung bình 0,28 0,26 0,09 0,42 0,26 Ghi chú: (-) khơng bắt gặp Năng suất khai thác trung bình cá phèn khoai vùng đánh cá chung vịnh Bắc Bộ có biến động tương đối rõ theo theo mùa gió thể qua tháng điều tra năm Năng suất khai thác trung bình cao chuyến điều tra tháng 10 (0,42 kg/h) tiếp đến chuyến điều tra tháng (0,26 kg/h), chuyến điều tra tháng (0,28 kg/h) thấp chuyến điều tra vào tháng (0,09 kg/h) Các chuyến điều tra vào tháng tháng 10 đại diện cho mùa gió Đơng Bắc chuyến điều tra vào tháng tháng đại diện cho mùa gió Tây Nam Ta thấy, mùa gió xu hướng biến động loài cá vùng đánh cá chung vịnh Bắc Bộ thời gian nghiên cứu có suy giảm đặc biệt năm 2015 (hình 3) Sự suy giảm biến động suất đánh bắt phèn khoai kiểm chứng qua việc phân tích phương sai (ANOVA) với mức ý nghĩa 95% độ tin cậy Tuy nhiên, khơng thấy có sai khác có ý nghĩa thống kê biến động suất khai thác năm có sai khác có ý nghĩa suất đánh tháng điều tra năm Hiện chưa có nghiên cứu đánh giá chuyên sâu biến động nguồn lợi loài cá vùng biển vịnh Bắc Bộ Một số kết nghiên cứu đưa thống kê sơ báo cáo chuyên đề chuyến điều tra khảo sát khác Do đó, số liệu để đánh giá biến động suất đánh bắt loài cá vùng biển nghiên cứu hạn chế Theo kết nghiên cứu Nguyến Bá Thông Mai Công Nhuận năm 2005 lồi cá vùng biển Đơng Nam Bộ có suất khai thác trung bình đạt khoảng 3,50 kg/h Tuy nhiên, vùng biển nghiên cứu ngư cụ khai thác khác số liệu mang tính chất tham khảo [7] Hình Biến động suất đánh bắt trung bình cá phèn khoai vùng đánh cá chung vịnh Bắc Bộ, giai đoạn 2013 - 2015 189 Mai Cơng Nhuận, Nguyễn Khắc Bát Hình Biến động suất đánh bắt cá phèn khoai vùng đánh chung vịnh Bắc Bộ theo mùa gió, giai đoạn 2013 - 2015 Mật độ phân bố - CPUA (kg/km2) Mật độ phân bố nguồn lợi cá phèn khoai vùng đánh cá chung vịnh Bắc Bộ từ năm 2013 đến năm 2015 thể xu hướng biến động tương ứng với biến động suất khai trung bình Biến động mật độ phân bố nguồn lợi loài số quan trọng việc đánh giá xu biến động sinh khối lồi theo thời gian Mật độ phân bố trung bình cá phèn khoai năm dao động từ 4,22 kg/km2 đến 12,49 kg/km2 Mật độ phân bố nguồn lợi theo tháng năm có biến động tương đối lớn dao động từ 1,38 kg/km2 (tháng 7) đến 18,73 kg/km2 (tháng 10) Sự biến động kiểm chứng qua việc phân tích phương sai (ANOVA) với mức ý nghĩa 0,05) có sai khác có ý nghĩa phân bố mật độ nguồn lợi tháng điều tra đại diện cho mùa năm (bảng 4) Mật độ phân bố lồi có liên quan chặt chẽ với yếu tố môi trường như: Nhiệt độ, độ muối, chất đáy, sinh vật phù du, yếu tố hải dương Mỗi lồi có đặc tính riêng phù hợp với điều kiện mơi trường thích hợp liên quan đến di cư phân bố loài Mối tương quan yếu tố môi trường, trường hải dương với suất phân bố cá nhà khoa học nghiên cứu từ sớm Tuy nhiên, kết mang lại nhiều hạn chế nhiều yếu tố khác đặc biệt Việt Nam Gần đây, trình độ khoa học với hỗ trợ trang thiết bị nghiên cứu 190 đại vấn đề quan tâm nghiên cứu mang lại kết tốt Theo thống kê PGS Phạm Thược (2003) công trình nghiên cứu nguồn lợi hải sản vịnh Bắc Bộ Từ 1960 - 1962: Hợp tác Việt Nam Liên Xô (1960 - 1961) Hợp tác Việt - Trung (1959 - 1962) tiến hành điều tra nguồn lợi cá vịnh Bắc Bộ vùng biển lân cận Lần thông tin quan trọng trữ lượng nguồn lợi hải sản đặc trưng điều kiện môi trường vịnh Bắc Bộ ước tính, xác lập sở khoa học tin cậy áp dụng điều kiện tối quan trọng đầu vào việc xây dựng phát triển nghề khai thác, chế biến hải sản quy mô công nghiệp tầm quốc gia Từ 1962-1988: Trạm nghiên cứu Cá biển (nay Viện Nghiên cứu Hải sản) triển khai số chương trình, đề tài điều tra tổng hợp nghiên cứu nguồn lợi cá môi trường vùng biển Việt Nam như: Điều tra nguồn lợi cá đáy vùng gần bờ phía tây vịnh Bắc Bộ (1962-1965); Điều tra nguồn lợi tôm ven bờ phía tây vịnh Bắc Bộ (1975-1976); phối hợp với Viện Nghề cá Na Uy triển khai đề tài: Điều tra nguồn lợi cá tầng tầng (19781980); Trong suốt giai đoạn này, có cố gắng định nghiên cứu tìm hiểu mối tương quan cá - môi trường, song không mang lại kết rõ rệt Lê Hồng Cầu (2008) báo cáo tổng kết đề tài “Nghiên cứu ảnh hưởng số yếu tố môi trường đến suất khai thác số lồi cá có giá trị kinh tế vịnh Bắc Bộ” nêu số hạn chế số lý chủ yếu sau đây: i) Hạn chế phương pháp luận khoa học phương tiện kỹ thuật để triển khai nghiên cứu; Biến động nguồn lợi số đặc điểm… ii) Nguồn số liệu gốc cá môi trường quản lý độc lập, phần số liệu mơi trường bị thất lạc; iii) Thiếu tính đồng bộ, đồng thời tính liên tục chuỗi liệu [9] Bảng Mật độ phân bố (kg/km2) cá phèn khoai vùng đánh cá chung vịnh Bắc Bộ qua tháng điều tra, giai đoạn 2013 - 2015 Năm Tháng Tháng Tháng Tháng 10 Trung bình 2013 2014 2015 Trung bình 0,00 7,15 0,08 2,41 13,68 16,60 3,49 11,26 0,50 2,90 0,73 1,38 20,30 23,32 12,57 18,73 8,62 12,49 4,22 8,45 Tháng 10-2013 Tháng 4-2013 Tháng 7-2013 Tháng 1-2014 Tháng 4-2014 Hình Phân bố CPUE (kg/h) tháng cá phèn khoai vùng đánh cá chung vịnh Bắc Bộ qua chuyến điều tra, giai đoạn 2013 - 2015 191 Mai Công Nhuận, Nguyễn Khắc Bát Tháng 7-2014 Tháng 10-2014 Tháng 1-2015 Tháng 4-2015 Tháng 7-2015 Tháng 10-2015 Hình Phân bố CPUE (kg/h) tháng cá phèn khoai vùng đánh cá chung vịnh Bắc Bộ qua chuyến điều tra, giai đoạn 2013 - 2015 (tiếp) Trong năm gần đây, số đề tài nghiên cứu dự báo ngư trường khai thác hải sản Viện nghiên cứu Hải sản thực phối hợp với số quan khác trường Đại học Khoa học tự nhiên, ĐHQG Hà Nội có nghiên cứu chuyên sâu mối quan hệ yếu tố môi trường, hải dương suất khai thác phân bố cá như: Đề tài khoa học trọng điểm cấp nhà nước KC.09/14; KC.09/18-15 nhằm phục 192 vụ nâng cao chất lượng dự báo ngư trường khai thác số yếu tố môi trường quan trọng như: Nhiệt độ, độ muối yếu tố quan trọng tương quan chặt chẽ đến phân bố cá Nghiên cứu Nguyễn Văn Hướng Đoàn Văn Bộ (2013) ảnh hưởng nhiệt độ đến phân bố suất khai thác vùng đánh cá chung vịnh Bắc Bộ có lồi cá phèn khoai rằng: Nhiệt độ nước biển tham số môi Biến động nguồn lợi số đặc điểm… trường có vai trị quan trọng liên quan đến tập tính cá Kết nghiên cứu tháng tháng 10 năm 2011 vùng đánh cá chung vịnh Bắc Bộ cho thấy nhiệt độ nước biển có ảnh hưởng rõ suất khai thác tập trung cá Trong tháng 4-2011, nhiệt độ nước biển tầng mặt phía bắc vùng đánh cá chung thấp, tương ứng suất khai thác thấp, khu vực có suất khai thác cao phân bố chủ yếu nửa phía nam, nơi có nhiệt độ nước biển tầng mặt 23,5oC Trong tháng 10-2011, nhiệt độ nước biển tương đối cao biến đổi (26,527,5oC) kéo theo mật độ cá suất khai thác tương đối cao phân bố đồng phạm vi toàn vùng đánh cá chung [8] Trữ lượng Trữ lượng nguồn lợi tức thời cá phèn khoai vùng đánh cá chung vịnh Bắc Bộ giai đoạn từ năm 2013 đến 2015 có biến động suy giảm rõ ràng thời gian nghiên cứu Trữ lượng nguồn lợi từ khoảng 3.360 khả khai thác 1.680 năm 2014 giảm xuống khoảng 50% 1.134 khả khai thác khoảng 567 năm 2015 Sự suy giảm nguồn lợi khoảng thời gian ngắn điều cho thấy nguồn lợi loài cá phèn khoai vùng đánh cá chung vịnh Bắc Bộ chịu tác động lớn gia tăng áp lực khai thác lên quần đàn loài vùng đánh cá chung Theo tháng điều tra, trữ lượng nguồn lợi loài cá cao tháng 10 (1.159 tấn) tiếp đến tháng (756,98 tấn), tháng (162,25 tấn) thấp tháng (92,62 tấn) (bảng 5) Kết nghiên cứu dự án Việt - Trung năm gần trữ lượng nguồn lợi hải sản vùng đánh cá chung nói chung trữ lượng nguồn lợi lồi cá phèn khoai nói riêng có suy giảm nghiêm trọng áp lực khai thác ngày tăng, công nghệ khai khai thác ngày đại nguồn lợi hải sản vô tận Các quan quản lý cần phải có giải pháp hợp lý để bảo vệ phát triển nguồn lợi loài cá vùng đánh cá chung theo hướng bền vững thời gian tới Bảng Trữ lượng ước tính tức thời cá phèn khoai vùng đánh cá chung vịnh Bắc Bộ qua chuyến điều tra, giai đoạn 2013 - 2015 Năm Tháng Tháng Tháng Tháng 10 Trữ lượng (tấn) KNKT (tấn) 2013 2014 2015 Trung bình 481,06 5,69 162,25 919,99 1.116,17 234,78 756,98 33,76 194,86 49,23 92,62 1.364,95 1.568,16 845,13 1.159,41 2.318,70 3.360,25 1.134,83 2.271,26 1.159,35 1.680,13 567,42 1.135,63 Ghi chú: (-) không bắt gặp Phân bố tần xuất chiều dài Biến động chiều dài trung bình cá phèn khoai vùng đánh cá chung vịnh Bắc Bộ giai đoạn từ năm 2013 đến 2015 qua tháng điều tra bắt gặp nhóm chiều dài từ - 18 cm (hình 5) Nhóm chiều dài chiếm ưu từ 10 14 cm (chiếm 82%) Trong tháng 1, bắt gặp nhóm chiều dài từ - 14 cm, tập chung chủ yếu nhóm chiều dài từ - 12 cm chiếm khoảng 78% Trong tháng 4, bắt gặp nhóm chiều dài từ - 17 cm, nhóm chiều dài chiếm ưu từ 11 - 14 cm (chiếm 85%) Trong tháng 7, bắt gặp nhóm chiều dài từ - 16 cm nhóm chiếm ưu bắt gặp từ 11 - 13 cm chiếm 87% Trong tháng 10, nhóm chiều dài bắt gặp từ - 18 cm nhóm chiều dài chiếm ưu bắt gặp từ - 13 cm (chiếm 91%) Theo kết nghiên cứu Nguyễn Bá Thơng Mai Cơng Nhuận (2005) nhóm chiều dài bắt gặp loài cá vùng biển Đông Nam Bộ từ - 12 cm (FL) Ở vùng Biển Đơng lồi cá có chiều dài kinh tế 15,50 cm (SL) chiều dài toàn thân bắt gặp lớn cá đực 20 cm (TL) [2, 7] Như vậy, phân bố chiều dài loài cá vùng đánh cá chung vịnh Bắc Bộ tương đối đa dạng so với vùng biển lân cận khu vực Qua kết phân tích số liệu thu thập từ chuyến điều tra Cá phèn khoai vùng đánh cá chung vịnh Bắc Bộ bị khai thác 193 Mai Công Nhuận, Nguyễn Khắc Bát nhiều thuộc nhóm chiều dài 11,50 cm Nhóm chiều dài có kích thước nhỏ (< 11,50 cm) chiếm tỷ lệ tương đối cao tổng sản lượng khai thác đặc biệt tháng 10 chiếm khoảng 34% Nhóm chiều dài có kích thước > 11,50 cm chiếm khoảng 16% Kích thước nhóm chiều dài bị khai thác nhiều thành phần sản lượng quan trọng phát triển loài tương lai Nhóm có kích thước nhỏ chưa đạt kích thước trưởng thành, kích thước sinh sản bị khai thác nhiều thành phần sản lượng ảnh hưởng đến lượng bổ sung vào quần đàn loài bị suy giảm ngược lại Hình Phân bố xuất chiều dài cá phèn khoai vùng đánh cá chung vịnh Bắc Bộ giai đoạn từ 2013 - 2015 Tương quan chiều dài khối lượng vùng đánh cá chung vịnh Bắc Bộ riêng cho giới tính Kết phân tích cho thấy hệ số đồng hóa b cho hai giống đực giống nhỏ 3,0 hế số di hóa a nhỏ chứng tỏ cá phèn khoai vùng đánh cá chung vịnh Bắc Bộ loài đồng sinh trưởng, cá có phát triển chiều dài tăng trưởng trọng lượng tương đối đồng Phương trình tương quan chiều dài khối lượng giống có dạng: W= 0,013 × 104*L2,10, R2 = 0,94 giống đực có dạng: W = 0,026×104*L2,07, R2 = 0,91 Chiều dài thành thục Lm50 Hình Phương trình tương quan chiều dài khối lượng cá phèn khoai vùng đánh cá chung vịnh Bắc Bộ, giai đoạn 2013 - 2015 Phân tích hồi quy tương quan chiều dài khối lượng xác định số a hệ số mũ b phương trình tương quan chiều dài khối lượng W=a*L b cá phèn khoai 194 Chiều dài Lm50 cá phèn khoai vùng đánh cá chung vịnh Bắc Bộ thời gian nghiên cứu phân tích riêng theo giới tính Số mẫu phân tích lớn, hệ số tương quan R cao cho thấy quan hệ hai biến chiều dài tỷ lệ thành thục chặt Kết phân tích cho thấy: Chiều dài Lm50 phèn khoai giống đực 11,60 cm thấp so với giống 12,60 cm (hình 7) Tuy nhiên, chiều dài thánh thục Lm50 giống giống Biến động nguồn lợi số đặc điểm… cao so với nhóm chiều dài bị khai thác nhiều (11,50 cm) phân tích Như vậy, nhóm chiều dài chưa đạt đến kích thước thành thục bị khai thác nhiều thành phần sản lượng, điều dẫn đến nguồn lợi đàn cá bố mẹ bị thiếu hụt suy giảm ảnh hưởng đến lượng bổ sung hệ bổ sung vào quần đàn lồi Kích thước mắt lưới khai thác áp lực khai thác ảnh hưởng sấu tới phát triển nguồn lợi loài cá vùng đánh cá chung vịnh Bắc Bộ thời gian tới khơng có biện pháp bảo vệ khai thác hợp lý nguồn lợi lồi cá Hình Chiều dài thành thục Lm50 cá phèn khoai vùng đánh cá chung vịnh Bắc Bộ, giai đoạn 2013 - 2015 Trữ lượng quần đàn cá bố mẹ Hình Biến động trữ lượng quần đàn cá bố mẹ cá phèn khoai vùng biển nghiên cứu Trên sở số liệu sinh học, độ chín muồi tuyến sinh dục, tần suất chiều dài, tương quan chiều dài - khối lượng, sinh khối quần đàn cá bố mẹ (SSB) loài cá phèn khoai vùng đánh cá chung vịnh Bắc Bộ dao động từ 160 đến 1.220 khoảng thời gian nghiên cứu Trữ lượng ước tính 980 chiếm khoảng 42% tổng trữ lượng loài (năm 2013); năm 2014 1.220 chiếm 36% tổng trữ lượng; Trữ lượng nguồn lợi lồi suy giảm mạnh năm 2015 xuống cịn 160 chiếm khoảng 14% tổng trữ lượng Song song với suy giảm nghiêm trọng trữ lượng nguồn lợi, trữ lượng quần đàn cá bố mẹ bị giảm 76% khoảng thời gian ngắn (2013 - 2015) Sự suy giảm nghiêm trọng ảnh hưởng xấu đến lượng bổ sung, khơi phục nguồn lợi lồi Đánh giá chung: Tổng trữ lượng loài, trữ lượng quần đàn cá bố mẹ (SSB) bị suy giảm nghiêm trọng thơi gian ngắn Đàn cá bố mẹ chiếm tỷ lệ thấp (%) quần đàn, nhóm có thước nhỏ chưa đạt đến kích thước thành thục bị khai thác nhiều thành phần sản lượng (L < Lm50) Tất phân tích đánh giá cho thấy, nguồn lợi cá phèn khoai vùng đánh cá chung vịnh Bắc Bộ bị khai thác mức Các quan quản lý cần có giải pháp hợp lý để bảo vệ khai thác nguồn lợi loài nói riêng nguồn lợi hải sản vùng đánh cá chung nói chung theo hướng phát triển bền vững 195 Mai Công Nhuận, Nguyễn Khắc Bát KẾT LUẬN VÀ KIẾN NGHỊ Kết luận Xu hướng biến động suất khai thác (kg/h) mật độ phân bố (tấn/km2) cá phèn khoai vùng đánh cá chung vịnh Bắc Bộ có biến động mạnh qua năm Năng suất khai thác trung bình lồi có có suy giảm hai mùa gió Đơng Bắc mùa gió Tây Nam Trữ lượng nguồn lợi cá phèn khoai vùng đánh cá chung vịnh Bắc Bộ bị suy giảm nghiêm trọng từ 3.360 năm 2014 xuống 1.134 năm 2015 Trữ lượng ước tính trung bình khoảng 2.271 khả khai thác cho phép khoảng 1.135 Trữ lượng quần đàn cá bố mẹ bị suy giảm nghiêm trọng chiếm tỷ lệ ngày thấp tổng trữ lượng quần đàn Cá phèn khoai loài đồng sinh trưởng, hệ số đồng hóa (b) xác định cho cá đực cá 2,10 2,07 Hệ số dị hóa (a) 0,013×104 với cá đực 0,026×104 với cá Chiều dài thành thục Lm50 cá phèn khoai vùng đánh cá chung vịnh Bắc Bộ ước tính cho cá 12,6 cm cm cá đực 11,6 cm Chiều dài thành thục Lm50 cá phèn khoai cá cao so với cá đực Kiến nghị Mở rộng phạm vị thu thập số không gian thời gian loài cá vùng biển vịnh Bắc Bộ Cá phèn khoai đối tượng quan trọng mặt kinh tế ý nghĩa hệ sinh thái vùng đánh cá chung vịnh Bắc Bộ Vì vậy, cần tiếp tục điều tra, thu phân tích bổ 196 xung thêm nguồn số liệu sinh học loài cá giúp cho việc đánh giá nguồn lợi lồi cá nói riêng nguồn lợi hải sản vùng đánh cá chung vịnh Bắc Bộ nói chung xác hiệu TÀI LIỆU THAM KHẢO Froese, R., 2004 FishBase World Wide Web electronic publication http://www fishbase org Fishbase.org, http://www.fishbase.org/summary/Upeneus japonicus Nikolsky, G V., 1963 The ecology of fishes London Pennington, M., 1983 Efficient estimators of abundance, for fish and plankton surveys Biometrics, 39(1), 281-286 Gulland, 1971 Estimation of maximum sustainable yield using surplus production models FAO Fisheries Department: pp 314-315 Jennings, S., Kaiser, M J., and Reynolds, J D., 2001 Marine fisheries ecology Nguyễn Bá Thông Mai Công Nhuận, 2006 Biến động nguồn lợi số đặc điểm sinh học chủ yếu cá phèn khoai (Upeneus japonicus Houttuyn, 1782) biển Đông Nam qua chuyến điều tra tàu giã đơn giai đoạn 2000-2005 Tuyển tập cơng trình nghiên cứu nghề cá biển Tập Nxb Nông nghiệp Nguyễn Văn Hường Đoàn Văn Bộ, 2013 Ảnh hưởng nhiệt độ nước biển đến suất khai thác cá vùng đánh cá chung vịnh Bắc Bộ Tạp chí Khoa học Đại học Quốc gia Hà Nội, Khoa học Tự nhiên Công nghệ Pp Biến động nguồn lợi số đặc điểm… VARIATION IN STOCK OF PACIFIC RUDDERFISH (Upeneus japonicus, Houttuyn, 1782) AT THE VIETNAM-CHINA SHARED ZONE FOR FISHING IN THE GULF OF TONKIN FROM 2013 TO 2015 Mai Cong Nhuan, Nguyen Khac Bat Research Institute for Marine Fisheries, Ministry of Agriculture and Rural Development ABSTRACT: Bensasi goatfish is one of the main target species of otter trawling and bottom gill net fishing at the Vietnam- China shared zone for fishing in the Gulf of Tonkin The catabolism and anabolism values (a, b) of the length-weight relationship equation were estimated to be from 0.013×104 - 0.026×104 and 2.07 - 2.10 respectively The length at first maturity of the female was 126 mm and that of the male was 116 mm The average of fishing productivity of this species varied from 0.12 kg/h to 0.40 kg/h during the investigation (2013 to 2015) The density, presented in biomass per area unit, ranged from 4.22 kg/km2 to 12.40 kg/km2 The total stock size was estimated to vary from 1,134 tonnes to 3,360 tonnes with an average of 2,271 tonnes The estimated allowable catch was 1,135 tonnes Keywords: Biomass, bottom gill net, density, fishing productivity, trawl 197 ... bắt trung bình cá phèn khoai vùng đánh cá chung vịnh Bắc Bộ, giai đoạn 2013 - 2015 189 Mai Cơng Nhuận, Nguyễn Khắc Bát Hình Biến động suất đánh bắt cá phèn khoai vùng đánh chung vịnh Bắc Bộ theo... loài cá vùng đánh cá chung vịnh Bắc Bộ thời gian tới biện pháp bảo vệ khai thác hợp lý nguồn lợi lồi cá Hình Chiều dài thành thục Lm50 cá phèn khoai vùng đánh cá chung vịnh Bắc Bộ, giai đoạn. .. trình nghiên cứu nguồn lợi hải sản vịnh Bắc Bộ Từ 1960 - 1962: Hợp tác Việt Nam Liên Xô (1960 - 1961) Hợp tác Việt - Trung (1959 - 1962) tiến hành điều tra nguồn lợi cá vịnh Bắc Bộ vùng biển lân