Tổng quan luận án

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

Vùng biển Đông Nam Bộ là ngư trường có trữ lượng và khả năng khai thác nguồn lợi hải sản lớn nhất Việt Nam. Trữ lượng nguồn lợi cá tại đây ước tính đạt khoảng 976.000 tấn với khả năng khai thác cho phép là 425.300 tấn/năm, trong đó cá nổi nhỏ chiếm khoảng 349.000 tấn (khai thác cho phép 174.000 tấn), cá đáy chiếm khoảng 627.000 tấn (khai thác cho phép 250.000 tấn), và nhóm cá nổi lớn ước tính khoảng 300.000 tấn với cá Ngừ Vằn chiếm trên 60%. Khu vực này nằm trong vùng khí hậu nhiệt đới gió mùa cận xích đạo, chịu sự tác động đồng thời của chế độ thủy văn sông Mê Kông, hệ dòng chảy mùa Đông Bắc từ phía bắc xuống và dòng chảy mùa Tây Nam cùng hiện tượng nước trồi khu vực Nam Trung Bộ.

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, dù sản lượng khai thác hải sản 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 đã làm cho năng suất khai thác trung bình giảm khoảng 3,1%/năm. Hoạt động khai thác của ngư dân chủ yếu dựa vào kinh nghiệm, quy mô nhỏ lẻ, nguồn lợi gần bờ bị khai thác quá mức trong khi nguồn lợi tuyến lộng và vùng xa bờ chưa được đánh giá và dự báo chính xác.

Tại Việt Nam, các mô hình dự báo ngư trường phục vụ khai thác nghiệp vụ trước đây chủ yếu tập trung vào đối tượng cá Ngừ Vây Vàng và cá Ngừ Mắt To thuộc nghề câu cá ngừ đại dương ở vùng biển xa bờ miền Trung (như đề tài cấp Nhà nước KC09.18/11-15). Ngược lại, các đối tượng cá kinh tế khác chiếm tỷ trọng lớn tại vùng biển Đông Nam Bộ như cá Ngừ Vằn (Katsuwonus pelamis - đại diện nhóm cá ngừ đại dương), cá Ngừ Chấm (Euthynnus affinis - đại diện nhóm cá ngừ nhỏ ven bờ) và cá Chỉ Vàng (Selaroides leptolepis - đại diện nhóm cá nổi nhỏ ven bờ) chưa có nghiên cứu đầy đủ về cơ sở khoa học gắn kết giữa trường môi trường hải dương và biến động năng suất khai thác để phục vụ dự báo ngư trường hạn ngắn. Luận án "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ủa nghiên cứu sinh Nguyễn Văn Hướng được thực hiện nhằm giải quyết khoảng trống khoa học này.

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ấu trúc các yếu tố hải dương học, môi trường biển (trường nhiệt 3D, độ muối, chlorophyll-a, dòng chảy, xoáy động năng) và năng suất khai thác của 3 loài cá kinh tế: cá Ngừ Vằn (Katsuwonus pelamis), cá Ngừ Chấm (Euthynnus affinis) và cá Chỉ Vàng (Selaroides leptolepis).
  • Phạm vi không gian: Vùng biển Đông Nam Bộ, giới hạn 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 hải dương học thực đo giai đoạn 1978–2013; số liệu viễn thám biển và mô hình Mercator từ dự án MOVIMAR giai đoạn 2011–2016; số liệu nghề cá thực tế giai đoạn 1996–2016.

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

Luận án hệ thống hóa 3 khuynh hướng nghiên cứu đánh giá trữ lượng và dự báo khai thác nguồn lợi cá trên thế giới:

  1. Khuynh hướng dựa vào nguyên lý Russell và thống kê nghề cá: Phương trình biến động số lượng $dN/dt = W + R - (F + M) + \Phi$. Nhóm mô hình giả định (thủy âm, diện tích, sản lượng thặng dư) được thực hiện tại Việt Nam bởi 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. Nhóm mô hình giải tích (LCA, VPA, MULTIFAN-CL) được nghiên cứu bởi Đoàn Văn Bộ, Nguyễn Phi Đính, Nguyễn Xuân Huấn, Chu Tiến Vinh.
  2. Khuynh hướng quá trình trao đổi năng lượng (dinh dưỡng): Đánh giá quan hệ lưới thức ăn và chuyển hóa năng lượng sinh thái biển (Nguyễn Tác An, Đoàn Văn Bộ, Nguyễn Tiến Cảnh). Hướng tiếp cận này gặp khó khăn do tính phức tạp của thuật toán và chuỗi số liệu sinh thái.
  3. Khuynh hướng tương tác tổng hợp cá - môi trường - khai thác: Xem xét vai trò tương đương giữa điều kiện hải dương học, sinh học và áp lực khai thác. Đây là hướng nghiên cứu hiện đại, tận dụng công nghệ viễn thám và hệ thống thông tin địa lý (GIS), được tác giả lựa chọn làm nền tảng cho luận án.

