Tổng quan về luận án

Nghiên cứu về số phận môi trường và độc tính của các chất ô nhiễm hữu cơ khó phân hủy (POPs), đặc biệt là nhóm Hydrocarbon thơm đa vòng (Polycyclic Aromatic Hydrocarbons – PAHs), đang là trọng tâm then chốt trong khoa học môi trường toàn cầu. Luận án tiến sĩ kỹ thuật của nghiên cứu sinh Đỗ Thị Lan Chi với đề tài "Nghiên cứu sự tồn lưu và rủi ro môi trường của các chất hữu cơ thơm đa vòng PAHs trong đất rừng ngập mặn xã Đồng Rui, huyện Tiên Yên, tỉnh Quảng Ninh" (Chuyên ngành: Môi trường Đất và Nước, Mã số: 62.44.03.01, Trường Đại học Thủy Lợi, 2018) là công trình nghiên cứu hệ thống, tiên phong tại Việt Nam tích hợp quan trắc hóa vi lượng phân đoạn không gian – thời gian với mô hình hóa nhiệt động lực học môi trường Fugacity cấp III và cấp IV.

Rừng ngập mặn (RNM) xã Đồng Rui là hệ sinh thái đất ngập nước ven biển tiêu biểu của miền Bắc Việt Nam, có diện tích tự nhiên 4.974,2 ha, nằm kẹp giữa sông Ba Chẽ và vịnh Tiên Yên với chế độ nhật triều thuần nhất có biên độ triều vào loại lớn nhất cả nước (3,5 – 4,0 m). Khu vực này chịu áp lực kép từ các nguồn phát thải pyrogenic quy mô lớn xung quanh (khai thác chế biến than Cẩm Phả – Cửa Ông cách 70 km, Nhà máy Nhiệt điện Mông Dương I & II cách 17 km, đốt sinh khối và giao thông thủy bộ) cùng quá trình bồi tụ trầm tích tự nhiên.

Khoảng trống nghiên cứu (Research Gap) lớn nhất được xác định là: Mặc dù các nghiên cứu trước đây tại Việt Nam (như Nghiêm Trung Dũng, 2005; Dương Thanh Nghị và cộng sự, 2011; Phạm Thị Kha, 2014; Đặng Hoài Nhơn và cộng sự, 2014) đã khảo sát PAHs trong không khí, nước, trầm tích cửa sông và sinh vật biển, nhưng hoàn toàn thiếu vắng các công trình điều tra định lượng có tính chu kỳ nhiều năm về sự tích lũy của 16 hợp chất PAHs ưu tiên (theo phân loại của US EPA) trong môi trường đất rừng ngập mặn theo cả lát cắt bề mặt, độ sâu phân tầng (0 – 20 cm) và chuỗi biến thiên thời gian thực nghiệm. Đồng thời, việc ứng dụng mô hình toán đa khoang Fugacity để mô phỏng sự chuyển vận liên pha và dự báo động thái tích lũy của Benzo(a)pyrene (BaP) trong đất ngập mặn tại Việt Nam chưa từng được triển khai bài bản.

Luận án đặt ra và giải quyết 3 câu hỏi nghiên cứu (Research Questions) và kiểm chứng 3 giả thuyết khoa học (Hypotheses) then chốt:

  • RQ1: Mức độ tồn lưu, đặc tính phân bố theo không gian (bề mặt, độ sâu 0–20 cm) và quy luật biến thiên theo thời gian (mùa mưa, mùa khô, qua các năm) của 16 hợp chất PAHs trong đất RNM Đồng Rui diễn ra như thế nào?
  • RQ2: Nguồn phát thải chính chi phối sự tích tụ PAHs tại RNM Đồng Rui có nguồn gốc từ đâu và mức độ rủi ro sinh thái (Risk Quotient – RQ), rủi ro ung thư (Cancer Risk – CR) đối với con người đạt mức nào?
  • RQ3: Sự chuyển vận khối lượng liên khoang giữa khí quyển, nước, đất và trầm tích cũng như xu thế tích lũy dài hạn của BaP trong đất RNM diễn biến ra sao dưới công cụ mô hình Fugacity cấp III và cấp IV?
  • H1: Sự phân bố của 16 PAHs trong đất RNM phụ thuộc chặt chẽ vào hàm lượng carbon hữu cơ tổng số (TOC), phản ứng đất (pH) và biến động theo hướng gió mùa cũng như chế độ thủy triều.
  • H2: Tỷ lệ đồng phân phân tử đặc trưng (diagnostic ratios) sẽ chỉ ra nguồn đốt cháy nhiên liệu hóa thạch/than đá chiếm ưu thế tuyệt đối so với nguồn dầu mỏ rò rỉ (petrogenic).
  • H3: Pha đất ngập mặn và trầm tích đóng vai trò là "bồn chứa" (sink) tích tụ chủ đạo đối với các hợp chất PAHs cao phân tử (như BaP) do dung lượng tập trung $Z$ lớn và khả năng giải phóng hạn chế.

Khung lý thuyết của luận án được xây dựng vững chắc dựa trên Lý thuyết Cân bằng Phân bố Đa môi trường của Donald Mackay (Fugacity Theory), Khung Đánh giá Rủi ro Sức khỏe Con người của Cơ quan Bảo vệ Môi trường Hoa Kỳ (US EPA Health Risk Assessment) và Học thuyết Địa hóa Hữu cơ về Động thái Hấp phụ Kỵ nước của Hợp chất Thơm. Luận án thực hiện quan trắc kéo dài từ năm 2013 đến tháng 11/2017 với 6 đợt lấy mẫu chuyên sâu (tháng 8/2014, 1/2015, 7/2015, 1/2016, 7/2016, 1/2017) trên 12 vị trí đại diện, thiết lập bộ số liệu chuẩn mực cho công tác bảo tồn đất ngập nước và quy hoạch môi trường tỉnh Quảng Ninh đến năm 2030.

Literature Review và Positioning

