Tổng quan về luận án

Hệ đầm phá Tam Giang - Cầu Hai (tỉnh Thừa Thiên Huế) là hệ thống đầ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 km², trải dài 68 km dọc bờ biển và là nơi cư ngụ của 889 loài sinh vật cùng hơn 1.000 loài động, thực vật thủy sinh có giá trị kinh tế cao, hàng năm cung cấp hơn 4.000 tấn thủy hải sản. Khu vực này đang tiếp nhận áp lực nhân sinh khổng lồ từ lưu vực rộng 2.648,76 km² với 615.946 nhân khẩu (năm 2016, mật độ trung bình 232 người/km²), cùng với 7.175 ha nuôi trồng thủy sản (NTTS), đàn gia súc gia cầm quy mô lớn gồm 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ải lượng ô nhiễm hữu cơ và dinh dưỡng phát sinh từ hoạt động sinh hoạt, du lịch, nông nghiệp và NTTS không ngừng đổ vào đầm phá, đe dọa trực tiếp đến tính toàn vẹn sinh thái và chức năng tự làm sạch tự nhiên của thủy vực.

Khoảng trống nghiên cứu (research gap) cốt lõi xuất phát từ thực tế: các công trình trước đây tại Tam Giang - Cầu Hai chủ yếu dừng lại ở việc quan trắc hiện trạng chất lượng nước đơn lẻ, đánh giá đa dạng sinh học hoặc mô phỏng thủy động lực học tĩnh, hoàn toàn thiếu vắng một công cụ định lượng tích hợp có khả năng mô hình hóa động học lan truyền đa chiều gắn với sức tải môi trường (Environmental Carrying Capacity - ECC) cho các thông số hữu cơ và dinh dưỡng (C, N, P) theo các kịch bản quy hoạch kinh tế - xã hội. Luận án giải quyết trực tiếp khoảng trống này thông qua 3 câu hỏi nghiên cứu và giả thuyết khoa học:

  • RQ1: Tải lượng các hợp chất C, N, P phát sinh và thực tế xâm nhập vào hệ đầm phá Tam Giang - Cầu Hai theo các phân vùng không gian và nguồn thải là bao nhiêu ở thời điểm cơ sở (2011) và xu thế biến động đến các mốc quy hoạch 2020, 2030? (Giả thuyết H1: Tải lượng xả thải gia tăng vượt cấp số nhân từ nguồn sinh hoạt và NTTS sẽ phá vỡ cân bằng dinh dưỡng tự nhiên).
  • RQ2: Cơ chế động lực học dòng chảy và các quá trình sinh địa hóa kiểm soát sự phân bố không gian - thời gian của DO, BOD₅, COD, N-NH₄⁺, N-NO₃⁻, P-PO₄³⁻ trong hệ đầm phá diễn ra như thế nào dưới tác động của chế độ bán nhật triều qua hai cửa Thuận An và Tư Hiền? (Giả thuyết H2: Tốc độ trao đổi nước hạn chế tại các vùng lõi như Đầm Sam - Thủy Tú tạo thành các túi tích tụ ô nhiễm cục bộ trong mùa khô).
  • RQ3: Ngưỡng sức tải môi trường tối đa và sức tải đề xuất cho các hợp chất C, N, P là bao nhiêu nhằm đảm bảo vừa đáp ứng Quy chuẩn Kỹ thuật Quốc gia (QCVN 08-MT:2015/BTNMT và QCVN 10-MT:2015/BTNMT) vừa duy trì sinh cảnh thủy sinh an toàn? (Giả thuyết H3: Sức tải môi trường thực tế thấp hơn đáng kể so với ngưỡng quy chuẩn lý thuyết do sự suy giảm tốc độ đồng hóa sinh học vào mùa khô).

Khung lý thuyết của công trình được thiết lập dựa trên lý thuyết sức tải môi trường của GESAMP (1986), lý thuyết cân bằng khối lượng chất ô nhiễm Vladimirov (1986) và cơ chế tự làm sạch tự nhiên của thủy vực (Chapra, 1997). Đóng góp đột phá của nghiên cứu là việc lượng hóa chính xác giới hạn chịu tải tối đa của hệ đầm phá lớn nhất Việt Nam bằng bộ công cụ mô hình hóa số trị tiên tiến Delft-3D (kết hợp Delft3D-FLOW và Delft3D-WAQ), thiết lập cơ sở dữ liệu định lượng phục vụ công tác quy hoạch phân vùng xả thải cho 5 huyện thị ven bờ.

