Tổng quan nghiên cứu

Tài nguyên đất đai đóng vai trò nền tảng không thể thay thế trong hệ sinh thái rừng, đặc biệt là tầng đất mặt từ 0 đến 20 cm nơi tập trung phần lớn hoạt động sinh học và dinh dưỡng khoáng. Tại khu rừng nghiên cứu thực nghiệm Núi Luốt thuộc Trường Đại học Lâm nghiệp (thị trấn Xuân Mai, huyện Chương Mỹ, Hà Nội; tọa độ 20°51’13” vĩ độ Bắc và 105°30’45” kinh độ Đông), việc nâng cao năng suất rừng trồng đang gặp nhiều trở ngại do độ phì đất biến động phức tạp và tốc độ sinh trưởng chiều cao của các loài cây chủ lực như Keo, Thông chỉ đạt mức trung bình khoảng 0,8 đến 0,9 m/năm. Vấn đề nghiên cứu trọng tâm là xác định cấu trúc phân bố không gian và mức độ biến thiên của các chỉ tiêu lý, hóa tính đất nhằm khắc phục những sai số lớn trong công tác điều tra thổ nhưỡng truyền thống.

Mục tiêu cụ thể của luận văn là định lượng hóa đặc điểm phân bố không gian của các tính chất vật lý (dung trọng, độ ẩm, thành phần cơ giới) và tính chất hóa học (hàm lượng chất hữu cơ, đạm dễ tiêu NH4+, lân dễ tiêu P2O5) trên diện tích lâm phần điển hình 6.000 m2 ở dải độ cao từ 76 đến 133 m. Đồng thời, đề tài đánh giá tính đại diện của các ô tiêu chuẩn 1.000 m2 và độ tin cậy của phương pháp lấy mẫu tổng hợp so với phương pháp lấy mẫu ngẫu nhiên độc lập. Ý nghĩa khoa học và thực tiễn của công trình thể hiện ở việc cung cấp cơ sở dữ liệu địa thống kê chuẩn xác, giúp giảm thiểu sai số lấy mẫu xuống dưới 5%, tối ưu hóa thiết kế mạng lưới điều tra lập địa và hỗ trợ lựa chọn loài cây trồng bản địa thích hợp nhằm phục hồi hệ sinh thái rừng bền vững.

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 hệ thống lý thuyết phát sinh học đất của Dokuchaev (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ẹ, sinh vật, khí hậu, địa hình và thời gian. Trong đó, thảm thực vật nhiệt đới đóng vai trò chủ đạo khi hàng năm bổ sung cho đất khoảng 2 đến 10 tấn chất hữu cơ trên mỗi hecta.

Song song với đó, nghiên cứu kế thừa lý thuyết mùn hóa hiện đại của Stevenson (1982) và Waksman, làm sáng tỏ các con đường chuyển hóa tàn tích sinh vật thành hợp chất mùn cao phân tử thông qua các phản ứng sinh hóa của vi sinh vật và các enzyme oxy hóa.

Ba khái niệm then chốt được chuẩn hóa trong nghiên cứu gồm:

  1. Dung trọng đất: Khối lượng của một đơn vị thể tích đất khô kiệt ở trạng thái tự nhiên (g/cm3), phản ánh độ chặt và độ xốp của cấu trúc đất.
  2. Dinh dưỡng khoáng dễ tiêu: Hàm lượng các ion nitơ amon (NH4+) và phosphat (P2O5) hòa tan mà rễ cây có thể hấp thu trực tiếp (mg/100g đất).
  3. Cấu trúc không gian địa thống kê (Variogram): Hàm số biểu diễn mối quan hệ phụ thuộc giữa phương sai sai khác của các tính chất đất với cự ly không gian giữa các điểm quan trắc.

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

Nghiên cứu thu thập dữ liệu thực nghiệm trên lâm phần 6.000 m2 phát triển trên nền đá mẹ Poocfiarit. Tổng cộng 108 điểm lấy mẫu đất tầng mặt (0 - 20 cm) được định vị theo mạng lưới hệ thống ngẫu nhiên bằng máy kinh vĩ, địa bàn ba chân và thước thép. Để so sánh phương pháp điều tra, 2 ô tiêu chuẩn diện tích 1.000 m2 được thiết lập để lấy mẫu riêng lẻ và mẫu tổng hợp theo hai sơ đồ: đường Ziczac (9 điểm, khoảng cách giữa các tuyến 13,3 m) và sơ đồ tán xạ từ tâm (9 điểm, gồm tâm ô và 8 điểm cách tâm 10 m).

Các mẫu đất được xử lý bằng cối đồng, rây qua lưới 1 mm và phân tích theo các phương pháp chuẩn quốc gia tại Trung tâm Thí nghiệm thực hành Khoa Lâm học:

  • Chất hữu cơ (OM%) xác định bằng phương pháp oxy hóa ướt Tyurin với hỗn hợp K2Cr2O7 và H2SO4 ở nhiệt độ 150°C.
  • Đạm dễ tiêu (NH4+) chiết bằng dung dịch KCl và so màu bằng thuốc thử Nessler trên máy quang phổ kế.
  • Lân dễ tiêu (P2O5) xác định bằng phương pháp so màu đo quang phổ.
  • Dung trọng xác định qua ống đóng dung tích 100 cm3; độ ẩm đo bằng phương pháp sấy khô kiệt ở 105°C; thành phần cơ giới phân tích theo phương pháp hút ống Pipet (phân cấp hạt sét <0,002 mm, thịt 0,002 - 0,02 mm, cát >0,02 mm).

Số liệu được phân tích bằng mô hình bán phương sai (Semivariogram), thuật toán nội suy khoảng cách nghịch đảo (IDW), kiểm định giả thuyết t-Student và phân tích phương sai F-Fisher ở mức ý nghĩa alpha = 0,05 trên phần mềm chuyên dụng.

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 108 mẫu đất tầng mặt trên toàn lâm phần ghi nhận các chỉ số định lượng cơ bản:

  1. Hàm lượng chất hữu cơ (OM%) trung bình đạt 2,16% (dao động từ 0,77% đến 3,39%), hệ số biến động Cv = 24,5%, sai tiêu chuẩn 0,53. Phân bố thực nghiệm có dạng gần đối xứng, trong đó hơn 75% số mẫu tập trung ở dải giá trị từ 1,6% đến 2,8%, xếp đất khu vực vào nhóm có độ phì hữu cơ trung bình.
  2. Dung trọng đất trung bình đạt 1,47 g/cm3 (dao động từ 0,93 đến 1,93 g/cm3), hệ số biến động Cv = 13,1%, sai tiêu chuẩn 0,19. Mật độ mẫu tập trung dày đặc nhất trong khoảng 1,35 đến 1,55 g/cm3, chiếm gần 90% dung lượng mẫu, cho thấy nền đất có độ chặt tương đối cao.
  3. Độ ẩm đất mặt trung bình đạt 0,30% (dao động từ 0,17% đến 0,46%), hệ số biến động Cv = 22,4%, sai tiêu chuẩn 0,07, với xu hướng ẩm giảm dần rõ rệt từ sườn Tây sang sườn Đông của lâm phần.
  4. Đất rất nghèo dinh dưỡng khoáng dễ tiêu trên 36 mẫu phân tích chuyên sâu: Đạm dễ tiêu NH4+ trung bình chỉ đạt 1,99 mg/100g đất (dao động từ 0,01 đến 5,10 mg/100g, Cv = 60,6%), có tới 75% vị trí lấy mẫu dưới 2,75 mg/100g. Hàm lượng lân dễ tiêu P2O5 trung bình đạt 1,40 mg/100g đất (dao động từ 0,50 đến 4,43 mg/100g, Cv = 55,6%), đa số tập trung dưới ngưỡng 1,25 mg/100g.
  5. Thành phần cơ giới có sự phân hóa: Cấp hạt cát chiếm trung bình 36,08%, cấp hạt sét chiếm 32,43%, và cấp hạt thịt chiếm 31,49%. Nhóm đất thịt pha sét (Clay loam) chiếm ưu thế tuyệt đối với 64% tổng số mẫu.

Thảo luận kết quả

Cấu trúc không gian địa thống kê qua đồ thị Variogram chỉ ra rằng phương sai của hàm lượng chất hữu cơ và dung trọng đất dao động rất mạnh ngay ở khoảng cách dưới 5 m. Điều này phản ánh tác động chi phối của các yếu tố ngẫu nhiên vi địa hình và sự phân bố rác rơi rụng cục bộ. Khoảng cách hữu hiệu để phương sai chất hữu cơ đạt 95% giá trị cực đại là 55 m (với phương sai đạt 0,29), trong khi khoảng cách hữu hiệu của độ ẩm đất là 50 m. Khi trực quan hóa dữ liệu qua biểu đồ phân bố tần suất và bản đồ mặt cầu nội suy 2D, các đốm suy thoái cục bộ hiện rõ, chứng minh các tính chất lý hóa không phân bố đồng nhất mà tạo thành các vi vũng sinh thái rải rác.

Khi so sánh giữa ô tiêu chuẩn 1.000 m2 và toàn bộ lâm phần 6.000 m2:

  • Hàm lượng chất hữu cơ trong ô tiêu chuẩn 1 (2,31%) và ô tiêu chuẩn 2 (2,32%) không có sự khác biệt mang ý nghĩa thống kê so với phần còn lại của lâm phần (2,10%) với giá trị kiểm định t tương ứng đạt p = 0,069 và p = 0,077 (>0,05). Phân tích kiểm định F về phương sai cũng cho kết quả tương đồng (p > 0,05), chứng minh ô tiêu chuẩn đại diện tốt cho hàm lượng chất hữu cơ và dung trọng.
  • Ngược lại, độ ẩm đất trung bình giữa ô tiêu chuẩn 1 (0,34%) và ô tiêu chuẩn 2 (0,26%) có sự khác biệt rất lớn so với lâm phần (p = 0,00003 < 0,05). Nguyên nhân do sự chênh lệch về độ tàn che tầng cây cao (0,7 ở ô 1 so với 0,5 ở ô 2) và độ che phủ thảm tươi (80% so với 65%), làm thay đổi trực tiếp cường độ bốc hơi vật lý bề mặt.
  • Đối với phương pháp lấy mẫu tổng hợp (Ziczac và tán xạ), kiểm định t cho thấy các trị số thu được đều sai khác đáng kể so với trung bình của các mẫu ngẫu nhiên độc lập ở mức ý nghĩa 0,05. Phương pháp lấy mẫu tổng hợp đã làm triệt tiêu tính biến thiên cực trị và che lấp các dị biệt không gian thực tế của đất rừng.

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

Nhằm nâng cao hiệu quả quản lý dinh dưỡng đất và quy hoạch lâm sinh, các giải pháp cụ thể được đề xuất như sau:

  1. Chuẩn hóa mạng lưới điều tra thổ nhưỡng: Áp dụng phương pháp bố trí điểm lấy mẫu ngẫu nhiên dạng lưới với khoảng cách tối thiểu giữa các điểm từ 50 đến 55 m nhằm bao quát trọn vẹn biến động không gian; loại bỏ phương pháp lấy mẫu tổng hợp đơn điểm trong các khảo sát lập địa chi tiết để kiểm soát sai số đo lường dưới 5%. Giải pháp này do các Trường Đại học và Viện Nghiên cứu Lâm nghiệp thực hiện trong giai đoạn 2026 - 2027.
  2. Bổ sung dinh dưỡng khoáng chuyên biệt: Tiến hành bón lót phân lân nung chảy kết hợp phân vi sinh cố định đạm với định mức 200 đến 300 kg/ha cho rừng trồng mới, nhằm mục tiêu nâng hàm lượng lân dễ tiêu P2O5 từ 1,40 mg/100g lên trên 3,50 mg/100g đất và cải thiện đạm dễ tiêu NH4+ trong 2 năm đầu chu kỳ canh tác. Chủ thể thực hiện là các Ban Quản lý rừng phòng hộ và doanh nghiệp trồng rừng.
  3. Chuyển đổi cơ cấu cây trồng theo hướng hỗn giao bản địa: Giảm mật độ trồng thuần loài Thông, Keo; tăng cường đưa các loài cây họ Đậu và cây gỗ lớn bản địa như Lim xanh, Lim xẹt, Sưa Bắc Bộ vào cấu trúc rừng với mật độ 200 đến 300 cây/ha. Mục tiêu duy trì độ tàn che tán trên 0,7, nâng độ che phủ thảm tươi trên 80% và giảm dung trọng đất xuống dưới 1,35 g/cm3 trong giai đoạn 5 đến 10 năm tới. Chi cục Kiểm lâm phối hợp cùng các chủ rừng chịu trách nhiệm tổ chức.
  4. Duy trì thảm mục và chống xói mòn tầng mặt: Tuyệt đối cấm các hoạt động đốt dọn thực bì toàn diện; áp dụng kỹ thuật làm đất cục bộ theo đường đồng mức để bổ sung từ 5 đến 8 tấn chất hữu cơ/ha/năm, giúp nâng tỷ lệ OM% từ mức 2,16% lên trên 3,0% vào năm 2030 dưới sự giám sát của lực lượng quản lý rừng cơ sở.

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

Nghiên cứu mang lại giá trị thực tiễn và học thuật sâu sắc cho 4 nhóm đối tượng chính:

  1. Giảng viên và nhà nghiên cứu khoa học đất lâm nghiệp: Khai thác bộ cơ sở dữ liệu thực nghiệm gồm 108 điểm đo chi tiết; ứng dụng mô hình Variogram và thuật toán nội suy không gian để phân tích quy luật thoái hóa đất nhiệt đới gió mùa.
  2. Kỹ sư quy hoạch lâm sinh và chuyên viên điều tra tài nguyên rừng: Tiếp cận các phân tích định lượng về sai số của phương pháp lấy mẫu tổng hợp; ứng dụng vào việc thiết kế mạng lưới ô tiêu chuẩn 1.000 m2 chuẩn xác cho các dự án kiểm kê tài nguyên đất quy mô lớn.
  3. Ban Quản lý các khu rừng thực nghiệm và Vườn Quốc gia: Nắm vững mối tương quan giữa độ tàn che (0,5 - 0,7) và lớp thảm tươi với độ ẩm đất mặt; xây dựng giải pháp kỹ thuật phục hồi rừng tự nhiên và bảo vệ lưu vực đầu nguồn hiệu quả.
  4. Học viên cao học và sinh viên ngành Lâm học, Quản lý tài nguyên rừng: Sử dụng luận văn như một tài liệu tham khảo chuẩn mực về phương pháp nghiên cứu khoa học, kỹ thuật xử lý thống kê nâng cao (kiểm định t, kiểm định F, semivariogram) và cách thức lập luận logic trong bảo vệ luận văn thạc sĩ.

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

  1. Tại sao hàm lượng chất hữu cơ trong đất tại núi Luốt chỉ đạt mức trung bình 2,16%? Nguyên nhân chủ yếu do tiền sử đồi trọc bị xói mòn rửa trôi kéo dài trong quá khứ, kết hợp với điều kiện nhiệt đới ẩm thúc đẩy quá trình khoáng hóa diễn ra nhanh hơn tích lũy mùn. Mặc dù rừng trồng đã phục hồi nhưng lượng thảm mục rơi rụng tích lũy chưa đủ lâu, dẫn đến biên độ dao động OM% từ 0,77% đến 3,39% giữa các vi địa hình.

  2. Phương pháp lấy mẫu đất tổng hợp bộc lộ những hạn chế kỹ thuật nào? Phương pháp trộn mẫu tổng hợp (dù theo đường Ziczac hay tán xạ từ tâm) làm bình quân hóa số liệu nhân tạo, che khuất hoàn toàn các giá trị cực trị và tính bất đồng nhất của lập địa. Kiểm định thống kê t cho thấy kết quả của mẫu tổng hợp luôn sai khác có ý nghĩa so với giá trị trung bình thực tế ở mức xác suất p < 0,05.

  3. Dung trọng đất trung bình 1,47 g/cm3 phản ánh tình trạng vật lý đất như thế nào? Trị số dung trọng 1,47 g/cm3 (với 64% mẫu thuộc nhóm thịt pha sét và hàm lượng sét đạt 32,43%) cho thấy đất có độ chặt tương đối lớn và độ xốp thấp. Hiện tượng này xuất phát từ quá trình nén chặt tự nhiên của đất Feralit sau giai đoạn mất rừng và sự xuất hiện của các lớp kết von nông từ độ sâu 20 cm.

  4. Mối quan hệ giữa độ tàn che tán rừng và độ ẩm tầng đất mặt diễn ra như thế nào? Độ tàn che đóng vai trò như một màng chắn bức xạ nhiệt 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 đạt 0,34%, cao hơn rõ rệt so với ô tiêu chuẩn có độ tàn che 0,5 với độ ẩm chỉ đạt 0,26% (p = 0,00003), do bốc hơi vật lý bề mặt bị hạn chế mạnh dưới tán rừng khép kín.

  5. Cần áp dụng chế độ bón phân gì cho cây rừng trên đất Feralit phát triển trên đá Poocfiarit? Do đất rất nghèo dinh dưỡng khoáng (đạm NH4+ chỉ đạt 1,99 mg/100g và lân P2O5 chỉ đạt 1,40 mg/100g), cần bón lót phân lân nung chảy để hạn chế sự cố định lân bởi oxit sắt, nhôm, đồng thời kết hợp bón phân hữu cơ vi sinh nhằm kích thích hoạt động phân giải dinh dưỡng của hệ sinh vật đất.

Kết luận

  • Đất Feralit tầng mặt (0 - 20 cm) tại khu rừng thực nghiệm Núi Luốt có hàm lượng chất hữu cơ trung bình 2,16%, dung trọng trung bình 1,47 g/cm3 và độ ẩm trung bình 0,30%.
  • Dinh dưỡng khoáng dễ tiêu ở mức nghèo rõ rệt với đạm NH4+ trung bình 1,99 mg/100g đất và lân P2O5 trung bình 1,40 mg/100g đất; thành phần cơ giới ưu thế là thịt pha sét chiếm 64%.
  • Phân bố không gian của các tính chất lý hóa đất mang tính không đồng nhất cao, với khoảng cách biến động bão hòa không gian đạt từ 50 đến 55 m.
  • Ô tiêu chuẩn 1.000 m2 có thể đại diện tốt cho toàn lâm phần 6.000 m2 về chất hữu cơ và dung trọng, nhưng không đại diện được cho chỉ tiêu độ ẩm đất.
  • Phương pháp lấy mẫu tổng hợp gây sai lệch thống kê có ý nghĩa so với phương pháp lấy mẫu ngẫu nhiên độc lập trong điều tra đất rừng.

Đóng góp lớn nhất của công trình là việc định lượng hóa chính xác cấu trúc không gian của đất Feralit dưới tán rừng thực nghiệm và thiết lập luận cứ khoa học để hoàn thiện quy trình điều tra lập địa lâm nghiệp. Trong giai đoạn 2026 - 2028, các nghiên cứu tiếp theo cần mở rộng quan trắc xuống các tầng đất sâu từ 20 đến 100 cm trên nhiều dạng lập địa dốc khác nhau. Các đơn vị quản lý tài nguyên và nghiên cứu lâm nghiệp cần ứng dụng ngay các phát hiện này vào thực tế để chuẩn hóa quy trình điều tra thổ nhưỡng và tối ưu hóa các mô hình phục hồi rừng bền vững.