Tổng quan về luận án

Hệ thống đầm phá Tam Giang – Cầu Hai thuộc tỉnh Thừa Thiên Huế là vùng đầm phá nước lợ ven biển lớn nhất khu vực Đông Nam Á với diện tích mặt nước xấp xỉ 22.000 ha, trải dài hơn 70 km từ cửa sông Ô Lâu đến đầm Cầu Hai và lưu thông với Biển Đông qua hai cửa biển huyết mạch là Thuận An và Tư Hiền (Nguyễn Đính và Phạm Thị Diệu My, 2005). Thủy vực này giữ vai trò trọng yếu trong việc điều hòa thủy văn, ngăn chặn xâm nhập mặn, bảo tồn đa dạng sinh học và là sinh kế cốt lõi của hơn 300.000 cư dân ven bờ (Nguyễn Văn Hợp và ctv., 2007). Sự phát triển nhanh chóng của hoạt động nuôi trồng thủy sản (NTTS) với diện tích canh tác đạt 4.215 ha vào năm 2017—chủ lực là tôm sú (Penaeus monodon), cua biển (Scylla paramamosain), cá dìa (Siganus guttatus) và cá kình (Siganus canaliculatus)—đang tạo sức ép nghiêm trọng lên khả năng tự làm sạch của môi trường đầm phá.

Mặc dù các chương trình nghiên cứu trước đây như Dự án hợp tác vùng Nord Pas de Calais (Pháp, 1998–2003) hay Dự án Quản lý tổng hợp vùng ven bờ CZMC/RIKZ (Hà Lan, 2001–2005) đã ghi nhận nguồn dữ liệu cơ sở về sinh thái, một khoảng trống nghiên cứu (research gap) cốt lõi vẫn tồn tại: các công trình quá khứ chưa thiết lập được mối liên kết động giữa cơ cấu hiện trạng nuôi trồng với sự phân hóa không gian - thời gian của chất lượng nước (CLN); đồng thời hoàn toàn thiếu vắng công cụ chỉ số CLN chuyên biệt cho đầm phá lợ mặn và mô hình toán thủy động lực học có khả năng dự báo ô nhiễm phục vụ quy hoạch nuôi trồng thủy sản bền vững. Các nghiên cứu quốc gia về chỉ số CLN (VN-WQI do Tổng cục Môi trường ban hành năm 2011) chỉ áp dụng cho nước mặt lục địa, không tương thích với các đặc tính hóa sinh phức tạp của hệ thống đầm phá ven biển.

Để giải quyết triệt để vấn đề trên, luận án tập trung vào 4 câu hỏi nghiên cứu (Research Questions - RQ) và các giả thuyết tương ứng:

  • RQ1: Hiện trạng kỹ thuật, cơ cấu mùa vụ và mức độ xử lý chất thải trong các mô hình NTTS ven đầm phá Tam Giang – Cầu Hai diễn biến ra sao?
  • RQ2: Sự biến động không gian và chu kỳ ngày đêm của các thông số lý - hóa - sinh môi trường nước chịu sự chi phối của những nhóm nhân tố cốt lõi nào?
  • RQ3: Làm thế nào để lượng hóa CLN đầm phá thành một chỉ số tích hợp (WQI) phản ánh chuẩn xác tính tương thích sinh thái đối với đối tượng tôm nuôi nước lợ?
  • RQ4: Mức độ lan truyền và diễn biến nồng độ các chất ô nhiễm hữu cơ và dinh dưỡng tại các vùng nuôi tập trung trong tương lai sẽ thay đổi như thế nào dưới tác động của thủy động lực học?

Khung lý thuyết của nghiên cứu được xây dựng trên sự giao thoa giữa Lý thuyết Thủy động lực học Môi trường đầm phá ven biển (Coastal Lagoon Hydrodynamics), Lý thuyết Chỉ số Chất lượng Nước (Water Quality Indexing Theory) và Lý thuyết Phân tích Không gian trong NTTS (Spatial Multi-Criteria Analysis). Phạm vi nghiên cứu bao quát 90 hộ NTTS tại 9 xã đại diện ven đầm phá, mạng lưới 44 trạm quan trắc không gian qua 5 đợt/năm (tháng 2, 5, 8, 10, 12) và 4 trạm quan trắc ngày đêm liên tục trong 7 ngày với tần suất 3 giờ/lần tại Tam Giang, Thuận An, Trường Hà và Tư Hiền.

