Tổng quan nghiên cứu

Vùng biển Nam Trung Bộ có diện tích khoảng 197.522 km2, trải dài từ vĩ độ 9°05'N đến 14°05'N và từ bờ biển Việt Nam đến kinh độ 111°05'E, với độ sâu đáy biến đổi mạnh mẽ từ vài mét ven bờ lên tới trên 3.000 m ngoài khơi. Đây là khu vực biển có ý nghĩa chiến lược đặc biệt về kinh tế và an ninh quốc phòng, đồng thời sở hữu hiện tượng nước trồi mùa hè độc đáo bậc nhất Biển Đông. Hoạt động của dòng nước trồi mang nguồn dinh dưỡng trù phú từ tầng sâu lên bề mặt, thúc đẩy năng suất sinh học sơ cấp đạt mức 60,7 đến 718 mgC/m3/ngày tại tâm trồi. Tuy nhiên, việc đánh giá định lượng và dự báo toàn diện các yếu tố sinh thái biển bằng phương pháp đo đạc thực địa truyền thống gặp rất nhiều rào cản do chi phí cao và hạn chế về độ bao phủ không gian lẫn thời gian.

Mục tiêu trọng tâm của nghiên cứu là khai thác và ứng dụng thành công hệ thống mô hình đại dương quy mô vùng (ROMS) kết hợp với mô-đun sinh thái 4 thành phần NPZD (Nitrate - Phytoplankton - Zooplankton - Detritus) nhằm mô phỏng chi tiết cấu trúc thủy động lực, trường nhiệt - muối và phân bố không gian ba chiều của các yếu tố sinh học sơ cấp. Phạm vi tính toán của mô hình bao phủ toàn bộ vùng biển Nam Trung Bộ và các khu vực lân cận trong chu kỳ mô phỏng liên tục 6 năm, phản ánh đầy đủ hai mùa gió Đông Bắc và Tây Nam. Nghiên cứu mang lại giá trị khoa học và thực tiễn sâu sắc, cung cấp cơ sở dữ liệu định lượng với hệ số tương quan ổn định đạt mức r = 0,9951, hỗ trợ đắc lực cho công tác quy hoạch khai thác thủy sản, bảo tồn đa dạng sinh học và phát triển kinh tế biển bền vững.

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

Khung lý thuyết áp dụng

Nghiên cứu được xây dựng trên nền tảng cơ học chất lưu địa vật lý và sinh thái học đại dương hiện đại. Khung lý thuyết thủy động lực học dựa trên hệ phương trình Navier-Stokes trung bình Reynolds (RANS) kết hợp với xấp xỉ Boussinesq và xấp xỉ thủy tĩnh. Hệ phương trình được giải trên hệ tọa độ trực giao phẳng ngang Đề-các (Arakawa-C) và tọa độ sigma biến đổi theo chiều sâu thẳng đứng gồm 20 tầng, cho phép mô tả chính xác sự tương tác phức tạp của dòng chảy đáy và mặt thoáng đại dương. Quá trình khép kín rối được thực hiện thông qua việc tham số hóa hệ số nhớt rối và hệ số khuếch tán rối theo phương pháp giải tích.

Về mặt sinh thái, mô hình NPZD do Fasham đề xuất được tích hợp để mô tả chu trình chuyển hóa vật chất giữa 4 biến trạng thái cốt lõi: muối dinh dưỡng Nitrat (N), sinh khối thực vật nổi (P), sinh khối động vật nổi (Z) và chất vẩn hữu cơ (D). Các quá trình sinh học chính bao gồm: quang hợp của thực vật nổi phụ thuộc vào bức xạ quang hợp và nhiệt độ với tỷ số C:N cố định là 6,625 mMol C / mMol N; động vật nổi tiêu thụ thực vật nổi với tốc độ ăn cực đại 0,9 ngày^-1; quá trình bài tiết và tử vong tạo thành chất vẩn; quá trình khoáng hóa chất vẩn tái tạo lại nguồn Nitrat với tốc độ bài tiết 0,10 ngày^-1. Chlorophyll-a được sử dụng làm chỉ thị định lượng then chốt để đánh giá sinh khối và sức sản xuất sơ cấp của thực vật phù du.

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

