Tổng quan về luận án

Luận án tiến sĩ Hải dương học "Nghiên cứu mối quan hệ giữa cấu trúc hải dương và năng suất khai thác một số loài cá kinh tế ở vùng biển Đông Nam Bộ" (mã số 62440228) của Nguyễn Văn Hướng, bảo vệ tại Trường Đại học Khoa học Tự nhiên, ĐHQGHN năm 2018, dưới sự hướng dẫn của TS Đoàn Văn Bộ và TS Nguyễn Khắc Bát, giải quyết một khoảng trống nghiên cứu rõ rệt: cho tới thời điểm thực hiện, Việt Nam "mới chỉ có duy nhất 01 mô hình và quy trình dự báo ngư trường khai thác cá ngừ đại dương ở vùng biển xa bờ Miền Trung do đề tài cấp Nhà nước KC09.18/11-15 nghiên cứu xây dựng", trong khi nhiều đối tượng cá kinh tế khác chưa hề có cơ sở khoa học cho dự báo nghiệp vụ.

Vùng biển Đông Nam Bộ (06°00'–11°30'N, 105°00'–113°00'E) là ngư trường có trữ lượng và khả năng khai thác lớn nhất cả nước, ước tính khoảng 976.000 tấn với khả năng khai thác cho phép 425.300 tấn, trong đó cá nổi nhỏ khoảng 349.000 tấn và cá đáy khoảng 627.000 tấn. Nghiên cứu chọn ba đối tượng đại diện cho ba nhóm sinh thái khác nhau: cá Ngừ Vằn (Katsuwonus pelamis, đại diện nhóm cá ngừ đại dương, chiếm trên 60% trữ lượng cá nổi lớn ~300.000 tấn), cá Ngừ Chấm (Euthynnus affinis, đại diện cá ngừ nhỏ ven bờ, trữ lượng khoảng 23,9 nghìn tấn) và cá Chỉ Vàng (Selaroides leptolepis, đại diện cá nổi nhỏ ven bờ, trữ lượng phía Nam ước 177–351 nghìn tấn).

Mục tiêu trung tâm được phát biểu rõ: "Xác định được mối quan hệ giữa năng suất khai thác một số loài cá kinh tế với các yếu tố hải dương học, môi trường ở vùng biển Đông Nam Bộ phục vụ xây dựng mô hình thực nghiệm dự báo ngư trường khai thác tại vùng biển này". Hai đóng góp mới được tác giả xác lập: (1) định lượng mối quan hệ giữa CPUE của ba loài cá với các yếu tố cấu trúc hải dương đặc trưng; (2) xây dựng bộ chỉ số thích ứng sinh thái (HSI) cho từng loài theo hai mùa gió, làm nền tảng nâng cao chất lượng bản dự báo ngư trường.

Literature Review và Positioning

Luận án hệ thống hóa ba khuynh hướng đánh giá nguồn lợi cá trên thế giới: (1) nguyên lý Russel và các cải tiến thống kê nghề cá (dN/dt = W + R – (F + M) + ε), với các mô hình giải tích LCA, VPA, MULTIFAN-CL; (2) tiếp cận trao đổi năng lượng dinh dưỡng qua chuỗi thức ăn; (3) tương tác tổng hợp cá–môi trường–khai thác. Tác giả chọn hướng thứ ba — hướng mà "hầu hết các công trình nghiên cứu dự báo biến động nguồn lợi cá trên thế giới" hiện theo đuổi.

Ở bình diện quốc tế, luận án đối sánh với nhiều nghiên cứu cụ thể. Silva và nnk (2000) chứng minh sản lượng cá thu, cá cơm, cá trích ở bắc Chile liên hệ mật thiết với SST và chlorophyll-a. Mansor và nnk (2001) xây dựng hệ thống dự báo ngư trường Biển Đông bằng GIS và phần mềm TroFFS cho vùng nhiệt đới. Mukti và nnk (2011) xác định khoảng thích ứng sinh thái của cá Ngừ Vằn với SST và chlorophyll-a tại vùng Bone Bay–Flores (Indonesia). Chen và nnk (2008, 2009) áp dụng hồi quy tuyến tính kết hợp mô hình HSI dự báo cá Ngừ Mắt To ở Ấn Độ Dương và cá thu Nhật Bản ở biển Đông Trung Hoa. Nhóm nghiên cứu Học viện Công nghệ Châu Á, Nhật Bản và Sri Lanka xác định chỉ số thích ứng của cá Ngừ Vây Vàng ở Ấn Độ Dương với sản lượng cao nhất tại SST 28–30°C, Chl0 0,1–0,4 mg/m³ và SSH 205–215 cm. Năm 2015, WCPFC vận hành mô hình SEAPODYM dự báo mật độ (g/m²) cá ngừ đại dương hạn 7 ngày ở Tây Thái Bình Dương.

