Tổng quan nghiên cứu

Khai thác khoáng sản là ngành kinh tế then chốt nhưng cũng tiềm ẩn rủi ro phát tán kim loại nặng vào môi trường đất, làm suy giảm đa dạng sinh học và đe dọa chuỗi thức ăn sinh thái. Tại thị trấn Bằng Lũng, huyện Chợ Đồn, tỉnh Bắc Kạn – khu vực sở hữu mỏ chì kẽm Lủng Váng và mỏ Bằng Lũng với trữ lượng quặng chì dao động từ 3,71% đến 9,5%, áp lực ô nhiễm đất đang ở mức báo động. Nhóm động vật thân mềm chân bụng trên cạn (Gastropoda) với đặc tính di chuyển chậm, gắn bó chặt chẽ với tầng thảm mục và đất mặt, trở thành đối tượng sinh vật chỉ thị lý tưởng cho các biến động hóa sinh thổ nhưỡng.

Nghiên cứu tập trung giải quyết bài toán định lượng mức độ phơi nhiễm kim loại nặng thông qua việc đánh giá cấu trúc quần xã sinh vật chỉ thị. Mục tiêu cụ thể gồm: xác định hàm lượng chì (Pb) và asen (As) tại 30 điểm thu mẫu đất thuộc 5 sinh cảnh; thống kê thành phần loài và các chỉ số đa dạng sinh học của ốc cạn; từ đó xây dựng các mô hình tương quan và phương trình hồi quy phản ánh mối quan hệ giữa độc chất môi trường và cấu trúc quần xã sinh vật.

Phạm vi thực địa được triển khai trong giai đoạn từ tháng 02 đến tháng 07 năm 2017 trên tổng diện tích tự nhiên 2.491 ha của thị trấn Bằng Lũng. Công trình mang lại ý nghĩa thực tiễn to lớn khi cung cấp bộ dữ liệu nền với 304 cá thể ốc cạn được định loại và phát hiện 83,33% số mẫu đất vượt ngưỡng chì quy định. Đây là cơ sở khoa học then chốt giúp các nhà quản lý môi trường thiết lập hệ thống quan trắc sinh học chi phí thấp, thay thế dần các phương pháp phân tích hóa lý đắt đỏ.

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 được xây dựng trên nền tảng lý thuyết sinh thái học quần xã và độc học môi trường, trọng tâm là lý thuyết tích lũy sinh học (Bioaccumulation) và chỉ thị sinh học (Bioindication). Lý thuyết sinh thái học quần xã khẳng định cấu trúc thành phần loài và mật độ cá thể phản ánh trực tiếp chất lượng sinh cảnh, trong đó các xáo trộn môi trường sẽ dẫn tới sự suy giảm loài nhạy cảm và gia tăng ưu thế của loài chống chịu.

Độc học môi trường giải thích cơ chế thẩm thấu trực tiếp qua biểu bì và chuỗi thức ăn, khi chì và asen xâm nhập vào cơ thể ốc làm rối loạn chuyển hóa tế bào, biến đổi mô gan và ức chế khả năng sinh sản. Khung nghiên cứu vận dụng 4 khái niệm cốt lõi:

  • Đa dạng sinh học thân mềm chân bụng: Biểu thị qua cấu trúc phân loại của lớp Gastropoda với 2 phân lớp chính là Mang trước (Prosobranchia) và Có phổi (Pulmonata).
  • Chỉ số đa dạng sinh học: Bao gồm độ phong phú loài Margalef (d), chỉ số đa dạng loài Shannon-Wiener (H') và độ đồng đều Pielou (J').
  • Ô nhiễm kim loại nặng thổ nhưỡng: Trạng thái tích tụ hàm lượng Pb, As vượt quá ngưỡng chịu tải tự nhiên và quy chuẩn kỹ thuật quốc gia.
  • Sinh vật chỉ thị môi trường đất: Các loài sinh vật phản ánh hiện trạng ô nhiễm thông qua sự biến đổi về mật độ, hành vi và cấu trúc quần xã.

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

Nguồn dữ liệu của luận văn được thu thập trực tiếp qua 2 đợt khảo sát thực địa chuyên sâu: Đợt 1 từ ngày 18/04/2017 đến 19/04/2017 và Đợt 2 từ ngày 02/06/2017 đến 03/06/2017. Cỡ mẫu nghiên cứu gồm 30 mẫu đất và 30 mẫu ốc cạn định lượng, phân bố đồng đều trên 5 dạng sinh cảnh: núi đá vôi gần nhà máy tuyển luyện Lủng Váng (5 điểm), núi đất đồi gần mỏ đá (6 điểm), rừng tự nhiên núi đất (6 điểm), sườn núi đá trung tâm (6 điểm) và rừng trồng núi đất gần mỏ Bằng Lũng (7 điểm).

