Tổng quan nghiên cứu

Nghiên cứu đặc điểm phân bố không gian các tính chất của đất là cơ sở khoa học then chốt trong công tác quy hoạch lâm nghiệp và phục hồi hệ sinh thái rừng bền vững. Đất đai là nguồn tài nguyên có hạn về không gian nhưng mang tính quyết định đến năng suất sinh khối cây trồng. Trong các hệ sinh thái rừng nhiệt đới, lớp đất mặt thường xuyên chịu tác động của quá trình rửa trôi, xói mòn và suy thoái độ phì. Nhằm giải quyết bài toán suy giảm dinh dưỡng và đánh giá sự sai khác thổ nhưỡng dưới các trạng thái rừng trồng, nghiên cứu được triển khai tại khu rừng thực nghiệm Núi Luốt thuộc Trường Đại học Lâm nghiệp tại thị trấn Xuân Mai, Hà Nội.

Mục tiêu trọng tâm của luận văn là xác định cấu trúc phân bố không gian của các chỉ tiêu lý hóa đất tầng mặt, đồng thời đánh giá độ tin cậy của các phương pháp lấy mẫu đất trong ô tiêu chuẩn so với toàn lâm phần. Phạm vi nghiên cứu được thực hiện trên một lâm phần thực nghiệm điển hình có diện tích 6.000 m2, đặc trưng bởi địa hình đồi thấp lượn sóng với độ dốc trung bình 15° đến 27°, độ cao tuyệt đối từ 76 m đến 133 m. Nghiên cứu tập trung phân tích 108 điểm mẫu đất tầng mặt ở độ sâu từ 0 đến 20 cm để đo lường các biến số về dung trọng, độ ẩm, thành phần cơ giới, hàm lượng chất hữu cơ, đạm dễ tiêu và lân dễ tiêu.

Kết quả nghiên cứu mang ý nghĩa thực tiễn to lớn trong việc cung cấp cơ sở dữ liệu định lượng, giúp các nhà quản lý lâm nghiệp tối ưu hóa mật độ trồng và lựa chọn loài cây trồng bản địa thích hợp nhằm cải thiện tốc độ sinh trưởng chiều cao đạt trên 0,8 đến 0,9 m/năm, đồng thời giảm thiểu sai số trong điều tra thổ nhưỡng định kỳ.

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

Khung lý thuyết áp dụng

Luận văn vận dụng nền tảng học thuyết phát sinh học đất của Vasily Dokuchaev năm 1879, khẳng định đất là một thể tự nhiên độc lập được hình thành dưới sự tác động tổng hợp của 5 yếu tố môi trường bao gồm đá mẹ, sinh vật, khí hậu, địa hình và thời gian. Trong đó, thảm thực vật rừng giữ vai trò chủ đạo trong việc cung cấp tàn tích hữu cơ và điều hòa vi khí hậu đất. Bên cạnh đó, nghiên cứu tích hợp học thuyết mùn hóa của Selman Waksman năm 1936 và Frank Stevenson năm 1982, giải thích các con đường sinh hóa chuyển hóa hợp chất polyphenol và lignin thành các chuỗi axit humic cao phân tử thông qua hoạt động của hệ vi sinh vật đất.

Mô hình nghiên cứu tập trung vào 4 khái niệm thổ nhưỡng học cốt lõi:

  • Phân bố không gian: Sự biến thiên liên tục hoặc gián đoạn của các thuộc tính đất trên bề mặt địa hình theo tọa độ hai chiều.
  • Dung trọng đất: Khối lượng của một đơn vị thể tích đất khô ở trạng thái cấu trúc tự nhiên, phản ánh độ chặt và độ xốp của đất.
  • Hàm lượng chất hữu cơ tầng mặt: Chỉ số phản ánh độ phì nhiêu, khả năng trao đổi cation và cấu trúc viên hạt của tầng đất 0 đến 20 cm.
  • Dinh dưỡng khoáng dễ tiêu: Hàm lượng ion amoni NH4+ và phosphat P2O5 hòa tan trong dung dịch đất mà rễ cây có thể hấp thụ trực tiếp.

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