Tổng quan tài liệu quốc tế cho thấy các nghiên cứu về PAHs trong đất ngập nước và trầm tích biển phát triển mạnh mẽ từ những năm 2000. S. K. Sarkar và cộng sự (2010) khi nghiên cứu vùng rừng ngập mặn Sundarban (Vịnh Bengal, Ấn Độ) ghi nhận tổng nồng độ của 16 PAHs dao động trong khoảng 32 – 2.938 ng/g (trung bình 634 ng/g), trong đó 5 hợp chất PAHs gây ung thư cao (BbF, BkF, BaP, Ind, DahA) chiếm từ 68% đến 73% tổng lượng tích lũy. Rivelino M. Cavalcante và cộng sự (2009) đánh giá ảnh hưởng của quá trình đô thị hóa đến rừng ngập mặn Fortaleza (Brazil) đã xác định hàm lượng của 17 hợp chất PAHs dao động từ 3,04 đến 2.234,76 $\mu\text{g/kg}$ tại sông Cocó và 3,34 đến 1.859,21 $\mu\text{g/kg}$ tại sông Ceará, vượt xa ngưỡng tích tụ ở nhiều vùng công nghiệp phát triển. Tại Hồng Kông, các công trình của Tam et al. (2001) và các cộng sự tại các dải rừng ngập mặn Ting Kok ($119 \pm 59\ \mu\text{g/kg}$), Sheung Pak Nai ($159 \pm 73\ \mu\text{g/kg}$), Yi O ($313 \pm 473\ \mu\text{g/kg}$) và đặc biệt là Ma Wan ($2.071 \pm 1.033\ \mu\text{g/kg}$) đã chứng minh tác động sâu sắc của nguồn thải công nghiệp lịch sử thông qua phân tích mẫu lõi theo độ sâu.

Trong bức tranh học thuật đó, tồn tại các cuộc tranh luận lý thuyết và thực nghiệm gay gắt:

  1. Mối tương quan giữa chất hữu cơ trong đất (SOM/TOC) và nồng độ PAHs: Trường phái truyền thống cho rằng các chất hữu cơ kỵ nước như PAHs luôn có hệ số tương quan thuận chặt chẽ với hàm lượng hữu cơ trong đất (Karickhoff et al., 1979). Tuy nhiên, các nghiên cứu của Bucheli và cộng sự (2004), Zhang và cộng sự (2006), Wan-Li Ma và cộng sự (2013) tại lưu vực sông Songhua (Trung Quốc, hàm lượng SOM từ 0,76% đến 10,6%), cũng như Xiang Bi và cộng sự (2016) tại dãy Himalayas ($r^2 = 0,001$) lại chỉ ra rằng không có mối tương quan tuyến tính có ý nghĩa giữa PAHs và TOC. Điều này đặt ra câu hỏi về vai trò chi phối của dạng carbon đen (black carbon/soot) và động lực thủy triều so với carbon hữu cơ tự nhiên.
  2. Quy luật phân bố theo độ sâu: Nghiên cứu tại Sundarban (Ấn Độ) và Miami (Mỹ, Kamaljit Banger, 2010) chỉ ra nồng độ PAHs cao nhất tập trung ở tầng mặt (0–5 cm hoặc 0–15 cm) do nguồn phát thải khí quyển hiện tại. Ngược lại, nghiên cứu tại Ma Wan (Hồng Kông) và sông Ceará (Brazil) ghi nhận nồng độ tăng vọt ở tầng sâu 10–20 cm hoặc giảm dần tùy thuộc vào lịch sử công nghiệp hóa trong quá khứ.

Luận án của Đỗ Thị Lan Chi định vị chính xác vào khoảng trống này: Đánh giá sự tương tác giữa TOC, pH và nồng độ 16 PAHs trong môi trường đất ngập mặn nhiệt đới gió mùa có biên độ triều cao, kết nối trực tiếp dữ liệu phân bố hóa học với mô hình toán Fugacity để làm sáng tỏ bản chất chuyển vận khối lượng.

+-----------------------------------------------------------------------------+
|                            ĐỊNH VỊ HỌC THUẬT                               |
+-----------------------------------------------------------------------------+
|  QUỐC TẾ:                                                                   |
|  - Sundarban (Ấn Độ): Sarkar et al. (2010) [32 - 2.938 ng/g]                |
|  - Fortaleza (Brazil): Cavalcante et al. (2009) [3,04 - 2.234,76 µg/kg]     |
|  - Ma Wan (Hồng Kông): Tam et al. (2001) [791 - 3.758 µg/kg]                |
|  - Mô hình Fugacity: Mackay (1979, 2001); Paterson et al. (1994)             |
+-----------------------------------------------------------------------------+
                                       |
                                       v
+-----------------------------------------------------------------------------+
|  VIỆT NAM (Các nghiên cứu trước):                                          |
|  - Khí quyển Hà Nội: Nghiêm Trung Dũng (2005)                               |
|  - Nước, trầm tích, sinh vật Vịnh Hạ Long: Dương Thanh Nghị (2011)          |
|  - Trầm tích ven bờ Hải Phòng - Quảng Ninh: Phạm Thị Kha (2014)             |
|  - Trầm tích bãi triều miền Bắc: Đặng Hoài Nhơn et al. (2014)               |
+-----------------------------------------------------------------------------+
                                       |
                                       v
+-----------------------------------------------------------------------------+
|  ĐÓNG GÓP KHẢO SÁT & MÔ HÌNH HÓA CỦA LUẬN ÁN (ĐỖ THỊ LAN CHI, 2018):       |
|  - Đất RNM Đồng Rui (Tiên Yên, Quảng Ninh): Chu kỳ 6 đợt (2014-2017)        |
|  - Động thái 3 chiều: Bề mặt, độ sâu (0-20 cm), chuỗi thời gian             |
|  - Tích hợp rủi ro kép (RQ sinh thái & CR sức khỏe theo US EPA)             |
|  - Mô phỏng cơ chế nhiệt động học Fugacity cấp III & IV cho BaP             |
+-----------------------------------------------------------------------------+

