Tổng quan nghiên cứu

Huyện Phú Lương thuộc tỉnh Thái Nguyên sở hữu tổng diện tích tự nhiên là 36.886,68 ha, trong đó địa hình đồi núi dốc chiếm trên 70% toàn huyện. Đáng chú ý, diện tích đất có độ dốc trên 3 độ chiếm tới 91,78%, tương đương 33.854,91 ha, phân bố chủ yếu trên nền địa hình chia cắt phức tạp bởi quá trình karst hóa mạnh mẽ. Sự kết hợp giữa lượng mưa trung bình năm lớn dao động từ 1.722,76 mm đến 2.005,08 mm cùng các hoạt động nhân sinh kéo dài như thâm canh nông nghiệp trên 4.300 ha chè và khai thác khoáng sản than, titan đã tạo áp lực nghiêm trọng lên tính chất lý hóa của đất. Đất dốc đối mặt với nguy cơ xói mòn, rửa trôi dinh dưỡng tầng mặt và tích tụ các chất độc hại thứ sinh tại các vùng chân dốc.

Mục tiêu nghiên cứu tập trung đánh giá thực trạng phân bố cấp độ dốc và tầng dày đất, phân tích mối quan hệ tương quan giữa các yếu tố tài nguyên sinh khí hậu với chỉ số thực vật NDVI, đồng thời lượng hóa các chỉ tiêu dinh dưỡng, độ chua và dư lượng kim loại nặng như chì, asen, cadimi. Trên cơ sở đó, đề tài xây dựng luận cứ khoa học để đề xuất hệ thống giải pháp công trình và phi công trình nhằm bảo vệ môi trường đất dốc bền vững.

Phạm vi nghiên cứu được triển khai thực địa tại 3 tiểu vùng đại diện gồm xã Yên Ninh, xã Yên Lạc và xã Động Đạt trong khung thời gian từ ngày 30/08/2016 đến ngày 30/07/2017. Ý nghĩa của đề tài thể hiện qua việc định lượng hóa dữ liệu thổ nhưỡng của 12 loại đất, hỗ trợ thiết lập giải pháp giảm thiểu từ 30% đến 40% lượng đất mất hàng năm do xói mòn, trực tiếp bảo vệ hơn 22.866,51 ha đất có tầng dày từ 70 cm đến 100 cm phục vụ phát triển kinh tế nông lâm nghiệp địa phương.

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 ứng dụng lý thuyết suy thoái đất do xói mòn thủy lực nhiệt đới kết hợp với lý thuyết sinh thái cảnh quan đồi núi để đánh giá sự tương tác đa chiều giữa điều kiện tự nhiên và áp lực nhân sinh. Khung nghiên cứu vận dụng 5 khái niệm chuyên ngành nền tảng:

Thứ nhất, môi trường đất là một hệ sinh thái chức năng hoàn chỉnh, phản ánh tổ hợp các đặc tính hóa học, vật lý và sinh học, đảm nhận vai trò duy trì chu trình sinh địa hóa và thanh lọc nguồn nước ngầm. Thứ hai, đất dốc được định nghĩa là bề mặt đất nằm nghiêng với địa hình mấp mô, phân thành các cấp dốc từ dốc nhẹ dưới 7 độ đến dốc rất mạnh trên 25 độ. Thứ ba, suy thoái môi trường đất biểu hiện qua sự sụt giảm dung lượng chất lượng, mất kết cấu hạt đất và giảm khả năng sản xuất sinh khối. Thứ tư, dung dịch đất và độ chua phản ánh trạng thái động của các cation kiềm trao đổi như canxi, magie bị thay thế bởi ion hydro và nhôm, quyết định chỉ số pH. Thứ năm, độ phì đất biểu thị hàm lượng chất hữu cơ, đạm tổng số, lân và kali dễ tiêu đáp ứng nhu cầu sinh trưởng cây trồng.

Khung pháp lý của nghiên cứu căn cứ vào Luật Bảo vệ Môi trường năm 2014, Nghị định số 19/2015/NĐ-CP, Quy chuẩn kỹ thuật quốc gia QCVN 03-MT:2015/BTNMT về giới hạn kim loại nặng trong đất và QCVN 15:2008/BTNMT về dư lượng hóa chất bảo vệ thực vật.

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

