Tổng quan luận án

Đề tài luận án tiến sĩ "Nghiên cứu động thái và khả năng di chuyển của một số kim loại nặng (Cu, Pb, Zn) trong đất lúa sử dụng nước tưới từ sông Nhuệ khu vực Thanh Trì, Hà Nội" được thực hiện bởi nghiên cứu sinh Chu Anh Đào, thuộc chuyên ngành Môi trường đất và nước (Mã số: 62440303), dưới sự hướng dẫn khoa học của PGS. Nguyễn Mạnh Khải và PGS. Nguyễn Ngọc Minh tại Trường Đại học Khoa học Tự nhiên, Đại học Quốc gia Hà Nội (năm 2016).

  • Tính cấp thiết và khoảng trống nghiên cứu: Ô nhiễm kim loại nặng (KLN) trong môi trường đất do các hoạt động nhân sinh đang là vấn đề môi trường nghiêm trọng, đe dọa trực tiếp đến chất lượng nông sản và sức khỏe con người qua chuỗi thức ăn. Quá trình đô thị hóa nhanh chóng làm gia tăng lượng nước thải đô thị trong khi hệ thống thoát nước bị quá tải và xuống cấp. Tại các vùng ven đô Hà Nội, do tình trạng khan hiếm nguồn nước tưới, nước thải từ hệ thống sông ngòi (như sông Tô Lịch, sông Nhuệ) được tái sử dụng để tưới cho đất canh tác lúa. Nguồn nước này chứa hàm lượng đáng kể các KLN và hợp chất ô nhiễm. Môi trường đất lúa có đặc thù ngập nước theo mùa, làm giảm sự trao đổi khí với khí quyển và tạo điều kiện khử chiếm ưu thế, dẫn đến sự chuyển hóa hóa học và động thái của KLN diễn ra rất phức tạp. Trước thời điểm nghiên cứu, các công trình trong nước chủ yếu khảo sát hàm lượng tổng số của KLN trong đất, trong khi những nghiên cứu mang tính hệ thống về động thái, sự chuyển hóa dạng liên kết và tốc độ di chuyển của Cu, Pb, Zn trong phẫu diện đất lúa dưới tác động của các yếu tố môi trường (Eh, pH, chất hữu cơ, các ion cạnh tranh) vẫn chưa được thực hiện đầy đủ.
  • Mục tiêu nghiên cứu:
    1. Đánh giá động thái, sự tích lũy, di chuyển, giải phóng và xu hướng chuyển hóa của các KLN (Cu, Pb, Zn) dưới ảnh hưởng của các điều kiện môi trường ($\text{pH}$, $\text{Eh}$, chất hữu cơ - $\text{CHC}$, $\text{SO}_4^{2-}$, $\text{PO}_4^{3-}$, $\text{Ca}^{2+}$, $\text{Mg}^{2+}$) trong điều kiện canh tác lúa sử dụng nước tưới từ sông Nhuệ.
    2. Góp phần cung cấp cơ sở khoa học trong việc tái sử dụng nước thải đô thị làm nước tưới cho nông nghiệp một cách bền vững.
  • Đối tượng và phạm vi nghiên cứu:
    • Đối tượng nghiên cứu: Nước tưới sông Nhuệ, đất trồng lúa nước, trầm tích sông và ba nguyên tố kim loại nặng chính gồm Đồng ($\text{Cu}$), Chì ($\text{Pb}$), Kẽm ($\text{Zn}$).
    • Phạm vi không gian: Khu vực huyện Thanh Trì, thành phố Hà Nội (tập trung tại lưu vực sông Nhuệ chảy qua huyện và các xã sản xuất nông nghiệp tiêu biểu như Đại Áng, Tả Thanh Oai, Hữu Hòa).
    • Phạm vi thời gian: Các số liệu và thí nghiệm được thực hiện và hoàn thành trong giai đoạn đến năm 2016.

Tổng quan tài liệu và vị trí của luận án

