Tổng quan về luận án

Nghiên cứu của tác giả Mai Trọng Hoàng (2021) với đề tài "Nghiên cứu sức chịu tải môi trường của sông Trường Giang, tỉnh Quảng Nam làm cơ sở khoa học cho việc sử dụng hợp lý" thuộc chuyên ngành Khoa học Môi trường (Mã số: 9440301.01), được thực hiện dưới sự hướng dẫn khoa học của PGS. Trần Văn Thụy và TS. Nguyễn Đăng Giáp tại Trường Đại học Khoa học Tự nhiên, Đại học Quốc gia Hà Nội. Công trình mang tính tiên phong trong việc giải quyết bài toán suy thoái môi trường nước tại một thủy vực sông ven biển có chế độ thủy động lực học lưỡng cực vô cùng đặc thù. Sông Trường Giang có chiều dài 67 km, chảy qua 18 xã thuộc 4 huyện, thành phố (Duy Xuyên, Thăng Bình, Tam Kỳ, Núi Thành), kết nối hai hệ thống sông lớn là Vu Gia - Thu Bồn ở phía Bắc (thông qua Cửa Đại) và Tam Kỳ ở phía Nam (thông qua Cửa Lở và Cửa An Hòa).

Vấn đề cốt lõi (research gap) mà luận án giải quyết xuất phát từ việc các nghiên cứu trước đây tại miền Trung (như Nguyễn Văn Cư, 2009; Nguyễn Thị Thảo Hương, 2011; Phạm Thị Ngọc Lan, 2013) hầu như chỉ tập trung quan trắc các chỉ tiêu chất lượng nước đơn lẻ hoặc khảo sát địa mạo, trầm tích mà chưa lượng hóa được khả năng tự làm sạch và sức chịu tải môi trường (environmental carrying capacity) dưới tác động đồng thời của dòng chảy triều hai chiều phức tạp và áp lực nhân sinh mãnh liệt.

