MỘT số đặc điểm SINH học cá CHỈ VÀNG (SELAROIDES LEPTOLEPIS, CUVIER 1833) ở VÙNG BIỂN tây NAM bộ

MỘT SỐ ĐẶC ĐIỂM SINH HỌC CÁ CHỈ VÀNG (SELAROIDES LEPTOLEPIS, CUVIER 1833) Ở VÙNG BIỂN TÂY NAM BỘ Vũ Thị HậufootnoteRef 1 và Phạm Quốc Huy, Nguyễn Viết Nghĩa 1 Phòng Nghiên cứu Nguồn lợi Hải sản, Vi. SELAROIDES LEPTOLEPIS, CUVIER 1833

MỘT SỐ ĐẶC ĐIỂM SINH HỌC CÁ CHỈ VÀNG (SELAROIDES LEPTOLEPIS, CUVIER 1833) Ở VÙNG BIỂN TÂY NAM BỘ Vũ Thị Hậu1 Phạm Quốc Huy, Nguyễn Viết Nghĩa Abstract SOME BIOLOGY CHARACTERS OF YELLOWSTRIP SCAD (SELAROIDES LEPTOLEPIS, CUVIER 1833) IN THE SOUTHWEST SEA OF VIETNAM Vu Thi Hau & Phạm Quốc Huy, Nguyen Viet Nghia Summary: Biology characters of Yellowstrip Scad (Selaroides leptolepis, Cuvier 1833) were analyzed by the fisheries biology data collected in landing sites from 2014-2015 in the Southwest Sea of Vietnam Results indicated that the main lengh of Yellowstrip scad in the Southwest sea estimated at 11,3 cm, between 4,5cm-15,4 cm The lengthweight relationship equations are: W= 0,000008L3,154 (Male), W= 0,000009L3,114 ( Female), W=0,00002L2,965 (Juvenie) Growth parameters (L∞ and K) of the von Bertalanffy equation were 16,3 cm and 1,2/year, respectively The length at first maturity of fish (Lm50) was 9,8cm Sex rate of fish was 1,1 The Yellowstrip scad breeded throughout the year and the spawning peaks are observed in Febuary and April The instantaneous total mortality Z=4,24, the instantaneous natural mortality due to predation, aging and other environmental causes M= 2,39, the instantaneous fishing mortality caused by fishing F=1,85 and exploitation rate E= 0,44 Keywords: Yellowstrip scad, Selaroides leptolepis, biology characters Tóm tắt Đặc điểm sinh học cá vàng (Selaroides leptolepis Cuvier, 1833) vùng biển Tây Nam Bộ phân tích dựa số liệu từ 14 chuyến thu mẫu sinh học Tiểu dự án I.9 “ Điều tra tổng thể trạng biến động nguồn lợi hải sản biển Việt Nam” Kết phân tích cho thấy, chiều dài trung bình cá vàng 11,3 cm, dao động từ 4,5-15,4 cm Sự tương quan chiều dài- khối lượng cá mơ tả theo phương trình W= 0,000008L 3,154 (cá đực), W= 0,000009L3,114 (cá cái), W=0,00002L2,965 (cá con).Chiều dài tối đa lý thuyết cá vàng tương đối thấp (16,3 cm), hệ số sinh trưởng K =1,2/ năm Cá tham gia sinh sản lần đầu chiều dài 9,8 cm Tỷ lệ đực/cái quần thể cá vàng 1,1 Cá đẻ rải rác quanh năm mùa vụ sinh sản từ tháng đến tháng Hệ số chết chung quần thể cá vàng xác định Z= 4,24, mức chết tự nhiên M=2,39, hệ số chết khai thác F= 1,85 Như vậy, hệ số khai thác E 0,44 Từ khóa: cá vàng, Selaroides leptolepis, sinh trưởng, sinh sản Phòng Nghiên cứu Nguồn lợi Hải sản, Viện Nghiên cứu Hải sản 224 Lê Lai, Quận Ngơ Quyền, Hải Phịng ĐT: 0313826986 Email: hauvuhp@gmail.com MỞ ĐẦU Cá vàng (Selaroides leptolepis) lồi cá nổi, có giá trị kinh tế cao, kích thước cá không lớn số lượng quần thể đông, có khả cho sản lượng khai thác lớn khai thác quanh năm (Võ Văn Phú, 2010) Theo FAO, cá vàng phân bố vùng biển Tây Ấn Độ Dương, vùng biển nhiệt đới từ phía Tây Thái Bình Dương tới Phillipine, giới hạn phía Bắc Nhật Bản phía Nam Australia Ở Việt Nam, cá phân bố rộng khắp từ vịnh Bắc đến Nam Bộ (Lê Hải Thiện, 2010) Các nghiên cứu cá vàng giới nghiên cứu vùng biển Malaysia (Mansor Mat Isa et al., 1994), Philipin (May R Guanco & Sheryll V Mesa, 2004), Thái Lan (Noopeth, 1984) Ở Việt Nam có nhiều nghiên cứu cá vàng, nhiên, nghiên cứu thường tập trung vào đánh giá trữ lương, nguồn lợi phân bố (Vũ Việt Hà, 2011), (Đặng Văn Thi & ctv, 2005) Các kết nghiên cứu đặc điểm sinh học sinh trưởng sinh sản cá vàng rời rạc, chưa cập nhật tập trung vài vùng biển nhỏ (Võ Văn Phú, Nguyễn Thị Hoàn, 2010), (Lê Hải Thiện, 2010) Từ nguồn số liệu thu thập chuyến thu mẫu sinh học nghề cá giai đoạn 2014-2015 đến nay, viết phân tích cung cấp thông tin số đặc điểm sinh học quần thể loài cá vàng vùng biển Tây Nam Bộ, góp phần bổ sung sở khoa học cho việc bảo vệ nguồn lợi phát triển nghề cá TÀI LIỆU VÀ PHƯƠNG PHÁP NGHIÊN CƯU 2.1 Tài liệu Tài liệu sử dụng viết thu thập từ tháng 5/2014-6/2015 Tiểu dự án I.9 “Điều tra tổng thể trạng biến động nguồn lợi hải sản biển Việt Nam” thuộc dự án “Điều tra tổng thể đa dạng sinh học nguồn lợi thủy, hải sản vùng biển Việt Nam; Quy hoạch xây dựng hệ thống khu bảo tồn biển phục vụ phát triển bền vững”.Tổng số mẫu thu 14 chuyến thu mẫu 2.314 mẫu 2.2 Phương pháp thu mẫu Mẫu sinh học thu ngẫu nhiên nhóm thương phẩm đội tàu khai thác vùng biển Tây Nam Bên cạnh thông tin sản lượng nhóm thương phẩm, thơng tin tàu thuyền ngư lưới cụ thu thập Với tần suất thu mẫu hàng tháng bến cá, lồi phân tích khoảng 100 -200 cá thể Các tiêu phân tích bao gồm: chiều dài đến chẽ vây theo ngun tắc làm trịn xuống (mm), cân khối lượng thể khối lượng tuyến sinh dục (gram), xác định độ chín muồi tuyến sinh dục theo tháng bậc Nikolski, 1963 Mẫu vật phân tích trường Phương pháp phân loại đặc điểm hình thái chủ yếu dựa vào tài liệu Vương Dĩ Khang; FAO; Nguyễn Văn Chung; Nguyễn Nhật Thi; Nguyễn Hữu Phụng; Phạm Thược; Tetsuji Nakabo; Fishbase… 2.3 Phương pháp phân tích số liệu Phân bố tần suất chiều dài chiều dài trung bình: Phân bố tần suất chiều dài chiều dài trung bình đối tượng nghiên cứu sản lượng khai thác phân tích phương pháp thống kê mơ tả thông thường theo hướng dẫn Fowler nnk, 1998: Trong đó: FL chiều dài đến chẽ vây trung bình cá; FLj chiều dài cá nhóm thứ j; fj số cá thể nhóm thứ j; n tổng số cá thể, m số nhóm chiều dài Tham số sinh trưởng chủng quần: Tham số sinh trưởng phương trình sinh trưởng von Bertallanfy (k, L∞) ước tính cho loài dựa số liệu tần suất chiều dài Lt = L∞ ∗ (1 − e − k ( t −t0 ) ) (Gayanilo, et al., 2002) Trong đó: Lt chiều dài cá; L∞ chiều dài cực đại cá; k hệ số sinh trưởng; t tuổi cá; t o chiều dài lý thuyết giả định tuổi Chiều dài Lm50: Chiều dài Lm50 chiều dài có 50% số cá thể tham gia vào sinh sản lần đầu, tính theo King, 1995;Sparre & Venema, 1998: P= 1 + e [ − r ( L − Lm 50 )] Trong đó: Lm50 chiều dài cá mà có 50% số cá thể tham gia vào sinh sản lần đầu; L chiều dài cá; r số Hệ số chết (hệ số chết chung Z, hệ số chế khai thác F, hệ số chế tự nhiên M): Hệ số chết chung (Z) tính tốn theo phương pháp Length-based phần mềm FiSAT II (Gayanilo, et al., 2002): ln(Ni/ti) = a + b * ti Trong đó: N số lượng cá thể nhóm thứ I, t i thời gian cần thiết để cá thể phát triển đến nhóm chiều dài thứ i, t i tuổi cá thể nhóm chiều dài i, a b tham số Hệ số chết tự nhiên (M) xác định theo công thức thực nghiệm Pauly (1980): log (M) = - 0,0066 - 0,279 log(L∞) + 0,6543 log (K) + 0,4634 log (T) Trong đó: L∞ chiều dài cực đại cá thể; K hệ số sinh trưởng; T nhiệt độ trung bình nơi sinh cư Hệ số chết khai thác (F) tính tốn theo công thức: F = Z – M Xác định hệ số thành thục: Để xác định mùa sinh sản đối tượng theo thời gian, cần xác định: tỉ lệ phát triển tuyến sinh dục theo thời gian, biến động hệ số GSI theo thời gian, tỉ lệ đực cái, thời điểm xuất cá (Juv.) cá phát triển tuyến sinh dục giai đoạn sớm Tỉ lệ đực/cái tỉ lệ giai đoạn phát triển tuyến sinh dục đối tượng phân tích theo tháng, tính tốn phương pháp thống kê mô tả thông thường theo hướng dẫn Fowler nnk, 1998 Mùa sinh sản loài cá dựa biến động hệ số thành thục sinh dục GSI (Gonado Somatic Index) biến động tỉ lệ giai đoạn phát triển tuyến sinh dục theo tháng GSI trung bình lồi cá xác định theo công thức West (1990): với GSI = wi*100/Wi Trong đó: wi khối lượng tuyến sinh dục cá thể i; Wi khối lượng cá 3 KẾT QUẢ VÀ THẢO LUẬN 3.1 Đặc điểm sinh học sinh trưởng 3.1.1 Phân bố tần suất chiều dài Cá vàng khai thác vùng biển Tây Nam có chiều dài đến chẽ vây dao động khoảng 4,5-15,4cm, chiều dài trung bình 11,3 cm Cá thể nhỏ 4,5 cm cá thể có chiều dài lớn 15,4 cm bắt gặp vào tháng năm 2014 Cấu trúc nhóm chiều dài bắt gặp, nhóm chiều dài ưu cá vàng khác theo tháng năm Đàn cá có kích thước nhỏ thường bắt gặp tháng 4-6 hàng năm với kích thước trung bình khoảng 10cm Đàn cá có kích thước lớn, trung bình khoảng 14cm bắt gặp nhiểu tháng 10,11,12 Các tháng 3,4,5,6 tháng 10 bắt gặp nhóm chiều dài Trong đó, nhóm thứ có chiều dài chiếm ưu cm (đây hệ sinh sau mùa sinh sản) nhóm thứ hai có chiều dài chiếm ưu 14 cm Hình Phân bố tần suất chiều dài cá vàng vùng biển Tây Nam 3.1.2 Tương quan chiều dài - khối lượng Phương trình tương quan chiều dài khối lượng quần thể cá vàng vùng biển Tây Nam xác định sở 2.314 cá thể đo chiều dài (L = 4,5÷15,4cm) cân khối lượng (W = 1÷66,7g) Kết tương quan chiều dài khối lượng cá vàng trình bày Hình cho thấy có biến động nhẹ đực, non Trong đó, tham số sinh trưởng b >3 cá đực cá cho thấy phát triển ưu trọng lượng lớn chiều dài ngược lại với non b < Kết tương đồng với kết nghiên cứu K K Tandon (1957) Ấn Độ với hệ số b cá đực (3,87) cá (3,52) lớn Selaroides le ptole pis Male W = 0.000008*L^3.154 Female W = 0.000009*L^3.114 Juv W = 0.00002*L^2.965 80 70 60 W (g) 50 40 30 20 10 40 60 80 100 120 140 160 M F J FL (mm) Hình Tương quan chiều dài - khối lượng cá vàng Tây Nam 3.1.3 Tham số sinh trưởng phương trình von Bertalanffy Sử dụng phương pháp ELEFAN I phân tích phân bố tần suất chiều dài hàng tháng xác định tham số sinh trưởng phương trình sinh trưởng chiều dài von Bertalanffy cá vàng vùng biển Tây Nam Bộ Đường cong biểu diễn sinh trưởng chiều dài loài cá trình bày Hình cho thấy chiều dài tối đa lý thuyết cá vàng L ∞ 16,3cm, hệ số sinh trưởng K 1,2/năm hệ số sinh trưởng trung bình ø’= 2,5 Kết nhỏ nhiều so với nghiên cứu Võ Văn Phú vùng biển Thừa Thiên Huế với chiều dài tối đa lý thuyết 22,8 cm nghiên cứu Noopeth (1984) vùng biển vịnh Thái Lan (19,2cm) Điều cho thấy áp lực khai thác tác động lên quần thể cá vàng vùng biển Tây Nam Bộ có dấu hiệu nghiêm trọng Hình Phương trình sinh trưởng cá vàng vùng biển Tây Nam 3.2 Đặc điểm sinh học sinh sản 3.2.1 Tỉ lệ đực/cái Từ 2.314 cá thể giải phẫu cho thấy tỷ lệ đực/cái quần thể lồi cá vàng xác định 1,1 (Hình 4) Cấu trúc giới tính biến động khác theo tháng, cá đực chiếm ưu so với cá vào hầu hết tháng năm Tuy nhiên vào tháng 4,5,6, tỷ lệ cá đực cá xấp sỉ có xuất thêm non; mùa sinh sản cá vàng tỷ lệ đực/cái xấp xỉ 50/50, đảm bảo trứng đẻ có tỷ lệ thụ tinh cao, cung cấp lượng bổ sung tốt vào quần đàn Ở thời điểm lại, cá đực chiếm ưu cấu trúc quần thể Hình Tỉ lệ giới tính cá vàng vùng biển Tây Nam 3.2.2 Chiều dài Lm50 Chiều dài thành thục lần đầu cá vàng xác định, sở phân tích từ số liệu sinh học nghề cá thu thập theo tháng bến cá từ tháng năm 2014 đến tháng năm 2015 (Hình 5) Kết nghiên cứu cho thấy, giá trị chiều dài thành thục loài cá Lm 50 = 9,8cm nhỏ so với kết nghiên cứu Novi Adriani (2015) Indonesia 13,3cm Chiều dài thành thục bị ảnh hưởng yếu tố khí hậu, đặc điểm sinh thái dinh dưỡng theo thời gian Ngoài ra, yếu tố áp lực khai thác tác động người ảnh hưởng đến giá trị Do mà nguồn lợi cá vàng bị ảnh hưởng, thường biến động theo chuỗi thời gian, có thời điểm hoạt động khai thác dồi có thời điểm xem Model: PSL=1/(1+exp(-r*(LSL-Lm))) y=1/(1+exp(-(.926463)*(x-(9.7 9904)))) 1.0 10 0.8 PSL 0.6 0.4 0.2 0.0 10 12 14 16 18 LSL Hình Phương trình Lm50 cá vàng vùng biển Tây Nam 3.2.3 Biến động tỉ lệ giai đoạn thành thục sinh dục Tỷ lệ thành thục cá vàng tổng hợp phân tích theo tháng sở liệu sinh học đối tượng thu thập từ tháng năm 2014 đến tháng năm 2015 (Hình 6) Hình Tỉ lệ giai đoạn chín muồi tuyến sinh dục cá vàng Xuất đàn cá bố mẹ có tuyến sinh dục thành thục (giai đoạn IV), tham gia sinh sản (giai đoạn V) sinh sản xong (giai đoạn VI) quanh năm Trong đó, tỷ lệ đàn cá thành thục sinh sản đạt tỷ lệ cao vào tháng 2-4 tương ứng với 90% tổng số cá thể bắt gặp Đàn cá có tuyến sinh dục giai đoạn phát triển (giai đoạn II III) bắt đầu xuất vào tháng 5-7 năm Cấu trúc quần thể có biến động mạnh tỷ lệ thành thục sinh dục biểu rõ Đàn cá non có tuyến sinh dục phát triển (giai đoạn Juv.) bắt gặp vào tháng năm Đồng thời với xu hướng giảm dần đàn cá bố mẹ tham gia sinh sản xong Khoảng thời gian sau đó, tỷ lệ đàn cá non đàn cá bố mẹ ln thay đổi Trên sở nhận định rằng, cá vàng biển Tây Nam đẻ rải rác quanh năm Nhận định tương đồng với nghiên cứu Lê Hải Thiện (2010) cho cá vàng đẻ từ tháng đến tháng 5, kéo dài đến tháng 9, tháng 10 3.2.4 Hệ số thành thục sinh dục Hệ số thành thục sinh dục biến động hệ số thành thục cá vàng xác định theo tháng thu mẫu trình bày Hình Hệ số thành thục tương đối đồng tháng Tuy nhiên, vào tháng 2,3,4, hệ số thành thục đạt giá trị cao Kết hợp với kết phân tích độ chín muồi tuyến sinh dục cho thấy cá vàng sinh sản rải rác quanh năm mùa vụ sinh sản từ tháng đến tháng Hình Hệ số GSI cá vàng theo thời gian vùng biển Tây Nam 3.2.5 Các hệ số tử vong hệ số khai thác Mức độ tử vong số lượng cá thể quần thể bị chết đơn vị thời gian Mức độ tử vong quần thể phụ thuộc trước hết vào tuổi thọ trung bình sinh vật điều kiện sống môi trường, biến đổi bất thường khí hậu, bệnh tật, lượng thức ăn có mơi trường, số lượng kẻ thù, mức độ khai thác người Mức chết quần thể cá vàng xác định dựa đường cong sản lượng tuyến tính hóa từ tần suất chiều dài với tham số đầu vào gồm: L ∞ = 16,3cm, K = 1,2/năm; T = 28oC Hệ số chết chung quần thể cá vàng xác định Z = 4,24 Theo cơng thức thực nghiệm Pauly (1980) ước tính mức chết tự nhiên quần thể loài cá M=2,39 Từ ta xác định hệ số chết khai thác là: F = 1,85 (Hình 8) Hình Ước tính hệ số chết từ phương trình đường cong sản lượng cá vàng Biến động mạnh mức chết chung tác động mức chết khai thác Trên sở mức chết khai thác mức chết tự nhiên xác định hệ số khai thác E, kết tương ứng E = F/Z = 0,44 Như vậy, kết chứng tỏ áp lực hoạt động khai thác lên quần thể cá vàng vùng biển Tây Nam chưa mức cho phép, có dấu hiệu suy giảm, cần có biện pháp để hạn chế ảnh hưởng áp lực khai thác tác động với quần thể loài cá KẾT LUẬN Chiều dài trung bình cá vàng vùng biển Tây Nam Bộ năm 2014-2015 11,3 cm, dao động từ 4,5-15,4 cm Phương trình tương quan chiều dài- khối lượng cá vàng thể hiện: W= 0,000008L3,154(cá đực), W= 0,000009L3,114( cá cái), W=0,00002L2,965 (cá con).Chiều dài tối đa lý thuyết cá vàng tương đối thấp (16,3 cm), hệ số sinh trưởng K =1,2/ năm Cá tham gia sinh sản lần đầu chiều dài 9,8 cm Tỷ lệ đực/cái quần thể cá vàng 1,1 Cá đẻ rải rác quanh năm mùa vụ sinh sản từ tháng đến tháng Hệ số chết chung quần thể cá vàng xác định Z= 4,24,mức chết tự nhiên M=2,39, hệ số chết khai thác F= 1,85 Như vậy, hệ số khai thác E 0,44 TÀI LIỆU THAM KHẢO Vũ Việt Hà (2011), Trữ lượng tức thời phân bố nguồn lợi cá vàng (Selaroides leptolepis) vùng biển phía nam Việt Nam dựa kết điều tra phương pháp thủy âm, Tạp chí Khoa học Công nghệ biển T11, Số Tr 85-96 Vương Dĩ Khang (1973) Ngư loại phân loại học, Nhà xuất Nông thôn Võ Văn Phú, Nguyễn Thị Hoàn (2010), Đặc điểm sinh trưởng cá vàng Selaroides 10 11 leptolepis (Cuvier, 1833) vùng biển Thừa Thiên Huế, Tạp chí Khoa học, Đại học Huế, Số 57, 2010, Trường Đại học Khoa học, Đại học Huế Lê Hải Thiện (2010), Một số đặc điểm sinh học cá Ngân (Alepes kneinii), cá Chỉ vàng (Selaroides leptolepis)và cá Tráo Mắt to (Selar crumenophthalmus) vùng biển Sóc Trăng, Bạc Liêu Luận văn tốt nghiệp, trường đại học Cần Thơ Vũ Việt Hà, Đặng Văn Thi & Nguyễn Duy Thành (2005), Tổng quan nguồn lợi hệ sinh thái vùng biển miền Trung Viện Nghiên cứu Hải sản, Hải Phòng Đặng Văn Thi, Nguyễn Bá Thông & Vũ Việt Hà (2005), Tổng quan nguồn lợi hệ sinh thái vùng biển vịnh Bắc Bộ Viện Nghiên cứu Hải sản, tr68 Đặng Văn Thi, Nguyễn Bá Thông & Vũ Việt Hà (2005), Tổng quan nguồn lợi hệ sinh thái vùng biển Đông Nam Bộ, Viện Nghiên cứu Hải sản Gayanilo, F.C.: P Sparre and D Pauly (1996), FISAT:FAO_ICLARM stock assessment tools, Ủ’s manual Rome, FAO, 126pp K K Tandon (1957), Biology and fishery of “ Choo parai”- Selaroides leptolepis (Cuvier and Valenciennes), Department of Zoology, Panjab University, Chandigarh Gayanilo, F.C.: P Sparre and D Pauly (1996), FISAT:FAO_ICLARM stock assessment tools, UR’s manual Rome, FAO, 126pp May R Guanco & Sheryll V Mesa (2004), Assessment of the commercial fisheries of western and central Visayan Sea, National Stock Assessment Project, Bureau of Fisheies and Aquatic Resources, Regional Field Office 6, Iloilo City 12 Mansor Mat Isa, Abd Haris Hilmi Ahmad Arshad & Samsudin Basir (1994),Distribution, Abundance and Biological Studies of Economically Important Fishes I the South China Sea, Area I: East Coast of Peninsular Malaysia, Marine fishery resources developmet and management department, Southeast Asian Fisheries development center, Chendering, 21080, Kuala Terengganu, Malaysia 13 Nikolsky (1963), The ecology of fishes London, acadimic Press, pp: 353 14 Noopeth, L (1984) Length frequency distribution, age and growth determination of yellow trevally Technical Report No 47 Pelagic Fisheries Unit, Marine Fisheries Division, DOF, Thailand 15 Novi Andriani (2015), Aspek Biologi Dan Tingkat Pemanfaatan Ikan Selar Kuning (Selaroides Leptolepis) Yang Tertangkap Jaring Cantrang Di Perairan Kabupaten Pemalang, Management Of Aquatic Resources Journal, Management Of Aquatic Resources Journal; Volume 4, Nomor 4, Tahun 2015; 24-32 16 Pauly and Martosubroto (1980) A selection of simple method for the assessent of tropical fish stock FAO, Fisheries Circular