Tổng quan luận án

Luận án tiến sĩ hải dương học với đề tài "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ộ" được thực hiện bởi nghiên cứu sinh Nguyễn Văn Hướng dưới sự hướng dẫn khoa học của PGS.TS. Đoàn Văn Bộ và TS. Nguyễn Khắc Bát, bảo vệ năm 2018 tại Trường Đại học Khoa học Tự nhiên – Đại học Quốc gia Hà Nội (Mã số chuyên ngành: 62440228).

                      +-------------------------------------------------------------+
                      |         LUẬN ÁN TIẾN SĨ HẢI DƯƠNG HỌC (NĂM 2018)            |
                      |            NCS. Nguyễn Văn Hướng - ĐHKHTN Hà Nội            |
                      +-------------------------------------------------------------+
                                                     |
             +---------------------------------------+---------------------------------------+
             |                                                                               |
+--------------------------+                                                    +--------------------------+
|  ĐỐI TƯỢNG HẢI DƯƠNG     |                                                    |     ĐỐI TƯỢNG NGHỀ CÁ    |
| - Nhiệt độ, Độ muối      |                                                    | - Cá Ngừ Vằn (khơi)      |
| - Dòng chảy, EKE         |                                                    | - Cá Ngừ Chấm (ven lộng) |
| - Chlorophyll-a, Alti    |                                                    | - Cá Chỉ Vàng (đáy/nổi)  |
+--------------------------+                                                    +--------------------------+
             \                                                                               /
              \                                                                             /
               +---------------------------------------------------------------------------+
               |                   PHÂN TÍCH TƯƠNG QUAN & MÔ HÌNH HÓA                      |
               | - Không gian: Ô lưới 0,5° x 0,5° tại VBĐNB (06°00'-11°30'N, 105°00'-113°00'E) |
               | - Thời gian: Trung bình mùa gió (Đông Bắc / Tây Nam), tháng, tức thời      |
               | - Công cụ: Hồi quy tuyến tính đa biến (AIC) & Chỉ số thích ứng sinh thái SI|
               +---------------------------------------------------------------------------+
                                                     |
                                                     v
                               +-------------------------------------------+
                               |     KẾT QUẢ VÀ ỨNG DỤNG THỬ NGHIỆM        |
                               | - Bộ chỉ số sinh thái tối ưu theo mùa vụ  |
                               | - Dự báo thử nghiệm cá Ngừ Vằn năm 2016   |
                               +-------------------------------------------+

Tính cấp thiết và khoảng trống nghiên cứu

Vùng biển Đông Nam Bộ (VBĐNB) là ngư trường có trữ lượng và khả năng khai thác hải sản lớn nhất Việt Nam, với tổng trữ lượng nguồn lợi cá ước tính khoảng 976.000 tấn và khả năng khai thác cho phép là 425.300 tấn/năm. Mặc dù tổng sản lượng khai thác hải sản hàng năm có xu hướng gia tăng, nhưng số lượng tàu thuyền và chi phí khai thác tăng nhanh đã khiến năng suất khai thác tính trên đầu tàu và công suất giảm trung bình 3,1%/năm (theo Đề á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). Hoạt động khai thác của ngư dân phần lớn còn mang tính tự phát, quy mô nhỏ lẻ, dựa trên kinh nghiệm truyền thống; nguồn lợi thủy sản ven bờ bị khai thác quá mức trong khi nguồn lợi và ngư trường vùng khơi chưa được đánh giá và dự báo chuẩn xác.

Về mặt hải dương học, VBĐNB nằm trong khu vực khí hậu nhiệt đới gió mùa cận xích đạo, chịu tác động của dòng chảy lạnh từ phía bắc tràn xuống trong mùa gió Đông Bắc và dòng chảy ngược hướng kèm hiện tượng nước trồi Nam Trung Bộ trong mùa gió Tây Nam, cùng với lượng nước ngọt lớn từ hệ thống sông Mê Kông đổ ra. Cấu trúc hải dương học tại đây biến động phức tạp theo mùa, chi phối trực tiếp đến môi trường sống và sự phân bố của các loài cá kinh tế.

Trước thời điểm luận án thực hiện, tại 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 bằng nghề câu vàng ở vùng biển xa bờ miền Trung (thuộc đề tài cấp Nhà nước KC09.18/11-15) được triển khai nghiệp vụ. Tuy nhiên, mô hình này chỉ xây dựng dựa trên cấu trúc trường nhiệt 3D mà chưa tích hợp các yếu tố môi trường quan trọng khác (độ muối, chlorophyll-a, dòng chảy), đồng thời chưa áp dụng cho các đối tượng cá kinh tế khác tại VBĐNB. Do đó, việc nghiên cứu mối quan hệ giữa cấu trúc hải dương đa yếu tố với năng suất khai thác nhằm cung cấp cơ sở khoa học để thiết lập mô hình dự báo cho vùng biển này là yêu cầu cấp thiết.

Mục tiêu nghiên cứu

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.