Phương pháp nghiên cứu chủ đạo là mô hình hóa số trị 3 chiều độ phân giải cao kết hợp phân tích chuỗi số liệu viễn thám và tái phân tích khí tượng hải văn toàn cầu. Cỡ mẫu tính toán được thiết lập trên lưới số trị có độ phân giải không gian 1/4 độ (khoảng 27 km), chia thành 20 tầng sigma theo phương đứng với tham số bề mặt bằng 6,0 nhằm tăng mật độ lưới ở lớp nước mặt 0 - 100 m và tham số đáy bằng 0,0. Lý do lựa chọn mô hình số trị ROMS-NPZD xuất phát từ tính ổn định cao, khả năng mô phỏng chi tiết các quá trình động lực học phi tuyến và sự tương tác phức tạp của hệ sinh thái ven biển mà các trạm quan trắc rời rạc không thể tái hiện được.

Dữ liệu đầu vào phục vụ mô hình được chuẩn hóa từ các nguồn dữ liệu tiêu chuẩn quốc tế:

  • Địa hình đáy biển trích xuất từ bộ dữ liệu ETOPO2 với độ phân giải 2 phút góc và đường bờ biển GSHHS.
  • Trường lực khí quyển bề mặt sử dụng bộ số liệu trung bình tháng COADS05 (độ phân giải 0,5 độ), bao gồm ứng suất gió, bức xạ sóng ngắn mặt trời (50 - 70 cal/m2.giờ), thông lượng nhiệt tịnh và thông lượng trao đổi nước ngọt (E-P).
  • Điều kiện ban đầu và điều kiện biên đại dương lấy từ tập dữ liệu hải văn WOA2009 (độ phân giải 1 độ) cho các trường nhiệt độ, độ muối, oxy hòa tan và nồng độ Nitrat.

Quy trình tính toán số trị được chia thành 2 giai đoạn: Giai đoạn 1 chạy mô hình thủy động lực học thuần túy trong 3 năm đầu (với bước thời gian tích phân dt = 2.400 giây) để trường dòng chảy và nhiệt - muối đạt trạng thái cân bằng động; Giai đoạn 2 kích hoạt mô-đun sinh thái NPZD chạy tiếp tục trong 3 năm tiếp theo để trích xuất trường sinh học ổn định.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Mô hình ROMS-NPZD đã tái hiện thành công bức tranh toàn cảnh về thủy động lực và sinh thái vùng biển Nam Trung Bộ với 3 phát hiện nổi bật:

Thứ nhất, hoàn lưu dòng chảy có sự phân hóa mùa rõ rệt. Trong mùa đông (tháng 1), dòng chảy mặt chịu sự chi phối của gió mùa Đông Bắc, chảy dọc bờ theo hướng Tây Nam với vận tốc trung bình từ 0,1 đến 0,4 m/s. Ngược lại, vào mùa hè (tháng 7), dòng chảy đổi hướng sang Đông Bắc với vận tốc tăng mạnh lên 0,5 đến trên 0,8 m/s ở dải ven bờ. Hiện tượng nước trồi xuất hiện mãnh liệt tại tọa độ 10,5°N - 11,5°N với vận tốc dòng dâng thẳng đứng vượt mức 8 x 10^-5 m/s, đưa khối nước lạnh và giàu dinh dưỡng từ độ sâu 75 - 200 m lên sát bề mặt.

Thứ hai, trường nhiệt - muối phản ánh chân thực ranh giới tâm nước trồi. Vào tháng 7, nhiệt độ bề mặt biển tại tâm trồi giảm mạnh xuống còn 25,6 - 27,0°C, thấp hơn từ 1,5 đến 2,0°C so với các vùng biển xung quanh (thường đạt trên 29,0°C). Đồng thời, độ muối tầng mặt tại tâm trồi tăng cao đạt 33,8 đến trên 34,4‰, tạo nên cấu trúc phân tầng thủy văn dị thường với các đường đẳng nhiệt và đẳng muối bị uốn cong đẩy ngược lên phía mặt thoáng.