Literature Review và Positioning

Lý thuyết sức tải môi trường đã phát triển qua 4 thập kỷ, khởi đầu từ định nghĩa nền tảng của Nhóm chuyên gia về các khía cạnh khoa học của ô nhiễm biển (GESAMP, 1986), trong đó sức tải môi trường được xác định: "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". GESAMP đã đề xuất công thức tính sức tải dựa trên thời gian lưu của nước: $EC_{max} = (C_{tiêu\ chuẩn} - C_{hiện\ tại}) \cdot V/\tau$, đặt nền móng cho trường phái quản lý hạn ngạch tải lượng. Tiếp đó, Pedro Duarte (2003) và Carrie Byron et al. (2011) đã mở rộng khái niệm này sang lĩnh vực NTTS với 4 cấp độ: sức tải vật lý, sản xuất, sinh thái và xã hội. Trường phái thứ hai tiếp cận qua mô hình quản lý tổng tải lượng ô nhiễm tối đa hàng ngày (Total Maximum Daily Load - TMDL) do Cục Bảo vệ Môi trường Hoa Kỳ (US EPA) phát triển và được triển khai mạnh mẽ tại châu Á bởi Chang Hee-Lee (2003) tại Hồ Shihwa (Hàn Quốc), Ayeon Lee et al. (2013) tại lưu vực Anyangcheon, và Li KeQiang et al. (2012, 2014) tại Vịnh Jiaozhou (Trung Quốc). Trường phái thứ ba sử dụng phương pháp phân tích thứ bậc (Analytic Hierarchy Process - AHP) kết hợp bộ chỉ thị đa mục tiêu để đánh giá sức tải tài nguyên nước đô thị (Meriem Naimi Ait-Aoudia et al., 2016; Shuping Wang, 2017; Han Mei et al., 2010).

Trong y văn tồn tại hai luồng quan điểm đối lập sâu sắc về việc xác định ngưỡng sức tải:

  • Quan điểm kỹ thuật - tiêu chuẩn (Normative Approach): Đại diện bởi các nghiên cứu của Dingjiang et al. (2012) và Ruibin Zhang et al. (2012), cho rằng sức tải môi trường chỉ đơn thuần là khoảng chênh lệch giữa tiêu chuẩn chất lượng nước do nhà nước ban hành ($C_{tiêu\ chuẩn}$) và nồng độ quan trắc nền ($C_{hiện\ tại}$). Điểm yếu của quan điểm này là bỏ qua tính năng động của hệ sinh thái và năng lực tự làm sạch thực tế của thủy vực.
  • Quan điểm sinh thái học thích ứng (Ecological Assimilation Approach): Đại diện bởi Paul Shin (2003), Ken Furuya (2003) và Bacher (2008), khẳng định sức tải phải được đo lường bằng năng lực đồng hóa sơ cấp (Primary Productivity) và cân bằng động học oxy hòa tan (DO kinetics) cùng nhu cầu oxy tầng đáy (Sediment Oxygen Demand - SOD). Việc duy trì nồng độ dưới ngưỡng quy chuẩn pháp lý chưa chắc đã bảo vệ được cấu trúc quần xã thủy sinh nhạy cảm nếu chu trình dinh dưỡng bị đảo lộn.

Luận án định vị chính xác tại điểm giao thoa giữa hai trường phái này: không chỉ tính toán sức tải dựa trên khung pháp lý Việt Nam (QCVN 08-MT:2015 và QCVN 10-MT:2015), mà còn thiết lập ngưỡng sức tải tối đa dựa trên khả năng đồng hóa tự nhiên và kịch bản bất lợi sinh thái.

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

  • So với nghiên cứu tại Hồ Fuxian (Trung Quốc) của Zhao et al. (2012, 2015) - nơi áp dụng mô hình 3 chiều tính toán TMDL đòi hỏi cắt giảm 66% TN, 68% TP và 57% COD để đạt chuẩn loại I CNWQS, luận án của NCS. Cao Thị Thu Trang xử lý một hệ thống mở ven biển phức tạp hơn gấp nhiều lần với sự chi phối đồng thời của dòng chảy sông Hương, sông Ô Lâu, sông Bồ và quá trình trao đổi nước hai chiều qua cửa biển Thuận An và Tư Hiền.
  • So với nghiên cứu Vịnh Goseung (Hàn Quốc) của Nakata et al. (2001) và Carrie Byron et al. (2011) sử dụng mô hình HEM-3D chỉ giới hạn đánh giá sức tải sinh học cho loài nhuyễn thể hai mảnh vỏ (mật độ thả tối ưu 16 cá thể/m³ cho 1.500 tấn sinh khối), nghiên cứu này giải quyết bài toán đa nguồn thải diện rộng tích hợp toàn diện các hợp chất C, N, P từ toàn bộ lưu vực kinh tế - xã hội.

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

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

