Tổng quan về luận án

Hệ đầm phá Tam Giang - Cầu Hai (tỉnh Thừa Thiên Huế) là hệ đầm phá ven bờ lớn nhất khu vực Đông Nam Á với diện tích mặt nước khoảng 216 $\text{km}^2$, chiều dài trải dài 68 km dọc bờ biển và là nơi cư trú của gần 1.000 loài động thực vật thủy sinh có giá trị kinh tế cao (Trần Đức Thạnh, 1985; Nguyễn Huy Yết, 2011). Tuy nhiên, áp lực nhân sinh từ 5 huyện/thị xã phụ cận (Phong Điền, Quảng Điền, Hương Trà, Phú Vang, Phú Lộc) với quy mô dân số 615.946 người (mật độ 232 người/$\text{km}^2$), sự bùng nổ diện tích nuôi trồng thủy sản (NTTS) từ 2.200 ha năm 2000 lên 7.175 ha năm 2016, cùng tổng đàn gia súc gia cầm đạt 22,4 nghìn con trâu, 33,6 nghìn con bò, 205,6 nghìn con lợn và 2,7 triệu con gia cầm đã tạo ra lưu lượng xả thải hữu cơ, dinh dưỡng khổng lồ đưa vào thủy vực.

Khoảng trống nghiên cứu (Research Gap) cốt lõi trước công trình này là sự thiếu vắng một công cụ định lượng tích hợp cơ chế thủy động lực học 3 chiều với các quá trình chuyển hóa sinh địa hóa (C, N, P) để xác định giới hạn chịu tải môi trường (Environmental Carrying Capacity - ECC) cho toàn bộ hệ thống đầm phá. Các nghiên cứu trước đây (Nguyễn Hữu Cử, 2006; IMOLA, 2008) chủ yếu dừng lại ở mức độ quan trắc hiện trạng tĩnh hoặc đánh giá cục bộ từng tiểu vùng sinh thái (như đầm Sam Chuồn) mà chưa thiết lập được ngưỡng xả thải tối đa theo các kịch bản động học đa thời gian.

Công trình khoa học của NCS. Cao Thị Thu Trang thuộc chuyên ngành Kỹ thuật Môi trường (Trường Đại học Bách khoa Hà Nội) dưới sự hướng dẫn của TS. Trịnh Thành và GS.TS. Trần Đức Thạnh đã giải quyết triệt để bài toán này. Nghiên cứu tập trung vào 3 câu hỏi nghiên cứu chính:

  1. Tổng tải lượng ô nhiễm các hợp chất hữu cơ và dinh dưỡng (BOD, COD, N, P) phát sinh từ các nguồn lục địa và nội vực đổ vào hệ đầm phá Tam Giang - Cầu Hai là bao nhiêu ở thời điểm cơ sở (2011-2012) và dự báo đến năm 2020, 2030?
  2. Cơ chế thủy động lực học, quá trình trao đổi nước qua hai cửa biển Thuận An - Tư Hiền và các quá trình tự làm sạch sinh - lý - hóa chi phối trường nồng độ các chất ô nhiễm diễn ra như thế nào trong các điều kiện thủy văn đối lập (mùa mưa và mùa khô)?
  3. Ngưỡng sức tải môi trường tối đa và sức tải đề xuất bền vững của hệ đầm phá đối với các thông số C, N, P là bao nhiêu nhằm bảo vệ hệ sinh thái thủy sinh và đáp ứng Quy chuẩn kỹ thuật Quốc gia (QCVN 08-MT:2015/BTNMT, QCVN 10-MT:2015/BTNMT)?

Hệ thống giả thuyết nghiên cứu gồm:

  • $H_1$: Khả năng tự làm sạch và sức chịu tải môi trường của hệ đầm phá Tam Giang - Cầu Hai phân hóa sâu sắc theo không gian (Tam Giang, Thủy Tú, Cầu Hai) và thời gian (mùa mưa lũ so với mùa kiệt), phụ thuộc chặt chẽ vào cấu trúc dòng chảy và chu kỳ triều tại các cửa mở.
  • $H_2$: Tải lượng ô nhiễm hữu cơ và dinh dưỡng dự báo theo kịch bản phát triển kinh tế - xã hội năm 2020 và 2030 sẽ vượt quá khả năng đồng hóa tự nhiên của hệ thống nếu không kiểm soát nguồn thải phân tán từ NTTS và sinh hoạt.
  • $H_3$: Việc áp dụng công cụ mô hình hóa số trị Delft-3D tích hợp mô-đun thủy động lực học (FLOW) và chất lượng nước (WAQ) cho phép xác định chính xác tổng tải lượng cho phép tiếp nhận hàng ngày (TMDL) theo từng phân vùng chức năng.

