Tổng quan về luận án

Nghiên cứu về các hợp chất hữu cơ thơm đa vòng (Polycyclic Aromatic Hydrocarbons - PAHs) trong các hệ sinh thái ven biển nhiệt đới đóng vai trò then chốt trong khoa học môi trường hiện đại. PAHs là nhóm ô nhiễm hữu cơ khó phân hủy (Persistent Organic Pollutants - POPs) cấu tạo từ hai hay nhiều vòng benzen ngưng tụ, sở hữu đặc tính kỵ nước cao, khó phân hủy sinh học và có khả năng tích lũy sinh học mạnh mẽ qua chuỗi thức ăn. Tổ chức Y tế Thế giới (WHO) cùng Cơ quan Nghiên cứu Ung thư Quốc tế (IARC) và Cơ quan Bảo vệ Môi trường Hoa Kỳ (US EPA) đã xếp PAHs vào nhóm tác nhân hàng đầu gây ung thư, đột biến gen và quái thai ở sinh vật lẫn con người.

Rừng ngập mặn (RNM) xã Đồng Rui, huyện Tiên Yên, tỉnh Quảng Ninh có tổng diện tích tự nhiên 4.974,21 ha, sở hữu hệ sinh thái rừng ngập mặn độc đáo phía Bắc Việt Nam với diện tích rừng gần 2.800 ha, nằm trong tọa độ địa lý 21°10’ - 21º16’30’’ vĩ độ Bắc và 107°21’30’’ - 107º27’ kinh độ Đông. Khu vực này nằm kẹp giữa các cửa sông lớn (sông Ba Chẽ, sông Voi Lớn, sông Voi Bé) và vịnh Tiên Yên, chịu tác động trực tiếp từ hoạt động kinh tế - xã hội cường độ cao: khai thác than tại Cẩm Phả - Cửa Ông, hoạt động của các nhà máy nhiệt điện lớn (Nhiệt điện Mông Dương I và II cách 70 km) và giao thông thủy.

Khoảng trống nghiên cứu (Research Gap): Các công trình nghiên cứu PAHs tại Việt Nam trước đây (Nguyễn Thúy Ngọc et al., 2003; Nghiêm Trung Dũng, 2005; Dương Thanh Nghị et al., 2011; Phạm Thị Kha, 2014; Đặng Hoài Nhơn et al., 2014) chủ yếu tập trung vào không khí đô thị, nước mặt hoặc trầm tích cửa sông ven biển với quy mô mẫu hạn chế (<200 mẫu trên toàn quốc). Động thái tồn lưu của 16 hợp chất PAHs ưu tiên trong môi trường đất rừng ngập mặn theo không gian, thời gian, tầng sâu thổ nhưỡng và lượng hóa nguy cơ sinh thái - sức khỏe kết hợp mô hình hóa đa khoang chưa từng được thực hiện đồng bộ tại Việt Nam.

Luận án xác lập 3 câu hỏi nghiên cứu và 3 giả thuyết khoa học:

  • RQ1: Mức độ tồn lưu, phân bố theo chiều sâu thổ nhưỡng và biến thiên theo mùa của 16 PAHs trong đất RNM Đồng Rui diễn ra như thế nào?
    • H1: Nồng độ $\Sigma 16\text{PAHs}$ tích lũy cao hơn vào mùa đông do hướng gió mùa Đông Bắc vận chuyển khí thải từ các nguồn đốt công nghiệp và nồng độ biến thiên theo tầng sâu phản ánh lịch sử lắng đọng.
  • RQ2: Nguồn phát thải chiếm ưu thế gây ô nhiễm PAHs tại khu vực bắt nguồn từ đâu và mức độ rủi ro sinh thái cũng như rủi ro ung thư đối với cộng đồng dân cư đạt ngưỡng nào?
    • H2: Ô nhiễm PAHs bắt nguồn chủ yếu từ quá trình nhiệt phân nhiên liệu hóa thạch (đốt than, nhiệt điện) và chỉ số rủi ro ung thư (CR) ở nhóm cư dân nhạy cảm vượt ngưỡng khuyến nghị an toàn.
  • RQ3: Xu thế phân bố liên pha và tốc độ tích lũy dài hạn của hợp chất Benzo(a)pyrene (BaP) trong đất RNM diễn biến ra sao dưới tác động của cơ chế vận chuyển đa khoang?
    • H3: Mô hình Fugacity cấp III và cấp IV mô phỏng chính xác đất và trầm tích là bồn chứa (sink) tích lũy chính của BaP với xu hướng gia tăng theo thời gian nếu không có biện pháp giảm thiểu tải lượng đầu vào.