Tổng quan tài liệu của luận án hệ thống hóa các kết quả nghiên cứu trong và ngoài nước về nguồn gốc, sự phân bố, các cơ chế chuyển hóa, con đường di chuyển của KLN trong đất, đặc tính đất lúa nước và hiện trạng sử dụng nước tưới:

  • Nguồn gốc và phân bố của kim loại nặng:

    • Nguồn tự nhiên: Phong hóa khoáng vật đá mẹ đóng góp một lượng nhỏ KLN tự nhiên (Levinson, 1974; Alloway, 1990a; Mitchell, 1964). Các khoáng vật như olivin, anorthit, augit, hornblend, biotit giải phóng $\text{Mn}$, $\text{Co}$, $\text{Ni}$, $\text{Cu}$, $\text{Zn}$, $\text{Pb}$ vào đất qua quá trình phong hóa. Campbell và nnk (1983) chỉ ra rằng hoạt động nhân sinh tạo ra lượng KLN gấp 13 lần đối với $\text{Cu}$, 15 lần đối với $\text{Cd}$, 21 lần đối với $\text{Zn}$ và 100 lần đối với $\text{Pb}$ so với nguồn tự nhiên.
    • Nguồn khí quyển: Khí thải công nghiệp và giao thông mang sol kim loại lắng đọng khô và lắng đọng ướt vào đất (Norton, 1986; Pacyna và nnk, 1984; Nriagu và Pacyna, 1988). Lindberg và Harriss (1989) ghi nhận lượng $\text{Pb}$ lắng đọng từ khí quyển đạt từ $3{,}1\text{--}31\text{ mg/m}^2\text{/năm}$ ở vùng nông thôn và $27\text{--}140\text{ mg/m}^2\text{/năm}$ ở vùng đô thị.
    • Nguồn nông nghiệp và công nghiệp: Bùn thải, phân bón photphat, phân chuồng, vôi và thuốc trừ sâu chứa hàm lượng đáng kể $\text{Cr}$, $\text{Cu}$, $\text{Zn}$, $\text{Cd}$, $\text{Pb}$ (Alloway và Fergusson, 1990; McGrath, 1987; Lane, 1989). Chất thải luyện kim và xỉ thải cũng chứa tỷ lệ lớn $\text{Zn}$ từ $1%$ đến $75%$ (Dương Văn Đảm, 2004).
    • Phân bố theo phẫu diện: Các KLN có nguồn gốc nhân tạo tích lũy tập trung ở tầng đất mặt ($0\text{--}25\text{ cm}$) và giảm dần theo chiều sâu (Brooks, 1987, 1990; Colbourn và Thornton, 1978).
  • Cơ chế chuyển hóa và các con đường di chuyển trong đất:

    • Quá trình lý hóa học: Hodgson (1963) xác lập sơ đồ chuyển hóa KLN qua các phản ứng hòa tan, kết tủa khoáng (cacbonat, hydroxit, photphat), hấp phụ trên bề mặt khoáng sét/oxit kim loại và tạo phức với chất hữu cơ hòa tan ($\text{DOM}$) (Stevenson, 1994; Saar và Weber, 1982). Hệ số hấp phụ $K_f$ của mô hình Freundlich thể hiện sự phân bố giữa pha rắn và pha lỏng (Buchter và nnk, 1989; Sauve và nnk, 2000, 2003).
    • Quá trình sinh học: Hoạt động của vi sinh vật và hệ rễ thực vật làm biến đổi điều kiện vùng rễ qua sự giải phóng axit hữu cơ, hô hấp kỵ khí làm hòa tan hoặc cố định KLN (Kurek, 2002; Clemens và nnk, 1999; Hall, 2002).
    • Con đường di chuyển: Luận án tổng hợp 5 con đường chính: (1) Khuếch tán và phân tán phân tử (Delolme và nnk, 2004; Wilcke và nnk, 1999); (2) Dòng chảy ưu thế qua khe nứt, lỗ rỗng sinh học (Richards và nnk, 1998; Jarvis và nnk, 1999; Sterckeman và nnk, 2000); (3) Di chuyển cùng hạt keo đất (Totsche và Kögel-Knabner, 2004; Buddemeier và Hunt, 1988; Keller và nnk, 2002; Karathanasis, 1999); (4) Di chuyển dưới dạng phức chất với chất hữu cơ hòa tan $\text{DOM}$ (Senesi, 1989, 1992; Karapanagiotis và nnk, 1991); (5) Rửa trôi, xói mòn và dòng chảy bề mặt (Barrel và nnk, 1993; Ongley và nnk, 2003). Quá trình bay hơi qua methyl hóa của $\text{As}$, $\text{Hg}$, $\text{Se}$ cũng được ghi nhận (Frankenberger và Benson, 1994; Bizily và nnk, 2002).
  • Đặc tính đất lúa và các yếu tố ảnh hưởng:

    • Điều kiện ngập nước làm giảm hệ số khuếch tán oxy (tỷ số $D_{\text{nước}}/D_{\text{không khí}} = 10^{-4}$), làm hạ thấp thế oxy hóa - khử ($\text{Eh}$), tăng $\text{pH}$ đất chua, tăng độ dẫn điện (Conner, 1918; Dennett, 1932; Nguyễn Vy và Trần Khải, 1978). Quá trình khử biến đổi $\text{Fe}^{3+}$ thành $\text{Fe}^{2+}$, $\text{Mn}^{4+}$ thành $\text{Mn}^{2+}$, khử sunfat thành sunfua ($\text{S}^{2-}$), tăng độ hòa tan photphat và tích lũy $\text{NH}_4^+$ (Harrison và Aiyer, 1913).
    • Sự linh động của KLN chịu sự kiểm soát của $\text{pH}$ (Lebourg và nnk, 1998; Maskall và Thornton, 1998), thế $\text{Eh}$ (Cantwell và nnk, 2002; Plant và Raiswell, 1983), chất hữu cơ và phytolith trong thân lá lúa (Wilding, 1967; Nguyễn Ngọc Minh và nnk, 2011), khoáng sét và tầng đế cày nén chặt (Eickhorst và Tippkötter, 2009).
  • Hiện trạng nghiên cứu nước tưới sông Nhuệ và vùng ven đô:

    • Các công trình của Nguyễn Mạnh Khải và nnk (2007a, 2008), Hồ Thị Lam Trà và Egashira (2000), Nguyễn Thị Lan Hương và nnk (2006, 2007a, 2010, 2014) chỉ ra rằng sông Tô Lịch và sông Nhuệ tiếp nhận hơn $400.000\text{ m}^3\text{/ngày đêm}$ nước thải đô thị, khiến nước và trầm tích bị ô nhiễm $\text{BOD}_5$, $\text{COD}$, $\text{Cu}$, $\text{Pb}$, $\text{Zn}$, làm gia tăng hàm lượng KLN trong đất lúa sử dụng nguồn nước tưới này.
    • Khoảng trống nghiên cứu được luận án giải quyết: Luận án đi sâu giải quyết khoảng trống về cơ chế và định lượng động thái giải phóng, khả năng dịch chuyển theo phương đứng, sự biến đổi các dạng liên kết hóa học của $\text{Cu}$, $\text{Pb}$, $\text{Zn}$ trong điều kiện khử - ngập nước liên tục tại phẫu diện đất lúa huyện Thanh Trì, kết hợp giữa mô phỏng thực nghiệm phòng thí nghiệm và mô hình hóa toán học/GIS.

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

