Tổng quan về luận án

Luận án tiến sĩ chuyên ngành Khoa học Môi trường với đề tài "Nghiên cứu phân vùng chất lượng nước vùng biển ven bờ Hải Phòng bằng chỉ số chất lượng nước" do nghiên cứu sinh Lê Văn Nam thực hiện (dưới sự hướng dẫn khoa học của TS. Lê Xuân Sinh và GS.TS Đặng Thị Kim Chi tại Học viện Khoa học và Công nghệ thuộc Viện Hàn lâm Khoa học và Công nghệ Việt Nam, 2021) là một công trình tiên phong trong lĩnh vực hải dương học môi trường và quản lý tài nguyên vùng bờ. Đặt trong bối cảnh dải ven biển Hải Phòng giữ vị trí chiến lược – là cửa ngõ giao thương hàng hải quốc tế của miền Bắc, nằm tại vùng chuyển tiếp giữa đới bờ châu thổ sông Hồng và đới bờ karst Đông Bắc – khu vực này đang đối mặt với xung đột gay gắt giữa phát triển kinh tế biển (cảng biển, công nghiệp hóa chất Đình Vũ, du lịch Cát Bà - Đồ Sơn, nuôi trồng thủy sản) và bảo tồn các hệ sinh thái đặc hữu.

Khoảng trống nghiên cứu (Research gap) được xác định rất rõ ràng: Phần lớn các mô hình tính toán chỉ số chất lượng nước (Water Quality Index - WQI) trên thế giới và tại Việt Nam đều được tối ưu hóa cho môi trường nước mặt lục địa dạng sông, hồ (như hướng dẫn VN_WQI theo Quyết định số 1460/QĐ-TCMT ngày 12/11/2019 của Tổng cục Môi trường). Đối với môi trường cửa sông ven biển – nơi chịu sự chi phối phức tạp của chế độ nhật triều, sự phân tầng nước, gradient độ mặn và quá trình trao đổi bùn cát – các hướng dẫn hiện hành vẫn chưa có bộ chỉ số chuẩn hóa phù hợp. Các nghiên cứu trước đây tại Việt Nam thường thiếu vắng sự kết hợp giữa mô hình thủy động lực học, thống kê đa biến và kỹ thuật GIS để lượng hóa bức tranh phân vùng chất lượng nước biển đa chiều.

Để giải quyết khoảng trống này, luận án thiết lập hệ thống câu hỏi nghiên cứu và giả thuyết khoa học chặt chẽ:

  • Câu hỏi nghiên cứu 1 (RQ1): Làm thế nào để xây dựng một bộ chỉ số WQI đặc thù cho vùng biển ven bờ ($WQI_{HP}$) có khả năng phản ánh chính xác hiện trạng ô nhiễm mà không bị rơi vào hiện tượng "che khuất" (eclipsing) hoặc "tính ảo" (ambiguity)?
  • Câu hỏi nghiên cứu 2 (RQ2): Quy luật phân bố không gian và biến động theo mùa/pha triều của các hợp phần ô nhiễm tại 4 tiểu vùng sinh thái trọng điểm ven bờ Hải Phòng diễn ra như thế nào dưới tác động của tải lượng phát thải lục địa?
  • Câu hỏi nghiên cứu 3 (RQ3): Phân vùng chất lượng nước bằng công cụ $WQI_{HP}$ kết hợp GIS có thể cung cấp luận cứ khoa học chuẩn xác nào phục vụ quy hoạch không gian biển và bảo vệ hệ sinh thái?
  • Giả thuyết nghiên cứu 1 (H1): Việc phân nhóm thông số dựa trên phân tích thành phần chính (PCA) và hiệu chỉnh hàm toán học tổng hợp dạng tích - tổng có trọng số sẽ loại trừ được sai số che khuất nồng độ cục bộ của các thông số độc hại/nhạy cảm.
  • Giả thuyết nghiên cứu 2 (H2): Phân vùng chất lượng nước biển theo WQI phản ánh mối tương quan chặt chẽ với động lực dòng chảy triều và tải lượng ô nhiễm từ lưu vực sông Bạch Đằng, phân định rõ rệt giữa đới cửa sông châu thổ phía Tây Nam và đới cửa sông hình phễu phía Đông Bắc.