Nguồn dữ liệu của đề tài tích hợp giữa dữ liệu thứ cấp về bản đồ thổ nhưỡng, số liệu khí tượng thủy văn của tỉnh Thái Nguyên cùng ảnh viễn thám vệ tinh Landsat 8 độ phân giải không gian 30x30 mét chụp vào tháng 6 năm 2016 và mô hình số độ cao kỹ thuật số DEM.

Cỡ mẫu nghiên cứu thực địa bao gồm 25 điểm mẫu đất hỗn hợp đại diện cho 5 loại hình thảm phủ chính gồm cây bụi tự nhiên, keo non, keo lâu năm, lúa nước và chè tại 3 xã Yên Ninh, Yên Lạc, Động Đạt. Phương pháp chọn mẫu tuân thủ nghiêm ngặt tiêu chuẩn TCVN 7538-2:2005 và TCVN 5297:1995, thu thập mẫu ở tầng mặt từ 0 cm đến 20 cm sau 2 ngày liên tục không mưa. Tại mỗi ô tiêu chuẩn, lấy từ 3 đến 5 điểm phụ theo đường chéo góc hoặc ziczac tạo mẫu hỗn hợp 1 kg, sau đó rút gọn đồng nhất còn 0,5 kg và định vị tọa độ bằng thiết bị GPS cầm tay. Đồng thời, trên hệ thống thông tin địa lý, sử dụng công cụ tạo điểm ngẫu nhiên để trích xuất 250 điểm cho mỗi tiểu vùng, tổng cộng 750 điểm dữ liệu không gian phục vụ phân tích tương quan sinh khí hậu.

Phương pháp phân tích trong phòng thí nghiệm áp dụng phương pháp đo điện cực thủy tinh xác định pHKCl theo TCVN 5979:2000, phương pháp oxy hóa kali bicromat xác định mùn tổng số theo TCVN 4050-1985, phương pháp Kendan cải biên xác định đạm tổng số theo TCVN 6498:1999, quang phổ phát xạ xác định kali dễ tiêu theo TCVN 8662-2011 và quang phổ hấp thụ nguyên tử AAS xác định kim loại nặng chì, asen, cadimi. Lý do lựa chọn tổ hợp phương pháp này nhằm đảm bảo độ nhạy cao và sai số cực tiểu đối với các mẫu đất miền núi.

Quy trình xử lý dữ liệu sử dụng phần mềm SPSS phiên bản 20.0 để phân tích tương quan Pearson, chạy ước lượng đường cong hồi quy phi tuyến Curve Estimation. Đồng thời, sử dụng phần mềm XLSTAT để phân tích thành phần chính PCA và phân tích cụm phân cấp AHC nhằm nhóm các vị trí có cùng mức độ tác động ô nhiễm. Việc mô hình hóa không gian được thực hiện trên ArcGIS phiên bản 10.2 qua phương pháp nội suy trọng số nghịch đảo khoảng cách IDW cho các chỉ tiêu dinh dưỡng và phương pháp Kriging địa thống kê cho kim loại nặng. Toàn bộ quá trình thu thập và phân tích hoàn thành trong 11 tháng nghiên cứu liên tục.

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

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

Thứ nhất, phân tích địa hình thổ nhưỡng cho thấy đất dốc cấp 3 đến 8 độ chiếm 12.184,40 ha (tương đương 33,03%) và cấp 8 đến 15 độ chiếm 11.672,10 ha (tương đương 31,64% tổng diện tích huyện). Toàn huyện có 12 nhóm loại đất, trong đó đất đỏ vàng trên đá sét và biến chất chiếm diện tích áp đảo với 20.906,03 ha (chiếm 56,68%), tập trung chủ yếu ở cấp dốc trên 25 độ với 11.409,56 ha. Về cấu trúc tầng dày, đất có tầng dày từ 70 cm đến 100 cm chiếm 22.866,51 ha (62,01%), đất có tầng dày từ 50 cm đến 70 cm đạt 6.942,54 ha (18,83%), còn lại là tầng dày trên 100 cm với 5.320,74 ha (14,43%).