Nguồn dữ liệu được thu thập từ khảo sát thực địa toàn diện trên lâm phần 6.000 m2 và 2 ô tiêu chuẩn có diện tích 1.000 m2 mỗi ô tại khu vực Núi Luốt trong năm 2017.

Phương pháp chọn mẫu và cỡ mẫu:

  • Trên toàn lâm phần 6.000 m2, thiết kế hệ thống 108 điểm lấy mẫu ngẫu nhiên hệ thống bằng máy kinh vĩ và địa bàn ba chân.
  • Trong mỗi ô tiêu chuẩn 1.000 m2, tiến hành thu thập mẫu đất theo hai phương pháp lấy mẫu tổng hợp: phương pháp Ziczac gồm 9 điểm chia đều trên 3 tuyến cách nhau 13,3 m và phương pháp Tán xạ gồm 8 điểm hướng tâm cách tâm ô 10 m cùng 1 điểm tại tâm.

Quy trình phân tích thí nghiệm:

  • Xác định chất hữu cơ bằng phương pháp oxy hóa ướt Tyurin với dung dịch K2Cr2O7 và H2SO4 đun ở nhiệt độ 150°C.
  • Xác định đạm dễ tiêu NH4+ bằng phương pháp chiết KCl, tạo màu với thuốc thử Nessler và đo trên máy quang phổ kế.
  • Xác định lân dễ tiêu P2O5 bằng phương pháp so màu quang phổ.
  • Xác định thành phần cơ giới bằng phương pháp hút pipet phân cấp hạt.
  • Xác định dung trọng bằng ống dung trọng thể tích chuẩn và độ ẩm bằng phương pháp sấy khô ở 105°C.

Lý do lựa chọn phương pháp phân tích: Việc kết hợp thống kê mô tả, kiểm định giả thuyết t-Student, phân tích phương sai F-test với mô hình hóa địa thống kê hàm bán phương sai Variogram và thuật toán nội suy nghịch đảo khoảng cách IDW cho phép mô tả chính xác tính dị hướng, tính biến thiên liên tục trong không gian và đánh giá khách quan sai số thống kê của các phương pháp trộn mẫu.

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

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

Thứ nhất, hàm lượng chất hữu cơ tầng mặt trong lâm phần đạt giá trị trung bình 2,16%, dao động trong khoảng từ 0,77% đến 3,39%. Hệ số biến động đạt 24,5% phản ánh mức độ không đồng nhất cao giữa các vị trí. Phân bố thực nghiệm có dạng gần đối xứng với 68% số mẫu tập trung trong khoảng 1,6% đến 2,8%. Mô hình bán phương sai dạng tuyến tính cho thấy khoảng cách hữu hiệu để phản ánh đầy đủ biến động mùn trong lâm phần lên tới 55 m.

Thứ hai, dung trọng đất trung bình của lâm phần đạt 1,47 g/cm3, dao động từ 0,93 g/cm3 đến 1,93 g/cm3 với hệ số biến động là 13,1%. Phân bố dung trọng có đỉnh tập trung rõ rệt ở mức 1,55 g/cm3, trong đó các mẫu đất có dung trọng dưới 1,25 g/cm3 hoặc trên 1,75 g/cm3 chỉ chiếm khoảng 10% tổng số mẫu. Trong khi đó, độ ẩm đất trung bình đạt 0,30%, dao động từ 0,17% đến 0,46% với hệ số biến động 22,4%, có xu hướng giảm dần từ sườn phía Tây sang phía Đông của lâm phần.