Literature Review và Positioning

Nghiên cứu về môi trường đầm phá ven biển trên thế giới được cấu thành từ ba dòng học thuật chính:

  1. Dòng nghiên cứu động thái thủy hóa và phú dưỡng đầm phá: Các công trình kinh điển của Abigail et al. (2009), Hungspreugs (1988) và Birt (1967) đã thiết lập nền tảng về cơ chế trao đổi nước và tích tụ vật chất tại các đầm phá bán khép kín. Các nghiên cứu thực nghiệm tại châu Phi của Eruola et al. (2011) tại Nigeria, Apau et al. (2012) tại Ghana (với nồng độ nitrat lên tới 2.905,71 mg/L và photphat đạt 487,14 mg/L) và Bawa et al. (2002) tại Togo đã chứng minh tác động tiêu cực của nước thải nông nghiệp và đô thị đến tải lượng hữu cơ đầm phá.
  2. Dòng nghiên cứu chỉ số chất lượng nước (WQI): Bắt nguồn từ chỉ số Horton (1965), tiếp nối bởi Brown et al. (1970) với NSF-WQI, mô hình O-WQI tại bang Oregon, F-WQI tại Florida (Hoa Kỳ), CCME-WQI (Canada, 2001) và chỉ số hàm nhạy của Bhargava (1983) tại Ấn Độ. Tại Việt Nam, các mô hình của Tôn Thất Lãng, Phạm Thị Minh Hạnh và Nguyễn Văn Hợp ctv. (2010) đã phát triển WQI cho các lưu vực sông nội địa.
  3. Dòng nghiên cứu ứng dụng GIS và mô hình toán thủy động lực học trong thủy sản: Khởi xướng từ các nghiên cứu tiên phong của Aguilar-Manjarrez & Ross (1995), Meaden & Thang (1996), Radiarta et al. (2008) tại Nhật Bản, Longdill et al. (2008) tại New Zealand và hệ thống mô hình MIKE 21 (DHI) ứng dụng trong mô phỏng tải lượng ô nhiễm.

Trong y văn tồn tại hai trường phái tranh luận lớn: Trường phái thứ nhất ủng hộ việc đánh giá CLN dựa trên từng thông số hóa lý riêng lẻ nhằm duy trì tính chính xác tuyệt đối của nồng độ chất ô nhiễm; trường phái thứ hai khẳng định việc tổng hợp các biến số thành một chỉ số tích hợp (WQI) là con đường duy nhất để lượng hóa hiện trạng môi trường phục vụ việc ra quyết định quy hoạch. Luận án định vị nghiên cứu tại điểm giao thoa của hai trường phái: kết hợp phân tích sâu từng thông số đơn lẻ với việc thiết lập bộ chỉ số WQI đặc thù và mô hình hóa lan truyền chất.

So sánh với các nghiên cứu quốc tế điển hình:

  • Tại đầm phá Rio Lagartos và Bojorquez (Mexico), Jorge & Morales-Ojeda (2010) ghi nhận thời gian lưu nước kéo dài 400–450 ngày dẫn đến tích tụ dinh dưỡng nghiêm trọng (P-PO43- đạt 1,38 µmol/L; N-NO3- đạt 8,96 µmol/L). Luận án của Trương Văn Đàn chứng minh rằng đầm phá Tam Giang – Cầu Hai có cơ chế lưu thông phức tạp hơn do chịu tác động đồng thời của chế độ bán nhật triều không đều qua 2 cửa biển và dòng chảy lũ sông Hương, tạo ra sự phân tầng và dị biệt hóa chất lượng nước sâu sắc hơn.
  • Tại đầm Tom Thumb (Australia), Josephine et al. (2006) ghi nhận hàm lượng P-PO43- vượt quy chuẩn từ 2 đến 12,2 lần và N-NH3 vượt từ 5,3 đến 8,6 lần. Tương tự, tại đầm Washdyke (New Zealand), Graeme (2015) ghi nhận nồng độ nitơ vô cơ hòa tan trung bình 3,7 mg/L, vượt ngưỡng an toàn sinh thái. Luận án đã làm rõ sự tương đồng về áp lực phú dưỡng hóa nhưng chỉ ra nét độc thù của Tam Giang – Cầu Hai: nguồn thải từ hoạt động NTTS và nông nghiệp đóng góp trên 98% vào tải lượng ô nhiễm hữu cơ và dinh dưỡng của toàn vùng đầm phá.

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

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