Khung lý thuyết của công trình được xây dựng dựa trên sự giao thoa của ba trụ cột: Lý thuyết Chỉ số Môi trường (Ott, 1978; Horton, 1965), Lý thuyết Quản lý biển theo không gian (Marine Spatial Planning - IOC/UNESCO, 2009) và Lý thuyết Tiếp cận hệ thống trong khoa học Trái đất (Von Bertalanffy). Về quy mô, nghiên cứu bao phủ toàn bộ vùng triều và biển ven bờ Hải Phòng ra đến quần đảo Cát Bà - Long Châu, tập trung vào 4 vùng trọng điểm: Tây Nam Đồ Sơn, Đông Bắc Đồ Sơn (cửa Nam Triệu), ven biển Cát Bà và vùng ngoài. Dữ liệu thực nghiệm bao gồm chuỗi quan trắc liên tục giai đoạn 2009–2019 và các đợt đo đạc chi tiết năm 2019 theo cả mùa mưa, mùa khô, tầng mặt, tầng đáy, pha nước lớn và pha nước ròng.

Literature Review và Positioning

Nghiên cứu về chỉ số chất lượng nước khởi nguồn từ công trình kinh điển của Horton (1965), người đầu tiên lượng hóa chất lượng nước thành một đại lượng vô thứ nguyên trên thang số. Tiếp đó, Brown et al. (1970) phát triển bộ chỉ số NSF-WQI cho Quỹ Vệ sinh Quốc gia Hoa Kỳ, đặt nền móng cho các hàm tổng hợp dạng tổng cộng hoặc nhân có trọng số. Ott (1978), Dunnette (1979) và Tebbutt (2002) đã hoàn thiện hệ thống tiêu chí tuyển chọn thông số đặc trưng. Tại Canada, Hội đồng Bộ trưởng Môi trường (CCME, 2001) đã ban hành chỉ số CCME-WQI dựa trên ba yếu tố: phạm vi (scope - $F_1$), tần suất (frequency - $F_2$) và biên độ (amplitude - $F_3$) của các thông số vượt ngưỡng quy chuẩn. Tại Châu Á, Bhargava (1983) tại Ấn Độ thiết lập WQI theo từng mục đích sử dụng nước cụ thể cho lưu vực sông Ganga; Liou et al. (2004) tại Đài Loan đề xuất phương pháp kết hợp hàm tổng và tích theo nhóm thông số nhằm giảm thiểu sai lệch toán học.

Trong y văn tồn tại hai luồng quan điểm đối nghịch sâu sắc về cấu trúc hàm WQI:

  1. Trường phái hàm tổng (Additive Formulations): Ưu điểm là tính toán trực quan, dễ áp dụng (như các biến thể của NSF-WQI), nhưng bộc lộ nhược điểm chí mạng là "tính che khuất" (eclipsing) – hiện tượng một thông số bị ô nhiễm cực nặng (như kim loại nặng hoặc coliform) bị bù trừ và che lấp bởi điểm số cao của các thông số khác, dẫn đến đánh giá sai lệch rằng nguồn nước vẫn an toàn.
  2. Trường phái hàm tích (Multiplicative Formulations): Khắc phục được tính che khuất khi một chỉ số phụ tiến về 0 sẽ kéo toàn bộ chỉ số tổng về 0, nhưng lại vấp phải "tính ảo" (ambiguity) và sự quá nhạy cảm khi có nhiều thông số tiệm cận ngưỡng dưới khiến chỉ số tổng bị sụt giảm quá mức so với thực tế.

Tại Việt Nam, các nghiên cứu của Lê Trình (2008, 2010), Nguyễn Lê Tú Quỳnh (2015), Tôn Thất Lãng (2014) và Phạm Ngọc Hồ (2011, 2015) đã nỗ lực bản địa hóa WQI nhưng hầu hết tập trung vào sông nội địa. Đối với môi trường biển, Cục Kiểm soát ô nhiễm (2010) đề xuất Sea Water Quality Index (SWQI) với 6 thông số (TSS, COD, $NH_4^+$, dầu mỡ, Pb, coliform) dạng trung bình cộng đơn giản, không có trọng số nên chưa phản ánh được mức độ độc tính sai biệt giữa các chất. Nghiên cứu của Nguyễn Thị Thế Nguyên (2014) tại Vịnh Hạ Long ($WQI_{HL}$) sử dụng hàm tích có trọng số cố định nhưng gặp khó khăn khi môi trường có sự biến động phức tạp và thiếu khả năng mở rộng thông số. Nghiên cứu của Trương Văn Đàn và nnk (2018) tại đầm phá Tam Giang - Cầu Hai hay Trần Thị Yên và Nguyễn Thị Thanh Thùy (2019) tại Quảng Bình vẫn dựa trên trọng số gán định tính, chưa gắn với phân tích cấu trúc dữ liệu thực tế.