Khung lý thuyết của luận án được xây dựng dựa trên sự tích hợp giữa lý thuyết cân bằng vật chất bảo toàn khối lượng trong thủy vực (Vlađimirov), nguyên lý sức tải môi trường biển của Nhóm chuyên gia GESAMP (1986) và lý thuyết động thái thủy động lực học - sinh thái đầm phá ven bờ (Estuarine and Coastal Hydrodynamics & Water Quality). Đóng góp đột phá của công trình là lượng hóa được sức tải môi trường cho 5 thông số ô nhiễm chính ($\text{BOD}_5$, $\text{COD}$, $\text{N-NH}_4^+$, $\text{N-NO}_3^-$, $\text{P-PO}_4^{3-}$), xác định chính xác tỷ lệ trao đổi nước ngày đêm giữa ba phân vùng Tam Giang - Thủy Tú - Cầu Hai, cung cấp cơ sở dữ liệu số trị đáng tin cậy phục vụ quy hoạch không gian biển và quản lý tích hợp vùng bờ tỉnh Thừa Thiên Huế.

Literature Review và Positioning

Nghiên cứu về sức tải môi trường (ECC) trên thế giới đã trải qua hơn ba thập kỷ tiến hóa với các luồng tư tưởng học thuật chính:

Luồng nghiên cứu thứ nhất tập trung vào mô hình sinh địa hóa - dinh dưỡng phục vụ xác định mật độ nuôi trồng tối đa trong NTTS. Điển hình là các công trình của Duarte và cs (2003), Paul Shin (2003), Brigolin và cs (2008), Ramón và cs (2010), Carrie Byron và cs (2011), Lotta và cs (2016). Các tác giả này ứng dụng mô hình tương tác giữa sinh khối loài hai mảnh vỏ, động học oxy hòa tan và chu trình vật chất để xác định sức tải sinh thái (Ecological carrying capacity). Tiêu biểu như nghiên cứu của Nakata và cs (1983) và mô hình HEM-3D tại vịnh Goseung (Hàn Quốc) xác định mật độ thả nghêu tối ưu là 16 cá thể/$\text{m}^3$ tương ứng sức tải 1.500 tấn thịt; hay nghiên cứu của vịnh Hiroshima (Nhật Bản) chứng minh để nuôi ngao bền vững thì tải lượng photpho tổng từ sông Ohta phải kiểm soát ở mức 0,5 tấn/ngày (Furuya, 2003).

Luồng nghiên cứu thứ hai tiếp cận theo phương pháp Tổng tải lượng xả thải tối đa hàng ngày (Total Maximum Daily Load - TMDL) kết hợp mô hình số trị thủy động lực học 3D. Hướng tiếp cận này khởi nguồn từ định nghĩa của GESAMP (1986): "Sức tải môi trường là một tài sản của môi trường và được xác định như khả năng chứa đựng các hoạt động hay là tốc độ các hoạt động... mà không có các tác động bất lợi". Nhóm tác giả Chang Hee-Lee (2003) tại hồ Shihwa (Hàn Quốc), Ayeon Lee và cs (2013) tại lưu vực Anyangcheon, Boyacioglu và Alpaslan (2008) tại Thổ Nhĩ Kỳ, và đặc biệt là chuỗi nghiên cứu tại Trung Quốc như Li KeQiang và cs (2012, 2014) tại vịnh Jiaozhou và vịnh Jinzhou, Zhao và cs (2012, 2015) tại hồ Fuxian và vịnh Quanzhou, Enhui và cs (2013) tại vịnh Hạ Môn đã chuẩn hóa quy trình phân bổ tải lượng ô nhiễm cho phép dựa trên mô hình mô phỏng lan truyền chất.

Luồng nghiên cứu thứ ba phát triển hệ thống chỉ số tích hợp (AHP - Analytic Hierarchy Process) đánh giá sức tải tổng hợp vùng ven biển (Yongquan Yin, 2010; Shuping Wang, 2017; Meriem Naimi Ait-Aoudia và cs, 2016). Han Mei và cs (2010) đã tổng hợp và phân loại sức tải môi trường nước thành 4 nhóm bản chất, nhấn mạnh mối quan hệ tương hỗ giữa khả năng tự làm sạch tự nhiên và các tiêu chuẩn môi trường mục tiêu.