Luận án mở rộng Lý thuyết Chỉ số Chất lượng Nước (WQI Theory) bằng cách giải quyết triệt để khiếm khuyết toán học trong các mô hình WQI truyền thống. Trong các mô hình của Tôn Thất Lãng hay Bhargava (1983), hàm tính chỉ số phụ $q_i$ bậc cao thường tạo ra giá trị biên bất thường (ngoài khoảng 0–100) khi nồng độ thông số tiến về 0 hoặc cực đại (ví dụ: mô hình của Tôn Thất Lãng cho $q_i = 150,23$ khi pH = 0). Luận án đã chuẩn hóa các hàm độ nhạy phi tuyến $q_i = f(x_i)$ đối với 10 thông số cốt lõi, ràng buộc giá trị $q_i$ tuyệt đối trong miền $[0, 100]$, phản ánh chính xác ngưỡng giới hạn sinh thái của tôm nuôi nước lợ.

Nghiên cứu tạo ra bước chuyển dịch quan niệm (paradigm shift) từ đánh giá tĩnh, đơn biến sang mô hình đánh giá động, đa chiều kết hợp cấu trúc phân tích thành phần chính (PCA) với mô hình động lực học chất lượng nước 2 chiều (2D Hydrodynamic-Water Quality Model).

[Nguồn thải NTTS & Nông nghiệp (>98%)] 
                   │
                   ▼
┌────────────────────────────────────────────────────────┐
│      Hệ thống Đầm phá Tam Giang - Cầu Hai              │
│  (Thủy động lực học: 2 cửa biển, gió mùa, triều)       │
└──────────────────────────┬─────────────────────────────┘
                           │
        ┌──────────────────┴──────────────────┐
        ▼                                     ▼
┌──────────────────────────────┐ ┌──────────────────────────────┐
│  Thành phần Chính 1 (PC1)    │ │  Thành phần Chính 2 (PC2)    │
│  60,7% Biến thiên CLN        │ │  18,7% Biến thiên CLN        │
│  [Hữu cơ - Dinh dưỡng - Khí độc]│ │  [Vật chất hòa tan - Nền]   │
│  (NH3, TC, BOD5, TSS,        │ │  (pH, TDS, Độ mặn, Độ kiềm)  │
│   N-NO3-, P-PO43-, DO)       │ │                              │
└──────────────┬───────────────┘ └──────────────┬───────────────┘
               └───────────────┬────────────────┘
                               ▼
┌────────────────────────────────────────────────────────┐
│           Chỉ số WQITGCH (10 thông số Delphi)          │
│   Mùa khô: 52 - 84 (TB - Tốt) | Mùa mưa: 44 - 79 (Xấu) │
└──────────────────────────────┬─────────────────────────┘
                               ▼
┌────────────────────────────────────────────────────────┐
│      Mô hình Dự báo MIKE 21 (HD & ECOLAB)              │
│   Dự báo 2020: Hotspots ô nhiễm tại Phú Mỹ, Lộc Điền   │
└────────────────────────────────────────────────────────┘

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