Luận án này định vị học thuật bằng cách kết hợp phân tích thành phần chính (PCA) trên chuỗi dữ liệu thực nghiệm 10 năm để chọn lọc khách quan 14 thông số đại diện, khắc phục hoàn toàn nhược điểm của phương pháp gán trọng số chủ quan kiểu Delphi truyền thống. Khi so sánh với các nghiên cứu quốc tế như mô hình đánh giá vịnh Maryland của Carruthers & Wazniak (2004) tại Hoa Kỳ hay phân vùng vịnh Gulluk của Yucel-Gier et al. (2013) tại Thổ Nhĩ Kỳ, luận án thể hiện bước tiến vượt bậc khi lồng ghép mô hình thủy động lực lan truyền chất ô nhiễm vào cấu trúc không gian của WQI, tạo ra công cụ quản lý thích ứng cao cho vùng ven biển nhiệt đới gió mùa.

Đó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à thách thức các lý thuyết chỉ số chất lượng nước truyền thống của Horton (1965) và Brown et al. (1970) thông qua việc phát triển công thức $WQI_{HP}$ đặc thù cho vùng biển ven bờ. Mô hình lý thuyết đề xuất đã giải quyết thành công "nghịch lý che khuất" (eclipsing paradox) bằng cách cấu trúc hóa quá trình tính toán theo 2 bậc: bậc 1 tính toán chỉ số phụ theo từng nhóm thông số chức năng bằng hàm nhân/hàm lũy thừa, bậc 2 tích hợp các nhóm thông số qua hàm tổng có trọng số chuẩn hóa vô thứ nguyên.

Bằng chứng toán học từ luận án chứng minh: Trong trường hợp giả định một chỉ số phụ cực trị đạt mức ô nhiễm nghiêm trọng ($q_i \in [0, 1]$) trong khi các chỉ số khác đạt trạng thái tối ưu ($q_j = 100$), công thức $WQI_{HP}$ cải tiến phản ánh ngay lập tức mức sụt giảm chất lượng nước tương ứng với cấp độ nguy hại, loại bỏ hoàn toàn tính ảo tưởng an toàn mà các hàm trung bình cộng mắc phải.

Mô hình Chỉ số Cấu trúc Toán học Ưu điểm cốt lõi Nhược điểm/Khoảng trống Khắc phục trong $WQI_{HP}$
Horton (1965) Tổng gia quyền tuyến tính Đơn giản, dễ tính toán Bị che khuất (Eclipsing) rất nặng Phân nhóm và dùng hàm phi tuyến
NSF-WQI (1970) Tích/Tổng có trọng số Đánh giá đa mục tiêu Trọng số Delphi chủ quan Xác lập trọng số qua ma trận PCA
CCME (2001) Vector 3 chiều ($F_1, F_2, F_3$) Không giới hạn thông số Không có trọng số mức độ độc Gán trọng số thứ nguyên hóa
SWQI (2010) Trung bình cộng 6 thông số Dễ áp dụng ven bờ Thiếu tính đại diện, không trọng số Mở rộng 14 thông số, gán trọng số
$WQI_{HP}$ (Luận án) Tích - Tổng 2 bậc phân nhóm Triệt tiêu che khuất, trọng số khách quan Đòi hỏi bộ dữ liệu đầu vào đầy đủ Tích hợp mô phỏng thủy động lực

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 3 lý thuyết nền tảng:

  1. Lý thuyết Hóa học Môi trường biển: Phân loại các thông số thành 5 nhóm chỉ thị cốt lõi: Nhóm thông số trạng thái hóa lý cơ bản (pH, DO, TSS, độ mặn); Nhóm chất hữu cơ tiêu hao oxy (BOD5, COD); Nhóm muối dinh dưỡng và phú dưỡng (Amoni $NH_4^+$, Nitrat $NO_3^-$, Phosphat $PO_4^{3-}$, Chlorophyll-a); Nhóm độc chất và kim loại nặng (Dầu mỡ khoáng, Fe, As, Cd, Pb...); Nhóm vi sinh vật chỉ thị (Coliform).
  2. Lý thuyết Thống kê Đa biến (PCA - Pearson, 1901): Sử dụng ma trận tương quan và ma trận hiệp phương sai để trích xuất các thành phần chính ($PC_j$), xác định hệ số tải ($b_{ij}$) và tính toán tính tương đồng ($S$) nhằm giảm chiều dữ liệu mà không làm mất thông tin gốc.
  3. Lý thuyết Quy hoạch Không gian Biển (MSP - IOC/UNESCO, 2009): Chuyển dịch từ quản lý đơn ngành sang quản lý không gian tích hợp, thiết lập thang phân loại 5 cấp tương ứng với 5 dải màu quy chuẩn: Rất tốt (Xanh dương đậm), Tốt (Xanh lá), Trung bình (Vàng), Kém (Cam), Rất kém (Đỏ).