Trong nước, tồn tại tranh luận phương pháp đáng chú ý: mô hình KC09 của Đoàn Văn Bộ và nnk (2010, 2015) — đang dự báo nghiệp vụ hạn tháng và hạn 7–10 ngày với độ chính xác trên 60% — tính cả 26 biến môi trường chỉ từ trường 3D nhiệt độ, trong khi các công trình quốc tế nêu trên nhấn mạnh vai trò độc lập của độ muối, chlorophyll-a và dòng chảy. Luận án định vị tại chính điểm giao này: mở rộng số biến môi trường (14 yếu tố cho cá ngừ, 10 yếu tố cho cá Chỉ Vàng) và mở rộng đối tượng từ nghề câu cá ngừ đại dương sang lưới rê và lưới kéo đáy.

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Nghiên cứu mở rộng khung lý thuyết sinh thái cá–môi trường theo ba trục. Thứ nhất, thay vì dùng SST như một biến đơn (truyền thống Froese và Pauly, 1996), luận án thao tác hóa "cấu trúc hải dương" thành hệ biến đa tầng: độ dày lớp đồng nhất trên (H0), độ sâu biên dưới lớp đột biến (H1), nhiệt độ biên dưới lớp đột biến (T1), độ sâu các mặt đẳng nhiệt 20°C và 24°C (H20, H24), gradien nhiệt ngang cực đại tầng mặt và tầng 50m (Grad0, Grad50). Thứ hai, tích hợp động lực học (dị thường độ cao động lực Alti, mật độ động năng rối EKE, tốc độ dòng chảy Spd_cur) — nhóm biến phản ánh front và nước trồi mà Lê Đức Tố và nnk (1995) đã chỉ ra vai trò nhưng chưa định lượng được do "sự thiếu hụt và tính không đồng bộ của nguồn dữ liệu thống kê nghề cá". Thứ ba, xác lập nguyên tắc điều kiện biên theo tập tính loài: bộ biến cho cá Ngừ Vằn/Ngừ Chấm phủ tầng 0–150m (nơi chúng "thường tập trung nhiều hơn ở lớp nước 50–80m"), còn cá Chỉ Vàng — với 62,3% mẻ lưới ở dải sâu 30–50m — chỉ dùng 10 biến tầng mặt đến 50m, thay Grad50 bằng Grad25.

Khung phân tích độc đáo

Khung phân tích tích hợp ba dòng lý thuyết: lý thuyết khoảng thích ứng sinh thái (ecological tolerance), mô hình Habitat Suitability Index (HSI) theo truyền thống Chen và nnk, và phân tích hồi quy tuyến tính nhiều biến với tiêu chuẩn thông tin Akaike (AIC) để chọn biến. Điểm mới về phương pháp là quy trình đồng bộ hóa cá–môi trường ba cấp thời gian (trung bình mùa gió, trung bình tháng nhiều năm, và tức thời theo ngày giai đoạn 2012–2015) trên lưới thống nhất 0,5×0,5 độ kinh vĩ — cho phép tách bạch quan hệ khí hậu học khỏi quan hệ synoptic phục vụ dự báo hạn ngắn.

Phương pháp nghiên cứu tiên tiến

Thiết kế nghiên cứu

Nghiên cứu theo triết lý thực chứng (positivism), thiết kế định lượng tương quan trên dữ liệu quan trắc dài hạn. Nguồn dữ liệu hải dương học gồm 4.325 trạm với 76.729 số liệu nhiệt độ, 2.776 trạm với 57.212 số liệu độ muối và 464 trạm với 22.835 số liệu chlorophyll-a từ các chuyến khảo sát của Viện Nghiên cứu Hải Sản giai đoạn 1978–2013, bổ sung 162.569 điểm dữ liệu viễn thám mỗi biến (SST, chlorophyll-a, dị thường mực biển, dòng chảy) từ dự án MOVIMAR 2011–2016 với độ phân giải 1/12 độ. Nguồn dữ liệu này đã được Bùi Thanh Hùng và nnk (2016) kiểm chứng với đo đạc CTD, sai lệch tối đa ±0,5°C với nhiệt độ và ±0,2‰ với độ muối.

Quy trình nghiên cứu rigorous