Đó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 biên độ ứng dụng của Lý thuyết Fugacity (Donald Mackay, 1979, 2001) trong môi trường đất ngập nước chịu tác động triều mạnh. Khái niệm độ đào thoát ($f$, đơn vị Pa) và năng lực dung nạp/độ tập trung ($Z$, đơn vị $\text{mol}/(\text{m}^3\cdot\text{Pa})$) được tham số hóa chi tiết cho 4 khoang chính (khí quyển, nước, đất ngập mặn, trầm tích) và các tiểu khoang (hạt bụi, nước lỗ rỗng, chất rắn hữu cơ). Mối quan hệ cơ bản:

$$C = Z \times f$$

đã được lượng hóa cụ thể, chứng minh rằng sự phân bố của hợp chất có độ tan trong nước cực thấp và hệ số phân bố octanol-nước cao ($\log K_{ow}$ của BaP = 6,50; BbF = 6,12; BkF = 6,84; BghiP = 7,10) bị chi phối áp đảo bởi thông số dung nạp $Z_{soil}$ và $Z_{sediment}$.

Đồng thời, luận án đóng góp vào Lý thuyết Tỷ lệ Đồng phân Phân tử Địa hóa (Molecular Diagnostic Ratios) để truy nguyên nguồn gốc hợp chất hữu cơ:

  • Proposition 1: Nếu $Ant/(Ant + Phe) > 0,1$ và $Flt/(Flt + Pyr) > 0,5$, nguồn tích tụ mang bản chất pyrogenic (đốt than đá, củi gỗ, sinh khối).
  • Proposition 2: Nếu $BaA/(BaA + Chr) > 0,35$ và $Ind/(Ind + BghiP) > 0,5$, khẳng định nguồn phát thải bắt nguồn từ các quá trình đốt cháy nhiệt độ cao công nghiệp (nhà máy nhiệt điện, luyện kim, chế biến than) kết hợp giao thông cơ giới.
  • Proposition 3: Các quá trình vận chuyển khối lượng đối lưu ($D_A$) và phản ứng phân hủy quang hóa/sinh học ($D_R$) xác lập trạng thái dừng phi cân bằng (Steady-state non-equilibrium) tại mô hình cấp III và tích lũy động học theo phương trình vi phân phi tuyến tại mô hình cấp IV.

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

Khung phân tích của luận án tích hợp đa tầng 3 trụ cột khoa học:

  1. Lý thuyết Động lực học Đa khoang Fugacity: Mô phỏng sự vận chuyển liên pha thông qua các hệ số tải lượng $D_{ij}$ (giữa khí – nước, khí – đất, nước – trầm tích qua các cơ chế lắng đọng khô, lắng đọng ướt, khuếch tán phân tử và sa lắng hạt).
  2. Khung Đánh giá Độc tính Tương đương (TEQ) và Rủi ro Ung thư (US EPA): Ứng dụng hệ số độc hại tương đương $TEF$ (lấy BaP làm chuẩn = 1,0; DahA = 1,0; BaA, BbF, BkF, Ind = 0,1; Chr = 0,01) để quy đổi tổng nồng độ độc học $TEQ = \sum (C_i \times TEF_i)$.
  3. Mô hình Hóa học Môi trường và Địa hóa Trầm tích: Phân tích ma trận tương quan Pearson giữa các biến ngẫu nhiên $\Sigma 16\text{PAHs}$, TOC, pH và các thông số cơ lý hóa của đất bồi tụ.

Điều kiện biên (Boundary Conditions): Mô hình áp dụng cho hệ sinh thái mở RNM Đồng Rui có tọa độ $21^\circ 10' - 21^\circ 16'30''\text{N}$ và $107^\circ 21'30'' - 107^\circ 27'\text{E}$; thể tích các khoang môi trường được chuẩn hóa theo diện tích xã đảo 4.974,2 ha; độ cao tầng khí quyển 1.000 m; độ dày tầng đất hoạt tính 0,2 m; tầng nước mặt dao động theo triều 1,5 – 4,0 m; tầng trầm tích đáy 0,05 m.

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

Thiết kế nghiên cứu

Nghiên cứu được thiết kế theo trường phái thực chứng luận (Positivism) với quy trình phân tích định lượng chuẩn mực cao. Thiết kế nghiên cứu đa cấp độ (Multi-level experimental design) bao gồm:

  • Tầng không gian: 12 vị trí lấy mẫu đất đại diện trên toàn bộ diện tích RNM Đồng Rui (từ ĐR1 đến ĐR12), kết hợp mạng lưới lấy mẫu không khí xung quanh, mẫu nước mặt bãi triều và mẫu trầm tích lõi.
  • Tầng độ sâu: Phân đoạn giải phẫu mẫu đất ở 4 độ sâu chi tiết: 0 – 5 cm, 5 – 10 cm, 10 – 15 cm và 15 – 20 cm (thực hiện chi tiết tại điểm nóng ĐR5).
  • Tầng thời gian: 6 đợt quan trắc thực địa cách nhau 6 tháng trải dài qua 4 năm (tháng 8/2014, 1/2015, 7/2015, 1/2016, 7/2016, 1/2017), phản ánh đầy đủ chu kỳ mùa mưa (tháng 7–8, gió Nam/Đông Nam) và mùa khô (tháng 1–2, gió Bắc/Đông Bắc).
+-----------------------------------------------------------------------------+
|                         SƠ ĐỒ THIẾT KẾ NGHIÊN CỨU                           |
+-----------------------------------------------------------------------------+
|  KHẢO SÁT & QUAN TRẮC THỰC ĐỊA (6 đợt: 2014 - 2017; 12 vị trí ĐR1 -> ĐR12)  |
|  - Mẫu đất (0-5, 5-10, 10-15, 15-20 cm)  - Mẫu khí quyển xung quanh         |
|  - Mẫu nước mặt thủy triều               - Mẫu trầm tích bãi triều          |
+-----------------------------------------------------------------------------+
                                       |
                                       v