Điều kiện biên (Boundary conditions) của khung phân tích được xác định rõ ràng: Áp dụng tối ưu cho vùng biển nông ven bờ, vịnh bán khép kín và vùng cửa sông chịu ảnh hưởng của chế độ nhật triều biên độ lớn (từ 2,5 m đến 4,0 m), có sự phân hóa mùa khô - mùa mưa rõ rệt dưới tác động của hoàn lưu cực đới và gió mùa Đông Nam Á.

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

Thiết kế nghiên cứu

Luận án tuân thủ triết lý nghiên cứu Thực chứng luận (Positivism) kết hợp với quan điểm Hiện thực phê phán (Critical Realism), xem xét chất lượng nước vừa là các thực thể hóa lý có thể đo lường định lượng chính xác, vừa là kết quả tương tác phức tạp giữa các quy luật tự nhiên (thủy hải văn) và cấu trúc kinh tế - xã hội nhân sinh (nguồn thải lưu vực).

Thiết kế nghiên cứu đa tầng (Multi-level mixed design) được triển khai qua 4 cấp độ:

  • Cấp độ 1 (Lưu vực và Nguồn thải): Kiểm kê, tính toán tải lượng ô nhiễm phát sinh và đưa vào vùng cửa sông Bạch Đằng từ 5 nguồn: Sinh hoạt, Công nghiệp, Chăn nuôi, Thủy sản và Rửa trôi đất.
  • Cấp độ 2 (Thủy động lực học): Thiết lập mô hình số trị mô phỏng trường dòng chảy triều, quá trình lan truyền chất ô nhiễm tầng mặt và tầng đáy.
  • Cấp độ 3 (Khảo sát thực địa và Phòng thí nghiệm): Mạng lưới quan trắc gồm hàng chục trạm chuẩn hóa, đo đạc đa thông số theo 2 pha triều (nước lớn - nước ròng) và 2 mùa (mưa - khô).
  • Cấp độ 4 (Không gian GIS): Số hóa bản đồ nền, nội suy không gian và biên tập bản đồ phân vùng chất lượng nước theo thang WQI chuẩn hóa.

Quy trình nghiên cứu rigorous

Quy trình thu thập và xử lý mẫu tuân thủ nghiêm ngặt các tiêu chuẩn quốc gia và quốc tế:

  • Thiết bị thu mẫu: Sử dụng thiết bị lấy mẫu nước chuyên dụng Niskin Van Dorn Sampler dung tích 2–5 lít, cho phép thu mẫu chính xác tại các tầng nước mặt (cách mặt nước 0,5 m) và tầng đáy (cách đáy 1,0 m).
  • Giao thức bảo quản mẫu: Mẫu phân tích $BOD_5$, COD, dinh dưỡng ($NH_4^+, NO_3^-, PO_4^{3-}$), kim loại nặng và coliform được bảo quản lạnh ở $4^\circ\text{C}$, cố định hóa học ngay tại hiện trường bằng axit chuyên dụng ($HNO_3, H_2SO_4$) theo đúng TCVN và Standard Methods for the Examination of Water and Wastewater (APHA).
  • Triangulation (Tam giác đạc): Kết hợp chặt chẽ giữa 4 nguồn dữ liệu: (1) Số liệu quan trắc thực địa tại trạm Đồ Sơn và các trạm khảo sát; (2) Chuỗi dữ liệu lịch sử 2009–2019 từ Trạm Quan trắc Phân tích Môi trường biển Miền Bắc; (3) Kết quả chạy mô hình mô phỏng lan truyền chất ô nhiễm; (4) Số liệu kiểm kê tải lượng phát thải từ Sở TN&MT và các đề tài cấp Nhà nước (KC.09).
  • Kiểm soát chất lượng (QA/QC): Độ tin cậy của phương pháp phân tích trong phòng thí nghiệm được kiểm soát qua mẫu trắng (blank), mẫu lặp (duplicate) và mẫu chuẩn thêm (spike) với độ thu hồi đạt 85–115%, hệ số biến thiên $CV < 5%$.

Data và phân tích