Cơ sở lý thuyết và khung phân tích

  • Lý thuyết chuyển hóa pha và động học dung dịch đất: Dựa trên mô hình chuyển hóa KLN của Hodgson (1963) và lý thuyết nhiệt động học cân bằng ion dung dịch. Sự phân bố ion kim loại tự do, phức chất hòa tan và dạng kết tủa phụ thuộc vào thế oxy hóa - khử ($\text{Eh}$), $\text{pH}$, nồng độ phối tử hữu cơ ($\text{DOM}$) và vô cơ ($\text{SO}_4^{2-}, \text{PO}_4^{3-}, \text{Cl}^-$).
  • Mô hình hấp phụ đẳng nhiệt Freundlich: $$S = K_f \cdot C^\beta$$ Trong đó $S$ là lượng kim loại hấp phụ trên pha rắn ($\text{mg/kg}$), $C$ là nồng độ kim loại cân bằng trong dung dịch ($\text{mg/L}$), $K_f$ là hệ số hấp phụ Freundlich thể hiện dung lượng hấp phụ, và $\beta$ là hệ số không thứ nguyên biểu thị cường độ hấp phụ.
  • Lý thuyết cấu trúc phẫu diện đất lúa: Phân tích sự phân tầng hình thái gồm: tầng canh tác ($A_c$, phân chia thành lớp oxy hóa mỏng phía trên có $\text{Eh} = 250\text{--}400\text{ mV}$ và lớp bùn nhão có $\text{Eh} < 200\text{ mV}$), tầng đế cày ($P$, độ sâu $10\text{--}20\text{ cm}$, bị nén chặt cơ học), tầng tích tụ ($B$, loang lổ oxit $\text{Fe-Mn}$) và tầng glây ($G$, khử mạnh, màu xám xanh).
  • Phương trình truyền dẫn khuếch tán khí: $$V = -a D \left(\frac{T}{T_0}\right)^2 \frac{dp}{dl}$$ Mô tả tốc độ khuếch tán khí qua khe rỗng đất dưới ảnh hưởng của độ bão hòa nước.