Trong bức tranh tổng quan đó, tranh luận học thuật nổi bật nằm ở mâu thuẫn giữa hai trường phái: Trường phái tiếp cận tĩnh dựa trên công thức dung tích lưu giữ của GESAMP ($EC_{max} = (C_{tiêu\ chuẩn} - C_{hiện\ tại}) \cdot V / \tau$) và trường phái tiếp cận động lực học số trị 3D. Nhược điểm lớn của cách tiếp cận GESAMP cổ điển là bỏ qua tính dị hướng của trường dòng chảy, sự biến thiên nồng độ theo thời gian thực và các phản ứng sinh hóa nội sinh phức tạp.

Luận án này định vị chính xác vào khoảng trống học thuật: Ứng dụng mô hình toán thủy động lực học phi tuyến Delft-3D kết hợp trường số liệu khảo sát thực địa chi tiết để mô phỏng sự lan truyền và phân rã các hợp chất $\text{C, N, P}$ trong một thủy vực đầm phá ven bờ nhiệt đới gió mùa có hình thái luồng lạch phức tạp và hai cửa biển trao đổi nước không đối xứng (cửa Thuận An và cửa Tư Hiền). So với các nghiên cứu quốc tế tại hồ Fuxian (Trung Quốc) hay hồ Shihwa (Hàn Quốc), hệ đầm phá Tam Giang - Cầu Hai có độ nông trung bình nhỏ (1,5 - 2m), cấu trúc bồi lắng luồng lạch phân nhánh mạnh và biên độ triều thấp (35 - 50 cm tại Thuận An), tạo ra thách thức tính toán thủy động lực học cao hơn đáng kể.

Đó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 và làm sâu sắc thêm lý thuyết bảo toàn khối lượng và cân bằng chất ô nhiễm trong thủy vực ven biển của Vlađimirov:

$$\Delta C = C_{vào} - C_{ra}$$

"Ở đây $\Delta C$ là gia lượng nồng độ chất ô nhiễm sau thời gian $t$, $C_{vào}$ là phát thải chất ô nhiễm vào thủy vực... và $C_{ra}$ là lượng chất ô nhiễm ra khỏi thủy vực trong thời gian $t$... Sự bảo tồn trạng thái bình thường của môi trường biển đòi hỏi phải thỏa mãn những điều kiện hạn chế: $0 < C \le \text{NĐTHCP}$ và $\Delta C < 0$" (Vlađimirov).

Luận án đóng góp vào lý thuyết sức tải môi trường qua các luận điểm then chốt:

  1. Lý thuyết phân hóa sức tải đa tầng: Luận án chứng minh rằng sức tải môi trường không phải là một giá trị hằng số cố định cho toàn bộ thủy vực, mà là một trường không gian biến thiên phụ thuộc vào "tỷ lệ trao đổi nước cục bộ" và "thời gian lưu giữ thủy lực". Khả năng đồng hóa chất ô nhiễm hữu cơ và dinh dưỡng của đầm Tam Giang cao gấp nhiều lần đầm Cầu Hai do tiếp giáp trực tiếp với cửa mở Thuận An và nguồn bổ cập từ sông Hương.
  2. Cơ chế liên kết phản ứng chuyển hóa C-N-P trong môi trường nước lợ ven bờ: Luận án làm sáng tỏ động học tiêu thụ oxy hòa tan ($\text{DO}$) dưới tác động đồng thời của quá trình oxy hóa sinh hóa ($\text{BOD}$), oxy hóa hóa học ($\text{COD}$), quá trình nitrat hóa ($\text{N-NH}_4^+ \to \text{N-NO}_2^- \to \text{N-NO}_3^-$) và nhu cầu oxy của trầm tích đáy ($\text{SOD}$).
  3. Mô hình định lượng 3 ngưỡng sức tải: Thiết lập khung lý thuyết so sánh giữa 3 mức sức tải: (i) Sức tải theo Quy chuẩn Quốc gia (QCVN), (ii) Sức tải tối đa đồng hóa (ngưỡng nguy hiểm - assimilative capacity limit), và (iii) Sức tải đề xuất an toàn cho hệ sinh thái thủy sinh.

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

Khung phân tích của luận án tích hợp liên ngành giữa:

  • Lý thuyết Thủy khí động lực học tính toán (Computational Fluid Dynamics): Mô tả trường vận tốc dòng chảy 2 chiều/3 chiều qua hệ phương trình Navier-Stokes trung bình Reynolds (RANS) kết hợp phương trình bảo toàn khối lượng liên tục.
  • Lý thuyết Sinh thái học hệ thống (Systems Ecology): Mô tả chu trình tuần hoàn vật chất và dòng năng lượng dinh dưỡng khoáng trong đầm phá ven bờ.
  • Lý thuyết Quản lý môi trường tích hợp vùng bờ (Integrated Coastal Zone Management - ICZM): Gắn kết sức tải môi trường với các kịch bản quy hoạch phát triển kinh tế - xã hội của tỉnh Thừa Thiên Huế đến năm 2020 và tầm nhìn 2030.