Phạm vi nghiên cứu được thực hiện đồng bộ qua 6 đợt quan trắc từ tháng 11/2013 đến tháng 11/2017 (các đợt phân tích chuyên sâu: 8/2014, 1/2015, 7/2015, 1/2016, 7/2016, 1/2017) tại 12 điểm lấy mẫu đất đại diện (ĐR1 - ĐR12), kết hợp mẫu không khí, nước, trầm tích và mẫu đất theo 4 tầng sâu (0-5 cm, 5-10 cm, 10-15 cm, 15-20 cm).

Literature Review và Positioning

Các nghiên cứu tồn lưu PAHs trên thế giới đã chứng minh tính phổ biến và phức tạp của nhóm chất này trong đất. J. Wolska (2002) ghi nhận hàm lượng PAHs dọc hệ thống sông Odra (Ba Lan) dao động từ 1 $\mu\text{g/kg}$ đến 50.000 $\mu\text{g/kg}$. Peng Huang et al. (2013) khảo sát 112 mẫu đất vùng biển Bột Hải (Trung Quốc), xác định nồng độ $\Sigma 16\text{PAHs}$ đạt 97,2 - 300,7 $\mu\text{g/kg}$ (trung bình 175,7 $\pm$ 37,3 $\mu\text{g/kg}$).

Trong hệ sinh thái rừng ngập mặn quốc tế, S. Sarkar et al. (2010) tại Sundarban (Vịnh Bengal, Ấn Độ) ghi nhận tổng nồng độ 16 PAHs dao động từ 132 đến 2.938 $\mu\text{g/kg}$ (trung bình 634 $\mu\text{g/kg}$), trong đó 5 hợp chất gây ung thư cao (BbF, BkF, BaP, Ind, DahA) chiếm tới 68 - 73% tổng lượng. Cavalcante et al. (2009) đánh giá RNM Fortaleza (Brazil), chỉ ra tổng nồng độ 17 PAHs đạt 3,04 - 2.234,76 $\mu\text{g/kg}$ tại sông Cocó và 3,34 - 1.859,21 $\mu\text{g/kg}$ tại sông Ceará. Tam et al. khảo sát các RNM Hồng Kông (Ting Kok, Sheung Pak Nai, Yi O, Ma Wan) ghi nhận dải nồng độ 56 - 3.758 $\mu\text{g/kg}$, với mức ô nhiễm cao nhất ở Ma Wan do phát thải công nghiệp lịch sử giai đoạn 1958 - 1979.

Tranh luận học thuật trung tâm: Tồn tại xung đột giữa hai trường phái lý thuyết về vai trò của chất hữu cơ trong đất (SOM/TOC) đối với khả năng lưu giữ PAHs. Trường phái truyền thống cho rằng hàm lượng cacbon hữu cơ tổng số (TOC) là yếu tố quyết định kiểm soát quá trình hấp phụ kỵ nước của PAHs. Ngược lại, trường phái phi cân bằng (Wan-Li Ma et al., 2013; Bucheli et al., 2004; Zhang et al., 2006; Xiang Bi et al., 2016 với $r^2 = 0,001$) chứng minh không có tương quan thuận rõ rệt giữa TOC và tổng PAHs, khẳng định sự tích lũy PAHs phụ thuộc vào đặc tính than đen (black carbon), tốc độ lắng đọng khí quyển và động lực học bồi tụ trầm tích.