Phương pháp nghiên cứu thực tế

  • Phương pháp lấy mẫu và phân tích cơ bản:
    • Mẫu nước: Đo nhanh tại hiện trường các thông số $\text{pH}$, $\text{Eh}$, nhiệt độ. Phân tích trong phòng thí nghiệm: $\text{COD}$, $\text{BOD}_5$, $\text{NH}_4^+$, $\text{PO}_4^{3-}$, $\text{SO}_4^{2-}$, các cation kiềm/kiềm thổ ($\text{Ca}^{2+}, \text{Mg}^{2+}$) và KLN tổng số ($\text{Cu}, \text{Pb}, \text{Zn}$) bằng phương pháp quang phổ hấp thụ nguyên tử (AAS) hoặc quang phổ phát xạ plasma (ICP-OES) theo Tiêu chuẩn Việt Nam (TCVN).
    • Mẫu đất và trầm tích: Lấy mẫu theo phẫu diện thẳng đứng tại các điểm đại diện ở huyện Thanh Trì (như xã Đại Áng). Phân tích thành phần cơ giới (TPCG), dung tích trao đổi cation ($\text{CEC}$), cacbon hữu cơ tổng số ($\text{CHC}$), $\text{pH}_{\text{H}2\text{O}}$, $\text{pH}{\text{KCl}}$, nồng độ cation trao đổi và hàm lượng KLN tổng số.
  • Phương pháp phân tích dạng liên kết (phân đoạn hóa học):
    • Áp dụng quy trình chiết liên hoàn nhằm xác định các dạng liên kết của KLN trong đất và trầm tích: dạng trao đổi, dạng liên kết cacbonat, dạng liên kết oxit $\text{Fe-Mn}$, dạng liên kết hữu cơ và dạng cặn dư silicat.
  • Phương pháp mô phỏng thực nghiệm di chuyển và động thái khử:
    • Mô hình cột đất (Lysimeter phòng thí nghiệm): Thiết lập hệ thống cột đất không xáo trộn hoặc xáo trộn có kiểm soát để mô phỏng dòng thấm dung dịch nước tưới chứa $\text{Cu}, \text{Pb}, \text{Zn}$, quan sát tốc độ di chuyển và lượng rò rỉ qua các tầng đất theo thời gian.
    • Thí nghiệm kiểm soát điều kiện khử: Ủ ngập nước mẫu đất trong bình kín có kiểm soát các công thức thí nghiệm (CT), theo dõi sự biến đổi theo chuỗi thời gian của các thông số $\text{Eh}$, $\text{pH}$, $\text{Fe}^{2+}$, $\text{Mn}^{2+}$, $\text{COD}$, $\text{PO}_4^{3-}$, $\text{SO}_4^{2-}$, $\text{Ca}^{2+}$, $\text{Mg}^{2+}$ và nồng độ hòa tan của $\text{Cu}$, $\text{Pb}$, $\text{Zn}$.
  • Phương pháp xử lý số liệu và mô phỏng không gian:
    • Phân tích thống kê mô tả, phân tích tương quan Pearson, thiết lập phương trình hồi quy đa biến giữa hàm lượng kim loại nặng hòa tan với các yếu tố môi trường.
    • Ứng dụng công nghệ GIS với các thuật toán nội suy không gian (Kriging, Inverse Distance Weighting - IDW, Spline) để thành lập bản đồ phân bố tích lũy không gian của $\text{Cu}$, $\text{Pb}$, $\text{Zn}$ trong khu vực nghiên cứu.
Tiêu chí Phương pháp / Công cụ Đối tượng áp dụng Chỉ tiêu xác định
Hóa học dung dịch Máy đo điện hóa, AAS, ICP-OES Nước tưới, dung dịch đất $\text{pH}$, $\text{Eh}$, $\text{COD}$, $\text{NH}_4^+$, $\text{PO}_4^{3-}$, $\text{SO}_4^{2-}$, $\text{Cu}$, $\text{Pb}$, $\text{Zn}$
Tính chất đất/trầm tích TCVN, phương pháp Walkley-Black, ống hút Robinson Đất phẫu diện, trầm tích sông TPCG, $\text{CHC}$, $\text{CEC}$, cation trao đổi, KLN tổng số
Phân đoạn KLN Chiết liên hoàn hóa học Đất phẫu diện, trầm tích Dạng trao đổi, cacbonat, oxit $\text{Fe-Mn}$, hữu cơ, cặn dư
Mô phỏng động thái Cột đất thí nghiệm, bình ủ ngập nước Đất lúa xã Đại Áng Tốc độ thấm, hệ số $K_f$, biến thiên ion theo điều kiện khử
Không gian & Thống kê Phần mềm GIS (IDW/Kriging), phân tích hồi quy Tập dữ liệu quan trắc Thanh Trì Bản đồ phân bố KLN, phương trình tương quan đa biến

Nội dung chính theo từng chương

Chương 1. Tổng quan

Chương 1 tổng hợp có hệ thống cơ sở lý luận về nguồn gốc, động thái và mô hình chuyển hóa của KLN trong hệ sinh thái đất - nước. Tác giả phân tích chi tiết sự đóng góp của các nguồn khoáng vật tự nhiên, lắng đọng khí quyển và hoạt động nông - công nghiệp.

Loại đá mẹ $\text{Cr}$ ($\text{mg/kg}$) $\text{Mn}$ ($\text{mg/kg}$) $\text{Cu}$ ($\text{mg/kg}$) $\text{Zn}$ ($\text{mg/kg}$) $\text{Pb}$ ($\text{mg/kg}$)
Đá siêu bazơ $2000\text{--}2980$ $1040\text{--}1300$ $15\text{--}15$ $50\text{--}58$ $0{,}1\text{--}14$
Đá bazơ $200$ $1500\text{--}2200$ $90\text{--}100$ $100$ $4\text{--}24$
Đá cát kết $4$ $400\text{--}500$ $4\text{--}15$ $40\text{--}52$ $20\text{--}24$
Đá phiến sét $42\text{--}68$ $620\text{--}1100$ $5{,}5\text{--}15$ $20\text{--}25$ $5{,}7\text{--}7$
Granit $35$ $220\text{--}700$ $30$ $16\text{--}30$ $42\text{--}65$
Đá vôi $90\text{--}100$ $850$ $39\text{--}50$ $100\text{--}120$ $20\text{--}23$