Thứ ba, hàm lượng dinh dưỡng khoáng dễ tiêu trong đất ở mức rất nghèo. Hàm lượng đạm dễ tiêu NH4+ trung bình chỉ đạt 1,99 mg/100g đất với hệ số biến động rất cao lên tới 60,6%, trong đó có 75% số điểm lấy mẫu đạt dưới mức 2,75 mg/100g đất. Hàm lượng lân dễ tiêu P2O5 trung bình đạt 1,40 mg/100g đất với hệ số biến động 55,6%, đa số mẫu phân tích tập trung dưới mức 1,25 mg/100g đất. Về thành phần cơ giới, 64% số mẫu (23 trên 36 mẫu) thuộc nhóm đất thịt pha sét, tỷ lệ cấp hạt cát chiếm 36,08%, cấp hạt sét chiếm 32,43% và cấp hạt thịt chiếm 31,49%.

Thứ tư, về tính đại diện của ô tiêu chuẩn và phương pháp lấy mẫu, kết quả kiểm định t-test cho thấy hàm lượng chất hữu cơ trong ô tiêu chuẩn 1 (2,31 g/100g) và ô tiêu chuẩn 2 (2,32 g/100g) không có sự sai khác mang ý nghĩa thống kê so với phần còn lại của lâm phần (2,10 g/100g) với p = 0,069 và p = 0,077 (p > 0,05). Tuy nhiên, độ ẩm đất lại có sự sai khác cực kỳ rõ rệt giữa ô tiêu chuẩn và lâm phần với p = 0,00003 (p < 0,05). Đặc biệt, phương pháp lấy mẫu tổng hợp theo đường Ziczac và Tán xạ đều cho kết quả sai lệch có ý nghĩa thống kê so với giá trị trung bình thực của các mẫu độc lập trong ô tiêu chuẩn.

Thảo luận kết quả

Sự phân bố không gian không đồng nhất của hàm lượng chất hữu cơ và dung trọng đất tại Núi Luốt bắt nguồn từ điều kiện địa chất và lịch sử tác động thảm thực vật. Đất tại khu vực là đất Feralit nâu vàng phát triển trên đá mẹ Poocfiarit trung tính. Ở tầng sâu 20 đến 60 cm, tỷ lệ kết von sắt đạt từ 20% đến 30%, chứng tỏ khu vực từng trải qua thời kỳ đồi trọc bị xói mòn mạnh trước khi được trồng rừng phục hồi.

Khả năng cố định lân của đất rất cao do sự tích lũy nhôm và sắt di động trong môi trường chua (pH dưới 7), dẫn đến hàm lượng P2O5 dễ tiêu cực kỳ thấp (1,40 mg/100g). Sự khác biệt về độ ẩm và mùn giữa các khu vực phụ thuộc trực tiếp vào độ tàn che tầng cây cao (ô tiêu chuẩn 1 đạt 0,7 so với ô tiêu chuẩn 2 đạt 0,5) và độ che phủ của thảm tươi (80% so với 65%). Kết quả này hoàn toàn thống nhất với các công trình nghiên cứu của các nhà khoa học lâm nghiệp tại Việt Nam về mối quan hệ giữa cấu trúc tán rừng và khả năng duy trì độ phì đất mặt.

Dữ liệu nghiên cứu được biểu diễn trực quan qua bản đồ đẳng trị nội suy không gian IDW, đồ thị bán phương sai Variogram và bảng tam giác thành phần cơ giới đất, giúp minh họa rõ nét các đốm biến động cục bộ mà các phương pháp thống kê cổ điển không thể hiện được.

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