Luận án đã đóng góp vào việc hoàn thiện và mở rộng lý thuyết cân bằng chất ô nhiễm trong môi trường biển của Vladimirov (1986), được mô tả bởi phương trình: $$\Delta C = C_{vào} - C_{ra}$$

Trong đó, luận án đã làm rõ cấu trúc vi mô của $\Delta C$ trong một đầm phá nhiệt đới gió mùa: $C_{vào}$ bao gồm phát thải trực tiếp từ bờ, tải lượng rửa trôi bề mặt sông mang ra, và chất ô nhiễm đưa vào do triều; trong khi $C_{ra}$ không chỉ phụ thuộc vào cơ chế pha loãng - trao đổi nước qua cửa biển ($V/\tau$) mà còn bị chi phối mạnh mẽ bởi tốc độ phân hủy sinh hóa hiếu khí, quá trình nitrat hóa - phản nitrat hóa và tốc độ hấp phụ/giải phóng photphat tại lớp trầm tích bề mặt (Adsorbed Phosphate - AAP).

Đồng thời, nghiên cứu làm sáng tỏ mâu thuẫn trong định nghĩa Luật Bảo vệ Môi trường Việt Nam (2005 và 2014). Nếu Luật BVMT 2014 định nghĩa: "Sức chịu tải của môi trường là giới hạn chịu đựng của môi trường đối với các nhân tố tác động để môi trường có thể tự phục hồi", luận án đã chứng minh bằng thực nghiệm mô hình rằng: khi môi trường đạt đến sát ngưỡng tự phục hồi vật lý - hóa học thì nồng độ các chất độc hại như amoni không phân ly ($NH_3$) và nitrit ($N-NO_2^-$) đã vượt xa ngưỡng tử vong của ấu trùng thủy sản. Do đó, luận án thiết lập mệnh đề lý thuyết mới: Sức tải môi trường của thủy vực đầm phá ven bờ phải được xác định dựa trên giới hạn sinh thái học của quần xã sinh vật chỉ thị nhạy cảm nhất, thay vì giới hạn phục hồi cơ học của hệ thống nước.

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

Khung phân tích của luận án tích hợp liên hoàn 3 trụ cột phương pháp:

  1. Lý thuyết phát thải đơn vị (Unit Emission Factor Theory): Kiểm kê và lượng hóa toàn diện tải lượng ô nhiễm phát sinh và đưa vào từ 4 nguồn chính (nước thải sinh hoạt, sản xuất công nghiệp, chất thải chăn nuôi, hoạt động NTTS và rửa trôi đất).
  2. Động lực học chất lưu môi trường (Environmental Fluid Dynamics): Giải hệ phương trình Navier-Stokes 2 chiều/3 chiều tích hợp bảo toàn khối lượng và phương trình truyền tải - khuếch tán chất ô nhiễm (Advection-Diffusion Equation).
  3. Động học chất lượng nước sinh địa hóa (Biogeochemical Water Quality Kinetics): Mô phỏng chu trình chuyển hóa các dạng hợp chất Carbon (BOD, COD), chu trình Nitơ ($NH_4^+ \rightarrow NO_2^- \rightarrow NO_3^-$), và chu trình Photpho ($PO_4^{3-}$, AAP) gắn liền với nhu cầu oxy đáy (SOD) và động lực hòa tan oxy (DO reaeration).