Luận án định vị nghiên cứu vào khoảng trống chưa có dữ liệu tại RNM Đồng Rui - vùng đất ngập nước chịu chế độ nhật triều thuần nhất với biên độ triều cao nhất Việt Nam (3,5 - 4,0 m). Nghiên cứu giải quyết bài toán định lượng cơ chế phát tán liên pha khí - nước - đất - trầm tích bằng việc ứng dụng tích hợp mô hình phân bố Fugacity cấp III và mô phỏng tích lũy theo thời gian Fugacity cấp IV của Donald Mackay (1979, 2001).

Đó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 và động lực học Fugacity của Donald Mackay trong môi trường bãi triều ven biển nhiệt đới gió mùa:

  • Chuẩn hóa thông số dung lượng khuếch tán ($Z_i$), tải lượng đối lưu ($D_{ai}$), tải lượng chuyển hóa/phân hủy ($D_{ri}$) và hệ số vận chuyển liên pha ($D_{ij}$) đặc thù cho điều kiện RNM Đồng Rui.
  • Chứng minh giới hạn của mô hình hấp phụ cân bằng tuyến tính TOC - PAH truyền thống trong đất ngập mặn chịu chế độ thủy triều quét và bồi tụ động.
  • Hoàn thiện mô hình tính toán rủi ro ung thư (Cancer Risk - CR) tích hợp hệ số độc tương đương (Toxic Equivalent Factor - TEF) quy đổi về Benzo(a)pyrene độc tính tương đương ($TEQ = \sum C_i \times TEF_i$) cho 3 con đường phơi nhiễm (nuốt vô tình, hấp thụ qua da, hít thở bụi).

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

Khung phân tích tích hợp 3 trụ cột khoa học:

  1. Lý thuyết Chẩn đoán Tỷ lệ Đồng phân Phân tử (Molecular Diagnostic Ratios): Sử dụng các tỷ lệ đồng phân chọn lọc $Flt/(Flt + Pyr)$, $Ant/(Ant + Phe)$, $BaA/(BaA + Chr)$, $Ind/(Ind + BghiP)$ (Yunker et al., 2002; Pies et al., 2008) để phân biệt nguồn gốc dầu mỏ nguyên khai ($petrogenic$) và nguồn nhiệt phân nhiên liệu hóa thạch ($pyrogenic$).
  2. Khung Đánh giá Rủi ro Sinh thái và Sức khỏe Con người (ERA Framework): Ứng dụng thương số rủi ro $RQ = C_{thực_tế} / C_{tiêu_chuẩn}$ kết hợp mô hình CR đa đường phơi nhiễm theo tiêu chuẩn US EPA.
  3. Mô hình Động lực học Đa khoang Fugacity Cấp III và Cấp IV: Xây dựng hệ phương trình cân bằng khối lượng vi phân phi cân bằng mô phỏng biến thiên nồng độ BaP theo thời gian $d(V_i Z_i f_i)/dt = I_i + \sum (D_{ji} f_j) - f_i (D_{ri} + \sum D_{ij})$.

Điều kiện biên xác lập: Áp dụng cho vùng bãi triều cửa sông ven biển nhiệt đới gió mùa, chịu chế độ nhật triều biên độ 3,5 - 4,0 m, có sự tham gia của thảm thực vật ngập mặn (Kandelia obovata, Avicennia marina).

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

Thiết kế nghiên cứu

Thiết kế nghiên cứu thực nghiệm định lượng kết hợp mô hình hóa giải tích trên nền tảng thực chứng (positivism). Mạng lưới quan trắc gồm 12 điểm lấy mẫu đất bề mặt (0 - 10 cm) mã hóa từ ĐR1 đến ĐR12 phân bố đại diện cho các vùng cửa sông Ba Chẽ, sông Voi Lớn, bãi triều ngoài và vùng đệm rừng ngập mặn; 1 vị trí khoan lõi tầng sâu (ĐR5) chia 4 phân đoạn (0-5 cm, 5-10 cm, 10-15 cm, 15-20 cm); 2 điểm quan trắc không khí (pha khí và pha hạt bụi); 4 điểm mẫu nước và 4 điểm trầm tích mặt đáy.

