Tổng quan nghiên cứu

Ô nhiễm kim loại nặng trong đất nông nghiệp đang là thách thức nghiêm trọng đối với an ninh lương thực và sức khỏe cộng đồng toàn cầu. Theo các thống kê địa hóa học, hoạt động nhân tạo đã phát thải lượng kim loại nặng vào môi trường vượt trội so với tự nhiên, cụ thể gấp khoảng 13 lần đối với đồng (Cu), 21 lần đối với kẽm (Zn) và tới 100 lần đối với chì (Pb). Tại vùng đồng bằng sông Hồng, việc thâm canh lúa nước kết hợp sử dụng nguồn nước tưới tiêu từ các hệ thống sông tiếp nhận nước thải như sông Nhuệ đã làm gia tăng nguy cơ tích lũy các độc chất này trong đất canh tác.

Đất lúa có chế độ ngập nước định kỳ, hình thành môi trường khử đặc thù khiến thế oxy hóa - khử (Eh) sụt giảm mạnh và biến đổi độ pH của môi trường. Điều này tác động trực tiếp đến sự hòa tan, chuyển hóa hóa trị và khả năng giải phóng của các cation kim loại. Luận văn Thạc sĩ Khoa học Môi trường của tác giả Phạm Vy Anh tại Trường Đại học Khoa học Tự nhiên - Đại học Quốc gia Hà Nội tập trung nghiên cứu ảnh hưởng của môi trường khử đến sự giải phóng các kim loại nặng gồm Cu, Pb và Zn trong đất trồng lúa tại xã Đại Áng, huyện Thanh Trì, Hà Nội.