Điều kiện biên và phạm vi áp dụng (boundary conditions) được xác định rõ ràng: Hệ đầm phá Tam Giang - Cầu Hai với diện tích tính toán bao phủ toàn bộ vùng nước mặt đầm phá và vùng phụ cận chịu ảnh hưởng triều, phân tách theo 2 mùa đặc trưng (mùa mưa tháng 11, mùa khô tháng 5) và các pha triều lớn - triều ròng tại 2 cửa Thuận An và Tư Hiền.

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

Thiết kế nghiên cứu

Nghiên cứu đứng trên lập trường nhận thức luận thực chứng (Positivism) kết hợp chủ nghĩa hiện thực phản biện (Critical Realism), sử dụng phương pháp luận định lượng nghiêm ngặt kết hợp giữa điều tra khảo sát thực địa đa thời gian và mô hình hóa toán học số trị chuyên sâu. Thiết kế nghiên cứu đa tầng (Multi-level design) bao gồm:

  • Tầng 1 (Lưu vực): Kiểm kê nguồn thải và ước tính tải lượng phát tán từ 5 huyện ven đầm và các phụ lưu thượng nguồn sông Hương, sông Bồ, sông Ô Lâu.
  • Tầng 2 (Thủy vực đầm phá): Khảo sát quan trắc thực địa trường thủy lực và các thông số hóa sinh tại 3 phân vùng chính: Phá Tam Giang, Đầm Sam - Thủy Tú và Đầm Cầu Hai.
  • Tầng 3 (Mô phỏng số trị): Thiết lập và giải bài toán thủy động lực học và lan truyền chất ô nhiễm trên lưới sai phân cong trực giao.

Quy trình nghiên cứu rigorous

Quy trình thu thập và phân tích dữ liệu được kiểm soát chặt chẽ:

  • Khảo sát thực địa và lấy mẫu: Triển khai các đợt khảo sát chi tiết trong 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 bổ sung năm 2016 và đợt tái kiểm chứng tháng 4/2017 tại Đầm Cầu Hai. Mẫu nước và trầm tích được thu thập, bảo quản theo TCVN và Standard Methods for the Examination of Water and Wastewater (APHA).
  • Phân tích phòng thí nghiệm: Xác định chính xác các chỉ tiêu: DO (máy đo điện cực quang học), BOD₅ (phương pháp ủ 5 ngày ở 20°C đo DO), COD (chuẩn độ pemanganat/bicromat), $N-NH_4^+$ (phương pháp Indophenol Blue), $N-NO_3^-$ (phương pháp so màu với axit sunfosalixylic), $N-NO_2^-$ (phương pháp Griess), $P-PO_4^{3-}$ (phương pháp Molipdat khử bằng axit ascorbic), TSS (phương pháp sấy khối lượng ở 105°C).
  • Tam giác đạc dữ liệu (Triangulation): Kết hợp đồng thời 4 nguồn dữ liệu: dữ liệu thống kê kinh tế - xã hội địa phương, ảnh viễn thám giải đoán biến động diện tích NTTS (2000-2016), chuỗi dữ liệu hải văn - khí tượng quốc gia (trạm Thuận An, trạm Huế), và dữ liệu quan trắc chuyên sâu của đề tài.

Data và phân tích mô hình hóa