Luận án tập trung giải quyết 3 câu hỏi nghiên cứu (Research Questions - RQ) và kiểm định các giả thuyết khoa học (Hypotheses - H) tương ứng:

  • RQ1: Tải lượng các chất ô nhiễm chính ($BOD_5$, $NH_4^+$, $NO_3^-$, $TSS$) phát sinh từ các nguồn thải trong lưu vực biến đổi như thế nào giữa năm cơ sở 2017 và dự báo đến năm 2030? $\rightarrow$ H1: Tải lượng từ nuôi trồng thủy sản (NTTS) và sinh hoạt chiếm tỷ trọng vượt trội, tạo áp lực ô nhiễm cục bộ nghiêm trọng trong mùa kiệt.
  • RQ2: Chế độ thủy động lực học chịu ảnh hưởng của triều lưỡng cực tác động như thế nào đến quá trình lan truyền, tích tụ chất ô nhiễm và biến đổi cấu trúc thủy sinh vật? $\rightarrow$ H2: Sự xuất hiện của điểm "0" lưu tốc tại khu vực giữa sông (khoảng Km27) làm triệt tiêu khả năng tự làm sạch, gia tăng ô nhiễm hữu cơ và làm suy giảm chỉ số đa dạng sinh học Shannon-Weiner ($H'$).
  • RQ3: Giới hạn sức chịu tải của từng đoạn sông theo các kịch bản phát triển là bao nhiêu để làm căn cứ quy hoạch sử dụng đất và bảo vệ nguồn nước? $\rightarrow$ H3: Việc kết hợp mô hình hóa 1 chiều với khung phân tích DPSIR và ma trận SWOT cho phép xác định chính xác ngưỡng tiếp nhận nước thải tối đa và lộ trình quản lý không gian thích ứng.

Khung lý thuyết của nghiên cứu tích hợp lý thuyết thủy động lực Saint-Venant, định luật bảo toàn khối lượng và khuếch tán Fick, lý thuyết chỉ thị sinh học môi trường và mô hình tải lượng tối đa hàng ngày (TMDL của US EPA, 2012; Don Butcher et al., 2011). Luận án được triển khai trên quy mô không gian 67 km sông, khảo sát qua 6 đợt thực địa liên tục trong giai đoạn 2016–2018 (tập trung mùa cạn từ tháng 4–6 và mùa mưa từ tháng 10–12), bao quát khu vực sinh sống của 153,3 nghìn dân với mật độ 540 người/$km^2$ (gấp 3,2 lần mức trung bình tỉnh Quảng Nam) và 1.489,3 ha mặt nước NTTS.

Literature Review và Positioning

Nghiên cứu chất lượng nước và sức chịu tải thủy vực trên thế giới đã phát triển qua nhiều giai đoạn với các trường phái tiếp cận khác nhau:

  • Đánh giá theo chỉ số chất lượng nước tổng hợp: Bắt đầu từ công trình nền tảng của Horton (1965) và Brown et al. (1970) với chỉ số NSFWQI, mở rộng sang chỉ số CCMEWQI của Canada (Khan et al., 2003; Debels et al., 2005) và chỉ số OWQI tại Anh (Abbasi, 2002; Kannel et al., 2007). Nhược điểm của trường phái này là che lấp sự biến thiên tức thời theo chu kỳ triều và không phản ánh được quá trình chuyển hóa động học của chất ô nhiễm trong không gian dòng chảy.
  • Mô hình hóa chất lượng nước và thủy động lực học: Sự phát triển của các công cụ mô phỏng 1D và 2D/3D như Aquatox (Lei et al., 2010), BLTM (Feaster & Conrads, 1999; Broshears et al., 2001), QUAL2K (Cristea & Pelletier, 2005; Kannel et al., 2007; Yali Zhang, 2016), WASP7 (Craig A. Stow et al., 2003), Delft 3D (David Taylor et al., 2005; Roy Van Weerdenburg, 2017) và hệ thống phần mềm MIKE của Viện Thủy lực Đan Mạch (DHI) (Kamel, 2008; Ahmed, 2010).

Các tranh luận học thuật nổi bật tồn tại giữa hai quan điểm:

  1. Quan điểm thủy lực tĩnh/bảo toàn tuyệt đối: Cho rằng việc tính toán sức chịu tải chỉ cần dựa trên lưu lượng dòng kiệt thiết kế $Q_{95%}$ và nồng độ giới hạn theo quy chuẩn (tương tự công thức của Ying Xie et al., 2014; Shu Ling Huang et al., 2011).
  2. Quan điểm thủy động lực động và cân bằng sinh thái: Cho rằng sức chịu tải là đại lượng biến đổi động phụ thuộc phi tuyến vào chu kỳ tương tác triều - sông, hiện tượng khuếch tán rối dọc và phản ứng sinh hóa của các bậc dinh dưỡng thủy sinh vật (GESAMP, 1986; Don Butcher et al., 2011; Roberto et al., 2016).

Luận án định vị nghiên cứu vào khoảng trống chưa được giải quyết: thiết lập mô hình tính toán sức chịu tải toàn diện cho hệ thống sông hở hai đầu lưu thông với biển, nơi chịu tác động của chế độ thủy văn phi đối xứng (mùa lũ chiếm 62,6–68,2% tổng lượng nước cả năm, trong khi mùa kiệt chỉ chiếm 31,8–37,4%, tháng IV chỉ đạt 2,1–2,6%).

Khi so sánh với các nghiên cứu quốc tế:

  • So với nghiên cứu của Lei Zhu et al. (2017) trên lưu vực sông Vị Hà (Trung Quốc) hay Eric Akomeah et al. (2015) trên sông Saskatchewan (Canada), nghiên cứu tại sông Trường Giang phức tạp hơn do không có trạm quan trắc thủy văn thường trực và dòng chảy thay đổi theo chu kỳ bán nhật triều/nhật triều không đều, đảo chiều liên tục giữa hai cửa sông.
  • So với mô hình quản lý lưu vực sông Fraser (Kent R. Gustavson et al., 1999) và lưu vực sông Tarim (Feng et al., 2001), nghiên cứu này đã tích hợp thành công chỉ thị sinh học thực nghiệm (chỉ số Shannon-Weiner $H'$ và Berger-Parker $d$) với mô hình mô phỏng MIKE 11 để định lượng hóa khả năng phục hồi sinh thái của thủy vực.

Đó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 cân bằng tải lượng môi trường (Mass Balance Theory) và lý thuyết tự làm sạch thủy vực (Self-purification Dynamics) trong điều kiện thủy văn biến động mạnh:

  1. Làm sáng tỏ cơ chế bẫy ô nhiễm thủy lực (Hydraulic Pollutant Trapping): Chứng minh rằng tại các dòng sông hẹp ven biển nối hai cửa vịnh/biển, điểm gặp nhau của hai sóng triều (điểm triệt tiêu vận tốc dòng chảy, $v \approx 0\text{ m/s}$) hoạt động như một "vùng bẫy cơ học", triệt tiêu khả năng phát tán chất lơ lửng và chất hữu cơ, biến đoạn giữa sông thành điểm nóng ô nhiễm dù tải lượng xả thải trực tiếp tại chỗ không phải là cao nhất.
  2. Khung lý thuyết tích hợp Động lực học - Sinh thái (Coupled Hydrodynamic-Ecological Framework): Mở rộng học thuyết chỉ thị sinh vật của May (1975) và Berger & Parker (1970), chứng minh sự suy giảm đơn điệu của chỉ số đa dạng Shannon-Weiner ($H'$) đồng biến với nồng độ các chất dinh dưỡng ($NH_4^+$, $NO_3^-$) và $BOD_5$, xác lập mối quan hệ định lượng giữa cấu trúc quần xã thực vật nổi (TVN), động vật nổi (DVN), động vật đáy (DVD) với sức chịu tải môi trường nước.
+-------------------------------------------------------------------------+
|                  KHUNG PHÂN TÍCH TỔNG HỢP CỦA LUẬN ÁN                   |
+-------------------------------------------------------------------------+
| [ĐỘNG LỰC NHÂN SINH]                                                    |
| - Dân số (153.300 người, mật độ 540 ng/km²)                             |
| - NTTS (1.489,3 ha tôm thẻ, tôm sú)                                     |
| - 3 KCN tập trung (Bắc Chu Lai, Tam Hiệp, Chu Lai - Trường Hải)         |
| - Chăn nuôi & Rửa trôi đất (99,627 km² đất nông nghiệp)                 |
+------------------------------------+------------------------------------+
                                     |
                                     v
+------------------------------------+------------------------------------+
| [PHÂN TÍCH TẢI LƯỢNG Ô NHIỄM (2017 & DỰ BÁO 2030)]                      |
| - Nguồn điểm (Point Sources): Công nghiệp, Sinh hoạt tập trung          |
| - Nguồn diện (Non-Point Sources): Rửa trôi đất, Chăn nuôi, Nông nghiệp  |
| - Nguồn NTTS phân tán ven sông                                          |
+------------------------------------+------------------------------------+
                                     |
                                     v
+------------------------------------+------------------------------------+
| [MÔ HÌNH HÓA THỦY ĐỘNG LỰC & CHẤT LƯỢNG NƯỚC (MIKE 11)]                 |
| - MIKE 11 HD: Giải hệ phương trình Saint-Venant 1 chiều                 |
| - MIKE 11 AD: Lan truyền - Khuếch tán Fickian 1 chiều                   |
| - MIKE 11 Ecolab: Quá trình biến đổi sinh hóa BOD, NH4+, NO3-, TSS      |
| - Hiệu chuẩn & Kiểm định với dữ liệu trắc đạc 6 đợt thực địa            |
+------------------------------------+------------------------------------+
                                     |
                                     v
+------------------------------------+------------------------------------+
| [ĐÁNH GIÁ SỨC CHỊU TẢI & PHỤC HỒI SINH THÁI]                           |
| - Tính toán khả năng tiếp nhận (L_tn) theo Thông tư 76/2017/TT-BTNMT    |
| - Tích hợp chỉ số sinh học: Shannon-Weiner (H') & Berger-Parker (d)    |
| - Đánh giá kịch bản mùa lũ / mùa kiệt (2017, 2030: A2, B1, B2)          |
+------------------------------------+------------------------------------+
                                     |
                                     v
+------------------------------------+------------------------------------+
| [ĐỊNH HƯỚNG QUY HOẠCH & QUẢN LÝ THÍCH ỨNG]                             |
| - Khung phân tích nguyên nhân - hệ quả: DPSIR                           |
| - Phân tích chiến lược đa mục tiêu: SWOT                                |
| - Giải pháp công trình (nạo vét, khơi thông) & phi công trình           |
+-------------------------------------------------------------------------+

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

Khung phân tích của luận án kết hợp đa tầng phương pháp luận:

  • Tầng 1 - Xác định tải lượng cơ sở: Áp dụng hệ số phát thải đơn vị chuẩn hóa cho nguồn thải sinh hoạt ($G_i$), công nghiệp ($D_i$), chăn nuôi và NTTS kết hợp hệ số rửa trôi mặt đất ($H_s$) và diện tích sử dụng đất ($A_i$) dựa trên các phương thức quản lý đất đai.
  • Tầng 2 - Mô phỏng số trị tương tác: Sử dụng bộ phần mềm MIKE (MIKE 11 HD, AD, Ecolab) tích hợp hệ thống thông tin địa lý (GIS) và mô hình số độ cao (DEM) nhằm tái hiện chính xác biến trình triều, trường vận tốc dòng chảy, sự biến thiên nồng độ chất ô nhiễm theo cả không gian và thời gian.
  • Tầng 3 - Tích hợp quản lý chiến lược: Sử dụng khung DPSIR (Drivers - Pressures - State - Impact - Responses) kết hợp ma trận SWOT nhằm chuyển hóa các kết quả mô hình toán định lượng thành các khuyến nghị chính sách sử dụng hợp lý lãnh thổ và tài nguyên nước.

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

Thiết kế nghiên cứu

Luận án áp dụng quan điểm thực chứng (positivism) kết hợp hệ thống luận môi trường (environmental systems theory). Thiết kế nghiên cứu là sự phối hợp chặt chẽ giữa quan trắc thực địa quy chuẩn, thí nghiệm sinh thái và mô phỏng số học (numerical modeling).

  • Hệ thống dữ liệu không gian: Toàn bộ mạng lưới sông 67 km được số hóa thành các mặt cắt ướt hình học chi tiết, tích hợp địa hình vùng cồn cát ven biển và hệ thống công trình vượt sông, đập Cổ Linh, các ao tôm lấn chiếm lòng sông.
  • Hệ thống phân đoạn nghiên cứu: Chia tách sông thành các phân đoạn chức năng dựa trên ngã ba hợp lưu với sông Thu Bồn, sông Tam Kỳ và các cửa xả tập trung của 3 KCN.

Quy trình nghiên cứu rigorous

Quy trình thu thập và phân tích dữ liệu tuân thủ nghiêm ngặt các tiêu chuẩn quốc gia (TCVN) và quốc tế (ISO, Standard Methods for the Examination of Water and Wastewater):

  1. Khảo sát thủy văn và địa hình: Thực hiện 6 đợt đo đạc đồng bộ tại các trạm MC1, MC2 trên sông Trường Giang và hệ thống sông lân cận; đo đạc lưu lượng, mực nước, vận tốc dòng chảy và mặt cắt địa hình lòng sông.
  2. Lấy mẫu và phân tích chất lượng nước: Lấy mẫu nước mặt tầng mặt tại các vị trí đại diện trong cả mùa cạn và mùa mưa; phân tích các thông số vật lý, hóa học ($BOD_5$, $NH_4^+$, $NO_3^-$, $TSS$, $DO$, độ mặn) theo tiêu chuẩn QCVN 08-MT:2015/BTNMT.
  3. Phân tích thủy sinh vật: Thu mẫu định tính và định lượng TVN (Phytoplankton), DVN (Zooplankton), DVD (Zoobenthos) tại các điểm thu mẫu; định danh phân loại học và tính toán các chỉ số đa dạng sinh học sinh thái:
    • Chỉ số Shannon-Weiner: $$H' = - \sum_{i=1}^{S} P_i \ln P_i$$
    • Chỉ số ưu thế Berger-Parker: $$d = \frac{N_{\max}}{N}$$ (Trong đó $N_{\max}$ là số lượng cá thể của loài chiếm ưu thế tuyệt đối; $N$ là tổng số cá thể của tất cả các loài trong mẫu thu).
+-------------------------------------------------------------------------+
|                  QUY TRÌNH THỰC NGHIỆM VÀ MÔ PHỎNG SỐ                   |
+-------------------------------------------------------------------------+
| [Đo đạc thực địa 6 đợt] --> [Phân tích Lab: Lý - Hóa - Sinh]             |
|          |                                   |                          |
|          v                                   v                          |
| [Hình học MC & Thủy văn]   [Chỉ số ĐDSH: H', d] & [Nồng độ: BOD, N, TSS]|
|          |                                   |                          |
|          +-----------------+-----------------+                          |
|                            |                                            |
|                            v                                            |
|             [Thiết lập mô hình MIKE 11]                                 |
|             - Module HD: Hiệu chuẩn & Kiểm định                         |
|             - Module AD: Hệ số khuếch tán E                             |
|             - Module Ecolab: Động học phân hủy k                        |
|                            |                                            |
|                            v                                            |
|             [Mô phỏng kịch bản tải lượng & Sức chịu tải]                |
+-------------------------------------------------------------------------+

Data và phân tích

Cơ sở toán học của các module mô phỏng trong luận án bao gồm:

1. Thủy động lực học 1 chiều (Module MIKE 11 HD): Giải hệ phương trình Saint-Venant:

  • Phương trình liên tục: $$\frac{\partial Q}{\partial s} + B \frac{\partial z}{\partial t} = q$$
  • Phương trình động lượng: $$\frac{\partial Q}{\partial t} + \frac{\partial (\alpha \frac{Q^2}{A})}{\partial s} + g A \frac{\partial z}{\partial s} + g \frac{Q |Q|}{K^2} - U_x q = 0$$ (Trong đó $K = A C \sqrt{R}$ là môđun lưu lượng, $C$ là hệ số Chézy, $R$ là bán kính thủy lực, $U_x$ là thành phần lưu tốc của dòng gia nhập $q$, $\alpha$ là hệ số hiệu chỉnh động năng).

2. Lan truyền và khuếch tán chất ô nhiễm (Module MIKE 11 AD): Dựa trên phương trình vi phân chuyển tải - khuếch tán 1 chiều: $$\frac{\partial (A C)}{\partial t} + \frac{\partial (Q C)}{\partial s} - \frac{\partial}{\partial s}\left(A E \frac{\partial C}{\partial s}\right) = - A K_{deg} C + C_q q$$ (Trong đó $C$ là nồng độ chất ô nhiễm, $E$ là hệ số khuếch tán dọc, $K_{deg}$ là hằng số suy giảm sinh hóa bậc một).

3. Đánh giá sức chịu tải môi trường nước: Thực hiện theo Thông tư 76/2017/TT-BTNMT và Quyết định 154/QĐ-TCMT: $$L_{tn} = (L_{td} - L_{nn} - L_{t}) \times F_s$$ (Trong đó: $L_{tn}$ là khả năng tiếp nhận nước thải/sức chịu tải (kg/ngày); $L_{td}$ là tải lượng tối đa cho phép của đoạn sông ứng với mục tiêu chất lượng nước; $L_{nn}$ là tải lượng chất ô nhiễm hiện có trong nguồn nước tiếp nhận; $L_t$ là tải lượng chất ô nhiễm đưa vào từ nguồn xả thải; $F_s$ là hệ số an toàn được lựa chọn nghiêm ngặt trong khoảng 0,3–0,7).

Mô hình đã được hiệu chỉnh (calibration) và kiểm định (validation) độc lập bằng cách so sánh chuỗi số liệu mực nước thực đo và mô phỏng tại trạm Tam Kỳ cùng các trạm phụ trợ trên sông Trường Giang, cho thấy độ khớp pha triều cao và sai số biên độ nằm trong giới hạn kiểm định thủy văn cho phép.

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

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

  1. Phát hiện điểm triệt tiêu vận tốc dòng chảy ("Điểm 0 thủy lực"): Kết quả mô phỏng trường dòng chảy chỉ ra sự tồn tại của điểm dừng triều tại vùng giữa sông (khoảng Km27, đoạn qua Thăng Bình). Vận tốc dòng chảy tại đây tiệm cận 0 m/s trong nhiều thời điểm chuyển triều. Hiện tượng này làm thời gian lưu nước tăng cao đột biến, ngăn cản quá trình pha loãng và biến đoạn sông này thành bể tích lũy bùn lắng hữu cơ và trầm tích phù sa mịn.
  2. Cấu trúc phát thải bất cân xứng: Năm 2017, hoạt động NTTS (với 1.489,3 ha) và nước thải sinh hoạt từ 153.300 dân cư là nguồn phát tán tải lượng chất hữu cơ ($BOD_5$) và amoni ($NH_4^+$) lớn nhất vào sông. Đất nông nghiệp (chiếm 47,37% diện tích lưu vực, tương đương 99,627 $km^2$) đóng góp lượng lớn tổng chất rắn lơ lửng ($TSS$) và chất dinh dưỡng rửa trôi trong các đợt mưa lớn mùa lũ (tháng X–XII).
  3. Sự thắt hẹp cục bộ làm tê liệt động lực dòng chảy: Lòng sông bị xâm lấn nghiêm trọng bởi các ao tôm tự phát (đoạn Bình Nam - Tam Thanh và Tam Thanh - Tam Tiến) và đập Cổ Linh (xã Bình Sa), khiến bề rộng mặt thoáng tại nhiều điểm co hẹp chỉ còn 2,5 m, gây cản trở dòng chảy kiệt, nghẽn lối thoát lũ và làm gia tăng nồng độ ô nhiễm tích tụ tại chỗ.
  4. Quy luật biến thiên nghịch đảo của chỉ thị sinh thái: Phân tích sinh học chỉ ra tương quan tuyến tính chặt chẽ giữa mức độ phú dưỡng hữu cơ và chỉ số ĐDSH. Tại các đoạn sông có $BOD_5$ và $NH_4^+$ tăng cao (vượt ngưỡng QCVN 08-MT:2015/BTNMT cột A2/B1), chỉ số đa dạng loài Shannon-Weiner ($H'$) của thực vật nổi và động vật đáy giảm mạnh, trong khi chỉ số ưu thế Berger-Parker ($d$) tiến sát về giá trị 1, phản ánh sự suy thoái nghiêm trọng của các loài nhạy cảm và sự thống trị của một vài loài chống chịu ô nhiễm cao.
  5. Dự báo sức chịu tải đến năm 2030: Dưới các kịch bản phát triển kinh tế (A2, B1, B2), nếu tỷ lệ lấp đầy các KCN (Bắc Chu Lai đạt 96,1%, Tam Hiệp 60%, Chu Lai - Trường Hải >80%) gia tăng cùng với việc tăng đàn gia súc (định hướng 800.000 trâu bò, 10 triệu gia cầm toàn tỉnh năm 2030) mà không có biện pháp xử lý triệt để, sức chịu tải đối với $BOD_5$ và $NH_4^+$ trên phần lớn các đoạn sông Trường Giang trong mùa cạn sẽ chuyển sang trạng thái âm (quá tải nghiêm trọng).
+-------------------------------------------------------------------------+
|                  BIẾN THIÊN NỒNG ĐỘ VÀ SỨC CHỊU TẢI SÔNG                |
+-------------------------------------------------------------------------+
| Đầu Bắc (Cửa Đại)          Giữa sông (Km27)         Đầu Nam (Cửa An Hòa)|
| [Triều Cửa Đại dồn vào] -> [ĐIỂM 0 LƯU TỐC] <- [Triều An Hòa đẩy lên]   |
| - Độ mặn: 1,3‰ -> 0,8‰     - Bồi lắng cực đại  - Độ mặn: 6,7‰ -> 2,0‰   |
| - Vận tốc: Dao động triều  - v ≈ 0 m/s         - Vận tốc: Dao động triều|
| - Sức chịu tải: Còn dư     - Quá tải BOD, NH4+ - Sức chịu tải: Tùy đoạn |
| - Chỉ số ĐDSH: H' cao      - H' thấp, d cao    - Chỉ số ĐDSH: H' khá    |
+-------------------------------------------------------------------------+

Implications đa chiều

  • Về mặt lý luận khoa học: Đặt nền móng cho việc ứng dụng mô hình động lực học liên kết để xác định sức chịu tải các hệ thống sông hở ven biển đặc thù chịu tác động triều hai chiều tại vùng duyên hải miền Trung Việt Nam.
  • Về mặt phương pháp luận: Thiết lập quy trình mẫu tích hợp giữa tính toán tải lượng phát thải từ sử dụng đất, mô phỏng số trị 1D (MIKE 11) và giải đoán biến động quần xã thủy sinh vật chỉ thị, có thể chuyển giao áp dụng cho các sông có điều kiện tương tự (như sông Cửa Lớn, Cà Mau).
  • Về mặt ứng dụng thực tiễn: Cung cấp bức tranh phân vùng tiếp nhận nước thải chi tiết cho từng đoạn sông, giúp các cơ quan quản lý biết chính xác đoạn sông nào còn dư địa tiếp nhận và đoạn nào bắt buộc phải ngừng cấp phép xả thải mới.
  • Về mặt chính sách quản lý: Khung DPSIR của luận án chỉ ra các giải pháp can thiệp có tính chiến lược:
    1. Giải pháp công trình: Nạo vét thông luồng toàn tuyến sông Trường Giang, phá bỏ các đăng đáy và công trình cản trở dòng chảy tại đập Cổ Linh để tăng cường lưu thông nước giữa hai cửa sông Cửa Đại và An Hòa, loại bỏ điểm 0 lưu tốc.
    2. Giải pháp phi công trình: Chuyển đổi mô hình nuôi tôm trên cát tự phát sang nuôi công nghệ cao có hệ sinh thái tuần hoàn nước; kiểm soát 100% nước thải công nghiệp từ KCN Bắc Chu Lai, Tam Hiệp và KCN Chu Lai - Trường Hải đạt chuẩn trước khi xả ra môi trường.

Limitations và Future Research

Luận án thẳng thắn thừa nhận các giới hạn khoa học mang tính khách quan:

  1. Giới hạn về không gian mô hình hóa 1 chiều: Việc sử dụng mô hình 1 chiều (MIKE 11) phản ánh rất tốt sự truyền triều và biến đổi nồng độ dọc trục sông chính, nhưng chưa mô phỏng chi tiết được sự phân tầng độ mặn và gradient nồng độ chất ô nhiễm theo phương ngang/thẳng đứng tại các khu vực mở rộng ngã ba sông và cửa biển.
  2. Giới hạn số liệu chuỗi thời gian thủy văn: Do sông Trường Giang không có trạm quan trắc thủy văn thường trực của mạng lưới quốc gia, số liệu hiệu chỉnh mô hình phụ thuộc chủ yếu vào các đợt đo đạc đồng bộ ngắn hạn (6 đợt trong các năm 2016–2018).
  3. Cơ chế phát thải nguồn phân tán: Tải lượng từ nguồn nông nghiệp và chăn nuôi chủ yếu được ước tính thông qua hệ số phát thải và số liệu niên giám thống kê, chưa bao quát hết sự biến động bất thường của các trận mưa cực đoan do biến đổi khí hậu.

Các hướng nghiên cứu tiếp theo được đề xuất:

  • Phát triển mô hình 2 chiều/3 chiều (như MIKE 21/3 hoặc Delft-3D) cho các vùng cửa sông nhạy cảm (Cửa Đại, Cửa An Hòa, Cửa Lở) để đánh giá chi tiết quá trình nêm mặn và tương tác bùn cát 3 chiều.
  • Nghiên cứu định lượng tác động của kịch bản biến đổi khí hậu toàn cầu và nước biển dâng kết hợp với quy hoạch vận hành xả lũ của các hồ chứa thủy điện bậc thang trên lưu vực sông Vu Gia - Thu Bồn đến sức chịu tải sông Trường Giang.
  • Ứng dụng công nghệ viễn thám và GIS phân giải cao để quan trắc tự động diện tích nuôi tôm và độ đục của nước sông theo thời gian thực.

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

  • Tác động học thuật: Đóng góp nguồn dữ liệu cơ sở toàn diện đầu tiên về thủy văn, môi trường và sinh thái cho sông Trường Giang; cung cấp tài liệu tham khảo chuẩn mực cho các nghiên cứu về động lực học sông ven biển và sức chịu tải sinh thái.
  • Chuyển đổi kinh tế - xã hội địa phương: Cung cấp luận cứ khoa học trực tiếp cho Ủy ban nhân dân tỉnh Quảng Nam trong việc tái cơ cấu ngành nuôi trồng thủy sản ven sông, hỗ trợ chuyển đổi phương thức sinh kế bền vững cho hàng vạn hộ dân thuộc 18 xã ven sông Trường Giang.
  • Tác động chính sách và quy hoạch: Đóng vai trò là cơ sở khoa học then chốt phục vụ công tác lập quy hoạch tỉnh Quảng Nam giai đoạn 2021–2030, tầm nhìn 2050; định hướng quy hoạch không gian thoát lũ, dự án nạo vét chỉnh trị sông Trường Giang và bảo vệ vùng đệm sinh thái rừng ngập mặn Cửa Lở - Cửa An Hòa.

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

  • Nghiên cứu sinh và Giảng viên ngành Môi trường/Thủy văn: Tiếp cận phương pháp luận mẫu về mô hình hóa tích hợp đa thông số (HD-AD-Ecolab-ĐDSH) và bộ dữ liệu thực địa trắc đạc phong phú.
  • Cơ quan Quản lý Nhà nước (Sở TN&MT, Sở NN&PTNT tỉnh Quảng Nam): Sử dụng các bản đồ phân vùng sức chịu tải để phê duyệt báo cáo đánh giá tác động môi trường (ĐTM), cấp phép xả thải và kiểm soát các nguồn xả vào sông.
  • Doanh nghiệp & Ban Quản lý các KCN (Bắc Chu Lai, Tam Hiệp, Chu Lai - Trường Hải): Xác định ngưỡng xả thải tối đa cho phép để đầu tư nâng cấp hệ thống xử lý nước thải tập trung đạt tiêu chuẩn kỹ thuật.
  • Cộng đồng cư dân địa phương: Hưởng lợi trực tiếp từ môi trường nước sạch hơn, giảm thiểu nguy cơ dịch bệnh trong nuôi trồng thủy sản và giảm thiểu rủi ro ngập lụt trong mùa bão lũ.

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ì?
Luận án đã mở rộng lý thuyết thủy động lực môi trường ven biển khi xác lập cơ chế "Bẫy ô nhiễm thủy lực" tại điểm dừng triều lưỡng cực (điểm "0" lưu tốc xấp xỉ Km27). Luận án chứng minh rằng tại các dòng sông nối hai cửa biển, hiện tượng suy thoái chất lượng nước nghiêm trọng nhất trong mùa kiệt không nhất thiết nằm ở cửa xả lớn nhất, mà nằm tại vùng hội tụ sóng triều nơi vận tốc dòng chảy bị triệt tiêu, làm mất khả năng khuếch tán và tự làm sạch tự nhiên.

2. Điểm đổi mới phương pháp luận so với các nghiên cứu trước đây là gì?
So với các nghiên cứu đánh giá đơn lẻ (chỉ dùng chỉ số WQI như Phạm Ngọc Hồ et al., 2009; Lê Văn Thăng et al., 2012) hoặc chỉ mô phỏng thủy lực thuần túy (Nguyễn Văn Tỉnh et al., 2015), luận án đã thiết lập quy trình giải tích 3 tầng khép kín: Lượng hóa tải lượng từ bản đồ sử dụng đất $\rightarrow$ Mô phỏng động học 1D bằng bộ công cụ MIKE 11 (HD-AD-Ecolab) $\rightarrow$ Đối sánh và kiểm chứng sinh thái bằng các chỉ số ĐDSH thực nghiệm ($H'$, $d$).

3. Phát hiện nào bất ngờ nhất từ kết quả số liệu thực nghiệm?
Phát hiện bất ngờ nhất là mặc dù lưu lượng tổng cộng của sông rất nhỏ trong mùa cạn (tháng IV chỉ chiếm 2,1–2,6% lượng nước cả năm), nhưng độ mặn ở hai đầu sông lại biến thiên dị thường: đầu Bắc (Duy Nghĩa) chỉ đạt 1,3‰ và giảm nhanh về 0,8‰ (do nước ngọt sông Thu Bồn đẩy ra), trong khi đầu Nam (Tam Tiến) bị mặn xâm nhập sâu tới 6,7‰ và lan tới Bình Sa (2‰), tạo ra sự phân hóa sinh thái thủy sinh vật sâu sắc giữa hai nửa Bắc - Nam của dòng sông.

4. Nghiên cứu có cung cấp quy trình tái lặp (replication protocol) không?
Luận án cung cấp đầy đủ thông số hình học mặt cắt, điều kiện biên thủy văn, vị trí tọa độ các trạm đo, nồng độ phát thải đơn vị cho từng loại hình sản xuất và quy trình thiết lập các module MIKE 11. Các nhà khoa học khác có thể áp dụng toàn bộ giao thức này cho bất kỳ dòng sông ven biển hở hai đầu nào tại khu vực nhiệt đới gió mùa.

5. Chương trình nghị sự nghiên cứu 10 năm tới được định hình ra sao?
Tập trung vào 3 hướng chính: (1) Tích hợp mô hình động lực học trầm tích 2D/3D để đánh giá hiện tượng xói lở bờ sông và bồi lấp luồng tàu sau các dự án nạo vét; (2) Đánh giá biến động sức chịu tải dưới tác động kép của biến đổi khí hậu và chuyển đổi sử dụng đất đô thị hóa dọc tuyến đường ven biển; (3) Thiết lập hệ thống giám sát tự động online kết hợp mô hình dự báo cảnh báo sớm chất lượng nước thời gian thực.

Kết luận

Luận án Tiến sĩ của Mai Trọng Hoàng (2021) là một công trình nghiên cứu công phu, chuẩn mực và có đóng góp học thuật đột phá cho chuyên ngành Khoa học Môi trường tại Việt Nam, thể hiện qua các điểm cốt lõi:

  1. Lượng hóa toàn diện tải lượng ô nhiễm: Xác định chính xác lượng phát thải $BOD_5$, $NH_4^+$, $NO_3^-$, $TSS$ từ các nguồn sinh hoạt, công nghiệp, NTTS, chăn nuôi và rửa trôi đất trên lưu vực sông Trường Giang giai đoạn 2017 và dự báo đến năm 2030.
  2. Làm rõ bản chất thủy động lực học lưỡng cực: Tái hiện thành công biến trình dòng chảy triều hai chiều và xác định điểm nghẽn lưu tốc làm gia tăng tích tụ ô nhiễm tại khu vực giữa sông.
  3. Xác lập giới hạn sức chịu tải sinh thái: Tính toán khả năng tiếp nhận nước thải chi tiết theo từng đoạn sông cho mùa lũ và mùa kiệt theo các kịch bản phát triển kinh tế - xã hội.
  4. Chứng minh quy luật chỉ thị sinh thái học: Lượng hóa mối quan hệ mật thiết giữa nồng độ các chất ô nhiễm hữu cơ/dinh dưỡng với sự suy giảm chỉ số đa dạng Shannon-Weiner ($H'$) và sự gia tăng chỉ số ưu thế Berger-Parker ($d$).
  5. Đề xuất hệ thống giải pháp quản lý khả thi: Chuyển hóa toàn bộ kết quả mô hình toán và đánh giá sinh thái thành khung hành động DPSIR và chiến lược SWOT, phục vụ trực tiếp cho quy hoạch phát triển bền vững và sử dụng hợp lý nguồn tài nguyên sông Trường Giang, tỉnh Quảng Nam.