Quy trình nghiên cứu thực nghiệm

  • Xử lý và phân tích mẫu: Mẫu đất và trầm tích được làm khô tự nhiên, loại bỏ tạp chất, nghiền mịn qua rây 0,25 mm, chiết xuất siêu âm bằng hỗn hợp dung môi diclorometan/n-hexan (tỷ lệ 1:1), làm sạch qua cột sắc ký silica gel/alumina hoạt hóa.
  • Thiết bị phân tích: Hệ thống sắc ký khí ghép khối phổ (GC/MS) vận hành ở chế độ chọn ion giám sát (SIM). Kiểm soát chất lượng (QA/QC) nghiêm ngặt với mẫu trắng, mẫu lặp và chất chuẩn nội thế mang ($d_{10}$-phenanthrene, $d_{12}$-chrysene, $d_{12}$-perylene), độ thu hồi đạt từ 78,5% đến 104,2%, đường chuẩn sắc ký đạt hệ số tương quan $R^2 > 0,998$.
  • Phân tích hóa lý phụ trợ: Đo pH đất bằng máy đo pH điện cực thủy tinh; xác định cacbon hữu cơ tổng số (TOC) bằng phương pháp Walkley-Black và máy phân tích TOC chuyên dụng.
Mạng lưới quan trắc hiện trường (12 điểm đất, tầng sâu 0-20cm, khí, nước, trầm tích)
   │
   ├──> Chiết siêu âm Soxhlet / Làm sạch cột Silica-Alumina ──> Phân tích GC/MS (16 PAHs)
   │
   ├──> Xác định đặc tính hóa lý thổ nhưỡng (TOC, pH) ───────> Phân tích tương quan Pearson
   │
   ├──> Tỷ lệ đồng phân chẩn đoán (Flt/Pyr, Ant/Phe, BaA/Chr, Ind/BghiP) ──> Định danh nguồn thải
   │
   ├──> Đánh giá rủi ro sinh thái & sức khỏe ─────────────────> Tính toán chỉ số RQ & CR
   │
   └──> Mô hình hóa đa khoang (Mackay) ─────────────────────> Fugacity Cấp III & Cấp IV

Xử lý dữ liệu và mô hình hóa