Thứ hai, chỉ số thực vật NDVI dao động trong khoảng từ -0,169 đến 0,650, thể hiện tính phân hóa sinh thái sâu sắc giữa các tiểu vùng. Tại tiểu vùng 1 gồm các xã Yên Ninh, Yên Trạch, Yên Đổ, chỉ số NDVI tương quan âm có ý nghĩa thống kê cao với lượng mưa năm (hệ số r = -0,485, p < 0,01, R² = 0,235) và tương quan dương với độ dốc (hệ số r = 0,339, p < 0,01). Tại tiểu vùng 2 gồm các xã Yên Lạc, Phú Đô, Tức Tranh, lượng mưa thể hiện tương quan nghịch rất mạnh với độ dốc (hệ số r = -0,706, p < 0,01) và độ cao (hệ số r = -0,739, p < 0,01) do thảm phủ cây chè chiếm hơn 4.300 ha chi phối mặt đệm. Tại tiểu vùng 3 gồm xã Động Đạt, Phấn Mễ, tương quan giữa NDVI với lượng mưa giảm xuống mức r = -0,273 và với độ dốc đạt r = 0,184 do mức độ đô thị hóa và hoạt động khai khoáng làm phân mảnh thảm phủ.

Thứ ba, đặc tính hóa lý đất dốc tại các điểm khảo sát thể hiện phản ứng chua đến rất chua với chỉ số pHKCl dao động từ 3,8 đến 4,8. Tỷ lệ mùn tổng số và đạm tổng số có xu hướng suy giảm rõ rệt ở các sườn dốc canh tác chè và keo non do thiếu thảm che phủ bảo vệ. Đối với kim loại nặng, hàm lượng chì, asen và cadimi tại các vị trí đất lúa nước và chân dốc xung quanh khu vực khai thác mỏ than Phấn Mễ và mỏ titan Động Đạt cao hơn từ 1,5 đến 2,8 lần so với khu vực rừng tự nhiên vùng cao, dù vẫn nằm trong ngưỡng giới hạn cho phép của QCVN 03-MT:2015/BTNMT.

Thảo luận kết quả

Sự suy giảm độ phì và chua hóa đất dốc tại huyện Phú Lương bắt nguồn từ tác động kép của các hiện tượng rửa trôi tự nhiên và phương thức canh tác bất hợp lý. Lượng mưa trung bình năm lớn từ 1.722,76 mm đến 2.005,08 mm trút xuống các sườn dốc đã nhanh chóng cuốn trôi các cation kiềm trao đổi như canxi, magie, giải phóng các ion hydro và nhôm tự do làm đất bị chua hóa mạnh. Hiện tượng này hoàn toàn trùng khớp với các công bố quốc tế của các nhà khoa học về xói mòn đất nhiệt đới, nơi ước tính hàng chục tỷ tấn đất mặt màu mỡ bị cuốn trôi mỗi năm khỏi các hệ thống canh tác đồi núi.

Bên cạnh đó, việc sử dụng phân bón vô cơ tại các nông hộ tăng 51,7% trong khi cây trồng chỉ hấp thụ được khoảng một phần ba đã thúc đẩy quá trình axit hóa đất và phát tán dư lượng hóa chất xuống các thung lũng tích tụ. Kết quả này tương đồng với nghiên cứu của các chuyên gia thổ nhưỡng trong nước, chỉ ra rằng đất canh tác dốc bị suy giảm tỷ lệ hạt kết đất từ 55,6% xuống chỉ còn 30% so với đất rừng tự nhiên.

Dữ liệu nghiên cứu được trực quan hóa sinh động thông qua ma trận tương quan Pearson và phương trình hồi quy phi tuyến bậc 3, mô tả chính xác sự suy giảm của chỉ số sinh khối thực vật khi lượng mưa và độ dốc thay đổi. Bản đồ phân bố không gian nội suy Kriging và sơ đồ nhánh phân cấp Dendrogram AHC thể hiện rõ rệt sự phân nhóm ô nhiễm kim loại nặng tập trung tại chân dốc, cung cấp bằng chứng trực quan về hướng lan truyền vật chất từ đỉnh đồi xuống các vùng sản xuất nông nghiệp trũng thấp.

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

Thứ nhất, quy hoạch và phân vùng sử dụng đất dốc theo khả năng thích ứng sinh thái. Ủy ban nhân dân huyện Phú Lương cùng Phòng Tài nguyên và Môi trường cần rà soát và chuyển đổi cơ cấu cây trồng trên 2.792,66 ha đất có độ dốc trên 25 độ và 2.740,83 ha đất dốc từ 20 đến 25 độ sang trồng rừng phòng hộ và khoanh nuôi tái sinh tự nhiên. Mục tiêu nâng tỷ lệ che phủ rừng đạt chỉ số NDVI trên 0,65, hoàn thành trong giai đoạn từ năm 2018 đến năm 2022.