(Nguồn tổng hợp trong luận án: Levinson, 1974 và Alloway, 1990a)

Trong chất thải công nghiệp, hàm lượng kẽm ($\text{Zn}$) chiếm tỷ trọng rất cao: xỉ pirit ($1\text{--}5%$), bụi lò cao luyện kim đen ($13{,}8%$), bụi lò nấu chảy đồng ($26\text{--}31%$), bụi lò nhà máy $\text{Pb-Zn}$ ($60\text{--}64%$), xỉ nổi mạ kẽm nóng ($75%$), bã thải sản xuất tơ nhân tạo ($60\text{--}67{,}5%$), và bã thải sản xuất natri hidrosunfit ($61{,}5\text{--}75%$) (Dương Văn Đảm, 2004).

Chương 1 cũng trình bày chi tiết về sự phân loại tính linh động tương đối của các kim loại trong đất theo nghiên cứu của Plant và Raiswell (1983):

Mức độ linh động tương đối Điều kiện trung tính - kiềm Điều kiện axit Điều kiện khử Điều kiện oxy hóa
Cao $\text{Zn}$ $\text{Zn}, \text{Cu}, \text{Co}, \text{Ni}, \text{Hg}, \text{Ag}, \text{Au}$ — —
Trung bình $\text{Cu}, \text{Co}, \text{Ni}, \text{Cd}$ $\text{Cd}$ — $\text{Hg}, \text{Ag}, \text{Au}$
Thấp $\text{Pb}$ $\text{Pb}$ $\text{Pb}$ —
Rất thấp – Cố định $\text{Fe}, \text{Mn}, \text{Al}, \text{Sn}, \text{Pt}, \text{Cr}, \text{Zr}$ $\text{Al}, \text{Sn}, \text{Pt}, \text{Cr}$ $\text{Al}, \text{Sn}, \text{Cr}, \text{Zn}, \text{Cu}, \text{Ni}, \text{Hg}$ $\text{Zn}, \text{Co}, \text{Cu}, \text{Ag}, \text{Ni}, \text{Hg}, \text{Ag}, \text{Au}, \text{Cd}, \text{Pb}$

Về thực trạng nông nghiệp tại huyện Thanh Trì, toàn huyện có $3.409\text{ ha}$ đất nông nghiệp; trong đó diện tích chuyên trồng rau là $1.177{,}6\text{ ha}$ (vùng bãi: $358{,}8\text{ ha}$ tại Yên Mỹ, Duyên Hà, Vạn Phúc; vùng trong đồng: $818{,}8\text{ ha}$ tại Thanh Liệt, Tam Hiệp, Liên Ninh, Vĩnh Quỳnh, Tứ Hiệp). Diện tích đất trồng lúa là $612\text{ ha}$, tập trung chủ yếu tại các xã Hữu Hòa, Tả Thanh Oai, Đại Áng. Các xã trong đồng phụ thuộc lớn vào nguồn nước dẫn từ sông Nhuệ và sông Tô Lịch qua hệ thống kênh mương và ao lắng $1\text{--}2\text{ ngày}$ trước khi bơm tưới.

Chương 2. Đối tượng, phạm vi, nội dung và phương pháp nghiên cứu

Chương 2 mô tả chi tiết địa bàn nghiên cứu tại huyện Thanh Trì, quy trình bố trí lấy mẫu và quy trình kỹ thuật phòng thí nghiệm:

  • Mô tả phẫu diện đất lúa tại các điểm lấy mẫu, ghi nhận đầy đủ các tầng hình thái $A_c, P, B, G$.
  • Quy trình phân tích tính chất lý hóa: TPCG xác định tỷ lệ cát, limon, sét; $\text{CEC}$ xác định bằng phương pháp amoni axetat ($\text{NH}_4\text{OAc}$); $\text{CHC}$ xác định bằng phương pháp Walkley-Black; các dạng anion $\text{PO}_4^{3-}, \text{SO}_4^{2-}$ và cation $\text{Ca}^{2+}, \text{Mg}^{2+}, \text{Fe}^{2+}, \text{Mn}^{2+}, \text{NH}_4^+$ được định lượng theo quy chuẩn.
  • Quy trình thực nghiệm mô phỏng cột đất trong phòng thí nghiệm nhằm xác định đường cong xuyên phá và khả năng lưu giữ của pha rắn đối với $\text{Cu}, \text{Pb}, \text{Zn}$.
  • Mô hình thí nghiệm khảo sát biến thiên thế oxy hóa - khử ($\text{Eh}$) và động thái hòa tan/kết tủa của các kim loại nặng theo thời gian ủ ngập nước trong các công thức thí nghiệm khác nhau (như CT2).