Nghiên cứu tập trung vào 3 phẫu diện đất phù sa có tầng glây (tuyển chọn từ 75 mẫu đất mặt khảo sát ở độ sâu 0 - 25 cm trên địa bàn huyện Thanh Trì). Phạm vi nghiên cứu đánh giá biến thiên của thế Eh từ -300 mV đến +700 mV, theo dõi động thái hòa tan của Fe(II), Mn(II) cùng các kim loại vết theo thời gian. Ý nghĩa cốt lõi của công trình là cung cấp luận cứ khoa học thực nghiệm giúp các nhà quản lý nông nghiệp kiểm soát quá trình giải phóng độc chất, hướng tới giảm thiểu từ 15% đến 30% nguy cơ tích lũy kim loại nặng vào chuỗi thức ăn thông qua chế độ tưới tiêu hợp lý.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Nghiên cứu vận dụng hệ thống lý thuyết sinh địa hóa đất ướt kết hợp các mô hình nhiệt động học và động học giải phóng vi lượng:

  • Mô hình nhiệt động học phản ứng oxy hóa - khử: Dựa trên phương trình Nernst ($Eh = E_0 - \frac{0.059}{n}\log\frac{[Rd]}{[Ox]} - 0.059\frac{m}{n}pH$) để định lượng mối liên hệ giữa thế oxy hóa - khử (Eh), độ hoạt động của các cặp redox ($O_2/H_2O, NO_3^-/N_2, Mn^{4+}/Mn^{2+}, Fe^{3+}/Fe^{2+}, SO_4^{2-}/S^{2-}$) và giá trị pH đất.
  • Mô hình động học giải phóng hai pha (Two-step kinetic release model): Cải biên từ lý thuyết của Cappuyns và Swennen (2008), phân chia quá trình giải phóng kim loại thành hai hệ đệm độc lập gồm hệ giải phóng nhanh ($RC_1$, hằng số tốc độ $r_1$) và hệ giải phóng chậm ($RC_2$, hằng số tốc độ $r_2$).
  • Phương trình hấp phụ đẳng nhiệt Freundlich: Mô tả sự tương tác phân bố kim loại giữa pha rắn và pha lỏng qua hệ số phân bố $K_d$ và tham số mũ $n$, phản ánh độ bão hòa bề mặt hấp phụ của khoáng sét và chất hữu cơ.
  • Các khái niệm nền tảng: Thế oxy hóa - khử (Eh), dung tích trao đổi cation (CEC), chất hữu cơ hòa tan (DOM), và quá trình cố định tại chỗ (in-situ immobilization).

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

  • Cỡ mẫu và phương pháp chọn mẫu: Thu thập 03 mẫu đất đại diện (ký hiệu ĐA1, ĐA2, ĐA3) được chọn lọc từ mạng lưới 75 điểm lấy mẫu ở độ sâu 0 - 25 cm trên toàn huyện Thanh Trì. Tiêu chí lựa chọn dựa trên 3 điều kiện nghiêm ngặt: vị trí có nồng độ kim loại nặng tích lũy cao, thuộc vùng trũng ngập nước quanh năm và sử dụng trực tiếp nguồn nước tưới từ sông Nhuệ.
  • Phương pháp phân tích lý hóa đất:
    • Đo $pH_{KCl}$ trong huyền phù tỷ lệ 1:5; xác định dung tích trao đổi cation (CEC) bằng phương pháp amoni axetat (Schachtschabel); định lượng chất hữu cơ (CHC) bằng phương pháp Walkley - Black.
    • Phân tích thành phần cấp hạt theo phương pháp ống hút Robinson (định luật Stokes).
    • Phân tích khoáng sét bằng nhiễu xạ tia X (XRD) trên máy Siemens D5005 với 4 trạng thái xử lý: bão hòa $K^+$, nung $550^\circ C$, bão hòa $Mg^{2+}$ và solvat hóa bằng etylen glycol nhằm phân biệt chính xác các khoáng sét nhóm 2:1 và 1:1.
  • Phương pháp chiết liên tiếp 5 bước (Tessier cải tiến): Tách phân đoạn kim loại thành 5 dạng liên kết: dạng trao đổi (F1), dạng liên kết cacbonat (F2), dạng liên kết oxit Fe-Mn (F3), dạng liên kết hữu cơ (F4) và dạng cặn dư (F5).
  • Thí nghiệm mô phỏng môi trường khử nhân tạo: Thiết lập tỷ lệ đất : nước = 2g : 40ml trong ống ly tâm kín, đặt trong buồng tối ở nhiệt độ phòng ($25^\circ C$). Quá trình lấy mẫu định kỳ được thực hiện nghiêm ngặt trong lồng khí trơ liên tục sục khí $N_2$ để ngăn ngừa oxy khí quyển xâm nhập. Đo thế Eh và pH bằng điện cực chuyên dụng Orion; định lượng nồng độ Cu, Pb, Zn, Fe, Mn hòa tan bằng phương pháp quang phổ hấp thụ nguyên tử ngọn lửa (AAS). Dữ liệu được xử lý thống kê ANOVA và tương quan Pearson trên phần mềm MINITAB 8.0 với độ tin cậy $p < 0.01$ trong timeline nghiên cứu kéo dài 12 tháng.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

  • Đặc tính lý hóa và thành phần khoáng: Mẫu đất nghiên cứu có phản ứng chua nhẹ đến gần trung tính ($pH_{KCl} = 5.82 - 6.12$), dung tích trao đổi cation cao ($CEC = 15.90 - 18.04\text{ cmol/kg}$), hàm lượng chất hữu cơ giàu ($CHC = 2.18% - 2.60%$). Cấp hạt thuộc loại đất thịt sét đến thịt pha limon với hàm lượng sét chiếm 22% - 43% và limon chiếm 50% - 61%. Nhiễu xạ tia X khẳng định khoáng sét chủ yếu gồm illit > kaolinit > clorit, là các khoáng không trương nở có khả năng giữ nước và hấp phụ cation ở mức trung bình đến thấp.
  • Biến thiên của thế Eh và sự hòa tan Fe, Mn: Khi ngập nước, thế Eh giảm sâu nhanh chóng từ mức oxy hóa ban đầu (+400 mV) xuống điều kiện khử mạnh (-150 mV đến -200 mV). Quá trình này kích hoạt vi sinh vật yếm khí sử dụng $Fe^{3+}$ và $Mn^{4+}$ làm chất nhận điện tử, làm nồng độ ion $Fe(II)$ và $Mn(II)$ hòa tan trong dung dịch tăng vọt gấp 5 đến 8 lần so với thời điểm ban đầu.
  • Động thái giải phóng kim loại nặng (Cu, Pb, Zn):
    • Đồng (Cu) và Chì (Pb) ban đầu được giải phóng nhẹ do sự hòa tan của màng oxit Fe-Mn (pha F3), nhưng sau đó giảm nhanh và duy trì ở mức rất thấp (dưới 5% tổng lượng giải phóng) do phản ứng kết tủa với ion sunfua hình thành $Cu_2S, CuS$ và $PbS$ khó tan.
    • Kẽm (Zn) thể hiện động thái giải phóng linh động hơn rõ rệt (chiếm khoảng 12% - 25% lượng trao đổi), nồng độ hòa tan duy trì kéo dài trong pha lỏng do ái lực tạo phức sunfua yếu hơn và tham số hấp phụ Freundlich $n < 1$.