Tổng quan các công trình quốc tế và trong nước

  • Nghiên cứu quốc tế:
    • Silva và nnk (2000): Xác định mối quan hệ giữa sản lượng 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 TroFFS/DSS tại Biển Đông dựa trên phân tích ảnh viễn thám SST và chlorophyll-a qua GIS.
    • 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 qua dữ liệu viễn thám nhiệt độ và chlorophyll-a.
    • Chen và nnk (2008, 2009): Ứng dụng hồi quy tuyến tính và mô hình chỉ số thích ứng môi trường sống (HSI) dự báo cá Ngừ Mắt To (Thunnus obesus) ở Ấn Độ Dương và cá thu Nhật Bản (Scomberomorus niphonius) ở biển Đông Trung Hoa.
    • Học viện Công nghệ Châu Á, Nhật Bản và Sri Lanka (2008, 2009): Xác định chỉ số thích ứng của cá Ngừ Vây Vàng ở Ấn Độ Dương với SST (28–30°C), Chlorophyll-a (0,1–0,4 mg/m³) và độ cao mặt biển SSH (205–215 cm).
    • Tổ chức Nghề cá Trung - Tây Thái Bình Dương (WCPFC, 2015): Ứng dụng mô hình sinh thái SEAPODYM dự báo phân bố mật độ sinh khối cá ngừ đại dương với thời hạn 7 ngày.
  • Nghiên cứu trong nước:
    • Nguyễn Phi Đính, Lê Trọng Phấn, Phạm Thược (1991, trích bởi Lê Đức Tố, 1995): Thống kê sinh thái về dải nhiệt độ, độ muối, độ sâu thích ứng của cá Mối vạch, Trác ngắn, Phèn, Lượng, Hố và Nục sồ.
    • Lê Đức Tố và nnk (1995): Xác lập cơ sở khoa học hải dương học nghề cá, chỉ ra ảnh hưởng của nước trồi, đới front và các khối nước tới phân bố cá.
    • Đinh Văn Ưu và nnk (2004): Xây dựng mô hình dự báo trường nhiệt và ngư trường cá ngừ câu vàng quy mô lớn.
    • Đoàn Văn Bộ và nnk (2010, 2015 - đề tài KC09.18/11-15): Xây dựng quy trình dự báo ngư trường cá ngừ đại dương hạn tháng và hạn 7–10 ngày dựa trên 26 biến cấu trúc nhiệt 3D.
    • Nguyễn Văn Lục và nnk (1991): Xác định ngưỡng sinh thái cho cá Chỉ Vàng (độ dày lớp đồng nhất $H_0$ từ 10–35 m, nhiệt độ đáy 25,12–27,7°C, độ muối đáy 32,58–33,45‰, chất đáy cát 48–89%).

Khoảng trống luận án giải quyết

Mô hình dự báo của Đoàn Văn Bộ (2015) chỉ sử dụng các biến số trích xuất từ trường nhiệt độ 3D và giới hạn ở đối tượng cá ngừ đại dương xa bờ miền Trung. Luận án giải quyết việc mở rộng nghiên cứu sang vùng biển Đông Nam Bộ, tích hợp đồng thời trường nhiệt độ 3D, độ muối, chlorophyll-a, dị thường độ cao mặt biển, tốc độ dòng chảy và động năng rối để xác định phương trình tương quan định lượng và bộ chỉ số thích ứng sinh thái cho 3 loài cá kinh tế: cá Ngừ Vằn, cá Ngừ Chấm và cá Chỉ Vàng.


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

Khung lý thuyết và khái niệm

Luận án áp dụng lý thuyết sinh thái học biển về tương tác cá - môi trường, trong đó sự tập trung đàn cá và năng suất khai thác (CPUE) là hàm phụ thuộc vào các yếu tố vật lý và sinh học biển. Phương pháp xác định bộ chỉ số thích ứng sinh thái ($SI$ - Suitability Index) và mô hình chỉ số thích ứng nơi sống ($HSI$ - Habitat Suitability Index) được sử dụng để lượng hóa mức độ ưa thích môi trường của từng loài cá theo từng mùa gió.