Nhằm nâng cao chất lượng lập địa và chuẩn hóa phương pháp điều tra thổ nhưỡng lâm nghiệp, 4 nhóm giải pháp cụ thể được đề xuất:

  1. Chuyển đổi mô hình lâm sinh sang rừng hỗn giao đa tầng tán: Thực hiện thay thế các lập địa rừng thuần loài sinh trưởng kém bằng các loài cây bản địa như Lim xanh, Lim xẹt, Đinh thối kết hợp cây họ Đậu cải tạo đất. Mục tiêu nâng độ che phủ thảm tươi lên trên 85% và bổ sung lượng vật chất hữu cơ rơi rụng đạt từ 4 đến 6 tấn/ha/năm. Thời gian thực hiện từ quý 1 năm 2026 đến hết năm 2028, do Ban Quản lý Rừng thực nghiệm Núi Luốt chủ trì.

  2. Cải tạo độ phì và giải tỏa lân cố định trong đất Feralit: Bón bổ sung vôi bột khử chua kết hợp phân lân nung chảy và phân vi sinh phân giải lân với liều lượng 200 đến 300 kg/ha cho các lâm phần nghèo dinh dưỡng. Mục tiêu nâng hàm lượng P2O5 dễ tiêu từ 1,40 mg/100g lên mức tối thiểu 3,50 mg/100g đất trong vòng 18 tháng, do Trung tâm Nghiên cứu Ứng dụng Lâm học triển khai.

  3. Chuẩn hóa quy trình bố trí điểm điều tra đất lâm nghiệp: Ban hành hướng dẫn kỹ thuật điều tra thổ nhưỡng yêu cầu khoảng cách lưới lấy mẫu tối thiểu đạt từ 50 m đến 55 m nhằm đảm bảo tính đại diện cho biến động không gian. Hạn chế sử dụng mẫu tổng hợp đơn lẻ để tránh sai số thống kê vượt quá 5%. Thời gian áp dụng bắt đầu từ quý 3 năm 2026, do Khoa Lâm học phê duyệt.

  4. Xây dựng công trình tiểu địa hình bảo tồn nước và chống xói mòn: Thiết lập hệ thống rãnh giữ nước và dải cây xanh chắn dòng chảy tràn theo đường đồng mức trên các sườn dốc từ 15° đến 27°. Mục tiêu duy trì độ ẩm đất mặt ổn định trên 0,32% trong suốt mùa khô kéo dài từ tháng 11 đến tháng 3 năm sau, do Đội Kỹ thuật Lâm sinh thực địa phụ trách.

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

Tài liệu luận văn mang lại giá trị học thuật và ứng dụng chuyên sâu cho 4 nhóm đối tượng:

  • Giảng viên và nghiên cứu sinh chuyên ngành Lâm học, Khoa học Đất: Cung cấp khung phương pháp luận hoàn chỉnh về ứng dụng địa thống kê (Geostatistics), mô hình bán phương sai Variogram và phương pháp nội suy IDW trong việc định lượng biến tính đất lâm nghiệp.
  • Kỹ sư lâm sinh và cán bộ quy hoạch vùng trồng rừng: Sử dụng số liệu thực chứng về dung trọng đất (1,47 g/cm3) và tỷ lệ cấp hạt thịt pha sét (64%) để thiết kế mật độ trồng rừng, cơ cấu cây trồng hỗn giao tối ưu cho vùng gò đồi Feralit.
  • Cán bộ kỹ thuật tại các Chi cục Kiểm lâm và Ban Quản lý Rừng phòng hộ: Tham khảo quy trình đánh giá độ phì, hướng dẫn khắc phục tình trạng thiếu hụt đạm (1,99 mg/100g) và lân dễ tiêu (1,40 mg/100g) dưới các tán rừng trồng phục hồi.
  • Học viên cao học và sinh viên các ngành Nông Lâm nghiệp: Làm tài liệu tham khảo chuẩn mực về phương pháp bố trí ô tiêu chuẩn 1.000 m2, kỹ thuật lấy mẫu theo tuyến Ziczac, Tán xạ và phương pháp phân tích phẫu diện đất Feralit tầng dày.

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

Đặc điểm thổ nhưỡng nổi bật nhất của đất tầng mặt tại khu vực Núi Luốt là gì?
Đất tại Núi Luốt là đất Feralit nâu vàng phát triển trên đá mẹ Poocfiarit trung tính. Đất có thành phần cơ giới chủ yếu là thịt pha sét chiếm 64% tổng số mẫu, kết cấu chặt với dung trọng trung bình 1,47 g/cm3, giàu nhôm sắt và nghèo dinh dưỡng khoáng dễ tiêu.