Thảo luận kết quả

Cơ chế điều khiển động thái kim loại trong đất lúa Đại Áng chịu sự chi phối tổng hợp của phản ứng chuyển pha khoáng và cân bằng thế Eh-pH. Khi môi trường chuyển sang trạng thái yếm khí sâu, việc khử các hợp chất sắt và mangan vô định hình làm mất đi các bề mặt hấp phụ tự nhiên của đất. Tuy nhiên, sự xuất hiện của phản ứng khử sunfat thành sunfua ($S^{2-}$) đã tạo ra cơ chế bẫy hóa học thứ cấp, cố định bền vững các cation $Pb^{2+}$ và $Cu^{2+}$.

Trong các báo cáo khoa học, các dữ liệu này được thể hiện tối ưu thông qua đồ thị tương quan Eh-pH (giản đồ Pourbaix) kết hợp với biểu đồ cột phân bố 5 phân đoạn chiết Tessier. Điều này làm nổi bật sự dịch chuyển của kim loại từ pha trao đổi sang pha liên kết sunfua bền vững. Kết quả này tương đồng với các công bố quốc tế của Kelderman & Osman (2007) trên trầm tích ngập nước và Ho et al. (2014) trên trầm tích sông Cam. Phát hiện khẳng định rằng việc duy trì đất ở trạng thái ngập nước liên tục sẽ hạn chế đáng kể sự hòa tan của Pb và Cu, nhưng nếu đất trải qua chu kỳ khô - ướt xen kẽ, quá trình tái oxy hóa sẽ kích thích giải phóng ồ ạt lượng lớn Zn và Cu vào dung dịch đất.

Đề xuất và khuyến nghị

  1. Điều tiết chế độ tưới tiêu ngập - khô xen kẽ (AWD) có kiểm soát: Ủy ban nhân dân xã Đại Áng và Hợp tác xã Nông nghiệp cần hướng dẫn nông dân duy trì mực nước ngập nông liên tục từ 3 - 5 cm trong suốt giai đoạn đẻ nhánh đến trổ bông để ổn định thế Eh ở mức khử nhẹ (-100 mV đến 0 mV), giúp giảm trên 35% nguy cơ hòa tan của Zn và Cd mà không làm phát sinh độc chất sunfua gây thối rễ lúa.
  2. Áp dụng kỹ thuật cố định kim loại nặng tại chỗ bằng vật liệu khoáng: Sở Nông nghiệp & Phát triển Nông thôn Hà Nội triển khai mô hình bổ sung vôi bột nông nghiệp ($CaCO_3$) với liều lượng 400 - 500 kg/ha kết hợp 200 kg/ha phân lân nung chảy giàu photphat (hydroxyapatit) trước mỗi vụ gieo cấy. Biện pháp này giúp nâng pH đất lên ngưỡng 6.5 - 7.0, giảm hơn 40% độ linh động của các ion kim loại tự do thông qua phản ứng đồng kết tủa pyromorphit.
  3. Quan trắc và xử lý sơ bộ nguồn nước tưới từ sông Nhuệ: Chi cục Bảo vệ Môi trường Hà Nội định kỳ lấy mẫu giám sát chất lượng nước tưới thủy nông 1 tháng/lần; xây dựng hệ thống kênh lắng sinh học có trồng thực vật thủy sinh (cỏ vetiver, bèo tây) tại các đầu mối trạm bơm nhằm giảm thiểu 50% hàm lượng cặn lơ lửng (TSS) và kim loại nặng trước khi đưa vào đồng ruộng.
  4. Tăng cường hàm lượng mùn hữu cơ bền vững cho đất: Khuyến cáo nông dân bổ sung 8 - 10 tấn phân chuồng ủ hoai mục kết hợp chế phẩm vi sinh mỗi năm giai đoạn 2026 - 2030, duy trì hàm lượng chất hữu cơ đất luôn trên mức 2.5%, nâng cao dung tích hấp phụ CEC của đất và hạn chế sự di chuyển của kim loại nặng xuống tầng nước ngầm.