Phân tích định lượng dựa trên bộ dữ liệu tải lượng và chất lượng nước chi tiết. Về tải lượng phát thải, luận án dẫn chứng các con số cụ thể: Năm 2012, tổng tải lượng ô nhiễm phát sinh hàng năm trong khu vực cửa sông Bạch Đằng là 86.547,35 tấn COD; 47.443,25 tấn $BOD_5$; 19.012,25 tấn N-T; 3.648,31 tấn P-T và 248.779,50 tấn TSS. Trong đó, lượng chất thải thực tế đưa vào vùng cửa sông Bạch Đằng tiếp nhận là 45.077,20 tấn COD; 8.086,80 tấn $BOD_5$; 13.747,40 tấn N-T; 3.120,20 tấn P-T và 105.801,40 tấn TSS. Nguồn sinh hoạt đóng vai trò chủ đạo chiếm 41,31%, chăn nuôi 30,27%, công nghiệp 13,50%, rửa trôi đất 12,65% và nuôi trồng thủy sản 2,27%.

Dự báo đến năm 2025, tổng tải lượng phát sinh tăng vọt lên 175.682,06 tấn COD (tăng hơn 100%) và 302.571,80 tấn TSS, tạo sức ép cực lớn lên môi trường tiếp nhận ven bờ.

Tải lượng phát sinh ô nhiễm khu vực cửa sông Bạch Đằng (tấn/năm):

Về mặt thống kê toán học, luận án thực hiện phân tích PCA trên phần mềm thống kê chuyên dụng. Các biến số được chuẩn hóa theo độ lệch chuẩn: $$x_1' = \frac{x_{1i} - \bar{x}1}{S{x1}}$$ Ma trận tương quan và các giá trị riêng ($A_j$), vector riêng ($e_{ij}$) được trích xuất để tính toán ma trận tải nhân tố $a_{ij} = e_{ij} \sqrt{A_j}$. Qua đó, xác định trọng số khách quan ($w_i$) cho từng nhóm thông số, đảm bảo tổng các trọng số $\sum w_i = 1$. Đồng thời, luận án áp dụng phương pháp Đánh giá Rủi ro Môi trường thông qua thương số rủi ro ($RQ = \frac{C_{thực tế}}{C_{quy chuẩn}}$) để kiểm chứng chéo mức độ nhạy cảm của các điểm quan trắc.

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

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

Nghiên cứu mang lại 4 phát hiện khoa học mang tính bước ngoặt về trạng thái môi trường biển ven bờ Hải Phòng:

  1. Sự phân hóa không gian chất lượng nước sâu sắc theo cấu trúc địa mạo bờ biển: Vùng bờ Tây Nam Đồ Sơn (cửa sông châu thổ Văn Úc - Thái Bình) chịu áp lực nặng nề của độ đục và tổng chất rắn lơ lửng (TSS thường xuyên vượt giới hạn cho phép do phù sa bồi tụ), trong khi vùng bờ Đông Bắc Đồ Sơn (cửa sông hình phễu Bạch Đằng - Nam Triệu) lại chịu áp lực chủ đạo từ chất hữu cơ ($COD, BOD_5$) và dinh dưỡng khoáng phát sinh từ nước thải đô thị và KCN Đình Vũ.
  2. Quy luật chi phối của động lực triều và mùa khí hậu đến trường phân bố WQI: Giá trị $WQI_{HP}$ đạt mức thấp nhất (ô nhiễm nặng nhất) tại tầng mặt trong pha nước ròng vào mùa mưa, khi lưu lượng dòng chảy lục địa mang theo tải lượng chất bẩn cực đại tràn ra biển. Ngược lại, vào pha nước lớn mùa khô, khối nước biển mặn từ ngoài khơi lấn sâu vào vùng cửa sông giúp quá trình xáo trộn và tự làm sạch diễn ra mạnh mẽ, nâng chỉ số WQI lên mức Tốt và Rất tốt.
  3. Nghịch lý ô nhiễm cục bộ tại vùng đệm di sản Quần đảo Cát Bà: Mặc dù vùng biển mở quanh quần đảo Cát Bà có chất lượng nước đạt mức Rất tốt ($WQI > 85$), nhưng tại các vụng kín, tùng áng (vùng nước karst đặc thù với 59 tùng, áng) và các khu vực tập trung bè nuôi thủy sản (Bến Bèo, Gia Luận), đã xuất hiện hiện tượng suy giảm DO cục bộ và gia tăng nồng độ Phosphat ($PO_4^{3-}$), Chlorophyll-a và Coliform, đe dọa trực tiếp đến 154 loài san hô và 538 loài động vật đáy quý hiếm.
  4. Hiệu năng vượt trội của công thức $WQI_{HP}$ trong việc triệt tiêu sai số toán học: Khi kiểm nghiệm thực nghiệm so sánh với công thức CCME của Canada và NSF-WQI của Mỹ, công thức $WQI_{HP}$ cho thấy độ nhạy cao hơn 28% trong việc phát hiện các điểm nóng ô nhiễm vi sinh và dầu mỡ khoáng tại luồng hàng hải Hải Phòng, hoàn toàn không bị ảnh hưởng bởi tính che khuất khi nồng độ kim loại nặng vượt ngưỡng quy chuẩn QCVN 10-MT:2015/BTNMT.
                    BIẾN ĐỘNG WQI THEO ĐỘNG LỰC TRIỀU VÀ MÙA
        Mùa Mưa / NR   Mùa Mưa / NL      Mùa Khô / NR      Mùa Khô / NL
        (Kém - Ô nhiễm) (Trung bình)      (Tốt)             (Rất Tốt)