Khung phân tích của luận án tích hợp ba trụ cột:

  1. Lý thuyết Đồng thuận Chuyên gia Delphi: Xác định 10 thông số ưu tiên và trọng số $w_i$ tương ứng: Nhiệt độ ($w=0,08$), pH ($w=0,11$), DO ($w=0,14$), Độ mặn ($w=0,09$), Độ kiềm ($w=0,09$), N-NO3- ($w=0,09$), P-PO43- ($w=0,09$), BOD5 ($w=0,11$), N-NH3 ($w=0,11$) và Tổng coliform ($w=0,09$).
  2. Kỹ thuật Phân tích Đa biến PCA: Phân rã ma trận tương quan môi trường thành các trục trực giao độc lập giải thích bản chất quá trình hóa sinh đầm phá.
  3. Mô hình Thủy động lực học MIKE 21 kết hợp Mô-đun Sinh thái ECOLAB: Mô phỏng tương tác giữa các trường vận tốc dòng chảy triều, quá trình khuếch tán bình lưu và các phản ứng phân hủy sinh hóa nội tại.

Điều kiện biên của khung phân tích được xác định rõ ràng: Áp dụng cho các thủy vực đầm phá ven biển nhiệt đới gió mùa có độ sâu nông (<3 m), biên độ triều nhỏ và chịu chi phối mạnh bởi chế độ thoát lũ theo mùa.

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

Thiết kế nghiên cứu

Nghiên cứu vận hành theo hệ hình thực chứng (positivism) với phương pháp tiếp cận định lượng đa tầng (multi-level design):

  • Tầng 1 (Kinh tế - Xã hội và Kỹ thuật canh tác): Khảo sát 90 hộ nuôi trồng thủy sản tại 9 xã trọng điểm ven đầm phá bằng bảng câu hỏi cấu trúc, đánh giá phương thức nuôi, mật độ giống, tỷ lệ kiểm dịch và tỷ lệ xử lý nước cấp/nước thải.
  • Tầng 2 (Quan trắc Không gian & Thời gian): 44 vị trí quan trắc môi trường nước trải đều trên toàn bộ không gian đầm phá (phá Tam Giang, đầm Sam Chuồn, đầm Thủy Tú, đầm Hà Trung, đầm Cầu Hai) qua 5 chu kỳ mùa vụ trong năm; 4 trạm quan trắc ngày đêm đo đạc liên tục 7 ngày/mùa với tần suất 3 giờ/lần.
  • Tầng 3 (Mô phỏng Toán học & GIS): Xây dựng lưới phi cấu trúc biểu diễn địa hình đáy đầm phá, thiết lập điều kiện biên thủy lực và tích hợp mô-đun ECOLAB trong phần mềm MIKE 21 để mô phỏng và dự báo CLN.

Quy trình nghiên cứu rigorous

Quy trình thu thập và phân tích mẫu tuân thủ nghiêm ngặt Standard Methods for the Examination of Water and Wastewater (APHA) và Tiêu chuẩn Quốc gia (TCVN):

  • Mẫu nước được thu bằng thiết bị lấy mẫu chuyên dụng ở tầng mặt và tầng giữa, bảo quản lạnh ở 4°C và vận chuyển về phòng thí nghiệm trong vòng 24 giờ.
  • Phương pháp phân tích 12 thông số: Nhiệt độ, pH, DO, Độ mặn, TDS được đo đạc trực tiếp tại hiện trường bằng máy đo đa chỉ tiêu điện tử đã hiệu chuẩn chuẩn độ. Độ kiềm được xác định bằng phương pháp chuẩn độ acid; TSS bằng phương pháp sấy khô ở 103–105°C; BOD5 bằng phương pháp ủ chai đo oxy hòa tan sau 5 ngày ở 20°C; N-NH3 xác định bằng phương pháp so màu Phenate; N-NO3- bằng phương pháp so màu khử Cadmium; P-PO43- bằng phương pháp so màu Ascorbic Acid; Tổng coliform xác định bằng phương pháp lên men nhiều ống (MPN).
  • Quy trình Delphi: Thực hiện 2 vòng khảo sát độc lập với 40 chuyên gia đầu ngành tại các trường đại học (ĐH Cần Thơ, ĐH Nông Lâm Huế) và viện nghiên cứu môi trường để xác lập độ tin cậy và trọng số thành phần của chỉ số $WQITGCH$.
  • Độ tin cậy của mô hình: Mô hình MIKE 21 được hiệu chỉnh (calibration) và kiểm định (validation) dựa trên chuỗi số liệu đo đạc thực tế về mực nước, lưu lượng và độ mặn tại các trạm hải văn Thuận An và Tư Hiền, đạt hệ số tương quan Nash-Sutcliffe $E > 0,75$.