Nội dung nghiên cứu

  1. Tổng quan về phương pháp dự báo ngư trường dựa trên quan hệ cá – môi trường, về điều kiện tự nhiên và một số đối tượng cá kinh tế ở vùng biển Đông Nam Bộ.
  2. Nghiên cứu, tính toán xác định cấu trúc các trường thủy động lực và môi trường biển (nhiệt độ, độ muối, chlorophyll-a, dòng chảy…) và biến động của chúng trong vùng biển nghiên cứu.
  3. Xác định khoảng thích ứng sinh thái các yếu tố môi trường biển đối với một số loài cá kinh tế (cá Ngừ Vằn, cá Ngừ Chấm và cá Chỉ Vàng) ở vùng biển Đông Nam Bộ.
  4. Nghiên cứu, đánh giá mối quan hệ giữa năng suất khai thác một số đối tượng cá kinh tế (cá Ngừ Vằn, cá Ngừ Chấm và cá Chỉ Vàng) với các yếu tố hải dương học, môi trường biển và ứng dụng trong dự báo thử nghiệm ngư trường ở vùng biển Đông Nam Bộ.

Đối tượng và phạm vi nghiên cứu

  • Đối tượng nghiên cứu: Các yếu tố cấu trúc hải dương học và môi trường biển (trường nhiệt, độ muối, dòng chảy, chlorophyll-a, mật độ động năng rối, độ cao mực nước biển) và 3 loài cá kinh tế đại diện:
    1. Cá Ngừ Vằn (Katsuwonus pelamis): Đại diện nhóm cá ngừ đại dương, cá nổi lớn vùng xa bờ.
    2. Cá Ngừ Chấm (Euthynnus affinis): Đại diện nhóm cá ngừ nhỏ ven bờ và tuyến lộng.
    3. Cá Chỉ Vàng (Selaroides leptolepis): Đại diện nhóm cá nổi nhỏ ven bờ, phân bố tầng đáy và tầng giữa.
  • Phạm vi không gian: Vùng biển Đông Nam Bộ, giới hạn trong tọa độ địa lý từ 06°00'N đến 11°30'N và từ 105°00'E đến 113°00'E, kéo dài từ bờ biển tỉnh Bình Thuận đến mũi Cà Mau.
  • Phạm vi thời gian: Chuỗi số liệu điều tra khảo sát thực tế giai đoạn 1978–2013; dữ liệu viễn thám và mô hình Mercator giai đoạn 2011–2016; dữ liệu năng suất khai thác nghề cá giai đoạn 1996–2016 (trong đó dữ liệu 1996–2015 dùng để xây dựng tương quan và bộ chỉ số sinh thái, dữ liệu năm 2016 dùng để kiểm chứng độc lập mô hình dự báo).

Tổng quan tài liệu và vị trí của luận án

Các hướng nghiên cứu dự báo ngư trường trên thế giới

Theo tổng quan của tác giả, lịch sử nghiên cứu đánh giá và dự báo biến động nguồn lợi cá trên thế giới tập trung vào 3 khuynh hướng chính:

  1. Khuynh hướng dựa vào nguyên lý Russell và thống kê nghề cá: Mô tả biến động số lượng cá thể đàn cá qua phương trình: $$\frac{dN}{dt} = W + R - (F + M) + \varepsilon$$ Trong đó $N$ là số lượng cá thể, $R$ là lượng bổ sung, $W$ là sinh trưởng, $F$ là hệ số chết do khai thác, $M$ là hệ số chết tự nhiên, và $\varepsilon$ là biến động ngẫu nhiên. Khuynh hướng này chia thành:
    • Mô hình giả định (đồng nhất sinh khối): Phương pháp thủy âm, diện tích, sản lượng thặng dư. Các tác giả Việt Nam áp dụng gồm: Nguyễn Khắc Bát, Bùi Đình Chung, Vũ Việt Hà, Nguyễn Viết Nghĩa, Phạm Thược, Đào Mạnh Sơn.
    • Mô hình giải tích: Xem xét cấu trúc tuổi và chiều dài như LCA, VPA, MULTIFAN-CL. Các tác giả Việt Nam áp dụng gồm: Đoàn Văn Bộ, Nguyễn Phi Đính, Nguyễn Xuân Huấn, Chu Tiến Vĩnh.
  2. Khuynh hướng dựa trên quá trình trao đổi năng lượng và dinh dưỡng: Nghiên cứu cân bằng sinh thái, lưới thức ăn và dòng năng lượng (ở Việt Nam có công trình của Nguyễn Tác An, Đoàn Văn Bộ, Nguyễn Tiến Cảnh). Hướng tiếp cận này chặt chẽ về mặt sinh thái nhưng gặp khó khăn khi mô hình hóa các thuật toán phức tạp.
  3. Khuynh hướng tương tác tổng hợp cá – môi trường – khai thác: Coi yếu tố hải dương học, sinh vật học và áp lực khai thác có vai trò bình đẳng dưới tác động không dừng của môi trường. Đây là hướng nghiên cứu hiện đại, được tác giả luận án lựa chọn kế thừa.

Các công trình quốc tế tiêu biểu được trích dẫn:

  • Silva và nnk (2000): Xác định mối quan hệ giữa sản lượng khai thác cá thu, cá cơm, cá trích ở bắc Chile với nhiệt độ bề mặt biển (SST) và chlorophyll-a.
  • Mansor và nnk (2001): Phát triển hệ thống hỗ trợ ra quyết định DSS (phần mềm TroFFS) dựa trên SST và chlorophyll-a từ viễn thám kết hợp GIS để tìm ngư trường cá tại Biển Đông.
  • Mukti và nnk (2011): Xác định khoảng thích ứng sinh thái của cá Ngừ Vằn tại Vịnh Bone-Flores với SST và chlorophyll-a viễn thám.
  • Nghiên cứu tại Đài Loan (2008): Hồi quy tuyến tính kết hợp phân tích thành phần chính (PCA) giữa CPUE câu vàng cá ngừ với cấu trúc nhiệt và chlorophyll-a ở Đại Tây Dương xích đạo.
  • Chen và nnk (2008, 2009): Ứng dụng hồi quy tuyến tính và mô hình thích ứng sinh thái HSI dự báo vùng đánh bắt cá Ngừ Mắt To (Thunnus obesus) ở Ấn Độ Dương và cá Thu Nhật Bản (Scomberomorus niphonius) ở biển Đông Trung Hoa.
  • WCPFC (2015): Ứng dụng mô hình sinh thái không gian SEAPODYM dự báo phân bố mật độ ($g/m^2$) cá ngừ đại dương hạn 7 ngày ở Tây Thái Bình Dương.