Dữ liệu nghề cá lấy từ CSDL VietFishBase (1996–2016) và số liệu WCPFC 2010–2016 thu tại cảng cá Quảng Ngãi, Bình Định, Phú Yên, Khánh Hòa: 4.116 mẻ lưới rê cá Ngừ Vằn (CPUE trung bình dao động 10,83–22,01 kg/km lưới theo tháng), 910 mẻ cá Ngừ Chấm và 247 mẻ lưới kéo đáy cá Chỉ Vàng. CPUE được chuẩn hóa: kg/km lưới cho lưới rê (CPUE=C/L), kg/giờ cho lưới kéo (CPUE=C/t). Kiểm soát chất lượng thực hiện hai lớp: loại sai số thô về thời gian/tọa độ và áp dụng luật "3 sigma" loại giá trị dị thường. Đáng chú ý về tính chặt chẽ kiểm định: dữ liệu 1996–2015 dùng xây dựng mô hình, còn toàn bộ số liệu năm 2016 được giữ lại làm tập kiểm chứng độc lập — thiết kế out-of-sample validation ít gặp trong các nghiên cứu nghề cá trong nước trước đó.

Data và phân tích

Các đặc trưng cấu trúc nhiệt được tính theo tiêu chí định lượng rõ ràng (H0 xác định theo ngưỡng Gz < 0,05°C/m; dòng chảy địa chuyển tính từ dị thường mực biển theo u = -(g/f)dz/dx, v = -(g/f)dz/dy). Phân tích gồm ma trận tương quan cá–môi trường theo từng mùa gió với ba loại chuỗi dữ liệu (tức thời, trung bình 2014–2015, trung bình nhiều năm), hồi quy tuyến tính nhiều biến, và mô hình HSI xác định chỉ số phù hợp SI theo tần suất xuất hiện CPUE trong các khoảng dao động của từng yếu tố. Bản đồ phân bố dựng bằng nội suy khoảng cách trong GIS.

Phát hiện đột phá và implications

Những phát hiện then chốt

Thứ nhất, nghiên cứu xác lập bức tranh cấu trúc hải dương hai mùa hoàn chỉnh cho VBĐNB: mùa gió Đông Bắc SST 25,0–29,0°C, dòng chảy ngoài khơi hướng Bắc–Nam đạt 20–30 cm/s; mùa gió Tây Nam SST 28,5–30,2°C, dòng chảy đảo hướng Tây Nam–Đông Bắc 10–15 cm/s, khu vực nước trồi có nhiệt tầng mặt thấp hơn xung quanh 1–2°C; độ muối biến thiên hẹp 29,0–34,0‰ giữa hai mùa.

Thứ hai, quan hệ cá–môi trường có tính mùa rõ rệt: các ma trận tương quan (bảng 4.1–4.10 của luận án) cho thấy tổ hợp yếu tố chi phối CPUE cá Ngừ Vằn ở mùa gió Đông Bắc khác biệt căn bản so với mùa gió Tây Nam — củng cố luận điểm lý thuyết rằng "ngay cả cùng một loài nhưng sống ở các điều kiện môi trường khác nhau (mùa vụ/khu vực) lại có những thích nghi và phản ứng riêng".

Thứ ba, ba bộ chỉ số thích ứng sinh thái theo mùa được xác lập lần đầu cho cả ba loài tại VBĐNB (bảng 4.11–4.12, 4.16–4.17, 4.21–4.22). Kết quả này mở rộng đáng kể so với công trình duy nhất trước đó của Nguyễn Văn Lục và nnk (1991) — vốn chỉ xác định được cho cá Chỉ Vàng các khoảng H0 10–35m, nhiệt độ tầng đáy 25,12–27,7°C, độ muối tầng đáy 32,58–33,45‰ với CPUE 0,2–45,0 kg/giờ bằng thống kê tần suất đơn thuần, chưa có quan hệ tổng hợp đa biến.

Thứ tư, dự báo thử nghiệm ngư trường cá Ngừ Vằn được triển khai cho các tháng 01, 07, 10 và 12/2016 và kiểm chứng độc lập bằng số liệu khai thác thực tế 2016 theo sai số tuyệt đối (bảng 4.23) — bước chuyển từ tương quan mô tả sang năng lực dự báo nghiệp vụ có đánh giá.