Chương 3. Kết quả nghiên cứu và thảo luận

Chương 3 trình bày toàn bộ kết quả phân tích định lượng, bao gồm:

  • Đặc tính nguồn nước tưới và trầm tích sông Nhuệ: Nước sông Nhuệ đoạn chảy qua Thanh Trì có mức độ ô nhiễm hữu cơ cao, hàm lượng $\text{COD}, \text{NH}_4^+, \text{PO}_4^{3-}$ tăng mạnh; hàm lượng $\text{Cu}, \text{Pb}, \text{Zn}$ trong nước và trầm tích đều ở mức cao, khẳng định đây là nguồn bổ sung kim loại nặng liên tục vào đất lúa. Trầm tích sông giàu chất hữu cơ và cấp hạt mịn, là nơi tích lũy kim loại dạng liên kết oxit và hữu cơ.
  • Sự tích lũy và phân bố không gian của kim loại nặng: Đất lúa Thanh Trì ghi nhận sự tích lũy rõ rệt của $\text{Cu}, \text{Pb}, \text{Zn}$ tại lớp đất mặt ($0\text{--}20\text{ cm}$). Kết quả mô phỏng bản đồ nội suy không gian bằng GIS thể hiện rõ gradient nồng độ giảm dần từ các kênh dẫn nước sông Nhuệ vào sâu trong nội đồng. Về dạng liên kết hóa học, phần lớn $\text{Pb}$ và $\text{Cu}$ liên kết chặt với chất hữu cơ và oxit $\text{Fe-Mn}$, trong khi $\text{Zn}$ có tỷ lệ dạng trao đổi và liên kết cacbonat cao hơn, làm tăng khả năng linh động của kẽm.
  • Khả năng cố định và tốc độ di chuyển trong đất:
    • Hằng số hấp phụ đẳng nhiệt Freundlich ($K_f$ và $\beta$) của các kim loại trên pha rắn đất lúa cho thấy ái lực hấp phụ tuân theo thứ tự: $\text{Pb} > \text{Cu} > \text{Zn}$.
    • Kết quả thử nghiệm trên mẫu đất xã Đại Áng chứng minh tầng canh tác giữ lại phần lớn lượng kim loại; tuy nhiên khi đạt trạng thái bão hòa hấp phụ, $\text{Zn}$ và $\text{Cu}$ bắt đầu dịch chuyển xuống tầng dưới qua khe nứt và dòng thấm.
    • Luận án đối chiếu giữa mô hình lý thuyết và kết quả thực nghiệm phòng thí nghiệm, đồng thời biện luận chi tiết các nguyên nhân gây sai số (sự phá hủy cấu trúc đoàn lạp khi đưa mẫu vào cột, hiện tượng dòng chảy thành ống trong điều kiện tự nhiên).
  • Ảnh hưởng của điều kiện khử đến động thái kim loại nặng:
    • Khi ngập nước kéo dài, thế $\text{Eh}$ đất giảm sâu (xuống dưới $-150\text{ mV}$), $\text{pH}$ dịch chuyển về khoảng trung tính ($6{,}5\text{--}7{,}2$).
    • Sự suy giảm $\text{Eh}$ thúc đẩy quá trình khử $\text{Fe}^{3+} \rightarrow \text{Fe}^{2+}$ và $\text{Mn}^{4+} \rightarrow \text{Mn}^{2+}$, làm nồng độ $\text{Fe}^{2+}$ và $\text{Mn}^{2+}$ hòa tan trong dung dịch đất tăng mạnh ở giai đoạn đầu, giải phóng các ion $\text{Cu}^{2+}, \text{Zn}^{2+}, \text{Pb}^{2+}$ từng bị hấp phụ trên oxit $\text{Fe-Mn}$.
    • Ở giai đoạn ngập nước sâu tiếp theo, vi khuẩn khử sunfat hoạt động tạo $\text{S}^{2-}$, dẫn đến sự hình thành kết tủa sunfua kim loại ($\text{CuS}, \text{PbS}, \text{ZnS}$) làm nồng độ kim loại hòa tan trong dung dịch đất giảm xuống mức rất thấp.
    • Luận án thiết lập ma trận hệ số tương quan Pearson và xây dựng các phương trình hồi quy tuyến tính mô tả mối quan hệ định lượng giữa nồng độ $\text{Cu}, \text{Pb}, \text{Zn}$ hòa tan với các thông số $\text{Eh}, \text{pH}, \text{Fe}^{2+}, \text{SO}_4^{2-}, \text{PO}_4^{3-}, \text{Ca}^{2+}, \text{Mg}^{2+}$.

