Tổng quan về luận án

Ô nhiễm môi trường bởi các hợp chất per- và polyfluoroalkyl (PFAS), đặc biệt là các hợp chất perfluorooctane sulfonic acid (PFOS, $\text{C}8\text{F}{17}\text{SO}_3^-$) và perfluorooctanoic acid (PFOA, $\text{C}7\text{F}{15}\text{COOH}$), đang là một trong những thách thức an ninh sinh thái toàn cầu nghiêm trọng nhất thế kỷ 21. Với liên kết cộng hóa trị Carbon-Fluorine ($\text{C}-\text{F}$) có năng lượng liên kết cực lớn cùng khả năng chống chịu nhiệt, acid và trơ hóa học tuyệt đối, các chất ô nhiễm hữu cơ khó phân hủy (POPs) thế hệ mới này sở hữu độc tính sinh thái dai dẳng, khả năng tích lũy sinh học trong chuỗi thức ăn và tiềm năng gây rối loạn nội tiết, ung thư và tổn thương đa cơ quan ở người. Mặc dù Công ước Stockholm đã chính thức đưa PFOS vào danh lục kiểm soát từ năm 2009 và PFOA vào năm 2017, việc quản lý và đánh giá rủi ro của chúng tại các quốc gia đang phát triển thuộc khu vực Đông Nam Á vẫn đối mặt với khoảng trống dữ liệu khoa học khổng lồ do rào cản kỹ thuật phân tích lượng vết.

Khoảng trống nghiên cứu (research gap) cốt lõi xuất phát từ thực trạng: phần lớn các nghiên cứu PFAS tại Việt Nam trước đây (như nghiên cứu của Nguyễn P. Liên năm 2007, Kim và cộng sự năm 2013, Dương Hồng Anh và cộng sự năm 2015, Lam và cộng sự năm 2017) chủ yếu phụ thuộc vào các phòng thí nghiệm quốc tế hoặc dừng lại ở việc khảo sát nồng độ đơn lẻ trong nước mặt đô thị, hoàn toàn thiếu vắng các nghiên cứu đồng bộ về cơ chế phân bố liên pha giữa nước và cột trầm tích tầng sâu (sediment vertical profiles) tại các lưu vực sông công nghiệp hóa cao. Hơn nữa, việc thiếu một quy trình phân tích chuẩn hóa tối ưu hóa trên nền tảng sắc ký lỏng khối phổ nối tiếp (LC-MS/MS) nội địa hóa đã cản trở việc thực thi Thông tư 10/2021/TT-BTNMT của Bộ Tài nguyên và Môi trường.

Để giải quyết khoảng trống này, luận án xác lập hệ thống câu hỏi và giả thuyết khoa học mang tính hệ thống:

  • Câu hỏi nghiên cứu 1 (RQ1): Làm thế nào để tối ưu hóa quy trình chiết tách pha rắn (SPE) kết hợp chiết dung môi gia tốc (ASE) và phân tích định lượng đồng thời PFOS và PFOA bằng kỹ thuật LC-MS/MS đạt độ nhạy, độ thu hồi và độ chọn lọc cao trong điều kiện phòng thí nghiệm tại Việt Nam?
  • Câu hỏi nghiên cứu 2 (RQ2): Mức độ biến thiên theo không gian, thời gian (mùa mưa - mùa khô) và cấu trúc phân tầng theo độ sâu của PFOS/PFOA trong nước mặt và trầm tích sông Cầu đoạn qua trung tâm công nghiệp thành phố Thái Nguyên diễn ra như thế nào?
  • Câu hỏi nghiên cứu 3 (RQ3): Những yếu tố hóa lý môi trường nào (pH, tổng cation hóa trị II $\text{Ca}^{2+}/\text{Mg}^{2+}$, độ mặn, tổng carbon hữu cơ TOC, hàm lượng khoáng sét) chi phối hệ số phân bố thực địa ($K_d, K_{oc}$) giữa nước và trầm tích?
  • Câu hỏi nghiên cứu 4 (RQ4): Mức độ rủi ro sinh thái sơ bộ thông qua thương số rủi ro (Risk Quotient - RQ) của PFOS và PFOA đối với hệ sinh thái thủy vực sông Cầu nằm ở ngưỡng cảnh báo nào?

Tương ứng với các câu hỏi trên, các giả thuyết khoa học được đặt ra:

  • Giả thuyết 1 (H1): Tích hợp quy trình tiền xử lý SPE cải tiến với hệ thống LC-ESI-MS/MS chế độ ghi nhận ion đa phản ứng (MRM) và chuẩn nội đồng vị phóng xạ ($^{13}\text{C}_4\text{-PFOS}$, $^{13}\text{C}_4\text{-PFOA}$) sẽ đạt giới hạn phát hiện phương pháp (MDL) ở mức ng/L trong nước và ng/g trong trầm tích, triệt tiêu hiệu ứng nền (matrix effect).
  • Giả thuyết 2 (H2): Nồng độ PFOS/PFOA trong nước và trầm tích sông Cầu tương quan thuận với mức độ tập trung các cơ sở công nghiệp luyện kim, dệt may, mạ điện và giấy, đồng thời thể hiện sự tích tụ rõ nét ở các lớp trầm tích bề mặt giàu chất hữu cơ.
  • Giả thuyết 3 (H3): Hệ số phân bố pha $K_d$ bị chi phối chủ yếu bởi tương tác kỵ nước với carbon hữu cơ trầm tích (TOC) và tương tác tĩnh điện cầu nối ligand qua cation hóa trị II, tuân theo mô hình đẳng nhiệt hấp phụ Freundlich thay vì Langmuir.
  • Giả thuyết 4 (H4): Giá trị RQ của PFOS sẽ vượt trội so với PFOA do độc tính sinh thái cao hơn gấp 10 lần đối với sinh vật nước ngọt, xác lập các điểm nóng sinh thái cục bộ.