Điều kiện biên lý thuyết (Boundary Conditions) được kiểm soát chặt chẽ: Thủy vực nghiên cứu được giới hạn bởi hệ đầm phá Tam Giang - Cầu Hai với hai biên mở thông ra Biển Đông tại cửa Thuận An và cửa Tư Hiền, các biên sông tiếp nhận nước ngọt từ hệ thống sông Ô Lâu, sông Hương, sông Truồi, sông Cầu Hai, cùng điều kiện khí tượng - thủy triều đại diện cho hai mùa đặc trưng (mùa mưa tháng 11 và mùa khô tháng 5).

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

Thiết kế nghiên cứu

Nghiên cứu được thiết kế theo trường phái thực chứng thực nghiệm kết hợp mô hình hóa số trị tiên tiến (Positivist Numerical Paradigm). Thiết kế bao gồm nghiên cứu đa mức (Multi-level empirical design):

  • Cấp độ vĩ mô: Toàn bộ lưu vực hứng nước đầm phá Tam Giang - Cầu Hai và dải ven bờ Thừa Thiên Huế.
  • Cấp độ trung mô: Ba phân vùng sinh thái chức năng gồm Phá Tam Giang, Đầm Sam - Thủy Tú và Đầm Cầu Hai.
  • Cấp độ vi mô: Các điểm xả thải cục bộ và các trạm đo đạc thông số thủy lý - thủy hóa, trầm tích học.

Khảo sát thực địa và lấy mẫu được triển khai theo quy trình chuẩn của Bộ Tài nguyên và Môi trường, tập trung vào hai đợt đo đạc liên tục đại diện cho hai chế độ thủy văn cực đoan: Mùa mưa (ngày 25-26/11/2011) và Mùa khô (ngày 19-20/5/2012), kết hợp dữ liệu kiểm định bổ sung vào tháng 4/2017 tại đầm Cầu Hai và chuỗi dữ liệu lịch sử giai đoạn 2004 - 2012.

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 các Tiêu chuẩn Việt Nam (TCVN) và Phương pháp chuẩn phân tích nước và nước thải (Standard Methods for the Examination of Water and Wastewater - APHA):

  • Thông số hiện trường: Nhiệt độ, $\text{pH}$, độ muối, oxy hòa tan ($\text{DO}$), độ đục được đo trực tiếp bằng máy đo đa chỉ tiêu tự động đã chuẩn hóa điện cực trước mỗi ca đo.
  • Thông số phòng thí nghiệm: $\text{BOD}_5$ (phương pháp ủ 5 ngày ở $20^\circ\text{C}$ đo DO), $\text{COD}$ (phương pháp chuẩn độ Bicromat), Amoni ($\text{N-NH}_4^+$ - phương pháp so màu Phenate), Nitrat ($\text{N-NO}_3^-$ - phương pháp so màu Axit Sunfosalixilic), Nitrit ($\text{N-NO}_2^-$), Phosphat ($\text{P-PO}_4^{3-}$ - phương pháp Axit Ascorbic), Tổng chất rắn lơ lửng ($\text{TSS}$ - phương pháp sấy khô ở $105^\circ\text{C}$).
  • Đo đạc hải văn - thủy động lực: Vận tốc và hướng dòng chảy được quan trắc tầng mặt và tầng đáy bằng máy đo dòng chảy chuyên dụng tại các mặt cắt then chốt (Thuận An, Tư Hiền, Tam Giang, Thủy Tú, Cầu Hai).
  • Kiểm định chất lượng dữ liệu (QA/QC): Phân tích mẫu lặp lại, mẫu chuẩn đối chứng, hệ số biến thiên $CV < 5%$, độ tin cậy của các phép phân tích đạt $\alpha = 0,05$.

Data và phân tích