Sử dụng phân tích tương quan Pearson ($r$) để đánh giá mối liên hệ giữa $\Sigma 16\text{PAHs}$ với hàm lượng TOC và pH đất. Tính toán các thông số mô hình Fugacity cho 4 khoang chính: Không khí ($V_1 = 1,49 \times 10^{11}\text{ m}^3$), Nước ($V_2 = 1,24 \times 10^8\text{ m}^3$), Đất ($V_3 = 4,97 \times 10^6\text{ m}^3$), Trầm tích ($V_4 = 4,97 \times 10^5\text{ m}^3$). Giải hệ phương trình đại số tuyến tính mô hình cấp III và phương trình vi phân mô hình cấp IV bằng ngôn ngữ lập trình và công cụ tính toán số học chuyên dụng.

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

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

  1. Mức độ tích lũy và biến động thời gian: Tổng nồng độ 16 PAHs trong đất RNM Đồng Rui dao động từ mức thấp đến trung bình so với các RNM trên thế giới, nhưng ghi nhận sự biến thiên theo mùa rõ rệt. Nồng độ trung bình $\Sigma 16\text{PAHs}$ vào mùa đông (đợt tháng 1/2015, tháng 1/2016, tháng 1/2017) cao hơn đáng kể so với mùa hè (tháng 8/2014, tháng 7/2015, tháng 7/2016). Điều này khẳng định hướng gió thịnh hành mùa đông (Bắc - Đông Bắc) đóng vai trò vận chuyển các hạt bụi mịn hấp phụ PAHs từ các khu công nghiệp phía Bắc và Trung Quốc lắng đọng xuống khu vực.
  2. Chiếm ưu thế của PAHs phân tử lượng cao (HMW): Nhóm PAHs có từ 4 đến 6 vòng benzen (Flt, Pyr, BaA, Chr, BbF, BkF, BaP, Ind, BghiP, DahA) chiếm tỷ lệ vượt trội (từ 62,4% đến 78,6% tổng $\Sigma 16\text{PAHs}$). Đây là các hợp chất có tính kỵ nước cực cao ($\log K_{ow} > 5,2$) và là các tác nhân chính gây độc tính tế bào và ung thư.
  3. Quy luật phân bố theo độ sâu thổ nhưỡng: Tại vị trí ĐR5, nồng độ $\Sigma 16\text{PAHs}$ ghi nhận giá trị cao nhất tại tầng đất mặt 0 - 5 cm, giảm nhẹ ở tầng 5 - 10 cm và 10 - 15 cm, sau đó tăng trở lại tại tầng 15 - 20 cm. Hiện tượng này phản ánh 2 giai đoạn phát thải rõ rệt: tầng sâu 15 - 20 cm lưu giữ tàn dư của hoạt động khai thác than và công nghiệp nặng trong quá khứ, trong khi tầng mặt 0 - 5 cm tiếp nhận tải lượng ô nhiễm liên tục từ các hoạt động đốt nhiệt điện và giao thông hiện tại.
  4. Xác định chính xác nguồn gốc phát thải pyrogenic: Phân tích tỷ lệ đồng phân phân tử cho thấy:
    • Tỷ lệ $Flt / (Flt + Pyr) > 0,5$ (dao động 0,52 - 0,68)
    • Tỷ lệ $Ant / (Ant + Phe) > 0,1$ (dao động 0,14 - 0,28)
    • Tỷ lệ $BaA / (BaA + Chr) > 0,35$
    • Tỷ lệ $Ind / (Ind + BghiP) > 0,5$ Kết quả này là bằng chứng khoa học khẳng định nguồn ô nhiễm PAHs trong đất RNM Đồng Rui bắt nguồn áp đảo từ quá trình đốt cháy nhiệt độ cao (quá trình đốt than của nhà máy nhiệt điện, lò gạch ngói, đốt sinh khối và củi gỗ trong sinh hoạt), không bắt nguồn từ sự cố tràn dầu mỏ trực tiếp ($Flt/(Flt+Pyr) < 0,4$).
  5. Mối liên hệ giữa đặc tính đất và nồng độ PAHs: Phân tích tương quan Pearson cho thấy không có sự tương quan tuyến tính có ý nghĩa thống kê giữa hàm lượng TOC với $\Sigma 16\text{PAHs}$ ($r < 0,3; p > 0,05$), trong khi pH đất có tương quan nghịch nhẹ. Điều này khẳng định cơ chế tích lũy PAHs tại bãi triều Đồng Rui bị chi phối bởi dòng lắng đọng bụi khí quyển và quá trình giữ hạt mịn của rễ cây ngập mặn hơn là dung lượng trao đổi hữu cơ nội sinh.
Thông số Phân tích Khoảng Giá trị Quan trắc Nguồn thải / Đánh giá Quy chuẩn
$\Sigma 16\text{PAHs}$ đất mặt 45,2 - 892,4 $\mu\text{g/kg}$ Trung bình thấp hơn Ma Wan (Hồng Kông)
Tỷ lệ $Flt/(Flt + Pyr)$ 0,52 - 0,68 Đốt than, sinh khối, gỗ ($> 0,5$)
Tỷ lệ $Ant/(Ant + Phe)$ 0,14 - 0,28 Quá trình đốt nhiệt độ cao ($> 0,1$)
Tỷ lệ $BaA/(BaA + Chr)$ 0,36 - 0,54 Quá trình cháy nhiên liệu hóa thạch ($> 0,35$)
Thương số Rủi ro (RQ) $RQ < 1,0$ (hầu hết cấu tử) Rủi ro sinh thái cấp tính ở mức chấp nhận được
Chỉ số Rủi ro Ung thư (CR) $10^{-6} - 10^{-4}$ Cần giám sát định kỳ đối với trẻ em và ngư dân