Tình hình nghiên cứu tại Việt Nam

  • Giai đoạn 1970–1980: Thiết lập dự báo ngư trường theo mùa vụ, theo quý phát trên Đài Tiếng nói Việt Nam (gián đoạn từ năm 1986 do thiếu kinh phí).
  • Nguyễn Phi Đính, Lê Trọng Phấn và Phạm Thược (1991) (trích bởi Lê Đức Tố, 1995): Tổng kết giới hạn sinh thái nhiệt độ, độ muối, độ sâu thích nghi của một số loài cá đáy và cá nục sò:
TT Loài cá Giới hạn nhiệt độ thích nghi / tối ưu (°C) Giới hạn độ muối thích nghi / tối ưu (‰) Độ sâu phân bố (m)
1 Mối vạch 15,0–30,0 / < 25,0 29,0–34,0 / 31,0–34,0 30–90
2 Trác ngắn 16,0–30,0 / 20,0–25,0 30,0–34,9 / 33,5–34,5 60–70
3 Hồng 15,0–32,0 / 21,0–23,0 30,0–34,9 / 33,0–34,0 > 30
4 Lượng 17,0–18,0 / - 31,0–33,0 / - > 20
5 Phèn 18,0–30,0 / 20,0–29,0 30,0–34,5 / 32,5–34,0 > 30
6 Nục sò - / 18,0–20,0 - / 31,0–34,5 Cá nổi
  • Lê Đức Tố và nnk (1995): Đặt nền móng luận cứ khoa học hải dương học nghề cá, nhấn mạnh vai trò của front, nước trồi và các khối nước.
  • Đinh Văn Ưu và nnk (2004): Mô hình hóa ngư trường cá ngừ đại dương quy mô lớn gắn với cấu trúc hải dương.
  • Đoàn Văn Bộ và nnk (2010, 2015): Xây dựng quy trình dự báo ngư trường câu vàng cá ngừ đại dương hạn tháng và hạn 7–10 ngày tại vùng biển xa bờ miền Trung sử dụng 26 biến cấu trúc nhiệt 3D.
  • Nguyễn Văn Lục và nnk (1991): Xác định khoảng thích ứng của cá Chỉ Vàng với độ dày lớp đồng nhất ($H_0 = 10-35m$), nhiệt độ đáy ($25,12-27,7^\circ C$), độ muối đáy ($32,58-33,45‰$) và chất đáy cát.

Khoảng trống nghiên cứu luận án giải quyết

Các công trình tại Việt Nam chủ yếu tập trung vào đặc điểm sinh học, trữ lượng chung hoặc dự báo nghề câu vàng cá ngừ đại dương dựa trên trường nhiệt. Luận án giải quyết khoảng trống về việc xác lập mối quan hệ định lượng đa biến đồng thời (bao gồm trường nhiệt, độ muối, dòng chảy, chlorophyll-a, động năng rối) và xây dựng bộ chỉ số thích ứng sinh thái chuyên biệt cho 3 loài cá kinh tế (Ngừ Vằn, Ngừ Chấm, Chỉ Vàng) tại vùng biển Đông Nam Bộ.

Cơ sở lý thuyết và phương pháp nghiên cứu

Cơ sở dữ liệu sử dụng trong nghiên cứu

Luận án khai thác và xử lý khối lượng dữ liệu thực nghiệm và dữ liệu số hóa quy mô lớn:

  1. Dữ liệu hải dương học thực đo (1978–2013): Thu thập từ các chuyến khảo sát của Viện Nghiên cứu Hải sản trên vùng biển Đông Nam Bộ.