Phương pháp chọn mẫu tuân thủ tiêu chuẩn TCVN 5297:1995 đối với đất và phương pháp ô tiêu chuẩn định lượng 1m² lặp lại 3 lần theo quy chuẩn quốc tế của Vermeulen và Maassen (2003) đối với ốc cạn. Việc lấy mẫu tầng đất mặt ở độ sâu 10 cm giúp phản ánh chính xác tầng tiếp xúc trực tiếp của thân mềm.

Phương pháp phân tích hóa lý thực hiện theo EPA Method 2008 trên thiết bị quang phổ hấp thụ nguyên tử AAS (Model AAnalyst 800 của hãng PerkinElmer, độ lặp lại sai số RSD dưới 2%), tuân thủ TCVN 6649:2000, TCVN 8246:2009 đối với chì và TCVN 8467:2010 đối với asen. Dữ liệu sinh thái được xử lý và phân loại theo hệ thống Schileyko (2011) và Poppe & Tagaro (2006), sau đó tính toán trên phần mềm sinh thái chuyên dụng Primer 6.18 kết hợp phân tích tương quan hồi quy đơn biến nhằm kiểm định độ tin cậy của các mối quan hệ sinh thái.

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

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

Kết quả phân tích 30 mẫu đất và quần xã thân mềm chân bụng tại thị trấn Bằng Lũng đã ghi nhận 4 phát hiện khoa học quan trọng:

  • Ô nhiễm chì ở mức nghiêm trọng cục bộ: Có tới 25/30 điểm khảo sát (tương đương 83,33%) có hàm lượng chì vượt giới hạn cho phép theo QCVN 03-MT:2015/BTNMT đối với đất lâm nghiệp (100 mg/kg). Hàm lượng Pb biến thiên rất lớn từ 57,05 mg/kg đến 8.025 mg/kg đất. Điểm khảo sát BL19 ghi nhận mức ô nhiễm cao nhất đạt 8.025 mg/kg, vượt 80,25 lần quy chuẩn; tiếp theo là BL010 đạt 7.734 mg/kg (vượt 77,34 lần) và BL17 đạt 7.085 mg/kg (vượt 70,85 lần).
  • Hàm lượng Asen nằm trong ngưỡng kiểm soát an toàn: Toàn bộ 30 mẫu đất khảo sát (100%) đều có hàm lượng Asen dao động trong khoảng từ 2,598 mg/kg (điểm BL08) đến 6,025 mg/kg (điểm BL02). Tất cả các giá trị này đều nằm sâu dưới ngưỡng tối đa cho phép là 20 mg/kg theo quy chuẩn quốc gia.
  • Cấu trúc đa dạng sinh học thân mềm chân bụng phong phú: Thu thập và định danh được 304 cá thể thuộc 31 loài, 23 giống, 10 họ, 2 bộ và 2 phân lớp. Trong đó, phân lớp Có phổi (Pulmonata) chiếm ưu thế áp đảo với 80,00% số họ, 69,57% số giống và 61,29% số loài, tiêu biểu là loài Bradybaena jourdyi. Phân lớp Mang trước (Prosobranchia) chiếm 20,00% số họ, 30,43% số giống và 38,71% số loài, nổi bật là loài Cyclotus discoideus.
  • Phân hóa sinh thái rõ rệt theo sinh cảnh: Sinh cảnh núi đá vôi ghi nhận chỉ số đa dạng sinh học cao nhất, đặc biệt tại điểm BL02 có số lượng loài và mật độ cá thể vượt trội so với các sinh cảnh rừng trồng và sườn núi đất thoái hóa.

Thảo luận kết quả

Mức độ ô nhiễm chì đặc biệt cao tại các điểm BL19, BL010, BL17 và BL16 xuất phát trực tiếp từ hoạt động khai thác quặng sunphua và quặng ôxít tại mỏ Bằng Lũng và mỏ Lủng Váng, kết hợp bụi phát tán từ các xe chuyên chở quặng dọc tuyến quốc lộ 254 và quốc lộ 257. Đối với Asen, nồng độ cao nhất tại BL02 (6,025 mg/kg) phản ánh sự tích tụ tự nhiên kết hợp với rửa trôi khoáng sản arsenopyrite gần khu vực tuyển quặng.