Đối tượng nên tham khảo luận văn

  • Cán bộ quản lý môi trường và nông nghiệp (Sở TN&MT, Sở NN&PTNT): Sử dụng các số liệu thực nghiệm về hệ số phân bố $K_d$ và biến thiên Eh-pH để xây dựng bản đồ phân vùng rủi ro ô nhiễm kim loại nặng và thiết lập quy chuẩn kỹ thuật canh tác an toàn cho các vùng đất trũng ven đô.
  • Giảng viên và nhà nghiên cứu chuyên ngành Khoa học Môi trường, Thổ nhưỡng học: Kế thừa phương pháp thiết lập hệ thí nghiệm khử kín khí sục $N_2$ và mô hình động học giải phóng hai pha làm tài liệu phục vụ giảng dạy sau đại học và phát triển các đề tài sinh địa hóa chuyên sâu.
  • Học viên cao học và sinh viên đại học khối ngành Môi trường, Hóa học, Nông học: Tham khảo quy trình phân tích chiết liên tiếp 5 bước Tessier chuẩn hóa, phương pháp định danh khoáng sét bằng nhiễu xạ tia X (XRD) và kỹ năng xử lý số liệu thống kê ANOVA trên phần mềm chuyên dụng.
  • Kỹ sư nông nghiệp và cán bộ khuyến nông cơ sở: Ứng dụng các quy luật chuyển hóa của kim loại trong điều kiện yếm khí để tối ưu hóa quy trình bón phân lân, vôi và lịch trình rút nước phơi ruộng phù hợp với từng thời kỳ sinh trưởng của cây lúa.

Câu hỏi thường gặp

Môi trường khử trong đất lúa ảnh hưởng như thế nào đến độ linh động của chì (Pb)?

Khi đất lúa ngập nước làm thế oxy hóa - khử (Eh) giảm sâu xuống dưới -100 mV, các vi sinh vật kỵ khí khử sunfat thành sunfua ($S^{2-}$). Ion $Pb^{2+}$ tự do sẽ phản ứng ngay với $S^{2-}$ để tạo thành khoáng chì sunfua ($PbS$ - galen) có độ tan cực kỳ thấp. Quá trình này giúp cố định trên 85% lượng chì trong đất, làm giảm đáng kể khả năng hấp thu của rễ lúa.

Tại sao kẽm (Zn) lại có khả năng giải phóng cao hơn đồng (Cu) và chì (Pb) trong điều kiện khử?

Kẽm có ái lực tạo phức với phối tử sunfua và chất hữu cơ yếu hơn nhiều so với đồng và chì. Đồng thời, tham số đẳng nhiệt Freundlich của Zn thường có $n < 1$, đồng nghĩa với việc ái lực hấp phụ của Zn trên bề mặt pha rắn sẽ giảm dần khi nồng độ dung dịch tăng. Khi ion $Fe(II)$ và $Mn(II)$ hòa tan cạnh tranh vị trí trao đổi, Zn dễ dàng bị đẩy ra dung dịch đất với tỷ lệ giải phóng đạt từ 12% đến 25%.