Tháng Số trạm nhiệt độ Số liệu nhiệt độ Số trạm độ muối Số liệu độ muối Số trạm Chlorophyll-a Số liệu Chlorophyll-a
1 562 5.282 186 2.155 30 932
2 305 2.440 128 1.184 1 19
3 309 16.631 168 8.691 67 3.596
4 251 7.543 178 6.751 18 4.063
5 376 8.214 199 6.796 26 760
6 473 7.837 417 7.239 92 5.690
7 363 4.943 315 4.548 36 1.835
8 229 2.632 220 2.589 - -
9 363 5.914 339 5.753 - -
10 388 6.340 264 5.019 16 1.074
11 359 6.204 260 5.594 165 4.440
12 347 2.749 102 893 13 426
Tổng 4.325 76.729 2.776 57.212 464 22.835

  1. Dữ liệu viễn thám và mô hình Mercator (Dự án MOVIMAR, 2011–2016): Độ phân giải $1/12^\circ$ kinh vĩ. Tổng cộng 162.569 bản ghi số liệu cho mỗi yếu tố: nhiệt độ mặt biển, chlorophyll-a, dị thường độ cao mực biển, dòng chảy.

  2. Dữ liệu nghề cá (VietFish-Base & WCPFC, 1996–2016):

    • Cá Ngừ Vằn: Khai thác bằng lưới rê, đơn vị CPUE là kg/km lưới ($CPUE = C/L$). Tổng số 4.116 bản ghi (1.680 từ sổ nhật ký, 2.055 từ khảo sát, 381 từ giám sát). CPUE trung bình tháng dao động từ 10,83 kg/km lưới (tháng 2) đến 22,01 kg/km lưới (tháng 11).
    • Cá Ngừ Chấm: Khai thác bằng lưới rê, đơn vị CPUE là kg/km lưới. Tổng số 910 mẻ lưới; CPUE trung bình tháng từ 6,7 kg/km lưới (tháng 4) đến 27,4 kg/km lưới (tháng 9).
    • Cá Chỉ Vàng: Khai thác bằng lưới kéo đáy, đơn vị CPUE là kg/giờ kéo ($CPUE = C/t$). Tổng số 247 mẻ lưới; CPUE trung bình tháng từ 0,5 kg/giờ (tháng 5, 8) đến 9,1 kg/giờ (tháng 12). Theo phân bố độ sâu khai thác: < 20m chiếm 8,1%; 20–30m chiếm 19,8%; 30–50m chiếm 62,3%; 50–100m chiếm 9,7%.

Các biến môi trường được lựa chọn

  • Đối với Cá Ngừ Vằn và Cá Ngừ Chấm (14 yếu tố trong tầng 0–150m):

    1. $T_0$: Nhiệt độ nước biển tầng mặt (°C).
    2. $Ano$: Dị thường nhiệt độ nước biển tầng mặt (°C).
    3. $H_0$: Độ dày lớp đồng nhất trên (m).
    4. $H_1$: Độ sâu biên dưới lớp đột biến nhiệt (m).
    5. $T_1$: Nhiệt độ biên dưới lớp đột biến nhiệt (°C).
    6. $H_{20}$: Độ sâu mặt đẳng nhiệt 20°C (m).
    7. $H_{24}$: Độ sâu mặt đẳng nhiệt 24°C (m).
    8. $Grad_0$: Gradien cực đại theo phương ngang nhiệt bề mặt (°C/10km).
    9. $Grad_{50}$: Gradien cực đại theo phương ngang nhiệt tầng 50m (°C/10km).
    10. $Sal_0$: Độ muối nước biển tầng mặt (‰).
    11. $Chlo$: Hàm lượng chlorophyll-a tầng mặt ($mg/m^3$).
    12. $Alti$: Dị thường độ cao động lực bề mặt biển (cm).
    13. $EKE$: Mật độ động năng rối ($cm^2/s^2$).
    14. $Spd_cur$: Tốc độ dòng chảy (cm/s).

  • Đối với Cá Chỉ Vàng (10 yếu tố trong tầng 0–50m): $T_0$, $Ano$, $H_0$, $Grad_0$, $Grad_{25}$ (Gradien nhiệt tầng 25m), $Sal_0$, $Chlo$, $Alti$, $EKE$, $Spd_cur$.

Phương pháp tính toán và xử lý số liệu

  • Tính toán cấu trúc nhiệt:
    • Lớp đồng nhất trên ($H_0$): Xác định theo tiêu chí gradien nhiệt thẳng đứng $G_z = \frac{|T_{H0} - T_0|}{H_0} < 0,05^\circ\text{C/m}$.
    • Lớp đột biến nhiệt ($H_1$): Tiếp giáp dưới lớp đồng nhất, thỏa mãn $G_z = \frac{|T_1 - T_{H0}|}{H_1 - H_0} > 0,05^\circ\text{C/m}$.
    • Mặt đẳng nhiệt ($H_{20}, H_{24}$): Nội suy tuyến tính theo phương thẳng đứng từ các tầng chuẩn.
    • Gradien nhiệt nằm ngang tầng H ($Grad_H$): $Grad_H = \sqrt{(G_{xH})^2 + (G_{yH})^2}$.
  • Tính toán trường động lực biển:
    • Vận tốc dòng chảy bề mặt $(u, v)$ tính từ dị thường mực biển: $$u = -\left(\frac{g}{f}\right)\frac{\partial z}{\partial y}, \quad v = \left(\frac{g}{f}\right)\frac{\partial z}{\partial x}$$ với $g = 980,\text{cm/s}^2$, $f = 2\Omega\sin\Phi$, $\Omega = 7,29 \times 10^{-5},\text{rad/s}$.
    • Mật độ động năng rối ($EKE$) xác định từ các thành phần vận tốc dòng chảy biến động.
  • Quy trình đồng bộ dữ liệu:
    • Kiểm tra sai số thô và áp dụng quy tắc $3\sigma$ loại trừ dị thường.
    • Không gian: Quy về tâm ô lưới $0,5^\circ \times 0,5^\circ$ kinh vĩ.
    • Thời gian: Đồng bộ theo 3 cấp độ: trung bình mùa gió (mùa Đông Bắc: tháng 10 – tháng 4 năm sau; mùa Tây Nam: tháng 4 – tháng 9), trung bình 12 tháng nhiều năm, và đồng bộ tức thời theo ngày giai đoạn 2012–2015.
  • Mô hình thống kê và mô hình sinh thái:
    • Hồi quy tuyến tính nhiều biến (Multiple Linear Regression): Kết hợp tiêu chuẩn thông tin Akaike (AIC) để chọn lọc tổ hợp biến môi trường tối ưu có tương quan cao nhất với CPUE.
    • Mô hình chỉ số thích ứng sinh thái HSI (Habitat Suitability Index): Thiết lập hàm Suitability Index ($SI \in [0, 1]$) cho từng yếu tố môi trường đơn lẻ dựa trên phân bố tần suất CPUE, từ đó tổng hợp thành chỉ số $HSI$ phục vụ khoanh vùng ngư trường tiềm năng.