Thứ ba, sự bùng nổ sinh thái và năng suất sinh học sơ cấp tập trung cao độ vào mùa hè. Nồng độ Chlorophyll-a tầng mặt tháng 7 tại dải nước trồi ven bờ đạt đỉnh từ 3,0 đến 3,5 mg/m3, cao hơn gấp 5 đến 7 lần so với mùa đông (chỉ đạt 0,1 - 0,6 mg/m3). Đi cùng với đó, sinh khối thực vật nổi tại tâm trồi tăng vọt lên mức 1,2 - 1,6 mmol N/m3 và sinh khối động vật nổi đạt 0,8 - 1,2 mmol N/m3, phân bố tập trung chủ yếu ở lớp nước chiếu sáng từ 0 đến 30 m.

Thảo luận kết quả

Cơ chế bùng phát năng suất sinh học tại vùng biển Nam Trung Bộ được giải thích bằng sự kết hợp giữa động lực gió mùa và quá trình cấp dưỡng chất. Khi gió Tây Nam thổi mạnh (với tần suất tốc độ gió trên 6 m/s chiếm hơn 60%), hiệu ứng phân kỳ dòng chảy tầng mặt kéo theo dòng nước giàu muối dinh dưỡng Nitrat từ đáy biển (nồng độ tầng đáy đạt trên 35 mmol N/m3) dâng lên tầng quang hợp. Dưới điều kiện bức xạ mặt trời dồi dào (đạt 130 - 180 Kcal/cm2.năm) và chiều dày lớp quang hợp đạt từ 80 đến 100 m, thực vật nổi hấp thụ Nitrat để quang hợp mạnh mẽ, tạo ra nguồn thức ăn dồi dào thúc đẩy sự gia tăng sinh khối của động vật nổi ngay sau đó.

Toàn bộ chuỗi dữ liệu mô phỏng được thể hiện trực quan qua các biểu đồ phân bố không gian 2D của trường vận tốc dòng chảy, bản đồ nhiệt độ bề mặt biển và các đồ thị mặt cắt thẳng đứng theo kinh tuyến 111°E cũng như mặt cắt ngang vuông góc với đường bờ. Các bảng tổng hợp chỉ số thống kê chỉ ra rằng lớp cực đại Chlorophyll-a ngầm (Subsurface Chlorophyll Maximum - SCM) thường tồn tại ổn định ở độ sâu 30 - 50 m vào mùa đông và dâng lên sát mặt ở độ sâu 0 - 20 m tại tâm nước trồi mùa hè.

So sánh với các nghiên cứu thực địa trước đây, giá trị tính toán từ mô hình có độ tương đồng rất cao với các công bố khoa học của những chuyên gia hải dương đầu ngành (nơi ghi nhận năng suất sơ cấp đạt 60,7 - 718 mgC/m3/ngày tại tâm trồi và hàm lượng phốt pho vô cơ đạt trên 8 mgP/m3). Điều này khẳng định mô hình ROMS-NPZD có độ tin cậy vượt trội trong việc lượng hóa chu trình sinh học đại dương.

Đề xuất và khuyến nghị