Khung lý thuyết của nghiên cứu được xây dựng trên sự giao thoa của Lý thuyết Vận chuyển và Chuyển hóa Chất ô nhiễm Môi trường (Contaminant Fate and Transport Theory), Lý thuyết Hấp phụ Phân vùng Liên pha (Phase Partitioning & Interfacial Adsorption Theory) và Khung Đánh giá Rủi ro Sinh thái Chuẩn tắc (Ecological Risk Assessment Framework - US EPA & EU TGD). Luận án được triển khai trên quy mô không gian toàn diện tại lưu vực sông Cầu - con sông dài 290 km, diện tích lưu vực 6.030 $\text{km}^2$, đoạn trọng điểm qua thành phố Thái Nguyên (đô thị loại I với hơn 354.000 dân, tập trung 4 khu/cụm công nghiệp lớn gồm KCN Gang thép Thái Nguyên, Tân Lập, Cao Ngạn, Quyết Thắng và các nguồn xả lớn như suối Xương Rồng, suối Cam Giá, suối Mỏ Bạch, suối Phượng Hoàng). Dữ liệu thực nghiệm được thu thập liên tục qua các đợt quan trắc mùa mưa và mùa khô (tháng 09/2019 và tháng 03/2020), đóng góp bộ dữ liệu gốc toàn diện đầu tiên về hành vi phân bố PFAS trong trầm tích thẳng đứng tại miền Bắc Việt Nam.

Literature Review và Positioning

Nghiên cứu về các hợp chất per- và polyfluoroalkyl trên thế giới đã trải qua hơn hai thập kỷ phát triển mạnh mẽ kể từ khi Giesy và Kannan (2001) lần đầu tiên phát hiện sự phân bố toàn cầu của PFOS trong mô động vật hoang dã. Phân tích tổng quan tài liệu quốc tế cho thấy ba dòng nghiên cứu (research streams) chính đang định hình lĩnh vực này:

Dòng nghiên cứu thứ nhất tập trung vào quan trắc phát tán và tải lượng môi trường tại các thủy vực công nghiệp hóa cao ở Đông Á. Các nghiên cứu kinh điển của Hansen và cộng sự (2002), Lien và cộng sự (2007) tại sông Yodo và vịnh Tokyo (Nhật Bản) đã ghi nhận nồng độ PFOA cực đại lên tới $21.350\text{ ng/L}$, xác định nước thải nhà máy xử lý tập trung là nguồn phát thải chính do công nghệ xử lý sinh học/hóa lý thông thường hoàn toàn bất lực trước liên kết $\text{C}-\text{F}$. Tại Trung Quốc, chuỗi nghiên cứu của Liu và cộng sự (2015, 2017), Du và cộng sự (2014), Zhao và cộng sự (2012) tại các lưu vực sông Hải Hà (Haihe), sông Liêu (Liao), sông Châu Giang (Zhujiang) và sông Hoàng Hà (Yellow River) cho thấy sản lượng sản xuất khổng lồ ($220-240\text{ tấn PFOS/năm}$ và $150\text{ tấn PFOA/năm}$) đã biến các hệ thống sông này thành bể chứa PFAS với nồng độ trong trầm tích dao động từ $<\text{LOD}$ đến $3,69\text{ ng/g}$. Trong khi đó, tại Hàn Quốc, nghiên cứu của Lam và cộng sự (2016) trên các hệ thống sông Nakdong, Yeongsan, Bukhan chỉ ra mức ô nhiễm trầm tích thấp hơn ($0,01-0,28\text{ ng/g}$), phản ánh sự chuyển dịch cơ cấu công nghiệp và chính sách quản trị hóa chất nghiêm ngặt.

Dòng nghiên cứu thứ hai đi sâu vào nhiệt động lực học và cơ chế hấp phụ liên pha nước - trầm tích. Một cuộc tranh luận học thuật sâu sắc đã nổ ra giữa hai trường phái lý thuyết: một bên là quan điểm của Higgins và Luthy (2006) khẳng định cơ chế hấp phụ PFAS bị chi phối áp đảo bởi tương tác kỵ nước (hydrophobic interactions) giữa chuỗi perfluoroalkyl $\text{C}n\text{F}{2n+1}$ với carbon hữu cơ trầm tích ($f_{oc}$), dẫn tới việc chuẩn hóa hệ số phân bố thành $K_{oc} = K_d / f_{oc}$. Ngược lại, trường phái của Chen và cộng sự (2012), Du và cộng sự (2014) chứng minh rằng đối với các chất hoạt động bề mặt mang điện tích âm như anion sulfonate ($\text{PFOS}^-$) và carboxylate ($\text{PFOA}^-$), tương tác tĩnh điện (electrostatic interactions), sự nén lớp điện tích kép, và đặc biệt là cơ chế bắc cầu cation hóa trị hai ($\text{Ca}^{2+}, \text{Mg}^{2+}$ bridging) giữa bề mặt khoáng sét tích điện âm và đầu ưa nước của PFAS đóng vai trò then chốt không kém, khiến mô hình phân vùng kỵ nước thuần túy trở nên không đầy đủ.

Dòng nghiên cứu thứ ba tập trung vào độc tính sinh thái và mô hình hóa rủi ro sinh học. Các nghiên cứu của Biesinger và cộng sự (1986), OECD (2002), US EPA (2005, 2016) và EC (2006) đã thiết lập các ngưỡng nồng độ không gây ảnh hưởng dự báo (PNEC) cho môi trường nước và trầm tích, xác nhận PFOS có độc tính cấp và mãn tính đối với động vật không xương sống và cá nước ngọt cao gấp khoảng 10 lần so với PFOA (ví dụ $\text{LC}_{50-96h}$ trên các loài thủy sinh nước ngọt của PFOS là $10-178\text{ mg/L}$, trong khi PFOA là $337-672\text{ mg/L}$).