Quá trình mô hình hóa sử dụng phần mềm số trị quốc tế Delft-3D do Viện Thủy lực Deltares (Hà Lan) phát triển:

  1. Thiết lập lưới tính toán và trường độ sâu: Xây dựng lưới cong trực giao (Curvilinear orthogonal grid) bám sát đường bờ phức tạp của hệ đầm phá, độ phân giải không gian cao tại các khu vực luồng lạch hẹp và các cửa biển. Trường độ sâu được nội suy từ bản đồ địa hình đáy biển tỷ lệ 1:25.000 và 1:50.000.
  2. Hiệu chỉnh mô hình thủy động lực (FLOW):
    • Hệ số nhám Manning ($n$): Hiệu chỉnh trong khoảng $0,020 - 0,028\ \text{s/m}^{1/3}$.
    • Hệ số nhớt xoáy ngang ($\nu_H$) và hệ số khuếch tán xoáy ngang ($D_H$): Hiệu chỉnh tối ưu để tái lập trường dòng chảy thực tế.
    • Kết quả kiểm định thống kê cho thấy sự tương thích xuất sắc giữa số liệu mô phỏng và thực đo: Vận tốc dòng chảy theo phương ngang ($U$) đạt $R^2 = 0,696$, sai số căn bậc hai trung bình $\text{RSE} = 0,0046$ trong mùa mưa; Vận tốc dòng chảy theo phương đứng ($V$) đạt $R^2 = 0,690$, $\text{RSE} = 0,0042$ trong mùa khô.
  3. Hiệu chỉnh mô hình chất lượng nước (WAQ):
    • Hiệu chỉnh các tham số sinh hóa: Tốc độ suy giảm $\text{BOD}$ ($k_{BOD} = 0,1 - 0,25\ \text{ngày}^{-1}$), tốc độ nitrat hóa ($k_{nit} = 0,05 - 0,15\ \text{ngày}^{-1}$), hệ số tái tạo khí qua bề mặt mặt nước, tốc độ lắng đọng photphat và nhu cầu oxy tầng trầm tích ($\text{SOD}$).
    • Kết quả so sánh giữa mô phỏng và quan trắc kiểm chứng tháng 4/2017 tại đầm Cầu Hai cho thấy sai số tuyệt đối của các thông số $\text{DO, BOD, COD, N, P}$ đều nằm trong khoảng sai số cho phép ($< 15%$), khẳng định độ tin cậy cao của bộ thông số mô hình.

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

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

  1. Sự suy giảm nghiêm trọng của khả năng tự làm sạch tự nhiên trong giai đoạn 2004 - 2012:

    • Hàm lượng oxy hòa tan ($\text{DO}$) có xu hướng giảm mạnh, cục bộ vào mùa khô tại đầm Thủy Tú ghi nhận giá trị thấp tới $4,7\ \text{mg/l}$ (dấu hiệu thiếu dưỡng khí đối với sinh vật thủy sinh đáy).
    • Ngược lại, nhu cầu oxy hóa học ($\text{COD}$) và $\text{BOD}_5$ tăng đột biến: Mùa mưa năm 2009, $\text{BOD}_5$ cao gấp 2,4 đến 4,4 lần so với năm 2005; $\text{COD}$ đạt đỉnh $19,69\ \text{mg/l}$ tại đầm Cầu Hai vào tháng 11/2011, vượt ngưỡng an toàn cho phép theo QCVN 10-MT:2015.
    • Nồng độ muối nitrat ($\text{N-NO}_3^-$) dao động trong khoảng $88,5 - 204,2\ \mu\text{g/l}$, vượt tiêu chuẩn chất lượng nước biển của ASEAN từ 1,4 đến 3,4 lần.
  2. Quy luật phân bố thủy động lực học và sự mất cân đối trao đổi nước:

    • Tại cửa Thuận An, dòng triều đóng vai trò thống trị: Pha triều lên có vận tốc trung bình tầng mặt đạt $51\ \text{cm/s}$ (cực đại $87\ \text{cm/s}$), tầng đáy đạt $47\ \text{cm/s}$ (cực đại $79\ \text{cm/s}$); pha triều xuống có vận tốc trung bình tầng mặt $42\ \text{cm/s}$ (cực đại $68\ \text{cm/s}$).
    • Ngược lại, đầm Cầu Hai là vùng "tù đọng động lực" với vận tốc dòng chảy rất nhỏ trong cả mùa mưa và mùa khô. Tỷ lệ trao đổi nước qua một ngày đêm tại Cầu Hai cực kỳ thấp so với Tam Giang và Thủy Tú, dẫn đến việc tích tụ các chất ô nhiễm hữu cơ và photphat hấp phụ ($\text{AAP}$) ở đáy đầm.
  3. Cấu trúc tải lượng ô nhiễm từ lục địa:

    • Nguồn phát thải sinh hoạt và nuôi trồng thủy sản đóng góp tỷ trọng lớn nhất vào tổng tải lượng các chất dinh dưỡng và chất hữu cơ đưa vào đầm phá.
    • Năm 2016, với 7.175 ha NTTS (chủ yếu là nuôi xen ghép và tôm cao triều/hạ triều), lượng thức ăn thừa và chất thải chưa qua xử lý là nguyên nhân cốt lõi gây ô nhiễm hữu cơ cục bộ.
  4. Kết quả mô phỏng các kịch bản tương lai (2020 và 2030):

    • Theo kịch bản phát triển cao và kịch bản đột xuất (xảy ra sự cố xả thải hoặc lũ cục bộ), nồng độ $\text{BOD}_5$ và $\text{COD}$ tại các khu vực nước nông, khu vực tập trung bè nuôi tôm đầm Sam Chuồn và đầm Cầu Hai sẽ vượt quy chuẩn QCVN 08/10 từ 1,5 đến 2,8 lần, đe dọa trực tiếp bãi đẻ và bãi giống tự nhiên.
  5. Lượng hóa Sức tải môi trường (ECC):

    • Luận án đã tính toán tường minh giới hạn sức chịu tải cho phép của hệ đầm phá Tam Giang - Cầu Hai đối với từng chất ô nhiễm ($\text{BOD}_5$, $\text{COD}$, $\text{N-NH}_4^+$, $\text{N-NO}_3^-$, $\text{P-PO}_4^{3-}$) theo ba kịch bản: Kịch bản QCVN, Kịch bản sức tải tối đa (đồng hóa cực hạn) và Kịch bản đề xuất bền vững sinh thái.