Phần mềm mô hình số trị chuyên dụng được sử dụng là bộ công cụ Delft-3D do Viện Deltares (Hà Lan) phát triển:

  • Mô-đun Delft3D-FLOW: Mô phỏng trường độ sâu, chế độ thủy triều, trường vận tốc dòng chảy $U$ (phương ngang) và $V$ (phương thẳng đứng). Kết quả hiệu chỉnh và kiểm định cho độ tương hợp cao với số liệu quan trắc thực tế:
    • Sai số tương đối bình phương trung bình vận tốc phương $U$ mùa mưa đạt $RSE = 0,0046$, hệ số xác định $R^2 = 0,696$.
    • Vận tốc phương $V$ mùa khô đạt $RSE = 0,0042$, hệ số xác định $R^2 = 0,690$.
    • Các tham số vật lý được tối ưu hóa: Hệ số nhám Manning ($n$), hệ số khuếch tán ngang ($D_H$), và độ nhớt xoáy theo phương ngang ($\nu_H$).
  • Mô-đun Delft3D-WAQ: Mô phỏng quá trình truyền tải và chuyển hóa sinh hóa của các biến chất lượng nước (DO, BOD₅, COD, $N-NH_4^+$, $N-NO_3^-$, $P-PO_4^{3-}$). Sai số giữa mô hình và quan trắc được kiểm chứng tại nhiều điểm đo trên toàn hệ thống và kiểm chứng lại tại Đầm Cầu Hai vào tháng 4/2017 cho thấy sai lệch nằm trong giới hạn cho phép ($< 15%$).

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

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

  1. Sự bùng nổ tải lượng ô nhiễm từ các nguồn lục địa: Kết quả kiểm kê và mô hình hóa cho thấy tổng tải lượng ô nhiễm đưa vào hệ đầm phá năm 2011 là rất lớn và có xu thế gia tăng vượt bậc đến năm 2020 và 2030:
    • Tải lượng BOD₅: Gia tăng mạnh từ nguồn sinh hoạt (chiếm trên 45%) và NTTS (chiếm 35%).
    • Tải lượng Nitrat ($N-NO_3^-$) trong nước đầm phá quan trắc được dao động từ 88,5 µg/L đến 204,2 µg/L, vượt tiêu chuẩn chất lượng nước ven biển của ASEAN từ 1,4 đến 3,4 lần.
    • Nhu cầu oxy sinh hóa (BOD₅) trong mùa mưa năm 2009 đạt 5,34 mg/L (gấp 2,4 - 4,4 lần so với năm 2005); COD dao động từ 2,50 đến 10,72 mg/L, tăng gấp đôi trong giai đoạn 2005 - 2009.
  2. Quy luật phân hóa thủy động lực và bẫy ô nhiễm không gian:
    • Tốc độ dòng chảy cực đại tập trung tại cửa Thuận An (đạt $105\text{ cm/s}$ trong pha triều rút mùa lũ và $51\text{ cm/s}$ trung bình pha triều lên) và cửa Tư Hiền ($62\text{ cm/s}$).
    • Ngược lại, khu vực Đầm Cầu Hai và Đầm Thủy Tú có tốc độ dòng chảy "rất nhỏ", tỷ lệ trao đổi nước kém. Do đó hình thành hiện tượng "nghịch lý trao đổi chất": Mặc dù mùa mưa có tải lượng xả thải từ sông đổ vào lớn nhất, nhưng nhờ lưu lượng dòng chảy mạnh nên thời gian lưu ngắn; trong khi vào mùa khô, lưu lượng sông giảm kiệt (chỉ chiếm 20 - 25% tổng lượng cả năm), biên độ triều tại Tam Giang chỉ đạt 35 - 50 cm, khiến Đầm Thủy Tú rơi vào trạng thái suy giảm DO nghiêm trọng (chỉ còn $4,7\text{ mg/L}$ vào tháng 4/2009), tiệm cận mức thiếu dưỡng khí sinh thái.
  3. Phát hiện hiện tượng phân tầng và phân dị nồng độ muối dinh dưỡng: Nồng độ photphat ($P-PO_4^{3-}$) dao động từ 3,6 đến 44,1 µg/L, gia tăng đột biến vào mùa mưa do rửa trôi phân bón nông nghiệp, trong khi amoni ($N-NH_4^+$) biến thiên từ 17,1 đến 112,5 µg/L.
  4. Lượng hóa ba mức sức tải môi trường (C, N, P) của hệ đầm phá:
    • Mức 1 (Theo QCVN 08 & 10-MT:2015): Hệ đầm phá có thể tiếp nhận tải lượng tối đa nhưng chất lượng nước ở nhiều vùng nhạy cảm sẽ bị đẩy lên sát ngưỡng giới hạn nguy hiểm.
    • Mức 2 (Sức tải tối đa / Kịch bản nguy hiểm): Xác định trần đồng hóa tối đa mà tại đó hệ sinh thái bắt đầu xuất hiện hiện tượng phú dưỡng hóa và suy thoái tầng đáy (SOD tăng cao).
    • Mức 3 (Sức tải đề xuất phát triển bền vững): Thiết lập mức tải lượng khống chế cho từng phân vùng (Phá Tam Giang, Đầm Sam - Thủy Tú, Đầm Cầu Hai) với hệ số an toàn sinh thái 25 - 30%, đảm bảo nồng độ DO luôn $\ge 5,5\text{ mg/L}$, BOD₅ $\le 4\text{ mg/L}$, và Nitrat $\le 100\text{ µg/L}$.