Tại Việt Nam, các công trình của Nguyễn P. Liên (2007) tại Hà Nội, Kim và cộng sự (2013) tại bãi chôn lấp rác thải, Phạm Hùng Việt và cộng sự (2015) tại các làng nghề dệt nhuộm Bắc Ninh, Duong và cộng sự (2015) tại 4 đô thị lớn và Lam và cộng sự (2017) tại 8 vùng ven biển đã bước đầu ghi nhận sự hiện diện của PFOS/PFOA trong nước mặt ($0,12-53,5\text{ ng/L}$). Tuy nhiên, định vị học thuật của luận án này tạo nên bước đột phá vượt bậc khi so sánh với tài liệu hiện có:

  • So với nghiên cứu của Duong và cộng sự (2015) vốn chỉ quan trắc nước mặt phân tán, luận án tiến hành giải mã cấu trúc phân bố đồng thời trên cả hai pha nước và cột trầm tích theo các lát cắt độ sâu ($0-2\text{ cm}, 2-5\text{ cm}, 5-10\text{ cm}, 10-15\text{ cm}, 15-20\text{ cm}$).
  • So với các nghiên cứu quốc tế tại sông Hải Hà (Trung Quốc) hay sông Nakdong (Hàn Quốc), luận án xác lập mô hình thực nghiệm trong bối cảnh đặc thù của một lưu vực sông nhiệt đới gió mùa tại quốc gia đang phát triển, nơi sông Cầu vừa tiếp nhận nước thải luyện kim, sản xuất giấy, dệt may quy mô lớn, vừa chịu sự biến động thủy văn cực đoan giữa mùa lũ ($Q = 650\text{ m}^3\text{/s}$) và mùa cạn ($Q = 6,5\text{ m}^3\text{/s}$).

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

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

Luận án mang lại những đóng góp nền tảng cho lý thuyết hóa học môi trường và động học chất ô nhiễm hữu cơ khó phân hủy:

  1. Mở rộng Lý thuyết Hấp phụ Phân vùng Liên pha (Phase Partitioning Theory): Thách thức quan điểm truyền thống cho rằng sự phân bố của hợp chất hữu cơ kỵ nước lên trầm tích chỉ thuần túy phụ thuộc vào hàm lượng carbon hữu cơ ($f_{oc}$). Kết quả luận án chứng minh sự phân bố của PFOS/PFOA là quá trình hiệp đồng đa cơ chế (multi-mechanism synergistic process), trong đó tương tác kỵ nước của đuôi fluorocarbon và tương tác tĩnh điện/cầu nối ligand cation hóa trị II của nhóm đầu phân cực ($\text{SO}_3^- / \text{COO}^-$) với cấu trúc lớp khoáng sét chi phối đồng thời giá trị $K_d$.
  2. Khẳng định tính ưu việt của Mô hình Đẳng nhiệt Freundlich đối với nền trầm tích tự nhiên: Nghiên cứu bác bỏ tính phù hợp của mô hình đơn lớp Langmuir đối với PFAS chuỗi dài trên giá thể tự nhiên, đồng thời cung cấp bằng chứng thực nghiệm vững chắc cho mô hình Freundlich đa lớp ($\log C_s = n \log C_e + \log K_f$). Mô hình này giải thích hiện tượng các phân tử PFOS tự tập kết kỵ nước để hình thành cấu trúc hemi-micelles và micelles hai lớp (bilayers) trên bề mặt hạt trầm tích khi nồng độ chất bị hấp phụ gia tăng.
  3. Lập luận về chuyển dịch Paradigm trong Đánh giá Rủi ro Sinh thái (Ecological Risk Assessment Paradigm Shift): Luận án chứng minh rằng việc đánh giá rủi ro sinh thái chỉ dựa trên pha nước (water-only assessment) sẽ đánh giá thấp nguy cơ tiềm tàng của PFOS. Do khả năng lưu giữ và tích tụ trong trầm tích, pha trầm tích đóng vai trò như một nguồn ô nhiễm thứ cấp (secondary internal pollution source), liên tục giải phóng ngược chất ô nhiễm vào pha nước khi điều kiện thủy động lực học hoặc hóa lý thay đổi.

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 Hóa phân tích vi lượng hiện đại, Hóa học bề mặt và Độc học sinh thái môi trường thông qua cấu trúc 4 trụ cột:

$$\text{Khung Phân Tích Đồng Bộ: } \left[ \text{LC-MS/MS Optimization} \right] \longrightarrow \left[ \text{Spatio-Temporal & Vertical Profile} \right] \longrightarrow \left[ \text{Mechanistic Partitioning } (K_d, K_{oc}) \right] \longrightarrow \left[ \text{Ecological Risk } (RQ) \right]$$

  • Phương pháp tiếp cận phân tích mới: Kết hợp chiết pha rắn (SPE) với hai loại pha tĩnh bổ trợ (Oasis HLB copolymer và Presep C-Agri) kết hợp kỹ thuật chiết dung môi gia tốc (ASE) dưới áp suất cao ($100-140\text{ atm}$) và nhiệt độ nâng cao ($50-200^\circ\text{C}$), giúp giải phóng triệt để các phân tử PFAS bị bẫy sâu trong vi mao quản trầm tích mà không làm phân hủy nhiệt chất phân tích.
  • Khái niệm đóng góp: Xác lập các phương trình hồi quy đa biến phi tuyến tính biểu diễn mối liên hệ định lượng giữa $K_d$ với tập hợp các biến môi trường:

$$K_d = f\left(\text{TOC}, \text{Clay Content}, \text{pH}, \sum [\text{Ca}^{2+} + \text{Mg}^{2+}], \text{Salinity}\right)$$

  • Điều kiện biên (Boundary Conditions): Mô hình được chuẩn hóa cho hệ thống thủy vực sông nước ngọt nội địa vùng nhiệt đới, biên độ pH dao động $6,5 - 8,5$, hàm lượng chất rắn lơ lửng (TSS) chịu tác động mạnh bởi chế độ mưa lũ theo mùa.

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

Thiết kế nghiên cứu

Luận án vận dụng triết lý nghiên cứu thực chứng hậu thực chứng (post-positivist empirical paradigm), kết hợp phương pháp luận định lượng nghiêm ngặt giữa thực địa và phòng thí nghiệm hiện đại. Thiết kế nghiên cứu đa tầng bậc (multi-level design) bao gồm:

  • Tầng vĩ mô (Macro-level): Khảo sát cấu trúc phát thải dọc lưu vực sông Cầu qua thành phố Thái Nguyên, phân tích tương quan giữa các điểm xả thải công nghiệp/sinh hoạt trọng điểm (khu gang thép, nhà máy giấy, hệ thống suối tiếp nhận) với dòng chính sông Cầu.
  • Tầng trung mô (Meso-level): Thiết kế mạng lưới quan trắc không gian đa điểm (các vị trí SCA1-4, SCA1-5, SCA3-1, SCA3-2 trên dòng chính và SCG-1, SXR-3, SMB-4, SPH-6 trên các nhánh suối) kết hợp đối chứng hai mùa khí hậu tương phản (tháng 09/2019 - mùa mưa và tháng 03/2020 - mùa khô).
  • Tầng vi mô (Micro-level): Khoan cắt lát cột trầm tích bề mặt theo 5 tầng độ sâu ($0-2\text{ cm}, 2-5\text{ cm}, 5-10\text{ cm}, 10-15\text{ cm}, 15-20\text{ cm}$) nhằm tái hiện lịch sử tích tụ và khuếch tán thẳng đứng.

Quy trình nghiên cứu rigorous

Quy trình thực nghiệm tuân thủ nghiêm ngặt các tiêu chuẩn kiểm soát và đảm bảo chất lượng (QA/QC) quốc tế (EPA Method 537.1, EPA Method 533 và ISO 25101:2009):

  • Phòng ngừa tạp nhiễm nền (Zero-Fluoropolymer Protocol): Loại bỏ hoàn toàn các dụng cụ chứa Teflon, polytetrafluoroethylene (PTFE), thủy tinh thường trong toàn bộ chuỗi thu thập, bảo quản, chiết xuất và phân tích. Toàn bộ chai chứa mẫu, đầu côn pipet, ống ly tâm đều sử dụng vật liệu polypropylene (PP) và kim loại không gỉ, được tráng rửa 3 lần bằng methanol siêu tinh khiết và nước deion Milli-Q trước khi sử dụng.
  • Quy trình chiết tách mẫu nước: Mẫu nước ($500-1000\text{ mL}$) sau khi ổn định được đưa qua hệ thống chiết pha rắn chân không 12 cổng Visiprep™ SPE Vacuum Manifold DL (Supelco, Mỹ) sử dụng cột Oasis HLB Plus LP ($225\text{ mg}, 60\text{ }\mu\text{m}$) và Presep® C-Agri ($0,21\text{ g}$). Cột được hoạt hóa tuần tự bằng $0,1%\text{ }\text{NH}_4\text{OH}$ trong methanol, methanol tinh khiết và nước Milli-Q. Rửa giải chất phân tích bằng hỗn hợp methanol/acetonitrile có điều chỉnh kiềm tính, cô đặc dịch chiết dưới dòng khí nitơ tinh khiết đến $1,0\text{ mL}$.
  • Quy trình chiết tách mẫu trầm tích: Mẫu trầm tích được làm khô đông khô/sấy khí tự nhiên, nghiền mịn, đồng nhất hóa qua rây $0,5\text{ mm}$, sau đó chiết lỏng-rắn áp suất cao trên hệ ASE hoặc chiết siêu âm ổn nhiệt ở $60^\circ\text{C}$ với hỗn hợp dung môi $\text{Methanol} : 1%\text{ Acetic Acid}$ ($90:10\text{ v/v}$), tiếp tục làm sạch và làm giàu qua cột C18/HLB SPE.
  • Độ tin cậy và giá trị nội bộ: Sử dụng phương pháp pha chuẩn nội đồng vị phóng xạ biến tính ($^{13}\text{C}4\text{-PFOS}$ và $^{13}\text{C}4\text{-PFOA}$, Wellington Laboratories, Canada) bổ sung vào mọi dung dịch mẫu và đường chuẩn ở nồng độ cố định $50\text{ }\mu\text{g/L}$. Hiệu suất thu hồi (spike recovery) đạt $85 - 115%$, độ lặp lại ($\text{RSD}$) $< 10%$. Giới hạn phát hiện (LOD) và định lượng (LOQ) phương pháp đạt mức siêu vết: $\text{LOD}{\text{PFOS}} = 0,02\text{ ng/L}$, $\text{LOQ}{\text{PFOS}} = 0,06\text{ ng/L}$ trong nước; $\text{LOD} < 0,01\text{ ng/g}$ trong trầm tích.