+-----------------------------------------------------------------------------+
|  PHÂN TÍCH HÓA LÝ & THỰC NGHIỆM CHUẨN (GC/MS, TOC, pH)                     |
|  - Chiết Soxhlet/Siêu âm DCM:n-Hexane (1:1) -> Cột sắc ký làm sạch Silica   |
|  - Đo sắc ký khí khối phổ GC/MS (Agilent) xác định 16 US EPA PAHs          |
|  - Phân tích TOC (Walkley-Black/Nhiệt phân), pH nước/đất                   |
+-----------------------------------------------------------------------------+
                                       |
                   +-------------------+-------------------+
                   |                                       |
                   v                                       v
+------------------------------------+  +------------------------------------+
| ĐÁNH GIÁ RỦI RO MÔI TRƯỜNG & Y TẾ  |  | MÔ HÌNH HÓA FUGACITY (MACKAY)      |
| - RQ (Risk Quotient sinh thái)     |  | - Cấp III: Cân bằng dừng đa khoang |
| - TEQ & Cancer Risk (CR theo EPA)  |  | - Cấp IV: Động thái tích lũy BaP   |
| - Đánh giá 3 đường phơi nhiễm      |  | - Tính toán thông số Z và D_ij     |
+------------------------------------+  +------------------------------------+

Quy trình nghiên cứu rigorous

Quy trình lấy mẫu, bảo quản và phân tích tuân thủ nghiêm ngặt các tiêu chuẩn của US EPA (Method 3540C, 3550B, 3630C và 8270D):

  • Xử lý mẫu đất và trầm tích: Mẫu đất sau khi thu thập bằng dụng cụ không gỉ được bọc giấy nhôm vô trùng, bảo quản lạnh ở $4^\circ\text{C}$, sấy khô tự nhiên trong bóng râm, loại bỏ sỏi vụn và rễ cây, nghiền mịn qua rây 0,15 mm (100 mesh).
  • Chiết tách và làm sạch: Chiết siêu âm lặp lại nhiều lần bằng hỗn hợp dung môi $n\text{-hexane} : \text{dichloromethane}$ ($1:1, v/v$). Dịch chiết được cô đặc bằng hệ thống cô quay chân không (Rotary Evaporator), sau đó làm sạch qua cột sắc ký hấp phụ chứa Silica gel/Alumina hoạt hóa để loại bỏ lipid và các hợp chất phân cực gây nhiễu.
  • Phân tích sắc ký khí – khối phổ (GC/MS): Thiết bị sắc ký khí khối phổ Agilent GC/MS với cột mao quản DB-5MS ($30\text{ m} \times 0,25\text{ mm} \times 0,25\ \mu\text{m}$). Chương trình nhiệt độ lò: $60^\circ\text{C}$ (giữ 1 phút), tăng $10^\circ\text{C/phút}$ lên $200^\circ\text{C}$, sau đó tăng $5^\circ\text{C/phút}$ lên $300^\circ\text{C}$ (giữ 10 phút). Khí mang Helium siêu tinh khiết ($99,999%$), tốc độ dòng $1,0\text{ mL/phút}$.
  • Kiểm soát chất lượng (QA/QC): Xây dựng đường chuẩn 5 điểm cho từng hợp chất PAH với hệ số tương quan $R^2 > 0,995$. Phân tích mẫu trắng (method blank), mẫu lặp (duplicate samples) và mẫu thêm chuẩn (matrix spike) với độ thu hồi (recovery) đạt từ 78,5% đến 106,2%, độ lệch chuẩn tương đối (RSD) $< 10%$. Giới hạn phát hiện (LOD) dao động từ 0,01 đến $0,05\ \mu\text{g/kg}$.

Data và phân tích toán học