Thứ hai, áp dụng các giải pháp canh tác nông lâm kết hợp chống xói mòn trên đất dốc. Trạm Khuyến nông huyện và chính quyền các xã Yên Ninh, Yên Lạc, Động Đạt cần hướng dẫn nông dân xây dựng mô hình ruộng bậc thang, trồng băng cỏ Vetiver chắn sóng và phủ thảm thực vật che phủ đất trên 11.672,10 ha đất dốc từ 8 đến 15 độ. Giải pháp này hướng tới mục tiêu cắt giảm từ 35% đến 45% lượng đất bị xói mòn tầng mặt hàng năm, duy trì triển khai liên tục từ năm 2018 đến năm 2025.

Thứ ba, chuẩn hóa kỹ thuật bón phân và quản lý hóa chất bảo vệ thực vật. Phòng Nông nghiệp và Phát triển nông thôn chỉ đạo các hợp tác xã chè tại vùng 2 thay thế 30% phân hóa học vô cơ bằng phân hữu cơ sinh học, áp dụng quy trình tưới tiêu khoa học trên 4.300 ha chè. Biện pháp này giúp cải thiện độ chua pHKCl đất lên ngưỡng trên 4,5, hạn chế phát tán khí nhà kính oxit nitơ và bảo vệ 22.866,51 ha đất có tầng dày 70 cm đến 100 cm.

Thứ tư, thiết lập mạng lưới quan trắc và giám sát môi trường định kỳ tại các điểm nóng khoáng sản. Sở Tài nguyên và Môi trường tỉnh Thái Nguyên phối hợp với các doanh nghiệp khai thác mỏ tại Động Đạt và Phấn Mễ xây dựng 15 trạm quan trắc định kỳ 6 tháng một lần đối với hàm lượng kim loại nặng chì, cadimi, asen trong đất và nước mặt. Đảm bảo 100% diện tích đất sản xuất nông nghiệp lân cận đáp ứng tiêu chuẩn an toàn sinh thái theo QCVN 03-MT:2015/BTNMT đến năm 2025.

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

Thứ nhất, các nhà quản lý và hoạch định chính sách tại các Sở, Phòng Tài nguyên và Môi trường, Nông nghiệp và Phát triển nông thôn. Nghiên cứu cung cấp hệ cơ sở dữ liệu số hóa GIS và bản đồ phân vùng thổ nhưỡng của 36.886,68 ha đất, phục vụ trực tiếp cho công tác lập quy hoạch sử dụng đất cấp huyện và thẩm định đánh giá tác động môi trường các dự án khai khoáng.

Thứ hai, cán bộ khuyến nông và kỹ sư lâm sinh thực địa. Luận văn cung cấp các chỉ dẫn kỹ thuật định lượng về tầng dày đất và độ dốc, giúp xây dựng các mô hình nông lâm kết hợp tối ưu trên 11.672,10 ha đất dốc 8 đến 15 độ, bảo vệ lớp đất mặt khỏi thoái hóa.

Thứ ba, học viên cao học, nghiên cứu sinh và giảng viên chuyên ngành Khoa học Môi trường, Quản lý Đất đai. Đề tài đóng vai trò là tài liệu tham khảo giá trị về phương pháp tích hợp viễn thám Landsat 8, mô hình số độ cao DEM với các thuật toán thống kê đa biến PCA, AHC và phép nội suy không gian địa thống kê Kriging.

Thứ tư, các doanh nghiệp khai khoáng và hợp tác xã sản xuất chè quy mô lớn. Tài liệu giúp các đơn vị nắm vững cơ chế lan truyền kim loại nặng và rửa trôi dinh dưỡng theo sườn dốc, từ đó xây dựng kế hoạch quản lý đất canh tác an toàn theo tiêu chuẩn VietGAP.

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

Yếu tố tự nhiên nào tác động mạnh nhất đến sự suy thoái đất dốc tại huyện Phú Lương? Độ dốc địa hình kết hợp với lượng mưa trung bình năm lớn từ 1.722,76 mm đến 2.005,08 mm là tác nhân tự nhiên hàng đầu gây xói mòn và rửa trôi đất. Hiện tượng này cuốn đi các cation kiềm trao đổi canxi, magie, khiến đất bị chua hóa nghiêm trọng với pHKCl giảm xuống dưới 4,5.