Dựa trên kết quả mô phỏng động lực và sinh thái học, nghiên cứu đưa ra 4 nhóm khuyến nghị hành động cụ thể:

  1. Thiết lập hệ thống dự báo ngư trường khai thác hải sản hạn ngắn: Ứng dụng mô hình ROMS-NPZD để phát bản tin dự báo vị trí tâm nước trồi và dải sinh khối Chlorophyll-a cao (trên 2,0 mg/m3) trước 7 đến 10 ngày. Mục tiêu nâng cao hiệu quả đánh bắt cá nổi nhỏ thêm 20 - 25% trong giai đoạn 2026 - 2028. Cơ quan chủ trì thực hiện: Cục Thủy sản phối hợp cùng Viện Hải dương học.

  2. Đầu tư mạng lưới phao biển quan trắc tự động chuyên sâu: Lắp đặt bổ sung 5 trạm phao hải văn - hóa sinh tự động đo đạc liên tục các thông số nhiệt - muối, oxy hòa tan và dòng chảy thẳng đứng tại vùng biển Ninh Thuận - Bình Thuận trong giai đoạn 2026 - 2027. Mục tiêu cung cấp chuỗi số liệu hiệu chỉnh theo thời gian thực nhằm giảm sai số mô hình xuống dưới 5%. Cơ quan thực hiện: Trung tâm Dự báo Khí tượng Thủy văn Quốc gia.

  3. Nâng cấp mô-đun sinh thái đa loài bậc cao: Mở rộng mô-đun NPZD hiện tại sang mô hình sinh thái đa kích thước (kết hợp thực vật nổi kích thước nano, micro và hai nhóm động vật nổi lớn - nhỏ) trước năm 2029 nhằm mô phỏng chính xác đường truyền thức ăn lên các loài hải sản thương mại. Cơ quan thực hiện: Trường Đại học Khoa học Tự nhiên - Đại học Quốc gia Hà Nội cùng các viện nghiên cứu chuyên ngành.

  4. Ban hành quy chế bảo tồn vùng sinh thái nước trồi trọng điểm: Áp dụng quy định cấm khai thác có thời hạn từ tháng 5 đến tháng 8 hàng năm tại vùng lõi nước trồi (vĩ độ 10,5°N đến 12,0°N), nhằm bảo vệ bãi đẻ và giai đoạn ấu trùng của các loài hải sản, hướng tới mục tiêu giảm 15% tốc độ suy thoái nguồn lợi thủy sản ven bờ vào năm 2030. Cơ quan thực hiện: Bộ Nông nghiệp và Phát triển Nông thôn cùng Ủy ban nhân dân các tỉnh Nam Trung Bộ.

Đối tượng nên tham khảo luận văn

Nội dung và phương pháp tiếp cận của luận văn mang lại giá trị thiết thực cho 4 nhóm đối tượng chính:

  1. Học viên cao học, nghiên cứu sinh và nhà khoa học ngành Hải dương học, Thủy văn, Biến đổi khí hậu: Nắm vững phương pháp thiết lập lưới tính sigma, kỹ thuật ghép nối mô-đun sinh thái NPZD vào mô hình ROMS và quy trình tiền - hậu xử lý dữ liệu NetCDF qua bộ công cụ ROMSTOOLS.

  2. Cơ quan quản lý nhà nước về Biển, Hải đảo và Thủy sản: Sử dụng các bản đồ phân vùng sinh thái và trữ lượng sinh vật sơ cấp làm luận cứ khoa học để quy hoạch không gian biển, xây dựng các khu bảo tồn biển và hoạch định chiến lược phát triển kinh tế biển xanh.

  3. Kỹ sư mô hình hóa môi trường và tư vấn đánh giá tác động môi trường biển: Tham khảo cấu trúc phương trình và bộ tham số hóa sinh học để ứng dụng mô phỏng sức tải môi trường, tính toán khả năng tự làm sạch của vực nước và dự báo hiện tượng phú dưỡng hay thủy triều đỏ.

  4. Doanh nghiệp khai thác thủy sản xa bờ và dịch vụ hậu cần nghề cá: Tận dụng quy luật biến động dòng chảy và vùng nước trồi để tối ưu hóa tuyến hành trình tàu cá, giảm 15 - 20% chi phí nhiên liệu tìm kiếm ngư trường và nâng cao sản lượng khai thác.

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

Mô hình ROMS kết hợp NPZD có ưu thế gì so với các mô hình thủy động lực học khác? Mô hình ROMS sử dụng hệ tọa độ sigma bám sát địa hình đáy kết hợp lưới ngang Arakawa-C, giúp giải quyết hoàn hảo bài toán dòng chảy ven bờ dốc đứng tại Nam Trung Bộ. Việc ghép nối trực tiếp mô-đun sinh thái NPZD cho phép mô phỏng tương tác 2 chiều giữa hoàn lưu đại dương và chu trình chuyển hóa sinh học liên tục trên diện tích gần 200.000 km2 mà các phép đo đạc đơn lẻ không thể bao quát được.