Kết quả và những đóng góp mới

  • Những đóng góp mới về mặt lý luận:
    • Làm sáng tỏ động thái biến đổi và mối tương quan định lượng giữa các dạng tồn tại, tính tan và sự giải phóng của $\text{Cu}$, $\text{Pb}$, $\text{Zn}$ với các yếu tố môi trường cốt lõi ($\text{Eh}$, $\text{pH}$, $\text{CHC}$, $\text{SO}_4^{2-}$, $\text{PO}_4^{3-}$, $\text{Ca}^{2+}$, $\text{Mg}^{2+}$) trong môi trường đất lúa ngập nước chịu tác động của nước thải đô thị.
    • Đánh giá một cách có hệ thống cơ chế phân bố pha và chuyển hóa dạng liên kết (trao đổi, liên kết oxit $\text{Fe-Mn}$, liên kết hữu cơ, cacbonat, cặn dư) của kim loại nặng trong đất lúa dưới sự biến thiên liên tục của thế oxy hóa - khử.
    • Xác lập các thông số hấp phụ thực nghiệm ($K_f, \beta$) và phương trình hồi quy mô tả động thái di chuyển của KLN trong phẫu diện đất lúa ven đô.
  • Những đóng góp mới về mặt thực tiễn:
    • Cung cấp bộ dữ liệu thực nghiệm toàn diện về hiện trạng chất lượng nước tưới, trầm tích sông Nhuệ và mức độ tích lũy kim loại nặng tại tầng đất mặt cũng như tốc độ thấm xuống các tầng sâu qua các chu kỳ canh tác lúa tại huyện Thanh Trì.
    • Cung cấp luận cứ khoa học thực chứng giúp các nhà quản lý môi trường và nông nghiệp đánh giá mức độ rủi ro, kiểm soát nguồn thải và xây dựng quy trình quản lý việc tái sử dụng nước sông Nhuệ phục vụ sản xuất nông nghiệp an toàn.
  • Kiến nghị và giải pháp đề xuất:
    • Cần kiểm soát chặt chẽ và xử lý sơ bộ nguồn nước sông Nhuệ trước khi dẫn vào hệ thống thủy nông nội đồng; duy trì thời gian lắng đọng tối thiểu tại các ao hồ trung chuyển để loại bỏ cặn lơ lửng chứa kim loại nặng.
    • Điều tiết chế độ ngập - rút nước hợp lý trong canh tác lúa nhằm kiểm soát thế oxy hóa - khử ($\text{Eh}$) và $\text{pH}$ đất, tránh hiện tượng giải phóng ồ ạt kim loại nặng linh động vào dung dịch đất trong các giai đoạn sinh trưởng nhạy cảm của cây trồng.
    • Tăng cường giám sát định kỳ hàm lượng kim loại nặng tích lũy tại tầng canh tác ở các xã trọng điểm (như Đại Áng, Tả Thanh Oai) để có biện pháp xử lý đất kịp thời khi nồng độ tiếp cận ngưỡng giới hạn quy chuẩn.

Hạn chế và hướng nghiên cứu tiếp

  • Hạn chế của luận án:
    • Phạm vi khảo nghiệm thực tế tập trung chủ yếu trên địa bàn huyện Thanh Trì, đại diện cho vùng đất phù sa ven đô sông Nhuệ, nên việc áp dụng các hệ số thực nghiệm cho các loại đất lúa có tính chất cơ lý khác biệt cần được hiệu chỉnh.
    • Mô hình mô phỏng cột đất trong phòng thí nghiệm chưa phản ánh hoàn toàn tính dị hướng và sự phức tạp của hệ thống lỗ rỗng sinh học (hang giun, khe nứt nẻ tự nhiên) ngoài đồng ruộng.
  • Hướng nghiên cứu tiếp:
    • Mở rộng nghiên cứu di chuyển kim loại nặng trên mô hình lysimeter quy mô hiện trường dài hạn qua nhiều vụ canh tác liên tiếp.
    • Khảo sát bổ sung các nguyên tố vết độc hại khác như Cadmi ($\text{Cd}$), Asen ($\text{As}$), Thủy ngân ($\text{Hg}$) và đánh giá định lượng sự tích lũy sinh học trong từng bộ phận của cây lúa (thân, lá, hạt gạo).