Nội dung chính theo từng chương

Chương 1: Tổng quan về phương pháp dự báo ngư trường dựa trên quan hệ cá – môi trường, điều kiện tự nhiên và một số đối tượng cá kinh tế vùng biển Đông Nam Bộ

Chương này tổng kết cơ sở lý luận về hải dương học nghề cá và bức tranh tự nhiên của khu vực nghiên cứu:

  • Điều kiện địa hình và trầm tích: Vùng biển có bờ biển kéo dài từ Bình Thuận đến mũi Cà Mau. Độ sâu 0–15m có dạng bờ dốc hẹp; vùng 15–50m trải rộng và thoải; vùng sâu > 50m có độ dốc tăng dần, đường đẳng sâu 100m mở rộng cách bờ tới 300 hải lý. Trầm tích bề mặt đáy phân hóa từ cát hạt lớn ven bờ Phan Thiết – Hàm Tân, cát hạt mịn chiếm ưu thế toàn vùng, bùn cát trước cửa sông Tiền, sông Hậu, đến trầm tích bùn sét và sét biển khơi.
  • Khí tượng – thủy văn: Khí hậu nhiệt đới gió mùa cận xích đạo. Mùa mưa (tháng 4/5 đến tháng 10/11) chiếm trên 90% lượng mưa cả năm. Chế độ gió gồm gió mùa Đông Bắc (tháng 11 – tháng 3, tốc độ gió trung bình trạm Phú Quý đạt 8,6 m/s vào tháng 1) và gió mùa Tây Nam (tháng 5 – tháng 10, tốc độ trung bình 2,5–3,5 m/s). Nhiệt độ không khí trung bình năm dao động 27,1–27,7°C.
  • Đặc tính sinh học 3 loài nghiên cứu:
    • Cá Ngừ Vằn: Chiếm > 60% sản lượng cá nổi lớn tại VBĐNB và miền Trung; kích thước khai thác phổ biến 40–80cm; nhiệt độ thích nghi chung 14,7–30,0°C, nhiệt độ đẻ trứng tối thiểu 25°C.
    • Cá Ngừ Chấm: Loài cá ngừ nhỏ ven bờ, kích thước khai thác phổ biến 24–45cm (chiếm ưu thế ở 36cm); nhiệt độ nước thích nghi 18,0–29,0°C.
    • Cá Chỉ Vàng: Phân bố tầng đáy và giữa ở độ sâu < 50m; kích thước khai thác trung bình 11–14cm; thức ăn chính là động vật phù du; bổ sung quần thể mạnh vào mùa gió Tây Nam.

Chương 2: Dữ liệu hải dương học, nghề cá và các phương pháp phân tích quan hệ cá-môi trường

Chương 2 mô tả chi tiết nguồn gốc, phương pháp thu thập, cấu trúc các trường dữ liệu và các bước tiền xử lý toán học:

  • Hệ thống hóa 76.729 số liệu nhiệt độ, 57.212 số liệu độ muối, 22.835 số liệu chlorophyll-a thực đo cùng 162.569 bản ghi viễn thám Movimar.
  • Quy chuẩn hóa cơ sở dữ liệu nghề cá VietFish-Base và WCPFC đối với 3 loài cá mục tiêu; chi tiết hóa phương pháp chuẩn hóa nỗ lực khai thác thành CPUE.
  • Thiết lập thuật toán phân tích gradien nhiệt, độ sâu lớp đồng nhất ($H_0$), biên dưới lớp đột biến ($H_1$), các mặt đẳng nhiệt ($H_{20}, H_{24}$), dòng chảy địa chuyển và năng lượng rối ($EKE$).
  • Xây dựng quy trình đồng bộ hóa dữ liệu cá – môi trường trên lưới không gian $0,5^\circ \times 0,5^\circ$ theo 3 thang thời gian (mùa, tháng, tức thời) và cấu trúc ma trận đồng bộ phục vụ phân tích tương quan và mô hình HSI.

Chương 3: Đặc điểm một số yếu tố hải dương học, môi trường ở vùng biển Đông Nam Bộ