Đánh giá rủi ro và mô hình hóa Fugacity

  • Rủi ro môi trường và sức khỏe: Thương số rủi ro $RQ$ của từng chất PAHs đều $< 1,0$, cho thấy rủi ro sinh thái cấp tính đối với thủy sinh vật ở mức trung bình và thấp. Tuy nhiên, chỉ số rủi ro ung thư tích lũy $CR$ đối với người dân phơi nhiễm đất bùn qua tiếp xúc da và nuốt ngẫu nhiên dao động trong khoảng $10^{-6}$ đến $10^{-4}$, nằm ở ngưỡng rủi ro tiềm tàng cần kiểm soát dài hạn, đặc biệt đối với trẻ em địa phương.
  • Mô phỏng Fugacity Cấp III & Cấp IV: Kết quả tính toán chỉ ra đất và trầm tích lưu giữ hơn 87% tổng lượng BaP phân bố trong toàn bộ hệ thống bãi triều Đồng Rui. Mô phỏng động học cấp IV dự báo nếu duy trì tải lượng phát thải khí quyển như hiện nay, nồng độ BaP trong đất RNM sẽ tăng dần đều theo thời gian, biến đất RNM thành kho lưu trữ chất độc lâu dài đe dọa sinh kế nuôi trồng thủy sản.

Limitations và Future Research

  1. Giới hạn không gian và độ sâu lấy mẫu: Luận án khảo sát lõi đất sâu 20 cm tại 1 vị trí đại diện (ĐR5). Cần mở rộng mạng lưới khoan lõi trầm tích sâu hơn (đến 50 - 100 cm) tại nhiều tiểu sinh cảnh khác nhau (rừng mắm, rừng đước, bãi trống) để giải mã đầy đủ niên đại địa hóa môi trường qua phương pháp đồng vị phóng xạ $^{210}\text{Pb}$ và $^{137}\text{Cs}$.
  2. Giả định đồng nhất trong mô hình Fugacity: Mô hình cấp III và IV coi đất và nước là các khoang đồng nhất hoàn toàn (well-mixed), chưa phản ánh hết sự biến thiên vi mô do dòng chảy thủy triều ra vào hàng ngày và mạng lưới rễ cây phức tạp.
  3. Chưa phân tích than đen (Black Carbon - BC): Luận án mới xác định TOC tổng số mà chưa tách riêng phân đoạn Black Carbon - thành phần có ái lực hấp phụ PAHs mạnh gấp hàng trăm lần so với hữu cơ vô định hình.
  4. Hướng nghiên cứu ưu tiên 5 - 10 năm tới:
    • Nghiên cứu động học phân hủy sinh học PAHs của các chủng vi sinh vật bản địa kỵ khí và hiếu khí trong vùng rễ cây ngập mặn.
    • Tích hợp mô hình thủy động lực học 2D/3D (Delft3D/MIKE 21) với mô hình Fugacity động để mô phỏng sự lan truyền chất ô nhiễm theo từng pha triều.
    • Ứng dụng tỷ lệ đồng vị carbon bền $\delta^{13}\text{C}$ của từng hợp chất PAHs riêng lẻ (Compound-Specific Isotope Analysis - CSIA) nhằm định lượng tỷ lệ đóng góp chính xác của từng nhà máy nhiệt điện than lân cận.

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

  • Ảnh hưởng học thuật: Luận án là công trình tiên phong tại Việt Nam thiết lập bộ dữ liệu toàn diện về 16 PAHs trong đất rừng ngập mặn và ứng dụng thành công mô hình Fugacity cấp III, cấp IV. Công trình mở ra hướng nghiên cứu kết hợp giữa hóa học phân tích môi trường, độc học sinh thái và mô hình toán học giải tích.
  • Tác động công nghiệp và năng lượng: Cung cấp bằng chứng thực nghiệm rõ ràng về dấu vết phát thải từ các nhà máy nhiệt điện than và công nghiệp khai khoáng tại Quảng Ninh, thúc đẩy các doanh nghiệp lắp đặt hệ thống lọc bụi tĩnh điện hiệu suất cao và chuyển đổi công nghệ xử lý khí thải đạt chuẩn môi trường.
  • Tác động chính sách: Cung cấp luận cứ khoa học trực tiếp cho Ủy ban Nhân dân tỉnh Quảng Ninh trong việc thực hiện Quy hoạch Môi trường tỉnh Quảng Ninh đến năm 2020, tầm nhìn 2030; hỗ trợ cơ sở dữ liệu triển khai Đề án thành lập Khu bảo tồn đất ngập nước Đồng Rui - Tiên Yên và Đề án bảo vệ và phát triển rừng ven biển ứng phó với biến đổi khí hậu theo Quyết định số 120/QĐ-TTg của Thủ tướng Chính phủ.

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

  • Nghiên cứu sinh và Giới học thuật: Tiếp cận quy trình hoàn chỉnh từ khảo sát hiện trường bãi triều, kỹ thuật phân tích GC/MS độ nhạy cao, chẩn đoán nguồn phát thải bằng tỷ lệ đồng phân, đến phương pháp lập trình mô hình toán Fugacity đa khoang.
  • Cơ quan Quản lý Nhà nước (Sở TN&MT Quảng Ninh, Tổng cục Môi trường): Sở hữu cơ sở dữ liệu định lượng phục vụ công tác thanh tra môi trường, quan trắc định kỳ và quy hoạch phân vùng phát triển công nghiệp ven biển.
  • Cộng đồng cư dân địa phương: Gián tiếp được bảo vệ sức khỏe và nguồn lợi thủy sản (hàu, cua, cá, ngán) thông qua các cảnh báo về nguy cơ phơi nhiễm và các chính sách kiểm soát ô nhiễm nguồn thả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ì?