Chỉ số thực vật NDVI có ý nghĩa như thế nào trong công tác bảo vệ môi trường đất dốc? Chỉ số NDVI chiết xuất từ ảnh vệ tinh phản ánh mật độ che phủ của thảm thực vật từ -0,169 đến 0,650. Thảm thực vật có NDVI cao trên 0,6 giúp tán xạ hạt mưa, giảm tốc độ dòng chảy mặt, qua đó giữ lại hàm lượng mùn và đạm tổng số trong đất, ngăn ngừa suy thoái tầng đất hữu hiệu 70 cm đến 100 cm.

Hàm lượng kim loại nặng trong đất nông nghiệp dốc tại Phú Lương hiện có vượt ngưỡng an toàn không? Kết quả kiểm nghiệm tại 25 điểm mẫu cho thấy hàm lượng chì, asen và cadimi vẫn nằm trong giới hạn cho phép của quy chuẩn QCVN 03-MT:2015/BTNMT. Tuy nhiên, tại các chân dốc gần mỏ titan Động Đạt và mỏ than Phấn Mễ, các kim loại này có xu hướng tích tụ cục bộ cao gấp 1,5 đến 2,8 lần so với vùng rừng tự nhiên.

Nhóm đất chiếm diện tích lớn nhất tại Phú Lương có đặc tính ra sao? Đất đỏ vàng trên đá sét và biến chất là nhóm đất chủ đạo với 20.906,03 ha, chiếm 56,68% tổng diện tích tự nhiên. Loại đất này phân bố nhiều ở cấp dốc trên 25 độ với 11.409,56 ha, có tầng dày từ 70 cm đến 100 cm chiếm tỷ lệ cao nhưng rất nhạy cảm với xói mòn khi bị mất thảm phủ thực vật.

Tại sao nghiên cứu lại phân chia huyện Phú Lương thành 3 tiểu vùng sinh thái riêng biệt? Việc phân chia 3 tiểu vùng dựa trên sự khác biệt về địa hình và tác động nhân sinh: Tiểu vùng 1 là núi cao chia cắt mạnh với rừng tự nhiên; Tiểu vùng 2 là vùng đồi bát úp tập trung hơn 4.300 ha chè; Tiểu vùng 3 là khu vực bằng phẳng chịu ảnh hưởng mạnh từ đô thị và khai thác mỏ.

Kết luận

  • Đất dốc trên 3 độ tại huyện Phú Lương chiếm tới 91,78% tổng diện tích, trong đó nhóm đất đỏ vàng trên đá sét và biến chất chiếm ưu thế tuyệt đối với 20.906,03 ha.
  • Cấu trúc tầng dày đất tương đối thuận lợi với 62,01% diện tích đạt độ sâu từ 70 cm đến 100 cm, song đối mặt với nguy cơ rửa trôi mạnh dưới lượng mưa năm từ 1.722,76 mm đến 2.005,08 mm.
  • Mối tương quan sinh khí hậu khẳng định vai trò cốt lõi của thảm thực vật NDVI trong việc giảm thiểu xói mòn đất và kiểm soát sự chua hóa môi trường đất.
  • Hàm lượng kim loại nặng chì, asen, cadimi dù nằm trong giới hạn cho phép nhưng đã xuất hiện dấu hiệu tích tụ thứ sinh tại các vùng chân dốc trũng thấp lân cận điểm mỏ khoáng sản.
  • Hệ thống giải pháp đề xuất kết hợp đồng bộ giữa quy hoạch rừng phòng hộ trên 5.533,49 ha đất dốc trên 20 độ, kỹ thuật canh tác nông lâm kết hợp và quan trắc môi trường định kỳ.

Đóng góp chính của luận văn là đã chuẩn hóa bộ dữ liệu không gian viễn thám kết hợp phân tích đa biến, cung cấp luận cứ thực tiễn phục vụ quản trị đất dốc bền vững. Kế hoạch tiếp theo cần tập trung số hóa mạng lưới giám sát xói mòn trực tuyến và nhân rộng mô hình thâm canh chè sinh thái giai đoạn 2018 đến 2025. Các cơ quan quản lý và cộng đồng địa phương hãy chủ động phối hợp hành động ngay hôm nay để bảo tồn tài nguyên đất dốc, bảo đảm an ninh sinh thái lâu dài cho tỉnh Thái Nguyên.