Implications đa chiều

  • Về mặt lý luận học thuật: Công trình chứng minh tính ưu việt của việc kết hợp mô hình số trị thủy động lực - chất lượng nước 3D trong việc giải bài toán sức tải môi trường cho các hệ đầm phá ven biển nhiệt đới có điều kiện địa hình phức tạp và bồi lấp cửa biển thường xuyên.
  • Về phương pháp luận ứng dụng: Thiết lập quy trình 4 bước chuẩn hóa từ điều tra nguồn thải, mô hình hóa số trị, hiệu chỉnh/kiểm định sai số đến phân bổ tải lượng xả thải tối đa, có thể chuyển giao và áp dụng trực tiếp cho các hệ đầm phá khác ở miền Trung Việt Nam (như đầm Lăng Cô, đầm An Khê, đầm Nước Ngọt, đầm Thị Nại).
  • Về chính sách và quản trị thực tiễn: Cung cấp luận cứ khoa học thực chứng để Sở Tài nguyên và Môi trường tỉnh Thừa Thiên Huế quy định hạn ngạch cấp phép xả thải, quy hoạch lại diện tích NTTS (cắt giảm hoặc chuyển đổi phương thức nuôi thâm canh sang nuôi sinh thái có kiểm soát nước thải), và thiết kế mạng lưới trạm quan trắc tự động liên tục tại các điểm nóng môi trường.

Limitations và Future Research

Mặc dù đạt được những kết quả đột phá, luận án thẳng thắn chỉ ra các hạn chế kỹ thuật và phạm vi:

  1. Dữ liệu biên khí tượng - hải văn: Chuỗi dữ liệu dòng chảy đo đạc liên tục ngoài hiện trường phục vụ hiệu chỉnh mô hình mới chỉ tập trung vào các đợt khảo sát ngắn ngày đại diện mùa mưa (tháng 11) và mùa khô (tháng 5), chưa bao quát hết các biến cố dị thường do biến đổi khí hậu như bão cực mạnh, lũ lịch sử hoặc hạn hán kéo dài nhiều tháng.
  2. Quá trình giải phóng dinh dưỡng từ trầm tích đáy: Tốc độ giải phóng nitơ và photpho từ lớp bùn đáy ($\text{SOD}$ và giải phóng $\text{PO}_4^{3-}$) chủ yếu được ngoại suy qua các tham số kinh nghiệm và hiệu chỉnh số trị, chưa thực hiện được các thí nghiệm vi mô trong phòng (Microcosm/Mesocosm cores) tại từng tiểu vùng trầm tích.
  3. Chưa tích hợp tương tác sinh thái bậc cao: Mô hình WAQ trong luận án tập trung vào các thông số thủy hóa chính ($\text{DO, BOD, COD, N, P}$) mà chưa ghép nối đầy đủ với mô-đun sinh khối chuỗi thức ăn (Phytoplankton - Zooplankton - Macrophytes - Benthos).