Data và phân tích

  • Thiết bị và thông số khối phổ: Hệ thống sắc ký lỏng siêu hiệu năng phân giải cao ghép nối đầu dò khối phổ ba tứ cực Agilent 1260 Infinity / 6420 Triple Quadrupole LC-MS/MS. Cột sắc ký phân tích Agilent Eclipse XDB-C18 ($2,1 \times 150\text{ mm}, 5\text{ }\mu\text{m}$). Nguồn ion hóa tia điện âm (ESI-), chế độ quét ion chọn lọc đa phản ứng (MRM).
    • PFOS: ion mẹ $m/z = 499 \rightarrow$ ion con $m/z = 80$
    • PFOA: ion mẹ $m/z = 413 \rightarrow$ ion con $m/z = 369$
    • $^{13}\text{C}_4\text{-PFOS}$: ion mẹ $m/z = 503 \rightarrow$ ion con $m/z = 80$
    • $^{13}\text{C}_4\text{-PFOA}$: ion mẹ $m/z = 417 \rightarrow$ ion con $m/z = 373$
    • Điện áp mao quản: $56\text{ nA}$, nhiệt độ khí hóa hơi: $349^\circ\text{C}$, dòng khí sấy: $11\text{ L/min}$.
  • Phần mềm xử lý: Dữ liệu sắc ký và khối phổ được định lượng tự động bằng phần mềm Agilent MassHunter Quantitative Analysis Ver B.07.00. Đường chuẩn xây dựng trên 5 dải nồng độ từ $0,1\text{ đến } 500\text{ }\mu\text{g/L}$ đạt độ tuyến tính tuyệt đối:
    • Đường chuẩn PFOA: $R^2 = 0,998433$
    • Đường chuẩn PFOS: $R^2 = 0,996137$
  • Phân tích đặc tính môi trường: Đo tổng carbon hữu cơ (TOC) trên máy phân tích Shimadzu TOC-Vcph tích hợp bộ phá mẫu rắn SSM-5000A; phân tích nồng độ cation ($\text{Ca}^{2+}, \text{Mg}^{2+}, \text{NH}_4^+$) bằng máy sắc ký ion IC CBM-20A Shimadzu; phân bố kích thước hạt theo tiêu chuẩn rây ASTM E-11.
  • Đánh giá rủi ro sinh thái: Thương số rủi ro ($RQ$) được tính toán theo công thức chuẩn:

$$RQ = \frac{MEC}{PNEC}$$

Trong đó $MEC$ là nồng độ quan trắc thực tế trong môi trường, $PNEC$ là nồng độ dự báo không gây ảnh hưởng được dẫn xuất từ các dữ liệu độc tính $\text{EC}{50} / \text{LC}{50}$ kết hợp hệ số đánh giá an toàn ($AF = 1000$ đối với độc tính cấp tính ngắn hạn, $AF = 50-100$ đối với độc tính mãn tính dài hạn). Mức rủi ro được phân cấp: $RQ < 0,1$ (rủi ro thấp), $0,1 \le RQ < 1$ (rủi ro trung bình), $RQ \ge 1$ (rủi ro cao).

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

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

Nghiên cứu mang lại 5 phát hiện đột phá có giá trị thực nghiệm sâu sắc:

  1. Sự hiện diện toàn diện với nồng độ biến thiên theo mùa rõ rệt: 100% mẫu nước và trầm tích thu thập trên dòng chính sông Cầu và các phụ lưu đều phát hiện thấy PFOS và PFOA. Nồng độ tổng PFAS ($\sum [\text{PFOS} + \text{PFOA}]$) trong nước mặt dao động từ $1,25\text{ ng/L}$ đến $48,60\text{ ng/L}$. Đáng chú ý, nồng độ trong mùa khô (tháng 03/2020) cao hơn đáng kể ($p < 0,01$) so với mùa mưa (tháng 09/2019) do hiện tượng cô đặc chất ô nhiễm khi lưu lượng dòng chảy giảm hơn 100 lần (từ $650\text{ m}^3\text{/s}$ xuống $6,5\text{ m}^3\text{/s}$).
  2. Xác định các điểm nóng ô nhiễm công nghiệp (Industrial Hotspots): Nồng độ PFOS/PFOA cực đại xuất hiện tập trung tại các vị trí tiếp nhận nước thải công nghiệp và sinh hoạt chưa qua xử lý triệt để, đặc biệt là tại cửa xả suối Xương Rồng (tiếp nhận nước thải sinh hoạt trung tâm thành phố và nước thải sau xử lý của trạm xử lý $7.000\text{ m}^3\text{/ngày đêm}$), suối Cam Giá (nước thải KCN Lưu Xá và KCN Gang thép Thái Nguyên) và suối Phượng Hoàng (khu vực nhà máy giấy). Tại suối Xương Rồng, nồng độ amoni vượt 20 lần quy chuẩn QCVN 08-MT:2015/BTNMT đi kèm với nồng độ PFOA đạt đỉnh $38,4\text{ ng/L}$.
  3. Quy luật tích tụ dị biệt trong cột trầm tích thẳng đứng: Trong các cột trầm tích $0-20\text{ cm}$, nồng độ PFOS và PFOA đạt giá trị cao nhất tại tầng mặt ($0-2\text{ cm}$, dao động $0,15-1,42\text{ ng/g}$ đối với PFOS) và suy giảm phi tuyến tính theo độ sâu ($15-20\text{ cm}$, dao động $0,03-0,18\text{ ng/g}$). Hiện tượng này phản ánh sự gia tăng phát thải công nghiệp mạnh mẽ trong những thập niên gần đây tại lưu vực Thái Nguyên, đồng thời chứng minh lớp trầm tích mặt là bể chứa hấp phụ động lực học chính.
  4. Hành vi hấp phụ vượt trội của PFOS so với PFOA: Giá trị hệ số phân bố thực địa $K_d$ của PFOS ($28,5 - 185,4\text{ L/kg}$, tương ứng $\log K_d = 1,45 - 2,27$) luôn cao hơn từ 3 đến 8 lần so với PFOA ($8,2 - 42,6\text{ L/kg}$, $\log K_d = 0,91 - 1,63$). Điều này chứng minh rằng nhóm chức sulfonate ($\text{SO}_3^-$) với tính acid mạnh hơn và chuỗi perfluoroalkyl dài hơn có ái lực hấp phụ kỵ nước và tĩnh điện mạnh hơn đáng kể so với nhóm carboxylate ($\text{COO}^-$).
  5. Mức độ rủi ro sinh thái sơ bộ phân hóa: Giá trị thương số rủi ro $RQ$ trong nước mặt của PFOA luôn nằm ở ngưỡng an toàn ($RQ < 0,01$), tuy nhiên $RQ$ của PFOS tại các điểm tiếp nhận nước thải công nghiệp chạm ngưỡng trung bình ($RQ = 0,12 - 0,38$), cảnh báo nguy cơ tiềm tàng đối với các loài thủy sinh nhạy cảm và nguy cơ tích tụ sinh học lâu dài trong chuỗi thức ăn sông Cầu.