Data và phân tích

Dữ liệu được xử lý thống kê bằng phần mềm SPSS và ngôn ngữ R. Phân tích thành phần chính (PCA) trên tập dữ liệu quan trắc đã phân tách rõ ràng 2 nhóm nhân tố chi phối chất lượng nước đầm phá:

  • Thành phần chính 1 (PC1 - Nhóm hữu cơ, vi sinh và khí độc): Bao gồm N-NH3, Tổng coliform, BOD5, TSS, N-NO3-, P-PO43- và DO, giải thích 60,7% tổng phương sai biến động môi trường.
  • Thành phần chính 2 (PC2 - Nhóm vật chất hòa tan và yếu tố nền): Bao gồm pH, TDS, độ mặn và độ kiềm, giải thích 18,7% tổng phương sai biến động. Tổng mức độ giải thích của hai thành phần chính đạt 79,4%, khẳng định tính đại diện tuyệt đối của tập thông số được lựa chọn.

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

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

  1. Sự phù hợp của nhóm thông số nền và sự suy thoái cục bộ theo mùa: Nhóm thông số cơ bản (nhiệt độ, pH, DO, độ mặn, độ kiềm) trên toàn bộ khối nước mở nhìn chung nằm trong ngưỡng thích hợp cho đời sống thủy sinh. Tuy nhiên, vào mùa mưa, giá trị pH và độ kiềm tại các khu vực ven bờ, vùng cửa sông tiếp nhận nước lũ và lân cận kênh xả thải NTTS/sinh hoạt bị suy giảm nghiêm trọng, rơi xuống dưới ngưỡng an toàn cho nuôi tôm sú nước lợ.
  2. Dị biệt hóa không gian giữa các tiểu vực đầm phá: Phân tích không gian chỉ ra rằng tiểu vực đầm Hà Trung – Thủy Tú có độ mặn và độ kiềm cao nhất trong toàn bộ hệ thống đầm phá, hàm lượng BOD5 thấp nhất nhưng nồng độ oxy hòa tan (DO) lại thấp hơn đáng kể so với phá Tam Giang và đầm Cầu Hai. Ngược lại, khu vực đầm Sam Chuồn là điểm nóng tích tụ dinh dưỡng với nồng độ N-NO3- và P-PO43- cao nhất toàn hệ thống.
  3. Ô nhiễm hữu cơ và dinh dưỡng tiệm cận và vượt ngưỡng quy chuẩn: Nồng độ BOD5 và P-PO43- tại hầu hết các trạm ven bờ đều vượt giới hạn cho phép của Quy chuẩn kỹ thuật Quốc gia đối với nước phục vụ NTTS. Khí độc N-NH3 tiệm cận mức nguy hiểm sinh thái, với nồng độ ban ngày cao hơn ban đêm và mùa mưa cao hơn mùa khô.
  4. Động thái biến thiên chỉ số $WQITGCH$: Chỉ số chất lượng nước đầm phá biến động từ 44 – 79 trong mùa mưa (chất lượng từ xấu đến trung bình) và từ 52 – 84 trong mùa khô (chất lượng từ trung bình đến tốt). Sự suy giảm chỉ số trong mùa mưa phản ánh áp lực trực tiếp từ dòng chảy tràn bề mặt mang theo chất ô nhiễm từ lưu vực nông nghiệp và khu dân cư đổ vào đầm phá.
  5. Dự báo ô nhiễm bằng mô hình MIKE 21 cho năm 2020: Kết quả mô phỏng lan truyền chất lượng nước trong các đợt tháng 2, tháng 5 và tháng 8 năm 2020 cảnh báo rõ ràng rằng các vùng mặt nước ven xã Phú Mỹ (huyện Phú Vang) và xã Lộc Điền (huyện Phú Lộc) đối mặt với nguy cơ ô nhiễm hữu cơ và phú dưỡng nghiêm trọng, thể hiện qua các dải nồng độ BOD5 và P-PO43- tăng cao đột biến.

