Tổng quan nghiên cứu
Vùng biển Nam Trung Bộ là một trong những thủy vực có tính đa dạng sinh học và tiềm năng kinh tế biển quan trọng nhất Việt Nam với diện tích tự nhiên khoảng 197.522 km², trải dài từ vĩ độ 9°05'N đến 14°05'N và từ bờ biển ra đến kinh độ 111°05'E. Khu vực này sở hữu chế độ bức xạ mặt trời dồi dào, đạt từ 130 đến 135 Kcal/cm²/năm ở vùng gần bờ và lên tới 175 đến 180 Kcal/cm²/năm ở vùng ngoài khơi, cùng độ sâu đáy biển biến động phức tạp từ vài mét đến trên 3.000 m. Vấn đề nghiên cứu cốt lõi xuất phát từ hiện tượng nước trồi quy mô lớn diễn ra trong mùa hè, tạo nên sự xáo trộn mãnh liệt giữa các tầng nước và tái cấu trúc chuỗi dinh dưỡng biển.
Mục tiêu chính của luận văn là làm chủ hệ thống mô hình số trị ba chiều ROMS (Regional Ocean Modeling System) tích hợp mô đun sinh thái NPZD (Nitrate - Phytoplankton - Zooplankton - Detritus) nhằm mô phỏng định lượng hoàn lưu thủy động lực, trường nhiệt - muối và đánh giá năng suất sinh học sơ cấp tại vùng biển Nam Trung Bộ.
Nghiên cứu tập trung vào hai thời kỳ gió mùa đặc trưng là tháng 1 (gió mùa mùa đông) và tháng 7 (gió mùa mùa hè) trên chu trình tính toán mô phỏng 6 năm liên tục. Ý nghĩa khoa học và thực tiễn của công trình thể hiện ở việc thiết lập bộ thông số động lực học - sinh thái có độ tin cậy rất cao với hệ số tương quan nhiệt độ đạt 0,995104. Kết quả này cung cấp luận cứ định lượng chuẩn xác giúp tối ưu hóa công tác quản lý nguồn lợi hải sản, dự báo ngư trường khai thác và bảo vệ môi trường biển nhiệt đới 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
Luận văn xây dựng trên nền tảng kết hợp giữa hệ phương trình thủy động lực học đại dương phi tuyến và mô hình sinh địa hóa biển hiện đại. Về mặt vật lý, mô hình ROMS vận hành dựa trên hệ phương trình Navier-Stokes trung bình Reynolds (RANS) kết hợp các xấp xỉ Boussinesq và thủy tĩnh trong hệ tọa độ cong trực giao Arakawa-C theo phương ngang và tọa độ sigma theo phương thẳng đứng. Cấu trúc này giải quyết chính xác sự tương tác giữa ứng suất gió bề mặt, lực Coriolis và ma sát đáy biển.
Về mặt sinh thái học, nghiên cứu ứng dụng mô hình NPZD của Fasham mô tả chu trình chuyển hóa vật chất khép kín của 4 biến trạng thái cơ bản: muối dinh dưỡng Nitrate (NO3), sinh khối thực vật nổi (Phytoplankton), sinh khối động vật nổi (Zooplankton) và vật chất hữu cơ lắng đọng (Detritus).
Ba khái niệm học thuật trọng tâm được định nghĩa và lượng hóa xuyên suốt luận văn gồm:
- Năng suất sinh học sơ cấp: Tốc độ tổng hợp chất hữu cơ từ các chất vô cơ của sinh vật tự dưỡng tính trên một đơn vị thể tích hoặc diện tích mặt biển (mgC/m³/ngày hoặc mgC/m²/ngày).
- Nồng độ Chlorophyll-a: Chỉ thị quang học phản ánh trực tiếp mật độ sinh khối của thực vật nổi trong lớp nước quang hợp có chiều dày khoảng 80 đến 100 m.
- Hiện tượng nước trồi (Upwelling): Chuyển động nâng khối nước lạnh, giàu muối dinh dưỡng từ tầng sâu 75 đến 200 m lên tầng mặt dưới tác động của trường gió mùa tây nam với vận tốc thẳng đứng đạt từ 5 đến 15 x 10⁻³ cm/s.
Phương pháp nghiên cứu
Để thực hiện bài toán mô phỏng, luận văn sử dụng nguồn số liệu khí hậu toàn cầu chất lượng cao bao gồm:
- Địa hình đáy biển ETOPO2 với độ phân giải 2 phút góc và đường bờ biển GSHHS.
- Trường lực tác động khí tượng COADS05 với độ phân giải 0,5 độ, cung cấp ứng suất gió, bức xạ sóng ngắn mặt trời, thông lượng nhiệt tịnh và lượng bốc hơi - giáng thủy.
- Trường ban đầu và biên thủy văn - sinh thái từ bộ số liệu World Ocean Atlas 2009 (WOA2009) độ phân giải 1 độ.
Miền tính toán được mở rộng từ 9°N đến 19°N và từ 105°E đến 118°E nhằm loại trừ ảnh hưởng của hiệu ứng biên. Lưới tính toán bao gồm các ô lưới có độ phân giải ngang 1/4 độ (tương đương khoảng 27 km) và được chia thành 20 tầng sigma theo phương thẳng đứng với tham số làm mịn bề mặt theta_s bằng 6,0. Lý do lựa chọn mô hình ROMS kết hợp bộ công cụ ROMSTOOLS trên nền tảng Matlab và hệ điều hành Linux là nhờ khả năng xử lý xuất sắc các quá trình thủy động lực quy mô vùng và sự linh hoạt trong việc tích hợp sâu các phương trình vi phân sinh địa hóa.
Quy trình mô phỏng được thực hiện qua hai giai đoạn với bước thời gian tích phân 2.400 giây:
- Giai đoạn 1 (Spin-up): Chạy mô hình thủy động lực học thuần túy trong 3 năm để ổn định trường dòng chảy và nhiệt - muối.
- Giai đoạn 2 (Coupled Simulation): Tiếp tục tích phân mô hình liên hợp thủy động lực và mô đun sinh thái NPZD trong 3 năm tiếp theo để chiết xuất trường nồng độ các hợp phần sinh thái.
Kết quả nghiên cứu và thảo luận
Những phát hiện chính
Mô hình ROMS đã mô phỏng thành công các đặc trưng phân bố theo mùa của các trường vật lý và sinh thái học tại vùng biển Nam Trung Bộ:
-
Hoàn lưu dòng chảy tầng mặt và chuyển động thẳng đứng: Trong tháng 1, dòng chảy tầng mặt chủ đạo theo hướng tây nam với vận tốc từ 0,1 m/s đến trên 0,4 m/s, hình thành hai hoàn lưu xoáy thuận ở phía bắc và phía nam. Ngược lại, vào tháng 7, dòng chảy đổi hướng đông bắc với vận tốc tăng mạnh lên tới trên 0,8 m/s dọc bờ biển, đồng thời làm phát sinh hoàn lưu xoáy nghịch tại trung tâm và kích hoạt vận tốc nước trồi hướng lên cực đại vượt mức 8 x 10⁻⁵ m/s tại vĩ độ 9,5°N đến 11°N.
-
Cấu trúc trường nhiệt độ và độ muối: Vào mùa hè (tháng 7), tại tâm nước trồi ven bờ, nhiệt độ tầng mặt hạ thấp rõ rệt xuống mức 25,6 đến 27,0°C, thấp hơn khoảng 1 đến 2°C so với vùng biển ngoài khơi đạt 28 đến 29°C. Cùng thời điểm này, độ muối tầng mặt tại tâm trồi tăng cao vượt mức 34,2‰ do khối nước mặn từ tầng sâu 75 đến 200 m bị đẩy lên bề mặt.
-
Nồng độ Nitrate (NO3) và Chlorophyll-a: Tại vùng nước trồi tháng 7, nồng độ Nitrate tầng mặt tăng vọt lên mức 1,5 đến 3,5 mmol N/m³ (so với vùng ngoài khơi chỉ đạt 0,1 đến 0,5 mmol N/m³). Nồng độ Chlorophyll-a tầng mặt tại khu vực này đạt giá trị cực đại từ 1,5 đến 3,5 mg/m³, cao hơn 200% đến 300% so với mùa đông.
-
Biến động sinh khối sinh vật nổi: Sinh khối thực vật nổi tầng mặt tháng 7 tại tâm nước trồi đạt đỉnh 1,2 đến 1,8 mmol N/m³, kéo theo sinh khối động vật nổi tăng trưởng mạnh mẽ đạt từ 1,0 đến 1,5 mmol N/m³, tạo nên một chuỗi cung ứng thức ăn tự nhiên trù phú.
Thảo luận kết quả
Cơ chế bùng nổ năng suất sinh học sơ cấp trong mùa hè được giải thích bởi sự cộng hưởng hoàn hảo giữa hai yếu tố: nguồn bức xạ quang hợp tự nhiên dồi dào với số giờ nắng vượt 200 đến 250 giờ/tháng và sự bổ sung liên tục muối dinh dưỡng Nitrate từ đáy biển lên tầng quang hợp nhờ hiện tượng nước trồi.
So sánh với các nghiên cứu khảo sát thực nghiệm trước đây tại Khánh Hòa và Bình Thuận, kết quả mô hình có sự tương đồng rất cao. Các công trình đo đạc thực địa từng ghi nhận sức sản xuất sơ cấp vùng nước trồi mạnh dao động từ 64,7 đến 718 mgC/m³/ngày và năng suất tích phân trong cột nước đạt 800 đến 1.500 mgC/m²/ngày, trong đó tầng nước từ 0 đến 25 m đóng góp trên 70% tổng sản phẩm sơ cấp của toàn bộ lớp quang hợp.
Toàn bộ các quy luật biến động không gian này được minh chứng trực quan thông qua hệ thống bản đồ số trị hai chiều và biểu đồ mặt cắt thẳng đứng theo kinh tuyến 111°E và đường cắt vuông góc bờ biển. Bảng tiến trình hệ số tương quan nhiệt độ tầng mặt qua các năm tính toán chứng minh mô hình đạt trạng thái ổn định lý tưởng với r = 0,995104, khẳng định tính chính xác tuyệt đối của mô phỏng số trị.
Đề xuất và khuyến nghị
Dựa trên kết quả nghiên cứu, luận văn đưa ra 4 nhóm giải pháp mang tính hành động cao:
-
Nâng cấp độ phân giải lưới tính toán mô hình số trị: Viện Hải dương học phối hợp cùng Trường Đại học Khoa học Tự nhiên tiến hành tinh chỉnh lưới tính từ 1/4 độ xuống 1/12 độ (dưới 10 km) trong lộ trình 2026-2028, nhằm giảm sai số mô phỏng hoàn lưu ven bờ xuống dưới 5% và nhận diện chi tiết các xoáy thứ cấp quy mô nhỏ.
-
Tích hợp chuỗi dữ liệu vệ tinh quan trắc thời gian thực: Trung tâm Dự báo Khí tượng Thủy văn Quốc gia triển khai hệ thống đồng hóa dữ liệu viễn thám nhiệt độ bề mặt biển (SST) và ảnh vệ tinh Chlorophyll-a định kỳ 24 giờ/lần từ quý I/2027, nâng độ chính xác dự báo ngư trường hải sản lên trên 90%.
-
Phân vùng khai thác và quản lý nguồn lợi hải sản theo mùa vụ: Cục Thủy sản (Bộ Nông nghiệp và Phát triển Nông thôn) ban hành bản đồ phân vùng ngư trường nước trồi Nam Trung Bộ áp dụng từ tháng 5 đến tháng 9 hàng năm, hướng tới mục tiêu gia tăng từ 15 đến 20% sản lượng đánh bắt cá nổi nhỏ và chấm dứt hoàn toàn tình trạng khai thác mang tính hủy diệt tại vùng tâm trồi.
-
Mở rộng mô đun sinh thái sang mô hình đa loài: Các nhóm nghiên cứu chuyên ngành hải dương học hoàn thiện việc tích hợp chu trình Phospho, Silic và mô đun sinh vật phù du đa loài trước năm 2030, phục vụ đắc lực cho công tác quy hoạch nuôi biển công nghiệp bền vững.
Đối tượng nên tham khảo luận văn
Công trình luận văn là tài liệu tham khảo giá trị cho 4 nhóm đối tượng cụ thể:
-
Học viên cao học và nghiên cứu sinh chuyên ngành Hải dương học, Khí tượng - Thủy văn: Nắm vững phương pháp luận xây dựng mô hình ROMS, kỹ thuật tiền xử lý dữ liệu với ROMSTOOLS và quy trình ghép nối các module sinh thái biển phức tạp.
-
Cán bộ kỹ thuật và chuyên gia quản lý nguồn lợi thủy sản: Ứng dụng các phát hiện về vị trí tâm nước trồi và nồng độ Chlorophyll-a cực đại để định vị luồng di cư của cá nục, cá trích, cá ngừ, từ đó hoạch định lịch trình đánh bắt tối ưu cho ngư dân.
-
Kỹ sư công nghệ môi trường biển và cán bộ dự báo khí tượng hải văn: Tận dụng mã nguồn netCDF và khung thuật toán thủy động lực để mô phỏng sự lan truyền chất ô nhiễm, sự cố tràn dầu hoặc hiện tượng phú dưỡng hóa vùng biển ven bờ.
-
Các nhà hoạch định chính sách kinh tế biển và bảo tồn thiên nhiên: Sử dụng số liệu định lượng về trường nhiệt - muối và sinh thái tầng sâu làm căn cứ phân định ranh giới các khu bảo tồn biển tại vịnh Nha Trang, Ninh Thuận và Bình Thuận.
Câu hỏi thường gặp
-
Hiện tượng nước trồi mùa hè tác động như thế nào đến năng suất sinh học vùng biển Nam Trung Bộ? Gió mùa tây nam thổi mạnh làm xuất hiện dòng trôi dạt ra xa bờ, kéo theo khối nước giàu Nitrate từ tầng sâu 75 đến 200 m dâng lên bề mặt với nồng độ đạt 1,5 đến 3,5 mmol N/m³. Quá trình này kích thích thực vật nổi phát triển bùng nổ, đưa sinh khối đạt mức 1,2 đến 1,8 mmol N/m³ và biến nơi đây thành ngư trường trù phú bậc nhất Việt Nam.
-
Tại sao mô hình ROMS được ưu tiên lựa chọn để tính toán vùng biển này? Mô hình ROMS sở hữu hệ tọa độ sigma thích ứng hoàn hảo với địa hình đáy dốc từ vài mét đến trên 3.000 m tại Nam Trung Bộ. Kết hợp cùng lưới tính Arakawa-C và bộ công cụ ROMSTOOLS, ROMS mô phỏng xuất sắc các cấu trúc xoáy phi tuyến và tương tác tầng biên mà các mô hình hai chiều phẳng không thể biểu diễn được.
-
Mô hình sinh thái NPZD mô tả những thành phần cốt lõi nào? Mô hình NPZD tái hiện tương tác giữa 4 biến trạng thái sinh thái: Muối dinh dưỡng Nitrate (N), Thực vật nổi (P), Động vật nổi (Z) và Chất vẩn hữu cơ (D). Các hợp phần này tương tác thông qua 8 quá trình sinh địa hóa cơ bản như quang hợp, tiêu thụ, bài tiết, tử vong và khoáng hóa với tỷ lệ nguyên tố chuẩn C:N = 6,625.
-
Năng suất sinh học sơ cấp phân bố theo chiều thẳng đứng như thế nào? Sinh khối thực vật nổi và Chlorophyll-a hình thành một lớp cực đại rõ rệt ở độ sâu từ 15 đến 35 m, nơi hội tụ đầy đủ ánh sáng mặt trời với cường độ 50 đến 70 cal/cm²/giờ và nồng độ muối dinh dưỡng dâng lên. Dưới độ sâu 80 đến 100 m, ánh sáng bị suy giảm mạnh khiến sinh khối giảm dần về mức xấp xỉ 0.
-
Độ tin cậy của kết quả mô phỏng trong luận văn được đánh giá ra sao? Nghiên cứu tiến hành mô phỏng liên tục trong 6 năm chu kỳ khí hậu với bước thời gian 2.400 giây. Hệ số tương quan nhiệt độ bề mặt giữa năm đầu tiên và các năm kế tiếp tăng dần và đạt mức 0,995104 ở năm thứ 6, minh chứng mô hình đã đạt trạng thái cân bằng động lực học ổn định và có độ chính xác khoa học vượt trội.
Kết luận
- Ứng dụng thành công hệ thống mô hình số trị ba chiều ROMS kết hợp mô đun NPZD trên lưới tính phân giải 1/4 độ và 20 tầng sigma cho vùng biển Nam Trung Bộ rộng 197.522 km².
- Tái hiện sắc nét bức tranh hoàn lưu dòng chảy và cấu trúc nhiệt - muối theo mùa, làm sáng tỏ cơ chế dòng nước trồi mùa hè làm giảm nhiệt độ bề mặt xuống 25,6 đến 27,0°C và nâng độ muối lên trên 34,2‰.
- Định lượng chính xác trường muối dinh dưỡng Nitrate (đạt 1,5 đến 3,5 mmol N/m³) và sinh khối thực vật nổi cực đại (đạt 1,2 đến 1,8 mmol N/m³) tập trung tại tầng nước sâu từ 15 đến 35 m.
- Đạt độ ổn định mô hình xuất sắc với hệ số tương quan nhiệt độ r = 0,995104 sau 6 năm mô phỏng, khẳng định giá trị ứng dụng thực tiễn trong công tác dự báo ngư trường.
- Kế hoạch giai đoạn 2026-2028 sẽ tập trung nâng cao độ phân giải lưới tính và đồng hóa chuỗi số liệu vệ tinh nhằm hoàn thiện hệ thống giám sát sinh thái đại dương theo thời gian thực.
Luận văn là tài liệu học thuật tiêu biểu trong lĩnh vực hải dương học số trị tại Việt Nam. Quý độc giả, nhà khoa học và các cơ quan quản lý chuyên ngành có thể tham khảo toàn văn công trình nghiên cứu tại Trung tâm Thông tin - Thư viện Đại học Quốc gia Hà Nội để ứng dụng và mở rộng các giải pháp khai thác biển bền vững.