Chương trình nghiên cứu tiếp nối 10 năm tới (Future Research Agenda):

  • Phát triển mô hình tích hợp biến đổi khí hậu - thủy văn lưu vực (SWAT) kết hợp Delft-3D để đánh giá tác động của nước biển dâng và các công trình thủy điện thượng nguồn sông Hương đến chế độ mặn và sức tải đầm phá.
  • Nghiên cứu cơ chế trao đổi chất ô nhiễm tại ranh giới nước - trầm tích (Sediment-Water Interface Exchange) bằng kỹ thuật buồng đo dòng đáy tại hiện trường (In-situ benthic chamber).
  • Xây dựng hệ thống số hỗ trợ ra quyết định thời gian thực (Real-time Decision Support System - DSS) tích hợp AI và mạng lưới cảm biến IoT cảnh báo sớm nguy cơ quá tải môi trường và phú dưỡng hóa.

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

  • Tác động học thuật: Đóng vai trò là công trình tiên phong tại Việt Nam lượng hóa sức tải môi trường đầm phá ven biển bằng mô hình số trị đa chiều. Luận án tạo tiền đề cho các công bố quốc tế uy tín trong danh mục ISI/Scopus về mô hình hóa môi trường nước ven bờ và quản lý lưu vực.
  • Tác động ngành và kinh tế: Định hướng tái cơ cấu ngành thủy sản tỉnh Thừa Thiên Huế, bảo vệ hơn 4.000 tấn thủy sản khai thác hàng năm và bảo tồn giá trị dịch vụ hệ sinh thái của gần 1.000 loài thủy sinh bản địa.
  • Tác động chính sách: Cung cấp công cụ khoa học trực tiếp phục vụ triển khai Luật Bảo vệ Môi trường và Đề án "Quản lý tổng hợp dải ven biển tỉnh Thừa Thiên Huế", hỗ trợ cơ quan quản lý xác định hạn ngạch xả thải cho các khu công nghiệp phụ cận và lưu vực sông Hương.
  • Ý nghĩa quốc tế: Hệ đầm phá Tam Giang - Cầu Hai là một trong những vùng đất ngập nước có tầm quan trọng quốc tế theo tiêu chí Ramsar. Kết quả nghiên cứu đóng góp trực tiếp vào nỗ lực toàn cầu về bảo tồn đa dạng sinh học và phát triển bền vững các hệ sinh thái ven bờ châu Á.

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

  • Nghiên cứu sinh và Giới học thuật: Tiếp cận một khung phương pháp luận hoàn chỉnh, bộ thông số thủy động lực và chất lượng nước đã được hiệu chỉnh chính xác bằng Delft-3D cho đầm phá nhiệt đới.
  • Cơ quan Hoạch định Chính sách (Sở TN&MT, Sở NN&PTNT, UBND tỉnh Thừa Thiên Huế): Nắm giữ bộ số liệu định lượng về sức chịu tải môi trường tối đa và hạn ngạch xả thải để phê duyệt các dự án đầu tư phát triển kinh tế, công nghiệp, du lịch và NTTS.
  • Doanh nghiệp và Cộng đồng nuôi trồng thủy sản: Nhận thức rõ ngưỡng an toàn mật độ nuôi thả và chế độ trao đổi nước, từ đó tối ưu hóa công nghệ nuôi, giảm thiểu rủi ro dịch bệnh và suy thoái nguồn nước.
  • Các tổ chức Bảo tồn Quốc tế (IUCN, WWF, Ramsar): Sở hữu cơ sở dữ liệu khoa học tin cậy để tài trợ và triển khai các dự án bảo tồn đất ngập nước và phục hồi hệ sinh thái cỏ biển, thảm thực vật ngập mặ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ì và đã mở rộng lý thuyết nào?
Luận án đã mở rộng lý thuyết cân bằng vật chất của Vlađimirov và khung sức tải môi trường của GESAMP bằng cách chuyển dịch từ mô hình cân bằng khối lượng tĩnh đơn giản ($EC_{max} = \Delta C \cdot V / \tau$) sang "Mô hình sức tải động học không gian 3 chiều". Đóng góp độc đáo nhất là việc thiết lập được mối quan hệ định lượng phi tuyến giữa trường dòng chảy triều, hệ số phân tán xoáy ngang ($D_H, \nu_H$) và tốc độ phân hủy sinh hóa của các hợp chất cacbon, nitơ, photpho theo từng tiểu vùng chức năng của hệ đầm phá ven bờ nhiệt đới.

2. Điểm cải tiến phương pháp luận so với các nghiên cứu tiền nhiệm là gì?
So với các nghiên cứu trước đây tại Việt Nam (Nguyễn Hữu Cử, 2006; Nguyễn Chí Công, 2007) vốn chỉ sử dụng mô hình 1D/2D đơn giản hoặc phân tích thống kê nồng độ tĩnh, luận án đã triển khai thành công mô hình số trị Delft-3D tích hợp đầy đủ mô-đun thủy động lực học (FLOW) và mô-đun sinh địa hóa chất lượng nước (WAQ). Phương pháp này đã được hiệu chuẩn và kiểm định đồng thời với số liệu dòng chảy (đạt $R^2 = 0,696$, $\text{RSE} = 0,0046$) và chất lượng nước thực đo tại 3 phân vùng đối lập trong cả mùa mưa, mùa khô và các pha triều lên - triều xuống.