Implications đa chiều

  • Về mặt Lý thuyết: Công trình đã bổ sung vào kho tàng khoa học môi trường một phương pháp luận hoàn chỉnh về thiết lập chỉ số WQI cho vùng biển ven bờ nhiệt đới chịu tác động triều, kết nối thành công giữa lý thuyết phân tích thống kê đa biến và lý thuyết thủy văn cửa sông.
  • Về mặt Phương pháp luận: Cung cấp quy trình 6 bước chuẩn hóa có khả năng chuyển giao và nhân rộng cho các vùng biển ven bờ khác tại Việt Nam như Quảng Ninh, Đà Nẵng, Bà Rịa - Vũng Tàu hay Kiên Giang.
  • Về mặt Thực tiễn và Chính sách:
    • Đóng vai trò là công cụ hỗ trợ ra quyết định (Decision Support System) trực quan cho UBND thành phố Hải Phòng trong việc lập Quy hoạch không gian biển (Marine Spatial Planning) đến năm 2030, tầm nhìn 2045.
    • Cung cấp cơ sở khoa học để xác lập ranh giới các bãi đổ bùn nạo vét luồng hàng hải, ngăn chặn bùn thải tràn vào khu bảo tồn biển Cát Bà.
    • Đưa ra lộ trình kiểm soát tải lượng phát thải từ đất liền, đặc biệt là đề xuất xây dựng trạm xử lý nước thải tập trung cho lưu vực sông Cấm, Lạch Tray trước khi đổ ra vịnh Bắc Bộ.

Limitations và Future Research

Nghiên cứu thẳng thắn thừa nhận 4 giới hạn mang tính khách quan:

  1. Độ phân giải thời gian của dữ liệu: Mặc dù đã kế thừa chuỗi số liệu 10 năm, nhưng các đợt khảo sát chi tiết chuyên sâu mới chỉ tập trung vào các tháng đại diện cho mùa khô và mùa mưa của năm 2019, chưa bao quát được các biến cố thời tiết cực đoan như siêu bão hoặc lũ lụt lịch sử.
  2. Cơ chế động lực trầm tích đáy: Luận án tập trung chủ yếu vào pha nước (nước mặt và nước đáy), chưa tích hợp đầy đủ quá trình trao đổi chất ô nhiễm tại ranh giới nước - trầm tích (benthic-pelagic coupling) và quá trình hấp phụ - giải hấp kim loại nặng của hạt keo phù sa.
  3. Các chất ô nhiễm mới nổi (Emerging Contaminants): Bộ chỉ số chưa bao gồm các thông số vi nhựa (microplastics), dư lượng dược phẩm và các chất gây rối loạn nội tiết (EDCs) – những yếu tố đang gia tăng nhanh chóng do hoạt động du lịch và đô thị hóa.

Chương trình nghiên cứu 10 năm tới (Future Research Agenda) được đề xuất:

  • Hướng 1: Tích hợp mạng lưới cảm biến IoT quan trắc tự động liên tục với thuật toán Trí tuệ Nhân tạo (Machine Learning/Deep Learning) để dự báo chỉ số $WQI_{HP}$ thời gian thực (Real-time WQI forecasting).
  • Hướng 2: Ứng dụng công nghệ viễn thám không gian (Satellite Remote Sensing) độ phân giải cao kết hợp WQI để giám sát diện rộng nồng độ Chlorophyll-a, TSS và nhiệt độ mặt biển ($SST$).
  • Hướng 3: Mở rộng nghiên cứu số phận sinh thái học độc chất (Ecotoxicological fate modeling) của các hợp chất hữu cơ khó phân hủy (POPs, PAHs) trong chuỗi thức ăn sinh vật biển từ sinh vật phù du đến cá đáy và thú biển tại Khu Dự trữ Sinh quyển Quần đảo Cát Bà.

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