Chương 3 trình bày các kết quả phân tích cấu trúc không gian và biến trình thời gian của các trường hải dương học:

  • Trường nhiệt: Nhiệt độ bề mặt biển ($T_0$) biến trình năm dạng 1 cực đại, 1 cực tiểu; thấp nhất vào tháng 1 (25,0–27,0°C ở ven bờ, tăng dần ra khơi) và cao nhất vào tháng 5 (29,0–30,5°C). Trong mùa gió Tây Nam (tháng 7–8), xuất hiện vùng nhiệt độ thấp cục bộ (27,5–28,5°C) ở ven bờ Bình Thuận do ảnh hưởng của dòng nước trồi Nam Trung Bộ lan tỏa xuống. Độ dày lớp đồng nhất ($H_0$) đạt giá trị lớn nhất vào mùa đông (tháng 12 – tháng 1, $H_0 = 40-70m$) do quá trình xáo trộn gió mạnh; đạt giá trị nhỏ nhất vào thời kỳ chuyển tiếp và mùa hè ($H_0 = 10-25m$).
  • Trường độ muối: Độ muối tầng mặt ($Sal_0$) biến động từ 29,0‰ đến 34,0‰. Mùa khô (tháng 12 đến tháng 5) độ muối cao và đồng nhất (32,0–34,0‰); mùa mưa (tháng 7–10) độ muối tầng mặt giảm mạnh ở vùng ven bờ cửa sông Mê Kông (xuống < 30,0‰) và hình thành gradien độ muối ngang rõ rệt.
  • Trường dòng chảy và EKE: Mùa gió Đông Bắc, dòng chảy tầng mặt có hướng chủ đạo Bắc – Nam và Đông Bắc – Tây Nam với tốc độ trung bình 20–30 cm/s ngoài khơi và 15–20 cm/s ven bờ. Mùa gió Tây Nam, dòng chảy đổi hướng sang Tây Nam – Đông Bắc với tốc độ trung bình 10–15 cm/s. Mật độ động năng rối ($EKE$) đạt cực đại tại các vùng giao thoa dòng chảy và khu vực thềm lục địa dốc phía ngoài khơi.
  • Hàm lượng Chlorophyll-a: Biến trình năm cho thấy chlorophyll-a tầng mặt đạt giá trị cao trong mùa gió Đông Bắc và tại vùng nước trồi ven bờ mùa hè (trung bình 0,3–1,5 $mg/m^3$, tại các cửa sông có thể > 2,0 $mg/m^3$). Vùng khơi duy trì giá trị thấp đặc trưng cho thủy vực nhiệt đới nghèo dinh dưỡng (0,08–0,25 $mg/m^3$). Phân bố thẳng đứng cho thấy cực đại chlorophyll-a dưới tầng mặt thường nằm ở độ sâu 30–60m.

Chương 4: 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ộ

Chương 4 là chương kết quả cốt lõi của luận án, xác lập các mô hình định lượng và thử nghiệm dự báo:

  • Đối với Cá Ngừ Vằn:
    • Ma trận tương quan chỉ ra CPUE cá Ngừ Vằn tương quan thuận chặt chẽ với độ dày lớp đồng nhất ($H_0$), độ sâu mặt đẳng nhiệt $H_{20}$, và tương quan nghịch với nhiệt độ bề mặt ($T_0$) trong mùa đông.
    • Bộ chỉ số thích ứng sinh thái ($SI$) mùa Đông Bắc: $T_0$ thích hợp từ 26,0–27,5°C; $H_0$ từ 40–60m; $Sal_0$ từ 33,0–34,0‰; chlorophyll-a từ 0,15–0,35 $mg/m^3$.
    • Mùa Tây Nam: $T_0$ thích hợp từ 28,5–29,8°C; $H_0$ từ 15–35m; $EKE$ từ 100–300 $cm^2/s^2$.
  • Đối với Cá Ngừ Chấm:
    • CPUE cá Ngừ Chấm có mối liên hệ mật thiết với gradien nhiệt độ mặt biển ($Grad_0$), độ muối tầng mặt ($Sal_0$) và hàm lượng chlorophyll-a tại vùng biển ven bờ và tuyến lộng.
    • Bộ chỉ số $SI$ xác định khoảng nhiệt độ thuận lợi trong mùa đông là 26,5–28,0°C, độ muối 31,5–33,2‰, chlorophyll-a từ 0,25–0,60 $mg/m^3$. Mùa hè, loài này tập trung tại vùng nước có nhiệt độ 28,0–29,5°C và các đới front thủy văn ven bờ.
  • Đối với Cá Chỉ Vàng:
    • Năng suất khai thác bằng lưới kéo đáy tương quan cao với cấu trúc nhiệt tầng đáy và tầng 25m ($Grad_{25}$), độ dày $H_0$ và hàm lượng chlorophyll-a tầng mặt.
    • Khoảng thích ứng sinh thái tối ưu mùa Đông Bắc: $T_0$ từ 26,0–27,5°C; $H_0$ từ 20–35m; $Sal_0$ từ 31,5–33,0‰. Mùa Tây Nam: $T_0$ từ 28,0–29,5°C; $H_0$ từ 10–25m.
  • Dự báo thử nghiệm ngư trường cá Ngừ Vằn:
    • Xây dựng bản đồ dự báo ngư trường tiềm năng cho các tháng 1, 7, 10, 12 năm 2016 tại VBĐNB dựa trên mô hình tích hợp tương quan và HSI.
    • Đánh giá kiểm chứng độc lập với số liệu khai thác thực tế năm 2016 bằng phương pháp sai số tuyệt đối, kết quả kiểm chứng cho thấy các vị trí dự báo bãi cá có mức độ trùng khớp tốt với các mẻ lưới có CPUE cao trong thực tế.

Kết quả và những đóng góp mới