3. Phát hiện bất ngờ nhất từ dữ liệu nghiên cứu là gì?
Phát hiện bất ngờ nhất là sự đối lập sâu sắc về khả năng tự làm sạch giữa Phá Tam Giang và Đầm Cầu Hai. Mặc dù đầm Cầu Hai có diện tích mặt nước rất lớn nhưng do đặc điểm hình thái đáy trũng, xa cửa mở Thuận An và cửa Tư Hiền thường xuyên bị bồi lắng thu hẹp, vận tốc dòng chảy trong đầm Cầu Hai ở mức "rất nhỏ" trong cả mùa khô và mùa mưa. Điều này biến đầm Cầu Hai thành một "bẫy lắng đọng ô nhiễm" với nồng độ $\text{COD}$ thực đo lên tới $19,69\ \text{mg/l}$, làm đảo lộn giả định ban đầu cho rằng vùng nước rộng sẽ có khả năng pha loãng và sức tải cao hơn.

4. Nghiên cứu có cung cấp quy trình nhân rộng (Replication Protocol) không?
Có. Luận án cung cấp một quy trình 4 bước chuẩn hóa với đầy đủ tham số kỹ thuật: (1) Phương pháp điều tra, kiểm kê và dự báo tải lượng ô nhiễm từ nguồn điểm và nguồn diện; (2) Thiết lập lưới tính toán cong trực giao, nội suy trường độ sâu và xác định điều kiện biên hải văn - khí tượng; (3) Quy trình hiệu chỉnh hệ số nhám Manning ($n$), độ nhớt ngang ($\nu_H$) và các hằng số phân rã sinh hóa; (4) Thuật toán tính toán sức tải và phân bổ tải lượng theo các kịch bản quy chuẩn QCVN và kịch bản đồng hóa tự nhiên.

5. Lộ trình nghiên cứu phát triển trong 10 năm tới được phác thảo ra sao?
Lộ trình 10 năm tập trung vào 3 hướng then chốt: (1) Mở rộng ghép nối mô hình thủy động lực biển ven bờ với mô hình khí hậu khu vực (RCMs) và mô hình dòng chảy tràn lưu vực (SWAT) dưới tác động của biến đổi khí hậu; (2) Đo đạc và mô phỏng chi tiết thông lượng trao đổi dinh dưỡng tại ranh giới bùn đáy - cột nước ($\text{Benthic Flux}$); (3) Xây dựng nền tảng bản sao số (Digital Twin) đầm phá Tam Giang - Cầu Hai phục vụ quản trị thông minh theo thời gian thực.

Kết luận

  1. Luận án đã giải quyết thành công bài toán khoa học cấp bách về đánh giá định lượng sức tải môi trường đối với các hợp chất $\text{C, N, P}$ cho hệ đầm phá Tam Giang - Cầu Hai, hệ đầm phá ven bờ lớn nhất Việt Nam và Đông Nam Á.
  2. Thiết lập thành công và làm chủ công nghệ mô hình hóa số trị tiên tiến Delft-3D (FLOW & WAQ), đạt độ chính xác cao trong mô phỏng trường dòng chảy ($R^2 \approx 0,70$, $\text{RSE} < 0,005$) và lan truyền các thông số chất lượng nước ($\text{DO, BOD, COD, N, P}$).
  3. Lượng hóa chi tiết tải lượng chất thải phát sinh từ 5 nguồn thải chính (sinh hoạt, công nghiệp, chăn nuôi, nuôi trồng thủy sản, rửa trôi đất) tại thời điểm 2011 và dự báo chính xác xu hướng gia tăng đến các năm 2020 và 2030.
  4. Phát hiện quy luật phân hóa không gian - thời gian của khả năng tự làm sạch tự nhiên: Phá Tam Giang có sức tải cao nhờ động lực triều cửa Thuận An, trong khi đầm Cầu Hai có sức tải rất thấp do chế độ lưu thông nước kém, tiềm ẩn nguy cơ phú dưỡng hóa và suy kiệt oxy hòa tan tầng đáy.
  5. Đề xuất hệ thống ngưỡng sức tải môi trường bền vững và các giải pháp quản lý kỹ thuật, chính sách phân vùng sinh thái, kiểm soát hạn ngạch xả thải, tạo cơ sở khoa học vững chắc cho chiến lược phát triển kinh tế - xã hội bền vững vùng đầm phá tỉnh Thừa Thiên Huế.
  6. Mở ra hướng nghiên cứu mới về mô hình hóa tích hợp đa tầng hệ sinh thái đầm phá nhiệt đới, đóng góp quan trọng vào kho tàng khoa học môi trường và hải dương học vùng ven bờ Việt Nam.