Thành phần khoáng sét trong đất lúa xã Đại Áng đóng vai trò gì đối với sự hấp phụ kim loại nặng?

Phân tích XRD chỉ ra khoáng sét tại khu vực nghiên cứu gồm illit, kaolinit và clorit với dung tích CEC thấp (8 - 30 cmol/kg). Do cấu trúc không trương nở và mật độ điện tích bề mặt hạn chế, pha khoáng sét đóng vai trò thứ yếu trong việc giữ giữ kim loại. Khả năng hấp phụ và giải phóng kim loại nặng của đất phụ thuộc chủ yếu vào hàm lượng chất hữu cơ (2.18% - 2.60%) và hệ màng oxit sắt - mangan vô định hình.

Nước tưới sông Nhuệ tác động ra sao đến hiện trạng ô nhiễm kim loại nặng tại khu vực nghiên cứu?

Nước sông Nhuệ tiếp nhận lượng lớn nước thải công nghiệp và sinh hoạt chưa qua xử lý, mang theo hàm lượng đáng kể các ion Cu, Pb, Zn và chất hữu cơ hòa tan. Việc sử dụng nguồn nước này tưới tiêu liên tục qua nhiều thập kỷ đã làm gia tăng nồng độ tích lũy kim loại nặng ở tầng đất mặt (0 - 25 cm) lên gấp 2 đến 4 lần so với các vùng đất lúa sử dụng nguồn nước sạch đối chứng.

Điểm ưu việt của phương pháp chiết liên tiếp Tessier so với phương pháp phân tích tổng số kim loại là gì?

Phân tích tổng số chỉ xác định được tổng lượng kim loại hiện diện mà không phản ánh được mức độ linh động hay độc tính sinh học. Phương pháp chiết liên tiếp Tessier bóc tách kim loại thành 5 phân đoạn hóa học riêng biệt (từ dạng trao đổi linh động F1 đến dạng cặn dư trơ F5), cho phép đánh giá chính xác phần kim loại có khả năng giải phóng gây độc cho cây trồng khi điều kiện môi trường (pH, Eh) thay đổi.

Kết luận

  • Luận văn đã xác lập định lượng mối quan hệ giữa thế oxy hóa - khử ($Eh = -200\text{ mV}$ đến $+400\text{ mV}$), pH ($5.82 - 6.12$) và động thái giải phóng của ba kim loại nặng Cu, Pb, Zn trong đất lúa phù sa sông Hồng.
  • Xác định cấu trúc khoáng sét illit - kaolinit - clorit có dung tích $CEC = 15.90 - 18.04\text{ cmol/kg}$, chứng minh hàm lượng chất hữu cơ và oxit Fe-Mn là hai yếu tố kiểm soát chính đến khả năng hấp phụ kim loại của đất Đại Áng.
  • Làm rõ cơ chế môi trường khử sâu giúp cố định Cu và Pb dưới dạng sunfua khó tan ($CuS, PbS$), nhưng làm tăng độ hòa tan của Fe, Mn và duy trì tính linh động cao của Zn.
  • Chuẩn hóa thành công quy trình thực nghiệm mô phỏng khử nhân tạo trong buồng khí trơ $N_2$ kết hợp phân tích chiết phân đoạn Tessier với độ tin cậy thống kê cao ($p < 0.01$).
  • Đóng góp luận cứ khoa học thực tiễn phục vụ quy hoạch canh tác lúa an toàn và kiểm soát ô nhiễm đất vùng tưới tiêu sông Nhuệ giai đoạn 2026 - 2030.

Trong giai đoạn tiếp theo (2026 - 2028), các nghiên cứu cần tiếp tục mở rộng đánh giá sự tích lũy sinh học của kim loại nặng trực tiếp trên hạt lúa và thử nghiệm các vật liệu cố định biochar quy mô đồng ruộng. Quý bạn đọc và các nhà nghiên cứu quan tâm có thể liên hệ thư viện chuyên ngành để khai thác trọn vẹn tài liệu khoa học giá trị này!