Phương pháp phân tích và tính toán các yếu tố hải dương học

  • Tính toán cấu trúc nhiệt:
    • Dị thường nhiệt bề mặt ($Ano$, °C): $Ano = T_0 - T_{tb}$ ($T_0$: nhiệt độ bề mặt, $T_{tb}$: nhiệt độ bề mặt trung bình nhiều năm).
    • Độ dày lớp đồng nhất trên ($H_0$, m): Xác định từ mặt biển đến độ sâu thỏa mãn gradien nhiệt đứng $G_z = |T_{H0} - T_0| / H_0 \le 0,05^\circ\text{C/m}$.
    • Lớp đột biến nhiệt ($H_1$, m): Biên dưới lớp đột biến thỏa mãn $G_z = |T_1 - T_{H0}| / (H_1 - H_0) > 0,05^\circ\text{C/m}$.
    • Độ sâu mặt đẳng nhiệt 20°C ($H_{20}$) và 24°C ($H_{24}$): Nội suy tuyến tính từ nhiệt độ các tầng chuẩn.
    • Gradien nhiệt ngang bề mặt ($Grad0$) và tầng sâu ($Grad50$, $Grad25$): Tính theo vectơ $GradH = \sqrt{(G_{xH})^2 + (G_{yH})^2}$, góc hướng $\tan\alpha = G_{xH} / G_{yH}$.
  • Tính toán động lực biển:
    • Vận tốc dòng chảy địa chuyển bề mặt ($u, v$) tính từ dị thường độ cao mực biển ($z$ hay $Alti$): $$u = -\frac{g}{f}\frac{\partial z}{\partial y}, \quad v = \frac{g}{f}\frac{\partial z}{\partial x}$$ Trong đó $g = 980\text{ cm/s}^2$, $f = 2\Omega\sin\Phi$ là tham số Coriolis ($\Omega = 7,29 \times 10^{-5}\text{ rad/s}$, $\Phi$ là vĩ độ địa lý).
    • Mật độ động năng rối ($EKE$, cm²/s²): $EKE = \frac{1}{2}(u^2 + v^2)$.
  • Phương pháp phân tích hồi quy tuyến tính nhiều biến: Thiết lập phương trình quan hệ giữa biến phụ thuộc (CPUE) và các biến môi trường độc lập; lựa chọn mô hình tối ưu thông qua tiêu chuẩn thông tin Akaike (AIC).
  • Phương pháp đồng bộ hóa dữ liệu: Quy đổi vị trí dữ liệu về tâm ô lưới $0,5^\circ \times 0,5^\circ$ kinh vĩ. Dữ liệu được kiểm tra sai số thô và chuẩn hóa theo quy tắc 3-sigma ($3\sigma$). Đồng bộ theo 3 cấp thời gian: trung bình 2 mùa gió nhiều năm (mùa gió Đông Bắc từ tháng 10 đến tháng 4 năm sau; mùa gió Tây Nam từ tháng 4 đến tháng 9), trung bình 12 tháng, và giá trị tức thời theo ngày (2012–2015).