Thông số Nồng độ quan trắc nền (2004 - 2012) Ngưỡng QCVN (08 & 10:2015) Sức tải theo QCVN Sức tải Đề xuất bền vững
DO $4,70 - 8,40\text{ mg/L}$ $\ge 5,0\text{ mg/L}$ Duy trì $\ge 5,0\text{ mg/L}$ $\ge 5,50\text{ mg/L}$
BOD₅ $0,82 - 5,34\text{ mg/L}$ $\le 6,0\text{ mg/L}$ Tiếp nhận mức biên $\le 4,00\text{ mg/L}$
COD $2,50 - 10,72\text{ mg/L}$ $\le 15,0\text{ mg/L}$ Dư địa hẹp mùa khô $\le 10,00\text{ mg/L}$
N-NO₃⁻ $88,5 - 204,2\text{ µg/L}$ Khuyến nghị ASEAN: $60\text{ µg/L}$ Vượt ngưỡng ASEAN $\le 100,0\text{ µg/L}$
P-PO₄³⁻ $3,60 - 44,10\text{ µg/L}$ Giới hạn ASEAN: $45\text{ µg/L}$ Tiếp cận ngưỡng $\le 30,0\text{ µg/L}$

Implications đa chiều

  • Về mặt lý thuyết: Chứng minh rằng mô hình động học phi tuyến tính (Delft3D) phản ánh chính xác tương tác phức tạp giữa chế độ thủy văn triều và các phản ứng sinh hóa hơn hẳn công thức tĩnh GESAMP.
  • Về mặt phương pháp luận: Chuẩn hóa quy trình 4 bước tích hợp: (1) Kiểm kê phát thải lưu vực $\rightarrow$ (2) Mô hình hóa động lực dòng chảy $\rightarrow$ (3) Mô phỏng lan truyền chất lượng nước $\rightarrow$ (4) Đảo ngược bài toán để tính toán sức tải tối ưu cho các thủy vực ven bờ nhiệt đới gió mùa.
  • Về mặt thực tiễn và chính sách:
    • Cung cấp luận cứ khoa học chính xác cho UBND tỉnh Thừa Thiên Huế và Sở Tài nguyên & Môi trường trong việc cấp phép xả thải công nghiệp, quy hoạch diện tích NTTS (khống chế không vượt quá 7.000 ha và chuyển dịch từ nuôi tôm thâm canh sang nuôi sinh thái kết hợp).
    • Đề xuất giải pháp công trình và phi công trình: nạo vét định kỳ thông luồng cửa Tư Hiền để gia tăng tốc độ trao đổi nước ($V/\tau$), xây dựng các trạm xử lý nước thải tập trung cho các khu dân cư ven đầm tại Hương Trà, Phú Vang và Phú Lộc.

Limitations và Future Research

  1. Hạn chế về chuỗi dữ liệu giám sát liên tục: Mặc dù luận án đã tổ chức các đợt đo đạc chi tiết vào mùa mưa 2011, mùa khô 2012 và khảo sát bổ sung 2016 - 2017, nhưng tần suất lấy mẫu chưa đạt mức quan trắc trực tuyến (online real-time monitoring) liên tục 365 ngày/năm để nắm bắt trọn vẹn các xung xả thải đột biến trong các trận bão lũ lịch sử.
  2. Giới hạn không gian của mô hình bùn cát và hấp phụ đáy: Mô-đun chất lượng nước chủ yếu tập trung vào các pha hòa tan của C, N, P trong cột nước; cơ chế trao đổi chất giữa lớp nước đáy và trầm tích (Sediment-Water Interface Exchange) mới được tham số hóa qua biến SOD và hệ số phân bố thực nghiệm mà chưa giải chi tiết mô hình động học bùn cát 3 chiều phi dừng.
  3. Biến động hình thái cửa biển khó lường: Cửa biển Tư Hiền có đặc tính bồi lấp và mở đột ngột sau các đợt lũ lớn (như từng bị lấp kín năm 1995 và mở lại), điều này làm biến đổi sâu sắc điều kiện biên thủy lực của mô hình trong các kịch bản dài hạn 2030.