Implications đa chiều

  • Ý nghĩa học thuật và lý thuyết: Xác lập cơ sở khoa học về động thái chất ô nhiễm trong đầm phá ven biển nhiệt đới; cung cấp công thức toán học chuẩn hóa cho chỉ số $WQITGCH$ và bộ chỉ số hiệu chỉnh khi thiếu số liệu ($WQITGCHhieuchinh$), mở đường cho các nghiên cứu đánh giá sinh thái đầm phá tại các quốc gia đang phát triển.
  • Đổi mới phương pháp luận: Thiết lập quy trình công nghệ kết hợp ba tầng: Viễn thám/GIS – Phân tích đa biến/WQI – Mô phỏng thủy động lực học MIKE 21, cho phép chuyển đổi từ quan trắc bị động sang quản lý và dự báo chủ động.
  • Ứng dụng thực tiễn cho ngành thủy sản: Cung cấp bản đồ phân vùng lấy nước cấp an toàn cho 4.215 ha ao nuôi ven đầm phá; khuyến cáo nông dân không lấy nước vào các thời điểm giao mùa hoặc sau các đợt mưa lớn khi độ kiềm và pH giảm thấp và độc tính khí độc N-NH3 tăng cao.
  • Khuyến nghị chính sách: Cung cấp cơ sở dữ liệu định lượng cho Chi cục Thủy sản và Sở Nông nghiệp & PTNT tỉnh Thừa Thiên Huế trong việc tái cơ cấu diện tích nuôi tôm, thu hẹp hình thức nuôi chắn sáo gây cản trở lưu thông thủy triều, và bắt buộc đầu tư hệ thống xử lý nước thải tập trung tại các xã Phú Mỹ và Lộc Điền.

Limitations và Future Research

Mặc dù đạt được những kết quả đột phá, luận án thẳng thắn thừa nhận 4 giới hạn nghiên cứu:

  1. Giới hạn mô hình hóa 2 chiều (2D): Mô hình MIKE 21 HD giả định khối nước đồng nhất theo chiều đứng, chưa phản ánh đầy đủ hiện tượng phân tầng mật độ (halocline và thermocline) xảy ra cục bộ tại các lạch sâu và luồng tàu gần cửa Thuận An trong mùa lũ.
  2. Độ phân giải nguồn thải lưu vực: Số liệu tải lượng ô nhiễm từ nông nghiệp và đô thị chủ yếu dựa trên ước tính hệ số phát thải và điều tra xã hội học, chưa có trạm đo tự động liên tục tại tất cả các cửa sông đổ vào đầm phá.
  3. Phạm vi thông số độc học chuyên biệt: Nghiên cứu chưa tích hợp chuỗi số liệu liên tục về dư lượng hóa chất bảo vệ thực vật và kim loại nặng dạng vết vào mô-đun mô phỏng động ECOLAB.