Đóng góp mới về mặt khoa học

  1. Xác định một cách có hệ thống và định lượng mối quan hệ giữa năng suất khai thác của 3 loài cá kinh tế đại diện (Cá Ngừ Vằn, Cá Ngừ Chấm, Cá Chỉ Vàng) với hệ thống 10–14 yếu tố cấu trúc hải dương học đặc trưng (cấu trúc nhiệt thẳng đứng $H_0, H_1, H_{20}, H_{24}$, gradien nhiệt ngang, độ muối, dòng chảy, chlorophyll-a, $EKE$, dị thường mực biển) tại vùng biển Đông Nam Bộ theo cả thang thời gian mùa gió, tháng và dữ liệu tức thời.
  2. Thiết lập được bộ chỉ số thích ứng sinh thái (Suitability Index – SI) của từng yếu tố môi trường đối với cá Ngừ Vằn, cá Ngừ Chấm và cá Chỉ Vàng phân theo hai mùa gió chính (Đông Bắc và Tây Nam), làm cơ sở khoa học cho việc ứng dụng mô hình HSI trong dự báo ngư trường.

Đóng góp mới về mặt thực tiễn

  1. Cung cấp quy trình công nghệ và thuật toán đồng bộ dữ liệu viễn thám – mô hình số trị – số liệu điều tra thực tế với dữ liệu nghề cá (CPUE), cho phép khai thác hiệu quả nguồn dữ liệu vệ tinh Movimar độ phân giải cao phục vụ nghề cá.
  2. Xây dựng và kiểm chứng thành công mô hình dự báo thử nghiệm ngư trường khai thác cá Ngừ Vằn tại vùng biển Đông Nam Bộ (năm 2016), góp phần hoàn thiện mạng lưới dự báo ngư trường khai thác hải sản trên toàn vùng biển Việt Nam.

Kiến nghị của luận án

  • Cần tiếp tục duy trì và mở rộng mạng lưới thu thập số liệu giám sát nghề cá, đặc biệt là sổ nhật ký khai thác đồng bộ tọa độ và thời gian mẻ lưới để cập nhật cơ sở dữ liệu VietFish-Base.
  • Tăng cường nghiên cứu ứng dụng dữ liệu viễn thám thời gian thực kết hợp mô hình dự báo hạn ngắn (7–10 ngày) phục vụ trực tiếp cho các đội tàu khai thác xa bờ tại vùng biển Đông Nam Bộ.

Hạn chế và hướng nghiên cứu tiếp

Hạn chế được ghi nhận từ văn bản

  • Nguồn số liệu nghề cá (CPUE) phân bố chưa thực sự đồng đều giữa các tháng trong năm; số lượng mẻ lưới quan trắc cá Chỉ Vàng (247 mẻ) và cá Ngừ Chấm (910 mẻ) còn ít hơn đáng kể so với cá Ngừ Vằn (4.116 mẻ), một số tháng mùa khô thiếu hụt mẫu quan sát cá Chỉ Vàng (tháng 3, tháng 4 không có số liệu trong CSDL).
  • Nguồn dữ liệu hải dương học viễn thám và mô hình Mercator dù có độ phủ rộng nhưng vẫn tồn tại sai số nhất định so với số liệu thực đo CTD (sai lệch nhiệt độ khoảng $\pm 0,5^\circ\text{C}$, độ muối khoảng $\pm 0,2‰$).
  • Thử nghiệm dự báo ngư trường nghiệp vụ và kiểm chứng độc lập mới chỉ thực hiện được cho đối tượng cá Ngừ Vằn trong năm 2016, chưa triển khai kiểm chứng thực địa đầy đủ cho cá Ngừ Chấm và cá Chỉ Vàng.

Hướng nghiên cứu tiếp

  • Bổ sung chuỗi số liệu điều tra nghề cá liên tục nhiều năm cho nhóm cá ngừ nhỏ ven bờ và cá nổi nhỏ tầng đáy để nâng cao độ chính xác của các ma trận tương quan.
  • Hoàn thiện phần mềm dự báo tự động và chuyển giao quy trình dự báo ngư trường hạn ngắn (hạn tuần, hạn 7–10 ngày) cho các đối tượng cá Ngừ Chấm và cá Chỉ Vàng tại vùng biển Nam Bộ.
  • Tích hợp thêm các yếu tố sinh học bậc cao và tương tác chuỗi thức ăn vào mô hình dự báo sinh thái không gian.

Giá trị tham khảo

Đối tượng thụ hưởng

  • Nghiên cứu sinh, học viên cao học và giảng viên: Thuộc các chuyên ngành Hải dương học, Quản lý nghề cá, Sinh học biển, Kỹ thuật khai thác thủy sản tại các trường đại học và viện nghiên cứu chuyên ngành.
  • Cơ quan quản lý và dự báo nghề cá: Viện Nghiên cứu Hải sản, Trung tâm Dự báo Ngư trường Khai thác Hải sản, Tổng cục Thủy sản, Chi cục Thủy sản các tỉnh ven biển Đông Nam Bộ (Bình Thuận, Bà Rịa – Vũng Tàu, Bến Tre, Trà Vinh, Sóc Trăng, Bạc Liêu, Cà Mau).
  • Ngư dân và các hợp tác xã khai thác hải sản xa bờ: Sử dụng thông tin ngư trường dự báo để tối ưu hóa hải trình chuyến biển, tiết kiệm nhiên liệu và nâng cao sản lượng đánh bắt.