Hướng nghiên cứu tiếp theo bao gồm:

  • Tích hợp mô hình biến đổi khí hậu toàn cầu (nước biển dâng, xâm nhập mặn cực đoan) vào điều kiện biên của Delft-3D.
  • Ứng dụng công nghệ đồng hóa dữ liệu vệ tinh (Sentinel-2, Landsat-8) và mạng lưới cảm biến IoT để cập nhật trường nhiệt độ, độ đục và diệp lục tố (Chlorophyll-a) theo thời gian thực.
  • Mở rộng đánh giá sức tải môi trường đối với các chất ô nhiễm mới nổi (Micropollutants) như dư lượng kháng sinh NTTS, hạt vi nhựa (Microplastics) và hóa chất bảo vệ thực vật hữu cơ khó phân hủy.

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

  • Tác động học thuật: Luận án mở ra hướng tiếp cận chuẩn mực cho các nghiên cứu đánh giá sức tải môi trường tại các hệ thống đầm phá, vũng vịnh ven bờ miền Trung Việt Nam (như Tam Giang - Cầu Hai, đầm Lăng Cô, đầm Thị Nại, đầm Cù Mông, vịnh Xuân Đài). Kết quả nghiên cứu là tài liệu tham khảo cốt lõi cho các chương trình đào tạo sau đại học ngành Khoa học Môi trường và Kỹ thuật Môi trường.
  • Chuyển đổi kinh tế - công nghiệp: Tái cấu trúc ngành nuôi trồng thủy sản tỉnh Thừa Thiên Huế theo hướng bền vững; giảm thiểu rủi ro dịch bệnh và hiện tượng cá chết hàng loạt do sốc thiếu oxy vào mùa khô, bảo vệ sinh kế cho hơn 300.000 cư dân sống phụ thuộc trực tiếp và gián tiếp vào nguồn lợi đầm phá.
  • Tác động chính sách: Trực tiếp đóng góp cơ sở dữ liệu số trị 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ế; định hình các quy định về phân vùng tiếp nhận nước thải và cấp phép xả thải lưu vực sông Hương - đầm Tam Giang.

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

  • Nghiên cứu sinh & Nhà khoa học: Tiếp cận bộ dữ liệu thủy văn - chất lượng nước độc chuẩn hóa, bộ thông số hiệu chỉnh mô hình Delft-3D (hệ số nhám Manning, hệ số khuếch tán, tốc độ phản ứng sinh hóa) đã được kiểm chứng thực tế tại thủy vực ven bờ nhiệt đới.
  • Cơ quan Quản lý Nhà nước (Bộ TN&MT, UBND tỉnh Thừa Thiên Huế): Sở hữu công cụ mô phỏng dự báo hỗ trợ ra quyết định (Decision Support System) để đánh giá tác động môi trường chiến lược (ĐMC) của các dự án công nghiệp, đô thị và du lịch ven biển đến năm 2030.
  • Doanh nghiệp & Hộ nuôi trồng thủy sản: Nhận được các khuyến cáo khoa học về lịch thời vụ xả đáy ao nuôi, mật độ thả giống tối ưu và giới hạn tải lượng thức ăn nhằm ngăn chặn hiện tượng phú dưỡng cục bộ.

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

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

Đó là việc mở rộng và chuẩn hóa lý thuyết cân bằng chất ô nhiễm Vladimirov (1986) cho thủy vực mở ven biển có trao đổi triều phức tạp, đồng thời thiết lập nguyên lý: Ngưỡng sức tải môi trường đầm phá phải được khống chế bởi giới hạn chịu đựng sinh thái của sinh vật thủy sinh nhạy cảm nhất chứ không thể neo theo các chỉ số tự phục hồi cơ học thuần túy của quy chuẩn xả thải thông thường.

2. Sự đổi mới về mặt phương pháp luận so với các nghiên cứu trước đây?

So với các nghiên cứu của Enhui et al. (2013) tại Vịnh Hạ Môn chỉ áp dụng mô hình 2D hộp đơn giản hoặc nghiên cứu của Nguyễn Hữu Cử (2006) chỉ dừng ở thống kê mô tả thực nghiệm, luận án đã tích hợp thành công mô hình động học 2D/3D liên hoàn Delft3D (FLOW + WAQ), hiệu chỉnh chính xác các tham số dòng chảy thực tế ($RSE = 0,0046$ và $RSE = 0,0042$; $R^2 \approx 0,70$) cho phép mô phỏng lan truyền chất ô nhiễm theo thời gian thực tương ứng với từng pha triều lên - xuống trong cả mùa khô và mùa mưa.