Giá trị tham khảo

  • Đối tượng hưởng lợi:
    • Nghiên cứu sinh, học viên cao học và sinh viên: Chuyên ngành Khoa học Môi trường, Môi trường Đất và Nước, Thổ nhưỡng học, Quản lý Tài nguyên Nước.
    • Giảng viên và nhà nghiên cứu: Sử dụng làm tài liệu tham khảo chuyên sâu về hóa học đất ngập nước, động thái oxy hóa - khử và phương pháp mô phỏng thực nghiệm di chuyển chất ô nhiễm.
    • Cán bộ quản lý môi trường và nông nghiệp: Các cơ quan quản lý tại Hà Nội và lưu vực sông Nhuệ - Đáy trong việc hoạch định chính sách sử dụng nước tưới và bảo vệ môi trường đất canh tác.
  • Các phần nội dung đáng tham khảo nhất:
    • Phần phương pháp bố trí thí nghiệm cột đất lysimeter và phương pháp ủ yếm khí kiểm soát $\text{Eh-pH}$ tại Chương 2.
    • Hệ thống bảng số liệu phân tích dạng liên kết kim loại nặng và các phương trình hồi quy tương quan đa biến tại Chương 3.
    • Tổng hợp lý thuyết về các quá trình hóa học, sinh học và các con đường di chuyển của kim loại nặng trong đất ngập nước tại Chương 1.

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

1. Yếu tố nào đóng vai trò chính kiểm soát khả năng hòa tan và di chuyển của $\text{Cu}, \text{Pb}, \text{Zn}$ trong đất lúa ngập nước theo luận án?
Thế oxy hóa - khử ($\text{Eh}$) và độ chua ($\text{pH}$) là hai thông số quyết định nhất. Khi đất ngập nước, $\text{Eh}$ giảm thúc đẩy quá trình khử sắt/mangan giải phóng kim loại bị hấp phụ, sau đó nếu $\text{Eh}$ giảm sâu, phản ứng khử sunfat tạo kết tủa sunfua ($\text{CuS}, \text{PbS}, \text{ZnS}$) lại làm giảm tính tan của kim loại.

2. Tầng đế cày ($P$) có vai trò gì đối với sự di chuyển của kim loại nặng trong phẫu diện đất lúa?
Tầng đế cày nằm ở độ sâu $10\text{--}20\text{ cm}$, được hình thành do canh tác lâu đời và nén chặt cơ học, có hàm lượng sét cao và độ thấm nước kém. Tầng này đóng vai trò như một lớp ngăn cản cơ học và hóa lý, làm chậm quá trình thấm và giữ lại phần lớn kim loại nặng ở tầng canh tác mặt ($A_c$).

3. Phytolith trong đất trồng lúa có tác động như thế nào đến tính linh động của kim loại nặng?
Cây lúa tích lũy silic tạo thành các thể phytolith trong thân lá. Phytolith có dung tích trao đổi cation ($\text{CEC}$) cao (xấp xỉ $160\text{ cmol(c+)/kg}$) và khả năng hấp phụ chất hữu cơ cao (axit humic đạt $15{,}25\text{ g/kg}$), giúp hấp phụ kép cation kim loại và anion hữu cơ, từ đó làm giảm tính linh động của kim loại nặng trong đất.

4. Nguồn nước sông Nhuệ được người dân huyện Thanh Trì khai thác và sử dụng tưới cho lúa theo quy trình nào?
Nước từ sông Nhuệ và sông Tô Lịch được các trạm bơm đưa vào các ao hồ chứa hoặc hệ thống kênh mương nội đồng, để lắng tự nhiên từ $1$ đến $2$ ngày nhằm giảm bớt cặn lơ lửng trước khi dẫn vào ruộng tưới phục vụ sản xuất.

5. Thứ tự ái lực hấp phụ của ba kim loại $\text{Cu}, \text{Pb}, \text{Zn}$ trên pha rắn đất lúa Thanh Trì diễn ra như thế nào?
Theo kết quả xác định hằng số hấp phụ Freundlich ($K_f$), ái lực hấp phụ của pha rắn đất đối với các kim loại tuân theo thứ tự giảm dần: $\text{Pb} > \text{Cu} > \text{Zn}$, trong đó $\text{Pb}$ liên kết bền vững nhất và $\text{Zn}$ có tính linh động cao nhất.

Kết luận

Luận án của Chu Anh Đào (2016) đã khảo sát và định lượng hóa một cách toàn diện động thái, sự phân bố dạng liên kết và khả năng di chuyển của $\text{Cu}$, $\text{Pb}$, $\text{Zn}$ trong đất lúa sử dụng nước tưới sông Nhuệ tại huyện Thanh Trì, Hà Nội. Bằng việc kết hợp giữa phân tích hóa học phẫu diện, mô hình hấp phụ Freundlich, mô phỏng cột đất thực nghiệm và mô hình hóa GIS, công trình đã xác lập mối quan hệ định lượng giữa nồng độ kim loại hòa tan với các biến số $\text{Eh}$, $\text{pH}$, $\text{CHC}$ và các ion trong điều kiện khử. Kết quả nghiên cứu cung cấp cơ sở dữ liệu thực nghiệm và luận cứ khoa học quan trọng cho công tác quản lý môi trường, kiểm soát ô nhiễm và định hướng sử dụng bền vững nguồn nước tưới ven đô.