Nội dung có giá trị tham khảo cao

  • Quy trình đồng bộ dữ liệu đa nguồn: Phương pháp kết nối ma trận dữ liệu viễn thám hàng ngày (Movimar) với dữ liệu tọa độ mẻ lưới trên lưới ô vuông $0,5^\circ \times 0,5^\circ$.
  • Bảng chỉ số thích ứng sinh thái ($SI$): Các ngưỡng nhiệt độ ($T_0$), độ dày lớp đồng nhất ($H_0$), độ muối ($Sal_0$), chlorophyll-a tối ưu theo mùa cho 3 loài cá kinh tế quan trọng.
  • Bộ công thức toán học tính toán cấu trúc nhiệt và động năng rối: Thuật toán xác định $H_0, H_1$, mặt đẳng nhiệt $H_{20}, H_{24}$, dòng chảy địa chuyển $(u, v)$ và $EKE$ từ dữ liệu độ cao mực nước biển.

Câu hỏi thường gặp

1. Luận án sử dụng những nguồn dữ liệu nào và khối lượng cụ thể ra sao?

Luận án sử dụng 3 nguồn dữ liệu chính:

  • Dữ liệu khảo sát hải dương học thực đo (1978–2013): 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; 464 trạm với 22.835 số liệu chlorophyll-a.
  • Dữ liệu viễn thám và mô hình Mercator (dự án MOVIMAR, 2011–2016): 162.569 bản ghi số liệu cho mỗi yếu tố (nhiệt độ, chlorophyll-a, dị thường mực biển, dòng chảy).
  • Dữ liệu nghề cá (VietFish-Base & WCPFC, 1996–2016): 4.116 bản ghi cá Ngừ Vằn (lưới rê), 910 bản ghi cá Ngừ Chấm (lưới rê) và 247 bản ghi cá Chỉ Vàng (lưới kéo đáy).

2. Ba loài cá được lựa chọn nghiên cứu đại diện cho những nhóm nguồn lợi sinh thái nào?

  • Cá Ngừ Vằn (Katsuwonus pelamis): Đại diện cho nhóm cá ngừ đại dương, cá nổi lớn di cư xa bờ.
  • Cá Ngừ Chấm (Euthynnus affinis): Đại diện cho nhóm cá ngừ nhỏ (neritic tuna) sống ở vùng ven bờ và tuyến lộng.
  • Cá Chỉ Vàng (Selaroides leptolepis): Đại diện cho nhóm cá nổi nhỏ sống ở vùng biển nông ven bờ (chủ yếu tập trung ở dải độ sâu 30–50m) và gắn liền với nền đáy mềm.

3. Những biến cấu trúc hải dương nào được chọn để phân tích tương quan với năng suất khai thác?

  • Đối với cá Ngừ Vằn và cá Ngừ Chấm: Phân tích 14 yếu tố trong tầng nước 0–150m ($T_0, Ano, H_0, H_1, T_1, H_{20}, H_{24}, Grad_0, Grad_{50}, Sal_0, Chlo, Alti, EKE, Spd_cur$).
  • Đối với cá Chỉ Vàng: Phân tích 10 yếu tố trong tầng nước 0–50m ($T_0, Ano, H_0, Grad_0, Grad_{25}, Sal_0, Chlo, Alti, EKE, Spd_cur$).

4. Tiêu chí xác định độ dày lớp đồng nhất nhiệt ($H_0$) và lớp đột biến nhiệt ($H_1$) trong luận án là gì?

  • Độ dày lớp đồng nhất trên ($H_0$): Lớp nước tiếp giáp mặt biển có gradien nhiệt thẳng đứng $G_z = \frac{|T_{H0} - T_0|}{H_0} < 0,05^\circ\text{C/m}$.
  • Biên dưới lớp đột biến ($H_1$): Lớp nước tiếp giáp dưới lớp đồng nhất có gradien nhiệt thẳng đứng $G_z = \frac{|T_1 - T_{H0}|}{H_1 - H_0} > 0,05^\circ\text{C/m}$.

5. Kết quả dự báo thử nghiệm ngư trường được tiến hành cho loài cá nào và vào thời gian nào?

Dự báo thử nghiệm ngư trường được xây dựng cho loài cá Ngừ Vằn (Katsuwonus pelamis) tại vùng biển Đông Nam Bộ cho các tháng 1, 7, 10 và 12 năm 2016, sau đó được đánh giá kiểm chứng độc lập bằng phương pháp tính sai số tuyệt đối so với số liệu khai thác thực tế thu nhận trong năm 2016.

Kết luận

Luận án tiến sĩ của tác giả Nguyễn Văn Hướng đã giải quyết bài toán khoa học về việc xác lập mối quan hệ định lượng giữa các yếu tố cấu trúc hải dương học đa biến với năng suất khai thác của ba loài cá kinh tế đại diện (Ngừ Vằn, Ngừ Chấm, Chỉ Vàng) tại vùng biển Đông Nam Bộ. Bằng việc tích hợp chuỗi số liệu khảo sát thực tế nhiều năm, dữ liệu viễn thám Movimar độ phân giải cao và cơ sở dữ liệu nghề cá chuẩn hóa, nghiên cứu đã xây dựng thành công bộ chỉ số thích ứng sinh thái ($SI$) theo mùa gió và mô hình thực nghiệm dự báo ngư trường. Kết quả thử nghiệm và kiểm chứng độc lập trên đối tượng cá Ngừ Vằn năm 2016 khẳng định độ tin cậy của phương pháp tiếp cận tương quan cá – môi trường, tạo tiền đề khoa học vững chắc để phát triển công nghệ dự báo ngư trường nghiệp vụ hạn ngắn phục vụ khai thác hải sản bền vững tại Việt Nam.