Luận án triển khai 3 công cụ mô hình hóa và thống kê chuyên sâu:

  1. Phương pháp Thương số Rủi ro Môi trường (Risk Quotient - RQ):

    $$RQ = \frac{C_i}{C_{QV}}$$

    Trong đó $C_i$ là nồng độ thực tế của hợp chất PAH thứ $i$ trong đất; $C_{QV}$ là giá trị tiêu chuẩn chất lượng môi trường định mức (dựa trên hướng dẫn chất lượng trầm tích/đất của EPA và tiêu chuẩn quốc tế). Phân loại: $RQ < 0,1$ (rủi ro thấp); $0,1 \le RQ \le 1,0$ (rủi ro trung bình); $RQ > 1,0$ (rủi ro sinh thái cao).

  2. Phương pháp Đánh giá Rủi ro Ung thư (Cancer Risk - CR): Xác định chỉ số rủi ro ung thư tích lũy cả đời qua 3 con đường: nuốt phải đất phơi nhiễm ($CR_{ing}$), hấp thụ qua da ($CR_{dermal}$) và hít thở bụi đất mang PAHs ($CR_{inh}$):

    $$CR_{ing} = \frac{TEQ \times IR_{soil} \times EF \times ED \times CFS_{oral} \times 10^{-6}}{BW \times AT}$$

    $$CR_{dermal} = \frac{TEQ \times SA \times AF \times ABS_d \times EF \times ED \times CFS_{dermal} \times 10^{-6}}{BW \times AT}$$

    $$CR_{inh} = \frac{TEQ \times InhR \times EF \times ED \times CFS_{inh}}{PEF \times BW \times AT}$$

    $$CR = CR_{ing} + CR_{dermal} + CR_{inh}$$

    Trong đó $TEQ = \sum (C_i \times TEF_i)$; các thông số sinh trắc học và phơi nhiễm ($BW$, $AT$, $EF$, $ED$, $SA$, $AF$, $ABS_d$, $PEF$, $CSF$) được chuẩn hóa theo khuyến cáo của US EPA cho cả nhóm người lớn và trẻ em.

  3. Mô hình toán Fugacity Cấp III và Cấp IV (Mackay):

    • Cấp III (Hệ mở, trạng thái dừng, có vận chuyển liên pha): Thiết lập hệ phương trình cân bằng khối lượng tuyến tính cho 4 khoang môi trường:

      $$I_i + \sum (D_{ji} \cdot f_j) = f_i \cdot \left( D_{Ri} + D_{Ai} + \sum D_{ij} \right)$$

      với $I_i$ là tải lượng phát thải trực tiếp vào khoang $i$; $D_{Ri}$ là thông số phân hủy hóa học/sinh học ($D_{Ri} = V_i \cdot Z_i \cdot k_i$); $D_{Ai}$ là thông số chuyển động đối lưu cuốn trôi; $D_{ij}$ là thông số vận chuyển liên pha giữa khoang $i$ và $j$.

    • Cấp IV (Hệ mở, trạng thái động phi dừng, biến thiên theo thời gian): Giải hệ phương trình vi phân đạo hàm bậc nhất theo thời gian $t$:

      $$\frac{d(V_i \cdot Z_i \cdot f_i)}{dt} = I_i(t) + \sum (D_{ji} \cdot f_j) - f_i \cdot \left( D_{Ri} + D_{Ai} + \sum D_{ij} \right)$$

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

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

  1. Hiện trạng tích lũy và đặc tính phân tử của 16 PAHs trong đất RNM Đồng Rui: Nồng độ tổng $\Sigma 16\text{PAHs}$ trong đất bề mặt dao động mạnh qua các đợt khảo sát, trung bình toàn khu vực đạt mức từ vài chục đến hàng trăm $\mu\text{g/kg}$ trọng lượng khô. Đợt lấy mẫu tháng 8/2014 ghi nhận sự biến thiên rõ nét giữa 12 vị trí lấy mẫu, trong đó các điểm ĐR1, ĐR5, ĐR9 có nồng độ vượt trội do nằm gần các luồng lạch đón nhận dòng chảy từ sông Ba Chẽ và vùng tiếp giáp hoạt động tàu thuyền. Nhóm PAHs cao phân tử từ 4 đến 6 vòng benzen (HPAHs gồm BaA, Chr, BbF, BkF, BaP, Ind, DahA, BghiP) chiếm tỷ lệ áp đảo (dao động từ 55% đến trên 75% tổng lượng PAHs), trong khi nhóm thấp phân tử 2–3 vòng (LPAHs như Nap, Acy, Ace, Flu, Phe, Ant) chiếm tỷ trọng nhỏ do dễ bị bay hơi và phân hủy quang hóa.
+-----------------------------------------------------------------------------+
|          BẢNG TỔNG HỢP NỒNG ĐỘ VÀ CHỈ SỐ RỦI RO QUA 6 ĐỢT QUAN TRẮC         |
+-----------------------------------------------------------------------------+
| Đợt khảo sát | Mùa khí hậu | Hướng gió chính | Đặc trưng nồng độ & rủi ro   |
+--------------+-------------+-----------------+------------------------------+
| Tháng 8/2014 | Mùa hè/mưa  | Nam - Đông Nam  | Σ16PAHs cao ở cửa sông;      |
|              |             |                 | HPAHs chiếm >65%; CR an toàn |
+--------------+-------------+-----------------+------------------------------+
| Tháng 1/2015 | Mùa đông/khô| Bắc - Đông Bắc  | Ảnh hưởng nhiệt điện, mỏ than|
|              |             |                 | Lắng đọng khí quyển tăng cao |
+--------------+-------------+-----------------+------------------------------+
| Tháng 7/2015 | Mùa hè/mưa  | Nam - Đông Nam  | Rửa trôi bề mặt, phù sa bồi  |
|              |             |                 | Tương quan Pearson TOC-PAHs  |
+--------------+-------------+-----------------+------------------------------+
| Tháng 1/2016 | Mùa đông/khô| Bắc - Đông Bắc  | Σ8PAHs gây ung thư tăng nhẹ; |
|              |             |                 | RQ trung bình của BaP > 0,1  |
+--------------+-------------+-----------------+------------------------------+
| Tháng 7/2016 | Mùa hè/mưa  | Nam - Đông Nam  | Biến thiên theo chu kỳ nhật  |
|              |             |                 | triều (biên độ 3,5 - 4,0 m)  |
+--------------+-------------+-----------------+------------------------------+
| Tháng 1/2017 | Mùa đông/khô| Bắc - Đông Bắc  | Dấu ấn pyrogenic rõ nét;     |
|              |             |                 | Mô hình Fugacity khớp thực đo|
+-----------------------------------------------------------------------------+

  1. Quy luật biến thiên theo chiều sâu và lịch sử phát thải: Kết quả phân tích mẫu lõi tại điểm nóng ĐR5 qua 4 tầng độ sâu (0–5 cm, 5–10 cm, 10–15 cm, 15–20 cm) cho thấy nồng độ $\Sigma 16\text{PAHs}$ có xu hướng tập trung cao ở lớp mặt (0–5 cm) và tầng đáy lõi (15–20 cm). Nồng độ cao ở tầng sâu 15–20 cm phản ánh dấu tích lịch sử ô nhiễm tích tụ từ các hoạt động khai thác than và vận tải thủy rầm rộ những thập kỷ trước. Lớp giữa (5–15 cm) có nồng độ thấp hơn do quá trình bồi tụ phù sa trung gian, trong khi lớp mặt (0–5 cm) tiếp tục gia tăng nồng độ do các nguồn phát thải công nghiệp hiện hữu từ Nhà máy Nhiệt điện Mông Dương và khói bụi đô thị.

  2. Mối tương quan giữa tính chất thổ nhưỡng (TOC, pH) với khả năng tồn lưu: Phân tích hệ số tương quan Pearson cho thấy mối quan hệ thuận có ý nghĩa thống kê giữa tổng hàm lượng $\Sigma 16\text{PAHs}$ và hàm lượng carbon hữu cơ tổng số ($TOC$) trong đất RNM ($p < 0,05$). Đất rừng ngập mặn với hàm lượng mùn hữu cơ cao hoạt động như một chất hấp phụ tự nhiên giữ chặt các hợp chất kỵ nước. Ngược lại, giá trị pH đất ít có mối tương quan trực tiếp đến hàm lượng PAHs nhưng ảnh hưởng gián tiếp đến hoạt tính của hệ vi sinh vật phân hủy kỵ khí/hiếu khí tại bãi triều.

  3. Nhận diện chính xác nguồn phát thải thông qua Tỷ lệ Chẩn đoán Đồng phân:

    • Tỷ lệ $Ant/(Ant + Phe)$ phân tích được đều $> 0,1$.
    • Tỷ lệ $Flt/(Flt + Pyr)$ phần lớn nằm trong khoảng $0,40 - 0,50$ và $> 0,50$.
    • Tỷ lệ $BaA/(BaA + Chr)$ dao động trong khoảng $> 0,35$.
    • Tỷ lệ $Ind/(Ind + BghiP)$ đạt giá trị $> 0,50$. Tập hợp các chỉ số này là bằng chứng địa hóa không thể bác bỏ, khẳng định nguồn PAHs trong đất RNM Đồng Rui bắt nguồn chủ yếu từ quá trình nhiệt phân (pyrogenic source) – bao gồm đốt than đá công nghiệp, khói thải nhiệt điện, đốt gỗ rác sinh hoạt và khí thải động cơ diesel của phương tiện vận tải thủy, hoàn toàn lấn át nguồn dầu mỏ rò rỉ nguyên thủy (petrogenic).

  4. Đánh giá Rủi ro Môi trường (RQ), Rủi ro Sức khỏe (CR) và Mô phỏng Fugacity:

    • Giá trị $RQ$ của đa số các PAHs đơn lẻ nằm trong ngưỡng an toàn đến rủi ro trung bình ($0,1 \le RQ \le 1,0$), tuy nhiên một số hợp chất độc tính cao như BaP, BbF, DahA tại các điểm nóng có thời điểm tiệm cận hoặc vượt ngưỡng 1,0.
    • Chỉ số rủi ro ung thư tổng hợp $CR$ đối với cộng đồng dân cư địa phương dao động trong khoảng an toàn cho phép của US EPA ($10^{-6}$ đến $10^{-4}$), trong đó con đường nuốt phải đất ngẫu nhiên ($CR_{ing}$) và hấp thụ qua da ($CR_{dermal}$) đóng góp hơn 99% tổng nguy cơ phơi nhiễm, trong khi hít thở ($CR_{inh}$) chiếm tỷ trọng không đáng kể.
    • Kết quả mô hình Fugacity Cấp III chỉ ra rằng trên 98% lượng BaP phát thải vào hệ sinh thái Đồng Rui sẽ lắng đọng và lưu giữ cố định trong khoang đất RNM và trầm tích đáy do dung tích tích tụ $Z$ của đất lớn gấp hàng nghìn lần so với nước và khí. Mô hình Cấp IV dự báo xu thế tích lũy BaP trong đất tiếp tục tăng dần đều nếu tải lượng phát thải công nghiệp khu vực không được kiểm soát chặt chẽ.