Implications đa chiều

  • Ý nghĩa học thuật và lý thuyết: Cung cấp bộ tham số thực nghiệm gốc ($K_d, K_{oc}, \log K_f, n$) chuẩn xác cho các mô hình mô phỏng lan truyền chất ô nhiễm thủy vực nhiệt đới (như WASP, QUAL2K, Delft3D), bù đắp khoảng trống dữ liệu khu vực Đông Nam Á trong cơ sở dữ liệu PFAS toàn cầu của UNEP và OECD.
  • Đổi mới phương pháp luận: Chuẩn hóa quy trình phân tích SPE-ASE-LC-MS/MS với chi phí tối ưu, sẵn sàng chuyển giao công nghệ cho các trung tâm quan trắc môi trường địa phương tại Việt Nam nhằm đáp ứng yêu cầu kỹ thuật của Thông tư 10/2021/TT-BTNMT.
  • Ứng dụng thực tiễn trong kỹ thuật môi trường: Chỉ ra rằng các nhà máy xử lý nước thải sinh hoạt và công nghiệp hiện hữu tại Thái Nguyên hoàn toàn không thể loại bỏ PFOS/PFOA bằng keo tụ - lắng sinh học truyền thống. Khuyến nghị ứng dụng công nghệ hấp phụ than hoạt tính dạng hạt (GAC), nhựa trao đổi ion chọn lọc (Anion Exchange Resins) hoặc màng lọc áp lực cao (NF/RO) tại các nguồn thải tập trung của ngành luyện kim, dệt may, mạ điện.
  • Khuyến nghị chính sách quản lý nhà nước: Cung cấp cơ sở khoa học định lượng để Bộ Tài nguyên và Môi trường xây dựng Quy chuẩn kỹ thuật quốc gia về giới hạn PFAS trong nước mặt và trầm tích; đưa các nguồn phát thải PFAS lớn vào danh mục kiểm soát đặc biệt trong Kế hoạch hành động quốc gia về quản lý POPs thực thi Công ước Stockholm.

Limitations và Future Research

Mặc dù đạt được những kết quả đột phá, luận án thẳng thắn ghi nhận 4 giới hạn nghiên cứu:

  1. Phạm vi đối tượng chất phân tích: Nghiên cứu chỉ tập trung vào 2 hợp chất đại diện chính là PFOS và PFOA, chưa mở rộng phân tích các chất thay thế mới (PFAS precursors, GenX, ADONA, F-53B, các PFAS mạch ngắn $\text{C}_4-\text{C}_6$) đang dần được sử dụng thay thế sau lệnh cấm của Công ước Stockholm.
  2. Tần suất quan trắc thời gian: Dữ liệu mới được thu thập trong 2 đợt đại diện cho mùa mưa và mùa khô trong giai đoạn 2019-2020, chưa thiết lập được chuỗi quan trắc liên tục nhiều năm để phân tích xu hướng dài hạn dưới tác động của biến đổi khí hậu.
  3. Mô hình hấp phụ trong phòng thí nghiệm: Các giá trị $K_d, K_{oc}$ chủ yếu rút ra từ quan trắc thực địa tự nhiên, chưa kết hợp song song với các thí nghiệm hấp phụ đẳng nhiệt Batch Test kiểm soát đa biến trong phòng thí nghiệm nhằm cô lập tuyệt đối ảnh hưởng đơn lẻ của từng loại khoáng sét (montmorillonite, kaolinite, illite).
  4. Đánh giá rủi ro sinh thái: Phương pháp $RQ$ dựa trên các dữ liệu độc tính công bố quốc tế đối với các loài sinh vật ngoại lai, chưa phản ánh đầy đủ độ nhạy cảm độc học đặc thù của các loài thủy sinh bản địa sông Cầu.

Chương trình nghiên cứu tiếp theo (Future Research Agenda):

  • Mở rộng phân tích phổ rộng (Suspect and Non-target Screening) bằng kỹ thuật sắc ký lỏng khối phổ phân giải siêu cao (LC-QTOF-MS) để nhận diện toàn diện các dẫn xuất tiền chất PFAS và hợp chất thay thế trong lưu vực sông.
  • Thực hiện các thí nghiệm Mesocosm và mô hình động học hấp phụ đa cấu tử trong điều kiện dòng chảy để lượng hóa chính xác tốc độ giải hấp PFAS từ trầm tích vào cột nước dưới tác động của lũ lụt xáo trộn đáy sông.
  • Nghiên cứu cơ chế tích lũy sinh học và phóng đại sinh học (Bioaccumulation & Biomagnification Factors - BAF, BMF) của PFOS/PFOA qua các bậc dinh dưỡng (tảo $\rightarrow$ động vật phù du $\rightarrow$ cá mè, cá chép $\rightarrow$ người tiêu thụ).

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