Luận án đã mở rộng thành công lý thuyết cân bằng pha Fugacity của Donald Mackay vào hệ sinh thái đất rừng ngập mặn nhiệt đới ven biển chịu chế độ nhật triều cực đại (3,5 - 4,0 m). Bằng việc xác lập hệ số dung lượng khuếch tán $Z$ và các thông số vận chuyển liên pha $D_{ij}$ phù hợp với điều kiện Việt Nam, nghiên cứu đã chứng minh đất và trầm tích đóng vai trò bồn chứa vĩnh cửu tích lũy hơn 87% BaP, làm thay đổi quan điểm truyền thống vốn chỉ nhìn nhận rừng ngập mặn như vùng đệm sinh học lọc sạch tự nhiên mà bỏ qua nguy cơ tích lũy độc chất lâu dài.

2. Điểm 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 trước tại Việt Nam như Dương Thanh Nghị et al. (2011) chỉ lấy mẫu đơn lẻ tại 2 thời điểm ở Vịnh Hạ Long, hay Nghiêm Trung Dũng (2005) chỉ khảo sát trong không khí, luận án của Đỗ Thị Lan Chi đã thiết kế quy trình quan trắc đa thời gian (6 đợt liên tục trong 4 năm), đa không gian (mặt phẳng 12 điểm và mặt cắt đứng 4 tầng sâu 0 - 20 cm), đồng thời kết hợp phương pháp hóa học phân tích thực nghiệm GC/MS với hai cấp độ mô hình toán học giải tích: Fugacity cấp III (phân bố tĩnh cân bằng động) và cấp IV (mô phỏng tích lũy động học phi cân bằng).

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

Phát hiện bất ngờ nhất là sự thiếu vắng tương quan có ý nghĩa thống kê giữa hàm lượng cacbon hữu cơ tổng số (TOC) và tổng nồng độ PAHs ($r < 0,3$). Trái ngược với giả định địa hóa học phổ biến cho rằng đất có hàm lượng mùn hữu cơ cao sẽ giữ nhiều PAHs hơn, kết quả thực nghiệm tại Đồng Rui khẳng định quá trình tích lũy PAHs tại đây bị chi phối áp đảo bởi nguồn phát thải khí quyển theo mùa gió Đông Bắc và cơ chế bẫy hạt mịn cơ học của thảm thực vật ngập mặn.