Một kết quả phản trực giác đáng chú ý: khác kỳ vọng phổ biến rằng nhiệt độ luôn là biến trội, phân tích khẳng định với một số quan hệ (như cá Ngừ Vây Vàng với cấu trúc nhiệt đại dương trong tài liệu tham chiếu) ảnh hưởng "không thể giải thích được bằng trực tiếp nhiệt độ, mà có thể thông qua chuỗi thức ăn" — biện minh cho việc đưa chlorophyll-a và EKE vào mô hình.

Implications đa chiều

Về lý thuyết, kết quả củng cố đồng thời hai truyền thống: lý thuyết ổ sinh thái nhiệt (khoảng thích ứng của cá Ngừ Vằn 14,7–30,0°C vùng nhiệt đới so với 18,0–28,0°C ở biển Nhật Bản theo Kawasaki, 1991) và lý thuyết front–năng suất sinh học. Về phương pháp, quy trình đồng bộ dữ liệu viễn thám–nghề cá theo ô lưới 0,5 độ có thể chuyển giao trực tiếp cho các vùng biển khác của Việt Nam. Về thực tiễn, sản phẩm "phục vụ trực tiếp công tác dự báo ngư trường tại vùng biển Đông Nam Bộ, góp phần tạo thành một hệ thống hoàn chỉnh về dự báo ngư trường trên toàn vùng biển Việt Nam" — trực tiếp phục vụ Đề án "Tái cơ cấu ngành thủy sản đến năm 2020" của Tổng cục Thủy sản, vốn chỉ rõ nguồn lợi gần bờ bị khai thác quá mức trong khi tuyến lộng và xa bờ "chưa được đánh giá, dự báo chính xác".

Limitations và Future Research

Luận án bộc lộ và thừa nhận một số giới hạn cụ thể. Một là, dữ liệu nghề cá phân bố không đều theo tháng (cá Chỉ Vàng hoàn toàn thiếu số liệu tháng 3–4; cá Ngừ Vằn tháng 1 chỉ có 45 mẻ so với 820 mẻ tháng 4), dù "cũng đủ đảm bảo tính khách quan cho mục đích nghiên cứu". Hai là, dự báo thử nghiệm và kiểm chứng mới thực hiện cho một loài (cá Ngừ Vằn) trong một năm (2016), chưa phủ cá Ngừ Chấm và cá Chỉ Vàng. Ba là, chuỗi chlorophyll-a thực đo mỏng hơn hẳn nhiệt–muối (464 trạm so với 4.325 trạm), khiến phân tích thẳng đứng của biến này phụ thuộc số liệu 2011–2013. Bốn là, mô hình hồi quy tuyến tính chưa xử lý được các quan hệ phi tuyến ngưỡng vốn phổ biến trong sinh thái.

Hướng nghiên cứu tương lai gồm: mở rộng dự báo nghiệp vụ và kiểm chứng đa năm cho cá Ngừ Chấm và cá Chỉ Vàng; tích hợp bộ chỉ số HSI vào quy trình dự báo hạn 7–10 ngày hiện hành của hệ KC09; thử nghiệm các mô hình phi tuyến (GAM, machine learning) trên cùng bộ dữ liệu đồng bộ; bổ sung biến sinh học (động–thực vật phù du) theo gợi ý của Patrick Lehodey và nnk (1998); và mở rộng khung phân tích sang các loài kinh tế khác cùng vùng biển.

Tác động và ảnh hưởng

Về học thuật, đây là công trình đầu tiên xác lập bộ chỉ số thích ứng sinh thái theo mùa cho ba loài cá kinh tế tại ngư trường lớn nhất Việt Nam, trở thành tài liệu tham chiếu bắt buộc cho các nghiên cứu dự báo ngư trường tiếp theo trong nước. Về ngành, kết quả đóng góp trực tiếp vào cơ sở dữ liệu hải dương học nghề cá của Viện Nghiên cứu Hải Sản và cải tiến công nghệ dự báo của Trung tâm Dự báo Ngư trường Khai thác Hải sản — đơn vị đang phát hành bản dự báo nghiệp vụ cho ngư dân. Về chính sách, sản phẩm phục vụ bài toán được luận án nêu ở tầm quốc gia: "đánh giá nguồn lợi và dự báo ngư trường chính xác nhằm đưa ra các chính sách về quản lý, tổ chức khai thác hợp lý, hiệu quả". Về quốc tế, phương pháp tương thích chuẩn WCPFC, tạo tiền đề tích hợp dữ liệu Việt Nam vào các đánh giá nguồn lợi cá ngừ Trung–Tây Thái Bình Dương, nơi cá Ngừ Vằn chiếm khoảng 50% tổng sản lượng cá ngừ toàn cầu ~4 triệu tấn/năm (FAO, 1996).