+-----------------------------------------------------------------------------+
|                CƠ CHẾ CHUYỂN VẬN KHỐI LƯỢNG THEO FUGACITY CẤP III           |
+-----------------------------------------------------------------------------+
|                                                                             |
|                           +-------------------+                             |
|                           |   KHÍ QUYỂN (Air) |                             |
|                           |   Z_air, f_air    |                             |
|                           +-------------------+                             |
|                              /      |      \                                |
|          Lắng đọng khô/ướt  /       |       \  Lắng đọng khô/ướt            |
|          Khuếch tán        /        |        \ Khuếch tán                   |
|                           v         |         v                             |
|       +---------------------+       |       +---------------------+         |
|       |   NƯỚC MẶT (Water)  |       |       |  ĐẤT RNM (Soil)     |         |
|       |   Z_water, f_water  |       |       |  Z_soil, f_soil     |         |
|       +---------------------+       |       |  (BỒN CHỨA CHÍNH)   |         |
|                  |                  |       +---------------------+         |
|     Sa lắng hạt  |   Hấp phụ mùn    |                  |                    |
|     Khuếch tán   v                  v                  | Rửa trôi bề mặt    |
|       +-------------------------------------+          |                    |
|       |        TRẦM TÍCH (Sediment)         |<---------+                    |
|       |        Z_sediment, f_sediment       |                               |
|       +-------------------------------------+                               |
|                                                                             |
+-----------------------------------------------------------------------------+

Implications đa chiều

  • Đóng góp lý luận học thuật: Bổ sung bộ dữ liệu gốc độc nhất về địa hóa hữu cơ đất ngập mặn nhiệt đới vào kho tàng dữ liệu khu vực Đông Nam Á; kiểm chứng thành công tính thích ứng của mô hình nhiệt động học Fugacity cấp III & IV trong điều kiện địa hình phức tạp có chế độ nhật triều biên độ cao.
  • Đổi mới phương pháp luận: Thiết lập quy trình phân tích tích hợp từ quan trắc thực nghiệm GC/MS đa phân đoạn, truy nguyên nguồn gốc bằng tỷ lệ đồng phân, đến đánh giá rủi ro kép (sinh thái $RQ$ & sức khỏe $CR$) và mô hình hóa dự báo động học, có thể nhân rộng cho các dải rừng ngập mặn ven biển khác (Cần Giờ, Cà Mau, Xuân Thủy).
  • Ứng dụng thực tiễn và 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 và Sở Tài nguyên và Môi trường trong việc thực thi Quyết định phê duyệt Quy hoạch Môi trường tỉnh Quảng Ninh giai đoạn 2020 – 2030, hỗ trợ Đề á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ệ rừng ven biển ứng phó biến đổi khí hậu theo Quyết định số 120/QĐ-TTg của Thủ tướng Chính phủ.