Mối tương quan giữa độc chất và các chỉ số sinh thái được lượng hóa chi tiết qua các biểu đồ phân tán và hàm hồi quy. Hàm lượng Asen trong đất thể hiện tương quan nghịch với số loài S (hệ số r = -0,367; phương trình hồi quy $y = 4,3852e^{-0,128x}$), mật độ cá thể N (r = -0,401; phương trình $y = 20,529e^{-0,27x}$), độ phong phú Margalef d (r = -0,182) và chỉ số Shannon H' (r = -0,236). Ngược lại, Asen có tương quan thuận mức trung bình với độ đồng đều Pielou J' (r = 0,478; phương trình $y = 0,176\ln(x) + 0,598$).

Đối với chì (Pb), nồng độ độc chất có tương quan nghịch nhẹ với số loài S (r = -0,071; phương trình $y = -0,00004x + 2,9145$) và mật độ N (r = -0,120; phương trình $y = -0,0005x + 10,768$), nhưng lại tương quan thuận rất mạnh với chỉ số đồng đều J' (r = 0,749; phương trình $y = 0,0558x + 0,4401$). Kết quả này hoàn toàn nhất quán với nghiên cứu của Hoàng Ngọc Khắc tại huyện Sơn Động (tỉnh Bắc Giang), chứng minh rằng khi ô nhiễm kim loại nặng gia tăng, các loài mẫn cảm bị tiêu diệt làm giảm tổng số loài và mật độ cá thể, trong khi số ít cá thể của các loài có sức chống chịu cao như Bradybaena jourdyi sống sót và phân bổ đồng đều hơn trong không gian sinh thái.

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

Dựa trên các bằng chứng khoa học thực nghiệm, nghiên cứu đề xuất 4 nhóm giải pháp kỹ thuật và quản lý môi trường cấp bách:

  • Thiết lập mạng lưới quan trắc sinh học định kỳ: Ủy ban nhân dân huyện Chợ Đồn và Sở Tài nguyên và Môi trường tỉnh Bắc Kạn cần áp dụng mô hình chỉ thị sinh học ốc cạn vào quy trình quan trắc môi trường đất hàng năm. Mục tiêu là kiểm soát độ biến thiên chỉ số Shannon H' và độ đồng đều J' tại 100% các khu vực tiếp giáp mỏ quặng, hoàn thành triển khai trong giai đoạn 2024 - 2025.
  • Kiểm soát nghiêm ngặt nguồn phát thải bụi chì và nước rỉ mỏ: Yêu cầu các doanh nghiệp khai thác tại mỏ Lủng Váng và mỏ Bằng Lũng đầu tư hệ thống phun sương dập bụi, che chắn 100% thùng xe vận chuyển trên tuyến quốc lộ 254 và 257. Target đặt ra là giảm 50% lượng bụi chì lắng đọng trên đất mặt trong vòng 18 tháng sau khi áp dụng chế tài kiểm tra.
  • Ứng dụng công nghệ thực vật xử lý đất ô nhiễm (Phytoremediation): Ban Quản lý rừng phòng hộ phối hợp với các viện nghiên cứu tiến hành trồng các loài cây siêu tích lũy kim loại kết hợp bón vôi nông nghiệp nhằm tăng pH đất, cố định ion chì tự do tại 25 điểm vượt quy chuẩn. Mục tiêu giảm nồng độ Pb trong đất mặt từ trên 5.000 mg/kg xuống dưới 300 mg/kg trong lộ trình 3 đến 5 năm.
  • Hoàn thiện quy hoạch phân vùng bảo tồn đa dạng sinh học núi đá vôi: Phòng Tài nguyên và Môi trường huyện Chợ Đồn cần xây dựng ranh giới hành lang đệm bảo vệ nghiêm ngặt các sinh cảnh núi đá vôi có độ đa dạng cao, nghiêm cấm khai thác đá vôi và san lấp tự phát, bảo tồn nguồn gen của 31 loài thân mềm bản địa trước năm 2026.

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

  • Cơ quan quản lý nhà nước về môi trường: Cán bộ Chi cục Bảo vệ Môi trường và Phòng Tài nguyên Môi trường cấp huyện có thể sử dụng các số liệu ô nhiễm chì và bộ phương trình tương quan để xây dựng định mức xả thải, cấp phép môi trường và lập kế hoạch phục hồi môi trường sau khai khoáng.
  • Các doanh nghiệp khai thác và chế biến khoáng sản: Lãnh đạo và kỹ sư môi trường tại các mỏ quặng chì kẽm có thể ứng dụng kết quả nghiên cứu để thiết kế chương trình giám sát phát tán độc chất cục bộ, đánh giá tác động môi trường thực tế từ khâu vận chuyển và tuyển khoáng.
  • Giảng viên, nghiên cứu sinh và học viên chuyên ngành Môi trường, Sinh học: Tài liệu cung cấp cơ sở lý luận, quy trình thu mẫu ốc cạn bằng ô tiêu chuẩn 1m² và phương pháp xử lý số liệu tương quan đa biến trên Primer 6.18 phục vụ đắc lực cho các đề tài nghiên cứu độc học sinh thái.
  • Các tổ chức bảo tồn thiên nhiên và phát triển bền vững: Các chuyên gia sinh thái có thể khai thác danh mục 31 loài ốc cạn và đặc tính phân bố sinh cảnh để đánh giá tính dễ bị tổn thương của hệ sinh thái karst trước áp lực công nghiệp hóa.

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