Chương trình nghiên cứu tương lai (Future Research Agenda) định hướng:

  • Nâng cấp mô hình lên không gian 3 chiều (MIKE 3 Flow Model) để mô phỏng chính xác quá trình phân tầng thủy văn và thiếu khí tầng đáy.
  • Tích hợp mạng lưới cảm biến IoT quan trắc tự động CLN thời gian thực (Real-time monitoring) truyền dữ liệu trực tiếp vào hệ thống hỗ trợ ra quyết định trực tuyến.
  • Mở rộng đánh giá sức tải môi trường (carrying capacity) của đầm phá dưới các kịch bản biến đổi khí hậu, nước biển dâng và đóng/mở các cửa biển định kỳ.

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

  • Tác động học thuật: Luận án cung cấp bộ dữ liệu quy chuẩn toàn diện nhất về môi trường đầm phá Tam Giang – Cầu Hai, đóng vai trò là tài liệu tham khảo nền tảng cho các nghiên cứu sau đại học về hải dương học hóa học, quản lý môi trường đới bờ và kỹ thuật nuôi trồng thủy sản nước lợ.
  • Tác động ngành và kinh tế: Tối ưu hóa hiệu quả kinh tế cho cộng đồng nuôi tôm ven đầm phá. Việc áp dụng bản đồ phân vùng CLN giúp giảm thiểu rủi ro dịch bệnh do suy giảm chất lượng nước, hạn chế thiệt hại hàng chục tỷ đồng mỗi năm cho người nuôi thủy sản tại Thừa Thiên Huế.
  • Tác động chính sách: Trực tiếp đóng góp luận cứ khoa học cho Quy hoạch tổng thể phát triển kinh tế - xã hội vùng đầm phá Tam Giang – Cầu Hai đến năm 2030 của UBND tỉnh Thừa Thiên Huế, định hướng chuyển đổi phương thức NTTS sang đa dạng sinh học và nuôi sinh thái bền vững.

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

  • Nghiên cứu sinh và Giới học thuật: Tiếp cận khung phương pháp luận kết hợp Delphi, PCA, WQI và MIKE 21 hoàn chỉnh để áp dụng cho 11 hệ thống đầm phá ven biển còn lại của miền Trung Việt Nam (đầm Lăng Cô, Trường Giang, Nước Mặn, Thị Nại, Đầm Nại...).
  • Cơ quan Quản lý Nhà nước (Tổng cục Thủy sản, Sở TN&MT, Chi cục Thủy sản): Sở hữu công cụ phần mềm trực tuyến tính toán nhanh chỉ số $WQITGCH$ và bản đồ dự báo ô nhiễm để chỉ đạo mùa vụ và cấp phép xả thải.
  • Doanh nghiệp và Hộ nuôi trồng Thủy sản: Nắm bắt chính xác lịch biến động chất lượng nước theo ngày đêm và theo mùa, chủ động kỹ thuật cấp/thoát nước và quản lý môi trường ao nuôi an toàn.

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

1. Đóng góp lý thuyết độc đáo nhất của luận án là gì?

Đóng góp lý thuyết độc đáo nhất là việc chuẩn hóa và tối ưu hóa toán học hàm độ nhạy phụ $q_i = f(x_i)$ cho 10 thông số môi trường nước lợ trong khuôn khổ Lý thuyết Chỉ số Chất lượng Nước (WQI), loại bỏ triệt để hiện tượng tràn giá trị biên thường gặp ở các mô hình WQI truyền thống (như của Tôn Thất Lãng hay Bhargava), đồng thời tích hợp thành công trọng số Delphi chuyên ngành nuôi tôm vào cấu trúc phân tích đa biến PCA.

2. Đổi mới phương pháp luận so với các nghiên cứu trước đây được thể hiện như thế nào?

So với các nghiên cứu đánh giá CLN đơn thuần bằng quan trắc thống kê mô tả (như Nguyễn Văn Hợp ctv., 2007 tại Tam Giang – Cầu Hai hay Nguyễn Hồng Thu ctv., 2013 tại Đầm Nại) hoặc các nghiên cứu ứng dụng GIS tĩnh (như Lê Công Tuấn & Lê Thị Hạnh, 2009), luận án đã tạo bước đột phá phương pháp luận khi ghép nối chuỗi dữ liệu đa tầng: từ điều tra xã hội học kỹ thuật NTTS $\rightarrow$ mạng lưới quan trắc không gian/ngày đêm 12 thông số $\rightarrow$ mô hình phân tích thành phần chính PCA $\rightarrow$ xây dựng chỉ số WQI trực tuyến $\rightarrow$ tích hợp mô hình thủy động lực học MIKE 21 2D để dự báo động thái chất lượng nước trong tương lai.

3. Phát hiện bất ngờ nhất (counter-intuitive result) từ dữ liệu thực nghiệm là gì?