3. Phát hiện nào gây bất ngờ nhất dưới góc nhìn dữ liệu?

Phát hiện bất ngờ nhất là: Mùa khô (tháng 5) mới là giai đoạn đầm phá bị tổn thương sinh thái nghiêm trọng nhất chứ không phải mùa mưa lũ (tháng 11). Mặc dù mùa mưa tiếp nhận tổng tải lượng ô nhiễm từ lục địa lớn gấp nhiều lần, nhưng lưu lượng xả lũ cực lớn qua cửa Thuận An ($V_{max} = 105\text{ cm/s}$) giúp xả rửa đầm phá nhanh chóng. Ngược lại, mùa khô thiếu dòng chảy bổ cập từ sông, hiện tượng phân tầng và tù đọng nước tại Đầm Thủy Tú đã kéo tụt nồng độ oxy hòa tan xuống mức báo động ($4,7\text{ mg/L}$), biến khu vực này thành bẫy ô nhiễm hữu cơ nghẹt thở đối với ấu trùng thủy sinh.

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 mô tả chi tiết từng bước: từ tọa độ và kỹ thuật lấy mẫu theo TCVN, phương pháp phân tích phòng thí nghiệm, bảng thông số phát thải đơn vị (kg/người/ngày, kg/con gia súc/năm, kg/ha NTTS/vụ), hệ số rửa trôi bề mặt theo ngày mưa, cho đến bộ tham số thiết lập lưới tính, điều kiện biên triều và hệ số nhám trong phần mềm Delft-3D, cho phép các nhà khoa học khác tái lập hoàn toàn mô hình cho các hệ đầm phá khác.

5. Chương trình nghị sự nghiên cứu 10 năm tới được phác thảo ra sao?

Luận án định hình lộ trình nghiên cứu 10 năm tập trung vào: (1) Tích hợp mô hình AI/Machine Learning với Delft-3D để dự báo chất lượng nước cực ngắn (Nowcasting); (2) Xây dựng bản đồ sức tải môi trường số thích ứng với các kịch bản nước biển dâng do biến đổi khí hậu; (3) Mở rộng tính toán sức tải cho các hợp chất hữu cơ khó phân hủy (POPs, PAHs) và kim loại nặng trong trầm tích đáy.

Kết luận

  1. Kiểm kê định lượng toàn diện: Xác lập bức tranh tổng thể về nguồn thải và tải lượng các hợp chất C, N, P đưa vào hệ đầm phá Tam Giang - Cầu Hai từ các hoạt động sinh hoạt, chăn nuôi, công nghiệp, NTTS và rửa trôi đất tại mốc cơ sở 2011 và dự báo chính xác đến năm 2020, 2030.
  2. Làm chủ công nghệ mô hình hóa tiên tiến: Ứng dụng thành công bộ mô hình số trị Delft-3D (FLOW và WAQ), hiệu chỉnh đạt độ tin cậy khoa học cao với sai số $RSE < 0,005$ và $R^2 \approx 0,70$, mô phỏng chi tiết trường nồng độ DO, BOD₅, COD, $N-NH_4^+$, $N-NO_3^-$, $P-PO_4^{3-}$ theo các pha triều và mùa khí văn.
  3. Xác lập 3 kịch bản sức tải khoa học: Cung cấp định lượng sức tải môi trường theo ngưỡng Quy chuẩn Việt Nam, ngưỡng sức tải tối đa (đồng hóa tự nhiên), và ngưỡng sức tải đề xuất phát triển bền vững cho từng phân vùng không gian của hệ đầm phá.
  4. Đổi mới tư duy quản lý môi trường lưu vực: Chuyển dịch từ phương thức "kiểm soát nồng độ cuối đường ống" thụ động sang "quản lý tổng tải lượng xả thải theo sức tải tiếp nhận sinh thái" chủ động.
  5. Mở ra các hướng nghiên cứu liên ngành: Đặt nền móng cho các nghiên cứu chuyên sâu về sinh thái học phục hồi đầm phá ven biển, quản trị lưu vực sông tích hợp (IWRM), và lượng giá kinh tế sinh thái dịch vụ môi trường nước tại Việt Nam.