Tại sao động vật thân mềm chân bụng trên cạn được chọn làm sinh vật chỉ thị kim loại nặng?

Ốc cạn di chuyển rất chậm, có vùng sống hẹp và tiếp xúc trực tiếp với đất qua chân bụng cũng như tầng thảm mục. Lớp vỏ vôi và mô mềm của ốc có khả năng tích lũy chì và asen qua chuỗi thức ăn, phản ánh chính xác mức độ ô nhiễm tích lũy theo thời gian tại một vị trí cố định.

Tình trạng ô nhiễm chì tại thị trấn Bằng Lũng nghiêm trọng ở mức độ nào?

Kết quả khảo sát thực tế cho thấy có 25/30 điểm lấy mẫu (83,33%) vượt quy chuẩn QCVN 03-MT:2015/BTNMT. Cá biệt tại điểm BL19, nồng độ chì đạt tới 8.025 mg/kg, vượt gấp 80 lần ngưỡng cho phép của đất lâm nghiệp, tập trung chủ yếu quanh các mỏ khai thác và tuyến đường vận chuyển quặng.

Hàm lượng Asen trong đất tại khu vực nghiên cứu có vượt quy chuẩn quốc gia không?

Không. Toàn bộ 30 điểm thu mẫu đều có hàm lượng Asen dao động từ 2,598 mg/kg đến 6,025 mg/kg, nằm trong giới hạn an toàn theo QCVN 03-MT:2015/BTNMT (ngưỡng tối đa cho phép là 20 mg/kg đất). Điểm cao nhất ghi nhận tại BL02 gần mỏ Lủng Váng với 6,025 mg/kg.

Mối quan hệ toán học giữa hàm lượng kim loại nặng và đa dạng loài ốc cạn là gì?

Hàm lượng Asen có tương quan nghịch với số loài S (r = -0,367; phương trình $y = 4,3852e^{-0,128x}$) và mật độ cá thể N (r = -0,401). Trong khi đó, chì có mối tương quan thuận mạnh với độ đồng đều Pielou J' (r = 0,749; phương trình $y = 0,0558x + 0,4401$).

Cần áp dụng biện pháp cấp bách nào để giảm thiểu rủi ro ô nhiễm chì cho cư dân địa phương?

Cần kiểm soát bụi phát tán bằng cách rửa xe và phủ bạt kín các phương tiện vận chuyển quặng trên quốc lộ 254 và 257, đồng thời thực hiện cô lập các bãi thải quặng đuôi và trồng cây cải tạo đất tại các điểm có hàm lượng chì trên 1.000 mg/kg.

Kết luận

  • Luận văn đã xác lập bức tranh toàn diện về ô nhiễm chì nghiêm trọng tại thị trấn Bằng Lũng với 83,33% mẫu đất vượt quy chuẩn, nồng độ cao nhất lên tới 8.025 mg/kg.
  • Khẳng định nồng độ Asen trong đất vẫn nằm trong ngưỡng kiểm soát an toàn tuyệt đối với 100% mẫu khảo sát đạt tiêu chuẩn quốc gia (dao động 2,598 - 6,025 mg/kg).
  • Xây dựng danh lục đa dạng sinh học chi tiết gồm 31 loài thân mềm chân bụng thuộc 23 giống và 10 họ, chỉ ra ưu thế vượt trội của phân lớp Pulmonata tại sinh cảnh núi đá vôi.
  • Lượng hóa thành công các mô hình tương quan và phương trình hồi quy sinh thái, chứng minh tính khả thi của việc sử dụng ốc cạn làm sinh vật chỉ thị kim loại nặng.
  • Đề xuất hệ thống giải pháp liên ngành từ quan trắc sinh học định kỳ, xử lý đất bằng thực vật đến kiểm soát giao thông vận tải khoáng sản với lộ trình 2024 - 2026.

Để bảo vệ môi trường sống và sức khỏe cộng đồng, các cơ quan chức năng và doanh nghiệp khai khoáng cần khẩn trương áp dụng ngay các khuyến nghị khoa học của công trình nghiên cứu này vào thực tiễn quản lý địa phương.