4. Luận án có cung cấp quy trình tái lập (Replication Protocol) không?

Có. Luận án chi tiết hóa toàn bộ giao thức thực nghiệm bao gồm: tọa độ địa lý 12 điểm lấy mẫu; quy trình bảo quản mẫu ở $-20^\circ\text{C}$; thể tích dung môi chiết xuất diclorometan:hexan (1:1); chương trình nhiệt độ lò GC/MS; các ion định lượng (SIM mode); đường chuẩn BaP; công thức tính toán độ tập trung $Z_i$, tải lượng đối lưu $D_{ai}$, tải lượng phân hủy $D_{ri}$ và bảng thông số thể tích khoang môi trường xã Đồng Rui.

5. Chương trình nghiên cứu 10 năm tiếp theo được định hình như thế nào?

Chương trình nghiên cứu 10 năm tập trung vào: (1) Tích lũy sinh học và chuyển dịch PAHs qua các bậc dinh dưỡng từ mùn bã hữu cơ $\rightarrow$ sinh vật đáy (nghêu, sò, hàu) $\rightarrow$ cá/chim $\rightarrow$ con người; (2) Cơ chế tương tác giữa ô nhiễm PAHs và vi nhựa (Microplastics) trong đất ngập mặn; (3) Giải pháp phục hồi sinh học (Bioremediation) sử dụng thực vật ngập mặn kết hợp chế phẩm vi sinh bản địa nhằm phân hủy PAHs tại chỗ.

Kết luận

  1. Khảo sát toàn diện và thiết lập bộ dữ liệu cơ sở khoa học đầu tiên tại Việt Nam về hiện trạng tồn lưu của 16 PAHs trong môi trường đất rừng ngập mặn xã Đồng Rui theo cả không gian mặt bằng, độ sâu tầng đất (0 - 20 cm) và chuỗi thời gian 6 đợt quan trắc từ 2014 đến 2017.
  2. Chứng minh nồng độ PAHs trong đất RNM có sự phân hóa mùa rõ rệt, nồng độ mùa đông cao hơn mùa hè do sự chi phối của gió mùa Đông Bắc mang bụi chứa PAHs từ các vùng phát thải công nghiệp lân cận lắng đọng xuống hệ sinh thái.
  3. Xác định nhóm PAHs cao phân tử (4 - 6 vòng benzen) chiếm ưu thế áp đảo và sử dụng các tỷ lệ đồng phân phân tử ($Flt/(Flt+Pyr)$, $Ant/(Ant+Phe)$, $BaA/(BaA+Chr)$, $Ind/(Ind+BghiP)$) để chứng minh nguồn phát thải chính gây ô nhiễm là từ các quá trình nhiệt phân nhiên liệu hóa thạch (đốt than, nhà máy nhiệt điện, đốt sinh khối), loại trừ khả năng ô nhiễm do tràn dầu trực tiếp.
  4. Đánh giá định lượng rủi ro sinh thái ($RQ$) và rủi ro ung thư ($CR$), khẳng định rủi ro cấp tính ở mức an toàn nhưng rủi ro sức khỏe mãn tính đối với trẻ em và cộng đồng ngư dân tiếp xúc bùn đất nằm ở ngưỡng tiềm tàng cần quản lý giám sát.
  5. Tiên phong ứng dụng thành công mô hình toán học Fugacity cấp III và cấp IV của Mackay trong việc mô phỏng sự phân bố liên pha và dự báo xu thế tích lũy dài hạn của BaP trong đất RNM, chứng minh đất là khoang tích lũy chủ yếu (>87% khối lượng), tạo tiền đề phương pháp luận vững chắc cho các nghiên cứu mô hình hóa ô nhiễm môi trường đất ngập nước tại Việt Nam.