Công trình luận án mang lại những tác động sâu rộng trên nhiều bình diện:

  • Tác động Học thuật: Được kỳ vọng trở thành tài liệu tham khảo chuẩn mực với tiềm năng trích dẫn cao trong các nghiên cứu về quản lý tổng hợp vùng bờ (ICZM), hải dương học môi trường và kỹ thuật môi trường biển tại khu vực Đông Nam Á.
  • Chuyển dịch Ngành và Công nghiệp: Tạo nền tảng để các doanh nghiệp cảng biển, khu công nghiệp ven biển (Deep C, Nam Đình Vũ) kiểm soát báo cáo xả thải định kỳ, giúp ngành nuôi trồng thủy sản Hải Phòng định vị an toàn các vùng nuôi lồng bè tuân thủ tiêu chuẩn VietGAP/ASC.
  • Tác động Chính sách cấp Trung ương và Địa phương: Là tài liệu khoa học nền tảng phục vụ thực thi Luật Tài nguyên, môi trường biển và hải đảo năm 2015, Nghị quyết 36-NQ/TW về Chiến lược phát triển bền vững kinh tế biển Việt Nam đến năm 2030, tầm nhìn 2045.
  • Bảo tồn Xã hội và Quốc tế: Đóng góp trực tiếp vào hồ sơ khoa học bảo vệ nguyên vẹn giá trị đa dạng sinh học của Khu Dự trữ Sinh quyển Thế giới Quần đảo Cát Bà và liên kết bảo tồn Di sản Thiên nhiên Thế giới Vịnh Hạ Long - Quần đảo Cát Bà được UNESCO công nhận.

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

  1. Nghiên cứu sinh và Giới học thuật: Tiếp cận phương pháp luận xây dựng chỉ số môi trường phức hợp, quy trình xử lý dữ liệu PCA đa biến và kỹ thuật lồng ghép mô hình toán học vào GIS.
  2. Cơ quan Quản lý Nhà nước: Tổng cục Biển và Hải đảo Việt Nam, Cục Kiểm soát ô nhiễm, Sở Tài nguyên và Môi trường Hải Phòng sở hữu bộ công cụ và bản đồ phân vùng chuẩn xác để cấp phép xả thải, quy hoạch vùng nuôi thủy sản và thiết lập khu bảo tồn biển.
  3. Doanh nghiệp Du lịch và Khai thác Thủy sản: Được cung cấp thông tin minh bạch, định lượng về chất lượng nước để lựa chọn vùng nuôi trồng an toàn sinh học và phát triển du lịch bãi biển bền vững.
  4. Cộng đồng Dân cư Ven biển: Hưởng lợi gián tiếp từ môi trường biển trong sạch, nguồn lợi thủy hải sản được phục hồi và sinh kế kinh tế biển được bảo đảm lâu dài.

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?

Đóng góp độc đáo nhất là việc tái cấu trúc hàm toán học $WQI_{HP}$ theo mô hình phân nhóm tích hợp 2 bậc kết hợp trọng số khách quan trích xuất từ ma trận hiệp phương sai PCA. Công trình đã mở rộng trực tiếp Lý thuyết Chỉ số Môi trường của Horton (1965) và Brown et al. (1970) sang môi trường biển nhiệt đới gió mùa, giải quyết triệt để sự xung đột giữa "tính che khuất" và "tính ảo" vốn làm tê liệt các mô hình WQI truyền thống.

2. Sự đổi mới về mặt phương pháp luận khi so sánh với ít nhất 2 nghiên cứu quốc tế?

So với mô hình WQI vịnh Maryland của Carruthers & Wazniak (2004, Hoa Kỳ) vốn chỉ dùng phép gán nhị phân đơn giản ($0$ và $1$) rồi lấy trung bình cộng làm mất mát độ nhạy cảm của nồng độ thực tế, và so với mô hình CCME (2001, Canada) không có cơ chế gán trọng số mức độ độc hại cho các thông số, phương pháp luận của luận án vượt trội nhờ: (1) Chuẩn hóa chỉ số phụ qua đường cong phi tuyến bám sát quy chuẩn QCVN 10-MT:2015/BTNMT; (2) Tích hợp mô phỏng thủy động lực lan truyền chất ô nhiễm để nội suy không gian thay vì chỉ dựa vào phép nội suy khoảng cách nghịch đảo (IDW) thuần túy.

3. Phát hiện bất ngờ nhất (counter-intuitive) có bằng chứng dữ liệu hỗ trợ là gì?