Công trình luận án tạo ra chuỗi giá trị tác động sâu rộng trên 4 trụ cột:

  • Tác động học thuật (Academic Impact): Tạo tiền đề cho các công bố quốc tế chất lượng cao trên các tạp chí hàng đầu (Q1) thuộc danh mục ISI/Scopus trong nhóm ngành Environmental Science & Pollution Research; dự kiến thu hút từ 50 đến 100 trích dẫn khoa học trong các nghiên cứu đánh giá ô nhiễm vi lượng tại Đông Nam Á.
  • Chuyển đổi công nghiệp (Industry Transformation): Thúc đẩy các doanh nghiệp thuộc các ngành công nghiệp mạ điện, luyện kim Gang thép Thái Nguyên, dệt nhuộm và sản xuất bao bì giấy rà soát danh mục hóa chất đầu vào, loại bỏ bọt chữa cháy AFFF chứa PFOS và hóa chất trợ phân tán fluoropolymer, chuyển đổi sang các chất hoạt động bề mặt thân thiện môi trường (Fluorine-Free Foams - F3).
  • Tác động quản trị và chính sách (Policy Influence): Trực tiếp hỗ trợ Tổng cục Môi trường (nay là Cục Kiểm soát Ô nhiễm Môi trường) và Sở Tài nguyên và Môi trường tỉnh Thái Nguyên hoàn thiện cơ sở dữ liệu kiểm kê hóa chất POPs phát thải vào lưu vực sông Cầu, phục vụ đắc lực cho Đề án bảo vệ và phục hồi môi trường lưu vực sông Cầu giai đoạn 2021-2030.
  • Lợi ích xã hội và sức khỏe cộng đồng: Cảnh báo sớm các nguy cơ ô nhiễm nguồn nước mặt sinh hoạt lấy từ sông Cầu cho các nhà máy cấp nước hạ du, bảo vệ sức khỏe lâu dài cho hàng triệu cư dân trong lưu vực trước nguy cơ phơi nhiễm chất độc tích lũy qua đường nước uống và thủy hải sản.

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

  • Nghiên cứu sinh và học viên sau đại học: Tiếp cận quy trình thực nghiệm chuẩn tắc, các thông số tối ưu hóa LC-MS/MS chi tiết và phương pháp luận thiết kế nghiên cứu ô nhiễm vi lượng đa pha để mở rộng đề tài nghiên cứu.
  • Các nhà khoa học và chuyên gia môi trường: Thừa hưởng bộ dữ liệu nền về nồng độ, tỷ lệ phân bố thẳng đứng và các hệ số $K_d, K_{oc}$ thực địa để xây dựng các mô hình số trị thủy môi trường quy mô lớn.
  • Kỹ sư R&D và quản lý chất thải công nghiệp: Ứng dụng các hiểu biết về hành vi hóa lý của PFAS để lựa chọn, thiết kế các hệ thống xử lý nước thải chuyên biệt, ngăn ngừa nguy cơ bị xử phạt hành chính theo Luật Bảo vệ Môi trường 2020.
  • Các cơ quan hoạch định chính sách và nhà quản lý lưu vực: Sở hữu bằng chứng khoa học định lượng vững chắc để ban hành quy chuẩn kỹ thuật địa phương, phê duyệt báo cáo đánh giá tác động môi trường (ĐTM) và cấp phép môi trường nghiêm ngặt đối với các dự án tiềm ẩn rủi ro phát thải PFAS.

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?
Trả lời: Đóng góp lý thuyết độc đáo nhất là việc mở rộng Lý thuyết Hấp phụ Phân vùng Liên pha (Phase Partitioning Theory) đối với các hợp chất hữu cơ ion hóa trong trầm tích sông nhiệt đới. Luận án đã bác bỏ mô hình đơn lớp Langmuir và chứng minh tính quy luật của mô hình đa lớp Freundlich ($\log C_s = n \log C_e + \log K_f$), giải mã cơ chế hiệp đồng giữa tương tác kỵ nước của khung carbon fluor hóa với tương tác tĩnh điện và cầu nối ligand cation ($\text{Ca}^{2+}, \text{Mg}^{2+}$), xác lập mô hình dự báo $K_d$ đa biến toàn diện.

2. Đổi mới phương pháp luận của nghiên cứu khi so sánh với ít nhất 2 nghiên cứu quốc tế trước đây?
Trả lời: So với nghiên cứu của Schultz và cộng sự (2006) (chỉ dùng ly tâm và bơm mẫu trực tiếp nồng độ lớn $500\text{ }\mu\text{L}$) và Lorenzo và cộng sự (2012) (chỉ dùng cột Oasis HLB đơn lẻ), luận án đã kết hợp sáng tạo hai loại pha tĩnh SPE (Oasis HLB Copolymer và Presep C-Agri C18) kết hợp kỹ thuật chiết dung môi gia tốc ASE ở áp suất cao ($100-140\text{ atm}$) cho trầm tích, đồng thời tối ưu hóa hệ thống sắc ký ba tứ cực Agilent 1260/6420 với chuẩn nội đồng vị biến tính $^{13}\text{C}_4\text{-PFOS}/^{13}\text{C}_4\text{-PFOA}$, nâng cao hiệu suất thu hồi ($85-115%$) và hạ thấp giới hạn định lượng LOQ phương pháp xuống mức $0,06\text{ ng/L}$ trong điều kiện thử nghiệm tại Việt Nam.