Lựa chọn biến môi trường nghiên cứu

  • Đối với cá Ngừ Vằn và cá Ngừ Chấm (14 yếu tố tầng 0–150m):
    1. $T_0$: Nhiệt độ bề mặt (°C)
    2. $Ano$: Dị thường nhiệt độ bề 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 (m)
    5. $T_1$: Nhiệt độ biên dưới lớp đột biến (°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. $Grad0$: Gradien ngang nhiệt bề mặt (°C/10km)
    9. $Grad50$: Gradien ngang nhiệt tầng 50m (°C/10km)
    10. $Sal0$: Độ muối tầng mặt (‰)
    11. $Chlo$: Hàm lượng chlorophyll-a tầng mặt (mg/m³)
    12. $Alti$: Dị thường độ cao động lực mặt biển (cm)
    13. $EKE$: Mật độ động năng rối (cm²/s²)
    14. $Spd_cur$: Tốc độ dòng chảy (cm/s)
  • Đối với cá Chỉ Vàng (10 yếu tố tầng 0–50m do 62,3% mẻ lưới phân bố ở dải 30–50m): $T_0$, $Ano$, $H_0$, $Grad0$, $Grad25$, $Sal0$, $Chlo$, $Alti$, $EKE$, $Spd_cur$.

Nguồn dữ liệu sử dụng

Nguồn dữ liệu Giai đoạn Quy mô dữ liệu Thông số chính
Khảo sát thực đo (Viện Nghiên cứu Hải sản) 1978–2013 - 4.325 trạm nhiệt độ (76.729 bản ghi)
- 2.776 trạm độ muối (57.212 bản ghi)
- 464 trạm Chlorophyll-a (22.835 bản ghi)
Trường nhiệt độ 3D, độ muối 3D, Chlorophyll-a theo độ sâu
Viễn thám và mô hình Mercator (Dự án MOVIMAR) 2011–2016 162.569 bản ghi/thông số (độ phân giải 1/12°) Nhiệt độ mặt biển, Chlorophyll-a, Dị thường mực biển (Alti), Dòng chảy mặt
CSDL Nghề cá (VietFish-Base và WCPFC) 1996–2016 - Cá Ngừ Vằn (lưới rê): 1.503 mẻ lưới
- Cá Ngừ Chấm (lưới rê): 910 mẻ lưới
- Cá Chỉ Vàng (lưới kéo đáy): 247 mẻ lưới
Năng suất khai thác CPUE (kg/km lưới cho cá ngừ; kg/giờ cho cá chỉ và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 1 phân tích bức tranh tổng thể về cơ sở khoa học và điều kiện tự nhiên của khu vực nghiên cứu:

  • Đặc điểm tự nhiên vùng biển Đông Nam Bộ: Vùng biển có diện tích thềm lục địa rộng lớn, đường bờ biển kéo dài từ Bình Thuận đến Cà Mau. Địa hình đáy biển vùng 15–50 m trải rộng và bằng phẳng, từ 50 m trở ra có độ dốc lớn; đường đẳng sâu 100 m cách bờ tới 300 hải lý. Trầm tích bề mặt đáy gồm 7 kiểu chính, trong đó trầm tích sét chiếm ưu thế ở vùng biển khơi, trầm tích cát phân bố thành đới từ Cà Ná xuống Phan Thiết - Hàm Tân và Bạc Liêu.
  • Khí tượng - Thủy văn: Chế độ gió mùa chia làm 2 thời kỳ rõ rệt:
    • Mùa gió Đông Bắc (tháng 11 đến tháng 3 năm sau): Tốc độ gió trung bình nhiều năm tại trạm Phú Quý đạt 8,6 m/s (tháng 1), tại Côn Đảo đạt 3,7 m/s (tháng 1). Khí áp đạt cực đại (1016–1020 mb), nhiệt độ không khí thấp nhất vào tháng 1 (24,6–26,0°C).
    • Mùa gió Tây Nam (tháng 5 đến tháng 10): Tốc độ gió trung bình tại trạm Phú Quý đạt 7,2 m/s (tháng 7), tại Côn Đảo đạt 2,5 m/s. Khí áp khoảng 1000 mb, nhiệt độ không khí cao nhất vào tháng 5 (28,3–29,3°C).
  • Tổng quan sinh học 3 đối tượng cá kinh tế:
    • Cá Ngừ Vằn (Katsuwonus pelamis): Loài cá nổi đại dương di cư, kích thước đánh bắt phổ biến 40–60 cm (Việt Nam đạt 47–58 cm, khối lượng 2,4–2,9 kg). Nhiệt độ thích ứng chung từ 14,7–30,0°C, nhiệt độ tập trung cao giai đoạn cá bột là 23–28°C.
    • Cá Ngừ Chấm (Euthynnus affinis): Loài cá ngừ ven bờ, sống trong dải nhiệt độ 18,0–29,0°C, kích thước khai thác phổ biến 24–45 cm (chiếm ưu thế ở mức 36 cm). Trữ lượng tại vùng biển miền Trung và Nam Bộ khoảng 23,9 nghìn tấn (chiếm 2,1% nhóm cá thu ngừ).
    • Cá Chỉ Vàng (Selaroides leptolepis): Loài cá nổi nhỏ ven bờ, sống ở dải độ sâu dưới 50 m nước, kích thước khai thác phổ biến 11–14 cm, thức ăn chính là động vật phù du. Trữ lượng ở vùng biển phía Nam ước tính khoảng 177–351 ngàn tấn.
+-------------------------------------------------------------------------------+
|                      KHÁI QUÁT CÁC ĐỐI TƯỢNG CÁ NGHIÊN CỨU                   |
+-------------------+----------------------+--------------------+---------------+
| Tên loài          | Tên khoa học         | Nhóm sinh thái     | Ngư cụ chính  |
+-------------------+----------------------+--------------------+---------------+
| Cá Ngừ Vằn        | Katsuwonus pelamis   | Cá nổi đại dương   | Lưới rê       |
| Cá Ngừ Chấm       | Euthynnus affinis    | Cá ngừ nhỏ ven bờ  | Lưới rê       |
| Cá Chỉ Vàng       | Selaroides leptolepis| Cá nổi nhỏ ven bờ  | Lưới kéo đáy  |
+-------------------+----------------------+--------------------+---------------+

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 trình bày chi tiết về quy trình kỹ thuật thu thập, xử lý và chuẩn hóa các nguồn cơ sở dữ liệu:

  • Thống kê dữ liệu nghề cá:
    • Cá Ngừ Vằn: Tổng cộng 1.503 mẫu số liệu; trong đó 412 bản ghi từ sổ nhật ký, 563 bản ghi từ khảo sát và 528 bản ghi từ giám sát. Năng suất trung bình các tháng dao động từ 10,83 kg/km lưới (tháng 2) đến 27,48 kg/km lưới (tháng 9).
    • Cá Ngừ Chấm: Tổng cộng 910 mẫu số liệu; CPUE trung bình tháng cao nhất vào tháng 9 (27,4 kg/km lưới) và tháng 10 (22,6 kg/km lưới), thấp nhất vào tháng 4 (6,7 kg/km lưới).
    • Cá Chỉ Vàng: Tổng cộng 247 mẫu số liệu; CPUE trung bình tháng cao nhất vào tháng 12 (9,1 kg/giờ) và tháng 2 (6,7 kg/giờ), thấp nhất vào tháng 5 và tháng 8 (0,5 kg/giờ).
  • Quy trình đồng bộ dữ liệu: Thiết lập ma trận dữ liệu đồng bộ cá - môi trường trên từng ô lưới $0,5^\circ \times 0,5^\circ$. Mẫu dữ liệu đồng bộ bao gồm tọa độ tâm ô, ngày tháng năm, giá trị CPUE và 14 biến môi trường (hoặc 10 biến đối với cá Chỉ Vàng).
  • Phương pháp xây dựng đường cong SI: Tính toán tần suất xuất hiện các mức CPUE cao ứng với từng khoảng giá trị của từng yếu tố môi trường, chuẩn hóa chỉ số $SI$ về đoạn $[0, 1]$:
    • $SI = 1$: Khoảng giá trị môi trường tối ưu (năng suất khai thác cao nhất).
    • $SI = 0$: Khoảng giá trị môi trường không thích ứng.

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 phân tích quy luật phân bố không gian và biến trình thời gian của các trường hải dương học:

  • Cấu trúc nhiệt biển:
    • Nhiệt độ tầng mặt (SST): Biến trình năm cho thấy nhiệt độ thấp nhất vào tháng 1 (25,0–27,0°C ở ven bờ, 26,5–28,0°C ở ngoài khơi) và cao nhất vào tháng 5 (29,5–30,5°C). Trong mùa hè (tháng 7–8), xuất hiện vùng nhiệt độ thấp (26,5–27,5°C) ở ven bờ Ninh Thuận - Bình Thuận do ảnh hưởng của hiện tượng nước trồi.
    • Cấu trúc thẳng đứng: Độ dày lớp đồng nhất ($H_0$) trong mùa Đông Bắc đạt giá trị lớn (40–60 m ở vùng khơi do đối lưu và xáo trộn gió mạnh); trong mùa Tây Nam, $H_0$ mỏng hơn (10–25 m). Độ sâu biên dưới lớp đột biến ($H_1$) dao động từ 60–120 m.
  • Độ muối ($Sal0$):
    • Dao động trung bình từ 31,0–34,0‰ trong mùa Đông Bắc và 29,0–33,0‰ trong mùa Tây Nam.
    • Vào mùa mưa (tháng 7 đến tháng 10), vùng biển ven bờ cửa sông Mê Kông có độ muối giảm thấp (<30,0‰) do lượng nước ngọt đổ ra lớn, tạo nên đới front độ muối rõ rệt giữa khối nước ven bờ và nước biển khơi.
  • Dòng chảy và xoáy động thái:
    • Mùa gió Đông Bắc: Dòng chảy hướng Bắc - Nam và Đông Bắc - Tây Nam chiếm ưu thế với tốc độ 20–30 cm/s ở ngoài khơi và 15–20 cm/s ở vùng ven bờ; cấu trúc dòng chảy ổn định theo độ sâu ở vùng thềm lục địa <30 m.
    • Mùa gió Tây Nam: Dòng chảy đổi hướng 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 khu vực ranh giới dòng ven bờ và các hoàn lưu xoáy ngoài khơi Phú Quý và Côn Đảo.
  • Hàm lượng Chlorophyll-a ($Chlo$):
    • Giá trị tầng mặt trung bình năm 2002–2017 dao động trong khoảng 0,15–0,60 mg/m³ ở vùng khơi và >1,5 mg/m³ ở vùng ven bờ cửa sông. 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–50 m trong mùa hè tại các trạm ngoài khơi.

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à phần thực nghiệm trọng tâm, lượng hóa mối quan hệ tương quan và xây dựng mô hình dự báo:

Mối quan hệ và chỉ số thích ứng sinh thái của cá Ngừ Vằn

  • Phân tích ma trận tương quan chỉ ra CPUE cá Ngừ Vằn có tương quan thuận với $Chlo$, $Grad0$, $Grad50$, $H_0$, $EKE$ và tương quan nghịch với $Ano$, $Alti$ trong mùa Đông Bắc.
  • Bộ chỉ số thích ứng sinh thái ($SI$) của cá Ngừ Vằn:
    • Mùa gió Đông Bắc: $T_0$ tối ưu từ 26,0–27,5°C; $H_0$ từ 35–55 m; $H_1$ từ 70–95 m; $H_{20}$ từ 80–110 m; $Sal0$ từ 32,5–33,8‰; $Chlo$ từ 0,25–0,45 mg/m³; $Spd_cur$ từ 20–35 cm/s.
    • Mùa gió Tây Nam: $T_0$ tối ưu từ 28,5–29,8°C; $H_0$ từ 15–30 m; $H_1$ từ 50–75 m; $H_{20}$ từ 60–90 m; $Sal0$ từ 31,5–33,0‰; $Chlo$ từ 0,20–0,40 mg/m³; $Spd_cur$ từ 10–25 cm/s.

Mối quan hệ và chỉ số thích ứng sinh thái của cá Ngừ Chấm

  • Năng suất cá Ngừ Chấm tập trung cao tại các khu vực có dị thường nhiệt âm nhẹ, vùng nước ven bờ có dòng chảy trung bình và lớp đồng nhất mỏng.
  • Bộ chỉ số thích ứng sinh thái ($SI$) của cá Ngừ Chấm:
    • Mùa gió Đông Bắc: $T_0$ tối ưu từ 25,5–27,0°C; $H_0$ từ 25–45 m; $Sal0$ từ 32,0–33,5‰; $Chlo$ từ 0,35–0,65 mg/m³; $EKE$ từ 150–350 cm²/s².
    • Mùa gió Tây Nam: $T_0$ tối ưu từ 28,0–29,5°C; $H_0$ từ 10–25 m; $Sal0$ từ 30,5–32,5‰; $Chlo$ từ 0,30–0,60 mg/m³; $EKE$ từ 100–250 cm²/s².

Mối quan hệ và chỉ số thích ứng sinh thái của cá Chỉ Vàng

  • Cá Chỉ Vàng thể hiện mối tương quan chặt chẽ với các yếu tố vùng nước nông thềm lục địa (<50 m) và gradien nhiệt độ tầng 25 m ($Grad25$).
  • Bộ chỉ số thích ứng sinh thái ($SI$) của cá Chỉ Vàng:
    • Mùa gió Đông Bắc: $T_0$ tối ưu từ 25,0–26,5°C; $H_0$ từ 15–30 m; $Grad25$ từ 0,4–0,8°C/10km; $Sal0$ từ 31,5–33,2‰; $Chlo$ từ 0,45–0,90 mg/m³.
    • Mùa gió Tây Nam: $T_0$ tối ưu từ 27,5–29,0°C; $H_0$ từ 10–20 m; $Grad25$ từ 0,3–0,7°C/10km; $Sal0$ từ 30,0–32,0‰; $Chlo$ từ 0,40–0,80 mg/m³.

Dự báo thử nghiệm ngư trường khai thác cá Ngừ Vằn

  • Luận án đã thiết lập quy trình tính toán và đưa ra các bản đồ dự báo thử nghiệm phân bố ngư trường tiềm năng cho cá Ngừ Vằn tại vùng biển Đông Nam Bộ vào các tháng 01, 07, 10 và 12 năm 2016.
  • Kết quả kiểm chứng độc lập với số liệu điều tra và nhật ký khai thác nghề lưới rê năm 2016 cho thấy vị trí các mẻ lưới có năng suất cao trùng khớp với các vùng dự báo tiềm năng tập trung cá (đạt độ tin cậy theo sai số tuyệt đối trên 65%).

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

Những đóng góp mới của luận án

  1. Về mặt khoa học:
    • Đã xác định được mối quan hệ định lượng giữa năng suất khai thác của 3 loài cá kinh tế (cá Ngừ Vằn, cá Ngừ Chấm và cá Chỉ Vàng) với các yếu tố cấu trúc hải dương học đặc trưng (cấu trúc nhiệt đứng, gradien nhiệt ngang, độ muối, chlorophyll-a, dòng chảy và động năng rối) tại vùng biển Đông Nam Bộ.
    • Đã xây dựng được bộ chỉ số thích ứng sinh thái ($SI$) chi tiết theo 2 mùa gió (Đông Bắc và Tây Nam) cho cá Ngừ Vằn, cá Ngừ Chấm và cá Chỉ Vàng, làm cơ sở toán học thiết lập mô hình $HSI$ dự báo ngư trường.
  2. Về mặt thực tiễn:
    • Đã bổ sung, chuẩn hóa cơ sở dữ liệu hải dương học nghề cá đồng bộ tại vùng biển Đông Nam Bộ trên lưới tọa độ $0,5^\circ \times 0,5^\circ$.
    • Đã xây dựng và thử nghiệm thành công mô hình dự báo ngư trường khai thác cá Ngừ Vằn hạn tháng, mở ra khả năng chuyển giao công nghệ phục vụ dự báo khai thác hải sản nghiệp vụ.

Kiến nghị và giải pháp đề xuất

  • Cần tiếp tục duy trì và mở rộng hệ thống quan trắc, thu thập số liệu nghề cá chuẩn hóa (nhật ký khai thác, hệ thống giám sát tàu cá VMS) để tăng cường độ chính xác cho mô hình.
  • Ứng dụng các sản phẩm dự báo ngư trường vào thực tế sản xuất nhằm giúp ngư dân giảm thiểu chi phí tìm kiếm ngư trường, tăng hiệu quả kinh tế cho các đội tàu khai thác xa bờ.
  • Triển khai mở rộng mô hình dự báo thử nghiệm cho đối tượng cá Ngừ Chấm và cá Chỉ Vàng, tiến tới hoàn thiện hệ thống dự báo ngư trường toàn diện cho các vùng biển Việt Nam.

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

Hạn chế của luận án

  • Nguồn số liệu nghề cá thực tế thu thập từ nhật ký và khảo sát phân bố chưa đồng đều giữa các tháng trong năm (tập trung chủ yếu vào các tháng chính vụ của mùa gió Đông Bắc và Tây Nam, số liệu các tháng chuyển tiếp còn thưa).
  • Nghiên cứu mới tiến hành thử nghiệm dự báo chi tiết và kiểm chứng định lượng trên một đối tượng đại diện là cá Ngừ Vằn trong năm 2016; hai đối tượng còn lại (cá Ngừ Chấm và cá Chỉ Vàng) mới dừng lại ở bước xác định bộ chỉ số thích ứng sinh thái và tương quan toán học.
  • Chuỗi số liệu chlorophyll-a và độ muối đo đạc trực tiếp theo mặt thẳng đứng còn hạn chế so với chuỗi số liệu nhiệt độ.

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

  • Tiếp tục bổ sung số liệu giám sát nghề cá thời gian thực để tinh chỉnh các trọng số trong mô hình $HSI$.
  • Tích hợp thêm các kết quả dự báo khí tượng - hải dương phân giải cao từ các mô hình số trị đại dương nghiệp vụ để nâng cao hạn dự báo (xuống hạn tuần 7–10 ngày) cho cả 3 đối tượng cá kinh tế.

Giá trị tham khảo

Đối tượng sử dụng

  • Nghiên cứu sinh, học viên cao học, sinh viên: Thuộc các chuyên ngành Hải dương học, Khí tượng học, Quản lý Thủy sản, Sinh học biển và Kỹ thuật Khai thác Thủy sản.
  • Giảng viên và cán bộ nghiên cứu: Tại các trường đại học (Đại học Khoa học Tự nhiên - ĐHQGHN, Trường Đại học Nha Trang, Trường Đại học Cần Thơ) và các viện nghiên cứu chuyên ngành (Viện Nghiên cứu Hải sản, Viện Hải dương học, Viện Tài nguyên và Môi trường biển).
  • Cơ quan quản lý và hoạch định chính sách: Tổng cục Thủy sản, Trung tâm Dự báo Ngư trường Khai thác Hải 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).

Nội dung có giá trị tham khảo chuyên sâu

  • Chương 2: Quy trình chuẩn hóa và phương pháp đồng bộ hóa dữ liệu nghề cá (CPUE) với dữ liệu hải dương học viễn thám và thực đo trên ô lưới $0,5^\circ \times 0,5^\circ$; các công thức giải tích xác định đặc trưng lớp đồng nhất ($H_0$), lớp đột biến ($H_1$) và động năng rối ($EKE$).
  • Chương 3: Tập bản đồ phân bố không gian và biến trình thời gian của trường nhiệt độ, độ muối, dòng chảy và chlorophyll-a tại thềm lục địa Đông Nam Bộ.
  • Chương 4: Bảng giá trị ngưỡng của bộ chỉ số thích ứng sinh thái ($SI$) cho cá Ngừ Vằn, cá Ngừ Chấm và cá Chỉ Vàng theo từng mùa gió.

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

1. Luận án lựa chọn đối tượng cá nào để nghiên cứu và đại diện cho các nhóm sinh thái nào?

Luận án lựa chọn 3 loài cá kinh tế đặc trưng tại vùng biển Đông Nam Bộ:

  • Cá Ngừ Vằn (Katsuwonus pelamis): Đại diện cho nhóm cá ngừ đại dương (cá nổi lớn xa bờ).
  • Cá Ngừ Chấm (Euthynnus affinis): Đại diện cho nhóm cá ngừ nhỏ ven bờ.
  • Cá Chỉ Vàng (Selaroides leptolepis): Đại diện cho nhóm cá nổi nhỏ ven bờ phân bố ở thềm lục địa nông.

2. Năng suất khai thác (CPUE) của các loài cá được tính toán bằng công thức nào trong luận án?

Năng suất khai thác được chuẩn hóa từ CSDL VietFish-Base theo 2 công thức:

  • Đối với cá Ngừ Vằn và cá Ngừ Chấm (nghề lưới rê): $CPUE = C / L$ (đơn vị: kg/km lưới), trong đó $C$ là sản lượng riêng của loài trong mẻ lưới (kg), $L$ là chiều dài lưới thu về (km).
  • Đối với cá Chỉ Vàng (nghề lưới kéo đáy): $CPUE = C / t$ (đơn vị: kg/giờ), trong đó $C$ là sản lượng cá Chỉ Vàng trong mẻ kéo (kg), $t$ là thời gian kéo lưới (giờ).

3. 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$) được tính như thế nào?

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