Phát hiện bất ngờ nhất là nồng độ một số kim loại nặng và chất dinh dưỡng tại khu vực cửa phễu Nam Triệu trong pha nước lớn không giảm tuyến tính theo khoảng cách ra biển mà lại tích tụ cục bộ tại vùng ranh giới nước lợ - mặn (độ mặn 15–20‰). Bằng chứng mô phỏng thủy động lực và phân tích hóa học xác nhận đây là vùng "bẫy bùn cát" (turbidity maximum) do quá trình keo tụ (flocculation) khi gặp môi trường điện ly cao, khiến các chất ô nhiễm bị giữ lại thay vì phát tán ra biển khơi.

4. Luận án có cung cấp giao thức tái lập (replication protocol) đầy đủ không?

Luận án cung cấp chi tiết toàn bộ quy trình tái lập trong 6 phụ lục: từ tọa độ trạm quan trắc, bảng thông số kỹ thuật thiết bị Niskin Van Dorn Sampler, công thức giải thuật PCA, các bước tính toán chỉ số phụ $q_i$, bảng ngưỡng phân loại $WQI_{HP}$ đến hướng dẫn thiết lập hệ tọa độ VN-2000 và biên tập lớp dữ liệu không gian trong GIS.

5. Khung chương trình nghiên cứu 10 năm tới được phác thảo như thế nào?

Luận án định hình lộ trình nghiên cứu 2021–2031 tập trung vào 3 trọng tâm: (1) Số hóa và tự động hóa hệ thống cảnh báo sớm WQI tích hợp AI; (2) Nghiên cứu cơ chế lan truyền và rủi ro sinh thái của vi nhựa và ô nhiễm dầu ngầm; (3) Thiết lập mô hình toán kinh tế sinh thái (Ecological-economic modeling) nhằm định giá dịch vụ hệ sinh thái nước biển ven bờ Hải Phòng.

Kết luận

Luận án tiến sĩ của NCS. Lê Văn Nam đã hoàn thành xuất sắc các mục tiêu đề ra với 6 đóng góp mang tính đột phá cụ thể:

  1. Hệ thống hóa toàn diện và nâng tầm phương pháp luận về chỉ số chất lượng nước áp dụng cho các vùng biển ven bờ, làm sáng tỏ các tương tác giữa quá trình thủy văn - hóa học - sinh thái vùng cửa sông nhiệt đới.
  2. Xây dựng thành công công thức chỉ số chất lượng nước đặc thù ($WQI_{HP}$) với cấu trúc toán học tối ưu, giải quyết triệt để hiện tượng che khuất nồng độ và thích ứng hoàn hảo với điều kiện tự nhiên của bờ biển Hải Phòng.
  3. Xác lập bộ trọng số khách quan cho các nhóm thông số chất lượng nước thông qua kỹ thuật phân tích thành phần chính (PCA) trên chuỗi dữ liệu thực nghiệm 10 năm, loại bỏ tính chủ quan của phương pháp chuyên gia truyền thống.
  4. Lượng hóa chi tiết tải lượng ô nhiễm phát sinh và đưa vào vùng cửa sông Bạch Đằng (với 45.077,20 tấn COD và 105.801,40 tấn TSS/năm năm 2012; dự báo tăng lên 175.682,06 tấn COD/năm vào 2025), chỉ rõ vai trò chi phối của nguồn thải sinh hoạt (41,31%) và công nghiệp.
  5. Thành lập bộ bản đồ phân vùng chất lượng nước chuẩn hóa theo 5 mức độ ô nhiễm cho toàn bộ dải ven bờ Hải Phòng ở các kịch bản: mùa khô, mùa mưa, nước lớn, nước ròng, tầng mặt và tầng đáy, bảo đảm tính khoa học và thẩm mỹ cao.
  6. Đề xuất hệ thống giải pháp đồng bộ từ quản lý quy hoạch không gian biển (MSP), kiểm soát phát thải lục địa liên tỉnh đến các giải pháp kỹ thuật công nghệ xử lý nước thải và bảo vệ 52.981 ha hệ sinh thái đáy mềm ven bờ.

Công trình đánh dấu bước chuyển dịch paradigm quan trọng trong quản lý môi trường biển tại Việt Nam – từ phương pháp quan trắc nồng độ đơn lẻ, thụ động sang phương pháp đánh giá chỉ số tích hợp, trực quan hóa không gian đa chiều. Luận án mở ra các nhánh nghiên cứu chuyên sâu mới về mô hình hóa chất lượng nước thời gian thực, quản trị biển thích ứng với biến đổi khí hậu và tạo lập di sản học thuật chuẩn mực, phục vụ sự phát triển kinh tế biển bền vững của thành phố Hải Phòng và cả nước.