Đối tượng hưởng lợi

Nghiên cứu sinh hải dương học và sinh thái nghề cá tìm thấy ở đây quy trình đồng bộ dữ liệu đa nguồn hoàn chỉnh cùng các khoảng trống đã được chỉ rõ (dự báo hạn ngắn theo loài, mô hình phi tuyến). Giới học thuật thâm niên có thể khai thác bộ ma trận tương quan mùa vụ để kiểm định lý thuyết ổ sinh thái động. Trung tâm dự báo ngư trường và doanh nghiệp khai thác xa bờ hưởng lợi trực tiếp từ bộ chỉ số HSI ba loài — đầu vào nâng chất lượng bản tin dự báo theo chuyến biển cho đội tàu lưới rê và lưới kéo. Nhà quản lý thủy sản có căn cứ định lượng để quy hoạch mùa vụ và vùng khai thác trên ngư trường chiếm tỷ trọng lớn nhất cả nước về khả năng khai thác cho phép (425.300 tấn/năm).

Câu hỏi chuyên sâu

1. Đóng góp lý thuyết độc đáo nhất là gì? Mở rộng khung HSI (Chen và nnk, 2008–2009) từ các biến bề mặt viễn thám sang hệ biến cấu trúc ba chiều (H0, H1, T1, H20, H24, Grad0, Grad50) cho vùng biển nhiệt đới gió mùa cận xích đạo, với điều kiện biên theo tầng sống của từng loài.

2. Đổi mới phương pháp so với nghiên cứu trước? So với Nguyễn Văn Lục và nnk (1991) chỉ dùng thống kê tần suất ba yếu tố, và so với mô hình KC09 chỉ tính từ trường 3D nhiệt độ, luận án kết hợp 14 biến đa nguồn (thực đo 1978–2013 + viễn thám MOVIMAR 1/12 độ) với hồi quy đa biến và HSI, kiểm chứng out-of-sample bằng dữ liệu 2016.

3. Phát hiện bất ngờ nhất? Độ muối VBĐNB biến thiên rất hẹp giữa hai mùa (chênh lệch mặt–đáy trung bình chỉ 1,48‰ mùa Đông Bắc), khiến vai trò phân dị ngư trường dồn về cấu trúc nhiệt và động lực (EKE, Alti) — trái với trực giác về vùng cửa sông Mê Kông.

4. Có quy trình tái lập không? Có: công thức tính từng biến, tiêu chí ngưỡng (Gz 0,05°C/m), định dạng ma trận đồng bộ (bảng 2.8–2.10) và quy tắc làm sạch dữ liệu được mô tả đủ chi tiết để tái lập trên vùng biển khác.

5. Chương trình nghiên cứu dài hạn? Kiến nghị của luận án mở lộ trình: hoàn thiện dự báo nghiệp vụ hạn ngắn cho cả ba loài, tích hợp vào hệ thống dự báo toàn vùng biển Việt Nam, và nâng cấp sang mô hình sinh thái đầy đủ kiểu SEAPODYM.

Kết luận

Luận án xác lập năm đóng góp cụ thể: (1) bộ dữ liệu đồng bộ cá–môi trường quy mô lớn nhất từng xây dựng cho VBĐNB (76.729 số liệu nhiệt độ, 162.569 điểm viễn thám/biến, 5.273 mẻ lưới ba loài); (2) bức tranh định lượng cấu trúc hải dương hai mùa gió của ngư trường lớn nhất cả nước; (3) các phương trình quan hệ CPUE–môi trường đa biến theo mùa cho cá Ngừ Vằn, Ngừ Chấm và Chỉ Vàng; (4) ba bộ chỉ số thích ứng sinh thái đầu tiên theo mùa cho ba loài này tại VBĐNB; (5) quy trình dự báo thử nghiệm có kiểm chứng độc lập năm 2016. Công trình đánh dấu bước chuyển của nghiên cứu nghề cá Việt Nam từ mô tả phân bố–sinh học sang dự báo định lượng dựa trên tương tác cá–môi trường, mở ít nhất ba dòng nghiên cứu mới (dự báo hạn ngắn theo loài ngoài cá ngừ đại dương, HSI đa biến cho cá tầng đáy ven bờ, tích hợp viễn thám độ phân giải cao vào nghiệp vụ), đồng thời kết nối nghề cá Việt Nam với chuẩn phương pháp của WCPFC và các nghiên cứu quốc tế từ Chile, Indonesia đến Ấn Độ Dương.