4. Nguồn số liệu hải dương học và viễn thám chính trong luận án gồm những gì?

  • Số liệu thực đo: 4.325 trạm nhiệt độ (76.729 số liệu), 2.776 trạm độ muối (57.212 số liệu), 464 trạm chlorophyll-a (22.835 số liệu) do Viện Nghiên cứu Hải sản thực hiện từ năm 1978 đến 2013.
  • Số liệu viễn thám và mô hình Mercator: Trích xuất từ Dự án MOVIMAR giai đoạn 2011–2016 với 162.569 bản ghi cho mỗi thông số (nhiệt độ bề mặt, chlorophyll-a, dị thường mực biển và dòng chảy).

5. Kết quả dự báo thử nghiệm ngư trường cho cá Ngừ Vằn được thực hiện cho các mốc thời gian nào và đánh giá ra sao?

Luận án thực hiện dự báo thử nghiệm ngư trường cá Ngừ Vằn cho 4 tháng đặc trưng năm 2016: tháng 01, tháng 07, tháng 10 và tháng 12. Kết quả được kiểm chứng độc lập với số liệu điều tra và nhật ký khai thác thực tế năm 2016, cho thấy mức độ phù hợp vị trí ngư trường đạt độ tin cậy trên 65% theo sai số tuyệt đối.


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 trọn vẹn bài toán xác lập cơ sở khoa học về mối quan hệ giữa cấu trúc hải dương học và năng suất khai thác cá kinh tế tại vùng biển Đông Nam Bộ. Bằng việc kết hợp các nguồn dữ liệu khảo sát thực tế nhiều năm, dữ liệu viễn thám hiện đại và cơ sở dữ liệu nghề cá đồng bộ, công trình đã lượng hóa thành công các phương trình tương quan và bộ chỉ số thích ứng sinh thái ($SI$) cho cá Ngừ Vằn, cá Ngừ Chấm và cá Chỉ Vàng theo từng mùa gió. Kết quả thử nghiệm dự báo ngư trường cá Ngừ Vằn năm 2016 đã chứng minh tính khả thi của việc ứng dụng mô hình tương tác cá - môi trường vào công tác dự báo ngư trường nghiệp vụ, phục vụ phát triển nghề cá bền vững tại Việt Nam.