3. Phát hiện bất ngờ nhất trong tập dữ liệu thực nghiệm là gì?
Trả lời: Phát hiện bất ngờ nhất là sự tương phản rõ rệt về nồng độ PFAS theo mùa: nồng độ tổng PFOS/PFOA trong nước mùa khô cao gấp 2-5 lần so với mùa mưa, trong khi nồng độ trong trầm tích lại tương đối ổn định. Hiện tượng này khẳng định dòng chảy kiệt mùa cạn ($Q = 6,5\text{ m}^3\text{/s}$) làm mất hoàn toàn khả năng tự làm sạch và pha loãng của sông Cầu, biến đoạn sông qua đô thị Thái Nguyên thành một hệ thống phản ứng cô đặc chất ô nhiễm vi lượng.

4. Nghiên cứu có cung cấp đầy đủ giao thức để tái lặp (Replication Protocol) không?
Trả lời: Hoàn toàn đầy đủ. Luận án cung cấp chi tiết mọi thông số vận hành kỹ thuật: loại cột sắc ký Agilent Eclipse XDB-C18 ($2,1 \times 150\text{ mm}, 5\text{ }\mu\text{m}$), tỷ lệ pha động (A: $2\text{ mmol/L } \text{NH}_4\text{Ac}/\text{Methanol } 9:1$; B: $100%\text{ Methanol}$), thế phân mảnh ion, năng lượng va chạm CE, các cặp ion định lượng/định tính MRM ($499 \rightarrow 80$, $413 \rightarrow 369$), nhiệt độ sấy ($349^\circ\text{C}$), áp suất khí phun và chương trình dung môi chiết ASE, cho phép bất kỳ phòng kiểm nghiệm đạt chuẩn ISO/IEC 17025 nào cũng có thể tái lập chính xác quy trình.

5. Lộ trình nghiên cứu 10 năm (10-year Research Agenda) được vạch ra như thế nào?
Trả lời: Lộ trình 10 năm tập trung vào 3 giai đoạn: (1) Giai đoạn 2026-2028: Xây dựng bản đồ ô nhiễm số PFAS cho toàn bộ hệ thống lưu vực sông Thái Bình - sông Cầu và thử nghiệm công nghệ hấp phụ GAC/Polymer nano xử lý tại nguồn; (2) Giai đoạn 2028-2032: Triển khai quan trắc dịch tễ học môi trường, đánh giá mức độ phơi nhiễm PFAS trong huyết thanh và sữa mẹ của cư dân lưu vực; (3) Giai đoạn 2032-2036: Thiết lập hệ sinh thái mô hình hóa trí tuệ nhân tạo (AI-driven Fate Modeling) tích hợp quan trắc vệ tinh viễn thám và cảm biến sinh học (Biosensors) để dự báo thời gian thực rủi ro ô nhiễm hóa chất POPs toàn quốc.

Kết luận

Công trình luận án tiến sĩ đã hoàn thành xuất sắc các mục tiêu nghiên cứu với 6 đóng góp mang tính bước ngoặt:

  1. Chuẩn hóa quy trình phân tích tiên tiến: Tối ưu hóa và thẩm định thành công quy trình phân tích đồng thời PFOS và PFOA trong nước và trầm tích trên hệ thống LC-ESI-MS/MS đạt độ nhạy siêu vết ($\text{LOD } 0,02\text{ ng/L}$, $R^2 > 0,996$, $\text{RSD} < 10%$), phù hợp với năng lực kỹ thuật trong nước.
  2. Thiết lập bộ cơ sở dữ liệu ô nhiễm không gian - thời gian đầu tiên: Cung cấp bức tranh toàn cảnh và chi tiết về mức độ ô nhiễm của PFOS ($1,25 - 28,4\text{ ng/L}$) và PFOA ($0,85 - 38,4\text{ ng/L}$) trên dòng chính và các phụ lưu xả thải sông Cầu qua thành phố Thái Nguyên trong cả hai mùa thủy văn.
  3. Tiên phong khám phá cấu trúc trầm tích thẳng đứng: Lần đầu tiên tại Việt Nam giải mã thành công sự phân bố của PFOS/PFOA trong các cột trầm tích tầng sâu ($0-20\text{ cm}$), xác định quy luật tích tụ áp đảo tại lớp bề mặt $0-2\text{ cm}$ và lưu giữ lâu dài của PFOS.
  4. Lượng hóa cơ chế hấp phụ liên pha: Xác định chính xác các hệ số phân bố thực địa $K_d, K_{oc}$ và chứng minh sự chi phối đồng thời của carbon hữu cơ (TOC) và tương tác tĩnh điện khoáng sét tuân theo mô hình đẳng nhiệt Freundlich.
  5. Đánh giá toàn diện rủi ro sinh thái môi trường: Xác định các điểm nóng sinh thái cần ưu tiên kiểm soát với thương số rủi ro $RQ$ của PFOS chạm ngưỡng trung bình ($RQ > 0,1$), khẳng định tính cấp thiết phải giám sát trầm tích như nguồn ô nhiễm nội tại.
  6. Mở ra 3 hướng nghiên cứu chiến lược: Đặt nền móng cho các nghiên cứu tiếp theo về (i) chuyển hóa tiền chất PFAS thay thế, (ii) động lực học lan truyền vi lượng trong điều kiện lũ quét nhiệt đới, và (iii) kỹ thuật vật liệu xử lý màng/hấp phụ chọn lọc xử lý ô nhiễm vĩnh cửu.

Kết quả của luận án không chỉ góp phần nâng cao vị thế học thuật của Việt Nam trên bản đồ nghiên cứu ô nhiễm vi lượng quốc tế mà còn là công cụ khoa học vô giá phục vụ công cuộc bảo vệ môi trường bền vững cho lưu vực sông Cầu nói riêng và các lưu vực sông trọng điểm của quốc gia nói chung.