Limitations và Future Research

Luận án thẳng thắn thừa nhận 4 giới hạn nghiên cứu:

  1. Phạm vi không gian đơn lẻ: Nghiên cứu tập trung tại xã đảo Đồng Rui, mặc dù có tính đại diện cao cho miền Bắc nhưng chưa thể bao quát toàn bộ các vùng rừng ngập mặn có đặc trưng thủy văn và thổ nhưỡng khác biệt ở miền Trung và miền Nam.
  2. Số lượng chất phân tích ưu tiên: Tập trung vào 16 hợp chất PAHs nguyên bản theo US EPA mà chưa mở rộng phân tích các dẫn xuất alkyl hóa (alkylated PAHs) hoặc các hợp chất thơm đa vòng chứa dị tố (N, S, O-PACs) có độc tính sinh học tiềm tàng.
  3. Giả định thông số môi trường trong mô hình Fugacity: Một số hệ số vận chuyển khối lượng, tốc độ phân hủy sinh học ($k$) và tải lượng phát thải đầu vào trong mô hình Fugacity phải kế thừa từ tài liệu chuẩn quốc tế hoặc tính toán gần đúng theo điều kiện khí tượng trung bình nhiều năm.
  4. Mẫu sinh vật đồng bộ: Chưa thực hiện phân tích chuỗi tích lũy sinh học (Bioaccumulation/Biomagnification) xuyên suốt các bậc dinh dưỡng từ thực vật ngập mặn (cây mắm, trang, sú) sang động vật không xương sống đáy (ngao, sò, ngao hai cùi, ốc) và cá bãi triều trong cả 6 đợt quan trắc.

Chương trình nghiên cứu tương lai (Future Research Agenda):

  • Mở rộng chuỗi tích lũy sinh học: Đánh giá chi tiết hệ số tích lũy sinh học ($BAF$) và hệ số chuyển hóa sinh học ($BSAF$) của PAHs từ đất ngập mặn sang các loài thủy sản thương mại chủ lực của Tiên Yên.
  • Tích hợp mô hình lan truyền thủy động lực học 3D: Ghép nối mô hình Fugacity với mô hình thủy lực Delft3D hoặc MIKE 21 để mô phỏng chính xác sự phát tán PAHs dưới tác động của dòng triều và sóng biển cực đoan.
  • Nghiên cứu công nghệ xử lý sinh học (Bioremediation/Rhizoremediation): Phân lập và tuyển chọn các chủng vi sinh vật nội sinh vùng rễ thực vật ngập mặn có khả năng phân hủy nhanh các PAHs cao phân tử để phục hồi các điểm nóng ô nhiễm.

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

  • Tác động học thuật: Công trình tạo nền móng vững chắc cho các nghiên cứu tiếp theo về hóa học môi trường đất ngập mặn tại Việt Nam. Phương pháp luận phân tích và bộ dữ liệu 6 đợt quan trắc thực nghiệm cung cấp tài liệu tham khảo chuẩn mực cho các chương trình đào tạo sau đại học ngành Khoa học Môi trường, Kỹ thuật Môi trường và Quản lý Đất đai.
  • Tác động chuyển đổi ngành công nghiệp: Kết quả nghiên cứu cung cấp cơ sở dữ liệu định lượng buộc các tập đoàn khai thác than (TKV) và các nhà máy nhiệt điện tại vùng Đông Bắc (Mông Dương, Cẩm Phả) phải nâng cấp hệ thống xử lý khí thải, tăng cường hiệu quả lọc bụi tĩnh điện và kiểm soát xả thải xỉ than nhằm giảm thiểu phát tán PAHs vào khí quyển.
  • Tác động quản lý nhà nước: Đóng góp trực tiếp vào cơ sở dữ liệu giám sát môi trường nền của Chi cục Bảo vệ Môi trường Quảng Ninh, hỗ trợ thiết lập hệ thống cảnh báo sớm rủi ro hóa chất độc hại tại các vùng đệm khu bảo tồn thiên nhiên.
  • Lợi ích xã hội và sinh kế: Bảo vệ môi trường đất rừng ngập mặn Đồng Rui đồng nghĩa với việc bảo vệ "vườn ươm" hải sản tự nhiên, duy trì nguồn lợi sinh kế bền vững cho hàng nghìn hộ dân địa phương sống bằng nghề khai thác nghêu, sò, cua, cá ven triều.

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

  • Nghiên cứu sinh & Giới học thuật: Tiếp cận một case study mẫu mực về ứng dụng sắc ký khí ghép khối phổ GC/MS kết hợp mô hình hóa nhiệt động học Fugacity cấp III & IV và tính toán rủi ro ung thư theo tiêu chuẩn US EPA.
  • Cơ quan Quản lý Môi trường (Bộ TN&MT, Sở TN&MT Quảng Ninh): Sở hữu luận cứ khoa học thực chứng để xây dựng quy chuẩn kỹ thuật địa phương về ngưỡng giới hạn PAHs trong đất ngập nước ven biển.
  • Doanh nghiệp Năng lượng & Khai khoáng: Nhận diện rõ trách nhiệm môi trường và dấu vết phát thải để chủ động đổi mới công nghệ sạch, chuyển đổi xanh trong sản xuất nhiệt điện và khai thác khoáng sản.
  • Cộng đồng Dân cư Ven biển Tiên Yên: Được đảm bảo môi trường sinh thái an toàn, giảm thiểu nguy cơ phơi nhiễm chất độc gây ung thư qua chuỗi thức ăn thủy hải sản.

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 mở rộng Lý thuyết Cân bằng Đa khoang Fugacity (Donald Mackay) vào một hệ sinh thái bãi triều rừng ngập mặn nhiệt đới có biến động nhật triều cực lớn (3,5 – 4,0 m). Luận án đã tham số hóa thành công độ tập trung dung nạp $Z$ và các thông số tải lượng vận chuyển liên pha $D_{ij}$ đặc thù cho vùng đất ngập nước ven biển miền Bắc Việt Nam, giải quyết bài toán định lượng phân bố chất ô nhiễm kỵ nước giữa 4 khoang môi trường mà các mô hình truyền thống chưa mô tả được.