Phát hiện bất ngờ nhất là đặc tính nghịch lý tại tiểu vực đầm Hà Trung – Thủy Tú: mặc dù đây là khu vực có độ mặn và độ kiềm cao nhất hệ thống đầm phá và hàm lượng chất hữu cơ BOD5 thấp nhất, nhưng nồng độ oxy hòa tan (DO) tại đây lại suy giảm nghiêm trọng hơn so với phá Tam Giang và đầm Cầu Hai. Nguyên nhân do hình thái lòng đầm hẹp, khả năng lưu thông trao đổi nước kém hơn và quá trình tiêu thụ oxy hóa học tầng đáy diễn ra mạnh mẽ.

4. Luận án có cung cấp quy trình tái lập nghiên cứu (replication protocol) không?

Có. Luận án cung cấp chi tiết toàn bộ bảng tọa độ 44 trạm quan trắc không gian, 4 trạm ngày đêm, quy chuẩn phân tích phòng thí nghiệm theo TCVN/APHA, ma trận trọng số Delphi ($w_i$), các phương trình hàm số $q_i$, các thông số thiết lập điều kiện biên và hệ số suy giảm sinh hóa trong mô-đun ECOLAB của mô hình MIKE 21, cho phép tái lập nghiên cứu hoàn toàn trên các hệ thống thủy vực đầm phá tương đồng.

5. Chương trình nghiên cứu 10 năm được vạch ra như thế nào?

Chương trình 10 năm tập trung vào 3 giai đoạn:

  • Năm 1–3: Xây dựng hệ thống trạm quan trắc trực tuyến IoT và chuyển đổi mô hình MIKE 21 sang mô hình 3D (MIKE 3) có xét đến phân tầng triều - lũ.
  • Năm 4–6: Mở rộng tích hợp bộ chỉ số $WQITGCH$ cho toàn bộ 12 đầm phá miền Trung Việt Nam, xây dựng nền tảng dữ liệu lớn (Big Data) cảnh báo dịch bệnh thủy sản do môi trường.
  • Năm 7–10: Nghiên cứu tác động tổng hợp của biến đổi khí hậu, axit hóa đại dương và suy giảm dòng chảy ngọt thượng nguồn lên sức tải sinh thái đầm phá phục vụ an ninh lương thực quốc gia.

Kết luận

Luận án tiến sĩ của Trương Văn Đàn đã hoàn thành xuất sắc mục tiêu nghiên cứu với 6 đóng góp cốt lõi:

  1. Xác lập bức tranh toàn diện và cập nhật nhất về hiện trạng NTTS trên diện tích 4.215 ha ven đầm phá Tam Giang – Cầu Hai, chỉ ra sự thiếu hụt nghiêm trọng trong hạ tầng xử lý nước cấp và nước thải.
  2. Phân rã và lượng hóa thành công các cơ chế chi phối chất lượng nước thành 2 nhóm thành phần chính: Nhóm hữu cơ - vi sinh - khí độc (chiếm 60,7% phương sai) và Nhóm vật chất hòa tan - yếu tố nền (chiếm 18,7% phương sai).
  3. Thiết lập thành công bộ chỉ số chất lượng nước đầm phá $WQITGCH$ gồm 10 thông số chuẩn hóa và bộ chỉ số hiệu chỉnh $WQITGCHhieuchinh$, cho phép lượng hóa chính xác mức độ an toàn sinh thái cho nuôi tôm nước lợ.
  4. Xây dựng phần mềm trực tuyến hỗ trợ tính toán nhanh chỉ số WQI và số hóa bản đồ phân vùng môi trường nước phục vụ quản lý trực quan.
  5. Hiệu chuẩn và ứng dụng thành công mô hình toán MIKE 21 (HD & ECOLAB) trong mô phỏng thủy động lực học và dự báo chính xác các điểm nóng ô nhiễm hữu cơ tại Phú Mỹ và Lộc Điền cho năm 2020.
  6. Mở ra phương pháp tiếp cận tích hợp công nghệ không gian - toán mô hình - chỉ số sinh thái trong quy hoạch không gian biển và đầm phá ven bờ tại Việt Nam.