Tại sao hàm lượng lân dễ tiêu trong đất lâm phần lại ở mức rất thấp (1,40 mg/100g đất)?
Đất Feralit tại khu vực có tính chất chua và tích lũy hàm lượng sắt nhôm tự do cao, tạo ra tầng kết von từ 20% đến 30% ở độ sâu 20 đến 60 cm. Các ion sắt nhôm này liên kết chặt với ion phosphat tạo thành các hợp chất khó tan, làm giảm nghiêm trọng lượng lân dễ tiêu cho cây.

Phương pháp lấy mẫu đất tổng hợp Ziczac và Tán xạ có thay thế hoàn toàn được mẫu ngẫu nhiên không?
Không. Kết quả kiểm định thống kê t-test trong luận văn chỉ ra rằng việc trộn mẫu tổng hợp theo đường Ziczac hay Tán xạ đều tạo ra sự sai lệch có ý nghĩa thống kê so với giá trị trung bình của các mẫu đất ngẫu nhiên độc lập trong ô tiêu chuẩn.

Độ tàn che của tán cây rừng ảnh hưởng như thế nào đến độ ẩm đất tầng mặt?
Độ tàn che đóng vai trò hạn chế bốc hơi bề mặt trực tiếp từ ánh nắng mặt trời. Tại ô tiêu chuẩn có độ tàn che 0,7 và che phủ thảm tươi 80%, độ ẩm đất mặt đạt 0,34%, cao hơn đáng kể so với khu vực có độ tàn che 0,5 với độ ẩm chỉ đạt 0,26%.

Khoảng cách bố trí điểm lấy mẫu tối ưu để phản ánh biến động hàm lượng mùn là bao nhiêu?
Dựa trên phân tích mô hình bán phương sai Variogram, khoảng cách hữu hiệu đạt 95% phương sai cực đại đối với hàm lượng chất hữu cơ trong lâm phần là 55 m. Do đó, cự ly lưới lấy mẫu cần được thiết kế ở khoảng cách này để đảm bảo độ chính xác.

Kết luận

  • Luận văn đã xác lập hệ thống dữ liệu định lượng về 5 chỉ tiêu lý hóa đất tầng mặt trên lâm phần thực nghiệm 6.000 m2 với 108 điểm mẫu tại Núi Luốt.
  • Phát hiện hàm lượng chất hữu cơ đạt mức trung bình 2,16%, dung trọng đạt 1,47 g/cm3, trong khi hàm lượng đạm NH4+ (1,99 mg/100g) và lân P2O5 (1,40 mg/100g) đều ở mức rất nghèo.
  • Chứng minh tính biến động không gian ngẫu nhiên cao của đất tầng mặt, xác định khoảng cách lấy mẫu hữu hiệu cho chất hữu cơ là 55 m và độ ẩm là 50 m.
  • Xác nhận ô tiêu chuẩn 1.000 m2 có tính đại diện tốt cho chất hữu cơ của lâm phần (p = 0,069) nhưng không đại diện cho độ ẩm (p = 0,00003).
  • Khẳng định phương pháp trộn mẫu tổng hợp Ziczac và Tán xạ làm phát sinh sai số thống kê so với giá trị trung bình thực của các mẫu độc lập.

Đóng góp khoa học then chốt của công trình là việc ứng dụng thành công kỹ thuật địa thống kê để làm sáng tỏ bản chất biến thiên thổ nhưỡng lâm nghiệp tầng mặt, tạo tiền đề cải tiến phương pháp điều tra rừng. Kế hoạch nghiên cứu tiếp theo giai đoạn 2026 đến 2028 sẽ mở rộng đánh giá phẫu diện đất ở các tầng sâu 60 đến 100 cm trên nhiều dạng địa hình phức tạp. Các nhà nghiên cứu và kỹ sư lâm sinh hãy khai thác và áp dụng ngay các bộ chỉ số này để nâng cao hiệu quả phục hồi rừng và quản lý đất bền vững.