Tại sao năng suất sinh học sơ cấp tại Nam Trung Bộ lại tăng đột biến vào tháng 7? Vào tháng 7, gió mùa Tây Nam hoạt động mạnh với tốc độ gió thường xuyên trên 6 m/s, tạo lực phân kỳ Ekman ven bờ và kích hoạt hiện tượng nước trồi mạnh ở độ sâu 75 - 200 m. Dòng dâng thẳng đứng với vận tốc hơn 8 x 10^-5 m/s đã bơm một lượng lớn Nitrat (trên 35 mmol N/m3) lên tầng nước mặt được chiếu sáng, kích thích thực vật nổi quang hợp và bùng nổ sinh khối.

Lớp cực đại Chlorophyll-a ngầm (SCM) được mô hình ghi nhận ở độ sâu bao nhiêu? Kết quả tính toán chỉ ra rằng lớp cực đại Chlorophyll-a ngầm thường phân bố ở độ sâu 30 - 50 m vào mùa đông và dâng lên sát mặt ở độ sâu 0 - 20 m tại tâm trồi mùa hè. Độ sâu này là nơi giao thoa lý tưởng giữa cường độ bức xạ quang hợp tối ưu (chiều dày lớp quang hợp đạt 80 - 100 m) và nguồn dinh dưỡng dồi dào khuếch tán từ dưới đáy lên.

Những bộ số liệu chuẩn quốc tế nào đã được tích hợp làm điều kiện biên cho mô hình? Nghiên cứu đã tích hợp 4 bộ dữ liệu toàn cầu có độ tin cậy cao: địa hình đáy ETOPO2 (độ phân giải 2 phút), đường bờ biển GSHHS, trường lực khí tượng bề mặt COADS05 (độ phân giải 0,5 độ) và trường dữ liệu hải văn - hóa sinh WOA2009 (độ phân giải 1 độ), bảo đảm tính chuẩn xác và đồng bộ cho toàn miền tính.

Độ tin cậy và tính ổn định của mô hình số trị trong luận văn được đánh giá ra sao? Mô hình được chạy liên tục trong 6 năm với bước thời gian tích phân 2.400 giây. Hệ số tương quan nhiệt độ bề mặt biển giữa các năm tính toán tăng dần và đạt mức r = 0,9951 ở năm thứ 6. Điều này minh chứng mô hình đã đạt trạng thái cân bằng động lực học rất cao, hạn chế tối đa sai số tích lũy và hoàn toàn đủ độ tin cậy để phục vụ các phân tích khoa học.

Kết luận

  • Ứng dụng thành công hệ thống mô hình số trị đại dương ba chiều ROMS kết hợp mô-đun sinh thái NPZD trên lưới tính phân giải 1/4 độ và 20 tầng sigma cho toàn vùng biển Nam Trung Bộ.
  • Làm sáng tỏ cấu trúc hoàn lưu dòng chảy và bản chất cơ chế của hiện tượng nước trồi mùa hè với vận tốc dòng dâng thẳng đứng cực đại vượt 8 x 10^-5 m/s tại tọa độ 10,5°N - 11,5°N.
  • Định lượng chính xác bức tranh phân bố nhiệt - muối và sự bùng nổ của các yếu tố sinh học sơ cấp, trong đó nồng độ Chlorophyll-a đạt đỉnh 3,5 mg/m3 và sinh khối thực vật nổi đạt 1,6 mmol N/m3 tại vùng nước trồi.
  • Xác lập độ tin cậy vượt bậc của mô hình số trị thông qua hệ số tương quan ổn định qua các năm đạt r = 0,9951, tương thích chặt chẽ với các chuỗi số liệu thực nghiệm lịch sử.
  • Cung cấp cơ sở khoa học then chốt phục vụ công tác quy hoạch không gian biển, dự báo ngư trường khai thác hải sản và ứng phó với biến đổi khí hậu giai đoạn 2026 - 2030.

Để phát huy tối đa giá trị thực tiễn của công trình, các nhóm nghiên cứu và cơ quan quản lý biển có thể tiếp tục kết nối, mở rộng mô hình với các mô-đun sinh thái bậc cao nhằm phục vụ đắc lực cho chiến lược phát triển bền vững kinh tế biển Việt Nam.