2. Sự đổi mới trong 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 trước đây tại Việt Nam (như Nghiêm Trung Dũng, 2005 chỉ khảo sát pha khí; Dương Thanh Nghị, 2011 khảo sát tĩnh pha nước – trầm tích – sinh vật; Phạm Thị Kha, 2014 khảo sát lát cắt trầm tích đơn lẻ), luận án đã đổi mới toàn diện:

  • Quan trắc chu kỳ 6 đợt liên tục trong 4 năm kết hợp phân tích lát cắt giải phẫu 4 tầng độ sâu đất (0–5, 5–10, 10–15, 15–20 cm).
  • Tích hợp đồng thời bộ công cụ hóa học phân tích GC/MS hiện đại với hệ số chẩn đoán nguồn phân tử, đánh giá rủi ro sinh thái ($RQ$), rủi ro sức khỏe ung thư ($CR$) và mô phỏng động học nhiệt động lực học Fugacity cấp III & IV.

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

Phát hiện bất ngờ nhất là sự tồn lưu nồng độ cao của PAHs tại tầng sâu 15–20 cm ở điểm nóng ĐR5, tương đương hoặc thậm chí cao hơn tầng mặt ở một số thời điểm. Điều này phản ánh rõ nét "dấu ấn địa chất môi trường" ghi lại giai đoạn phát triển công nghiệp khai thác than và vận tải thủy cực thịnh trong quá khứ trước khi các chính sách bảo vệ rừng ngập mặn được siết chặt. Đồng thời, tỷ lệ chẩn đoán đồng phân khẳng định nguồn ô nhiễm không phải do tràn dầu tại chỗ mà do bụi lắng đốt than/nhiệt điện từ các khu vực lân cận chuyển vận qua khí quyển và dòng triều.

4. Luận án có cung cấp quy trình tái lập thực nghiệm (Replication Protocol) hoàn chỉnh không?

Hoàn toàn có. Luận án mô tả chi tiết, minh bạch từng bước quy trình: từ tọa độ mạng lưới 12 điểm lấy mẫu GPS, kỹ thuật bảo quản mẫu đất/khí/nước/trầm tích, hóa chất chuẩn nội, quy trình chiết siêu âm và cột làm sạch Silica/Alumina, chương trình cài đặt nhiệt độ và thông số ion định lượng trên GC/MS Agilent, đến toàn bộ hệ thống công thức giải tích toán học tính toán $TEQ$, $CR$, hệ phương trình ma trận phân bố $Z$ và $D$ trong mô hình Fugacity cấp III và vi phân cấp IV.

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

Chương trình 10 năm hướng tới: (1) Xây dựng bản đồ GIS số hóa 3D động học về rủi ro ô nhiễm hóa chất hữu cơ bền cho toàn bộ dải ven biển vịnh Bắc Bộ; (2) Tích hợp mô hình Fugacity với các mô hình biến đổi khí hậu để dự báo tác động của nước biển dâng và hiện tượng thời tiết cực đoan đến sự giải phóng PAHs từ "bồn chứa" đất ngập mặn; (3) Triển khai các giải pháp xử lý sinh thái dựa vào tự nhiên (Nature-based Solutions – NbS) bằng thảm thực vật ngập mặn phục hồi.

Kết luận

  1. Thiết lập bộ số liệu chuẩn mực quốc gia: Luận án đã hoàn thành cuộc điều tra toàn diện, kéo dài 4 năm (2013–2017) với 6 đợt quan trắc thực địa, cung cấp bộ dữ liệu hoàn chỉnh đầu tiên về hiện trạng 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 Đồng Rui, Tiên Yên, Quảng Ninh.
  2. Khẳng định quy luật phân bố không gian – thời gian và cơ chế hấp phụ: Xác định rõ xu thế phân bố PAHs theo chiều sâu (0–20 cm) và theo mùa khí hậu; chứng minh mối tương quan thuận có ý nghĩa thống kê giữa hàm lượng carbon hữu cơ tổng số ($TOC$) với khả năng giữ PAHs trong đất ngập mặn.
  3. Làm sáng tỏ nguồn gốc ô nhiễm pyrogenic: Sử dụng thành công hệ thống tỷ lệ đồng phân phân tử [$Ant/(Ant+Phe) > 0,1$; $Flt/(Flt+Pyr) > 0,4-0,5$; $BaA/(BaA+Chr) > 0,35$; $Ind/(Ind+BghiP) > 0,5$] để chứng minh nguồn phát thải chi phối bắt nguồn từ các quá trình đốt than, nhiệt điện và phương tiện cơ giới.
  4. Định lượng toàn diện rủi ro sinh thái và sức khỏe con người: Xác lập bức tranh rủi ro sinh thái ($RQ$) và chứng minh chỉ số rủi ro ung thư ($CR$) nằm trong giới hạn kiểm soát của US EPA, trong đó phơi nhiễm qua đường tiêu hóa và da chiếm vai trò chủ đạo.
  5. Tiên phong ứng dụng mô hình toán Fugacity cấp III và cấp IV: Mô phỏng thành công sự phân bố đa khoang và dự báo xu thế tích lũy dài hạn của hợp chất chỉ thị BaP, khẳng định đất ngập mặn và trầm tích là bồn chứa lưu giữ hơn 98% tổng lượng ô nhiễm trong hệ sinh thái.
  6. Giá trị ứng dụng thực tiễn cao: Công trình là căn cứ khoa học trực tiếp phục vụ công tác quy hoạch bảo tồn đất ngập nước Đồng Rui – Tiên Yên và hoạch định chính sách kiểm soát môi trường bền vững cho tỉnh Quảng Ninh cũng như dải ven biển Việt Nam.