Tổng quan nghiên cứu

Tại Việt Nam, diện tích đất trống đồi núi trọc và thoái hóa từng chiếm khoảng 6,16 triệu ha (tương đương 18,59% tổng diện tích tự nhiên theo thống kê lâm nghiệp năm 2008). Trong bối cảnh độ che phủ rừng cả nước đạt 38,7%, việc đẩy mạnh trồng rừng sản xuất bằng các loài cây mọc nhanh như Keo tai tượng (Acacia mangium Willd.) giữ vai trò then chốt nhằm đáp ứng nhu cầu nguyên liệu gỗ công nghiệp và phục hồi thảm che phủ. Tuy nhiên, phương thức canh tác thuần loài luân kỳ ngắn tiềm ẩn nhiều nguy cơ làm suy giảm độ phì nhiêu của đất nếu không có chiến lược quản lý lập địa khoa học.

Mục tiêu trọng tâm của nghiên cứu là xác định quy luật tương tác hai chiều giữa các chỉ tiêu sinh trưởng sinh khối rừng trồng Keo tai tượng thuần loài và một số tính chất lý hóa học của đất. Đề tài khảo sát chuyên sâu tại 12 điểm nghiên cứu thuộc hai vùng sinh thái trọng điểm gồm tỉnh Tuyên Quang (đại diện vùng Trung tâm Bắc Bộ) và tỉnh Thừa Thiên Huế (đại diện vùng Bắc Trung Bộ), với các lâm phần có độ tuổi từ 5 đến 16 năm.

Ý nghĩa khoa học và thực tiễn của công trình thể hiện ở việc định lượng chính xác các chỉ số lý hóa học đất then chốt quyết định đến năng suất rừng, đồng thời đánh giá mức độ biến đổi của lập địa sau các giai đoạn sinh trưởng so với trạng thái rừng thứ sinh liền kề. Nghiên cứu cung cấp hệ thống cơ sở dữ liệu thực nghiệm chuẩn xác, hỗ trợ đắc lực cho công tác quy hoạch lập địa, dự báo sản lượng và xây dựng giải pháp kỹ thuật lâm sinh bền vững cho hơn 320.000 ha rừng keo trên toàn quốc.

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 tuần hoàn dinh dưỡng sinh học trong hệ sinh thái rừng trồng và lý thuyết phân hạng lập địa lâm nghiệp. Cơ chế cộng sinh cố định đạm khí quyển giữa rễ cây họ Đậu (Fabaceae) với vi khuẩn nốt sần Rhizobium và nấm rễ Mycorrhizal đóng vai trò trung tâm trong quá trình cải tạo độ phì nhiêu đất.

Các khái niệm chuyên ngành then chốt được lượng hóa gồm:

  • Dung trọng đất (Bulk Density - BD, tính bằng g/cm³): Chỉ số phản ánh độ chặt, độ xốp và độ thoáng khí của tầng đất mặt.
  • Tổng cacbon hữu cơ đất (Soil Organic Carbon - SOC, tính bằng %): Thước đo hàm lượng mùn và khả năng lưu trữ cacbon của lập địa.
  • Tổng cation bazơ trao đổi (Total Exchangeable Base Cations - TEC, tính bằng cmol/kg đất): Biểu thị dung lượng hấp phụ khoáng kiềm trao đổi như canxi, magie, kali.
  • Độ chua hoạt tính (pH nước) và hàm lượng Lân dễ tiêu (P2O5, tính bằng mg/kg đất): Các yếu tố dinh dưỡng giới hạn trực tiếp đến quá trình sinh trưởng của thực vật.

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

Nghiên cứu thiết lập hệ thống mẫu điều tra tại 12 lâm phần rừng trồng Keo tai tượng thuần loài (7 điểm tại Tuyên Quang và 5 điểm tại Thừa Thiên Huế). Đối với nội dung đánh giá biến đổi đất theo kiểu sử dụng, 6 cặp lập địa (3 cặp tại Tuyên Quang và 3 cặp tại Thừa Thiên Huế) được chọn ngẫu nhiên có rừng trồng keo và rừng thứ sinh liền kề cùng lịch sử canh tác.

Cỡ mẫu và phương pháp thu thập dữ liệu cụ thể:

  • Mỗi lâm phần bố trí ngẫu nhiên 3 ô tiêu chuẩn (OTC) độc lập với kích thước 20m x 25m (diện tích 500 m²/ô), tổng cộng thực hiện trên 36 OTC rừng trồng và 18 OTC rừng thứ sinh đối chứng.
  • Trong từng OTC, tiến hành đo đếm toàn bộ đường kính ngang ngực (D1.3) và chiều cao vút ngọn (Hvn) của tất cả các cây. Sinh khối cây cá lẻ được tính theo phương trình sinh khối của Phan Minh Sáng (2008) cho Keo tai tượng: Sinh khối = 0,00373 * D^2,9211; và phương trình của Brown (1997) cho rừng thứ sinh.
  • Lấy mẫu đất tại tầng mặt 0 - 10 cm theo sơ đồ 5 điểm đào phẫu diện ngẫu nhiên (1 điểm trung tâm, 4 điểm phụ tại 4 hướng) rồi trộn đều thành 1 mẫu hỗn hợp đại diện cho mỗi OTC. Mẫu đất được sấy khô ở 65°C và phân tích bằng các thiết bị hiện đại như máy phân tích nguyên tố dò hồng ngoại CNS-2000 (LECO) và phương pháp chuẩn hóa quốc tế.

Lý do lựa chọn phương pháp phân tích: Nghiên cứu áp dụng Mô hình ảnh hưởng hỗn hợp tuyến tính (Linear Mixed-effects Model - LME) trên phần mềm thống kê R (gói nlme, HH, multcomp). Phương pháp LME là chuẩn mực phân tích dữ liệu phân cấp hiện đại, cho phép tách bạch rõ ràng giữa biến cố định (Kiểu sử dụng đất) và biến ngẫu nhiên (Lập địa - Site), loại trừ triệt để sai số do hiện tượng giả lặp (pseudoreplication) mà các mô hình phân tích phương sai ANCOVA truyền thống thường gặp phải.

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

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

Quá trình phân tích thực nghiệm trên toàn bộ hệ thống ô tiêu chuẩn tại hai tỉnh Tuyên Quang và Thừa Thiên Huế đã xác định 4 phát hiện quan trọng:

  1. Biến động sinh khối và tốc độ tăng trưởng: Sinh khối rừng trồng Keo tai tượng tăng nhanh theo tuổi. Tại Tuyên Quang, sinh khối tích lũy dao động từ 47,58 ± 3,36 tấn/ha ở rừng 7 tuổi lên đến 144,70 ± 3,65 tấn/ha ở rừng 16 tuổi (trung bình toàn vùng đạt 100,87 ± 2,84 tấn/ha). Tại Thừa Thiên Huế, sinh khối biến thiên từ 24,90 ± 4,99 tấn/ha ở tuổi 5 đến 115,30 ± 3,21 tấn/ha ở tuổi 9 (trung bình đạt 56,11 ± 2,65 tấn/ha). Tốc độ tăng trưởng sinh khối bình quân hàng năm giữa hai vùng có sự tương đồng lớn, đạt 8,48 tấn/ha/năm tại Tuyên Quang và 8,24 tấn/ha/năm tại Thừa Thiên Huế.

  2. Cải thiện tính chất vật lý đất: Trị số dung trọng đất (BD) có xu hướng giảm tỷ lệ nghịch với tuổi rừng, giúp tầng đất mặt ngày càng tơi xốp. Ở Tuyên Quang, dung trọng đất duy trì ở mức tối ưu từ 1,02 đến 1,21 g/cm³. Tại Thừa Thiên Huế, dung trọng biến động từ 1,11 đến 1,51 g/cm³, trong đó các lâm phần keo 5 - 6 tuổi có dung trọng cao hơn (1,28 - 1,51 g/cm³) gây độ chặt bí nhất định cho rễ cây hô hấp.

  3. Tích lũy hàm lượng Cacbon và Nitơ: Tổng cacbon hữu cơ trong đất (SOC) tầng mặt dưới rừng keo tại Tuyên Quang đạt trung bình 2,52 ± 0,32% (dao động từ 2,30% đến 2,79%), cao hơn đáng kể so với mức trung bình 1,85 ± 0,13% (dao động từ 1,45% đến 2,54%) tại Thừa Thiên Huế. Rừng keo lớn tuổi trả lại lượng mùn và đạm hữu cơ dồi dào qua thảm mục rơi rụng hàng năm.

  4. Hiện tượng chua hóa đất dưới tán keo: Độ pH đất dao động từ 3,73 ± 0,09 đến 4,97 ± 0,11 (đất từ rất chua đến chua vừa). Đất dưới các lâm phần keo thuần loài 5 - 6 tuổi thể hiện phản ứng chua mạnh nhất (như lập địa TQ6 đạt pH 3,73; TTH5 đạt pH 4,23), phản ánh rõ nét sự suy giảm tổng cation bazơ trao đổi (TEC) ở các giai đoạn cây sinh trưởng sung mãn.

Thảo luận kết quả

Mức tăng trưởng sinh khối bình quân 8,24 - 8,48 tấn/ha/năm ghi nhận trong nghiên cứu phản ánh đúng tiềm năng sinh học của Keo tai tượng trên các lập địa đồi núi dốc vùng nhiệt đới gió mùa. So sánh với các công bố của Ngô Đình Quế (11,14 tấn/ha/năm) và Vũ Tấn Phương (15,53 tấn/ha/năm), năng suất trong nghiên cứu này thấp hơn do điều kiện lập địa nghiên cứu đa phần là đất thoái hóa nghèo dinh dưỡng và mức độ quản lý lâm sinh còn hạn chế.

Dữ liệu được thể hiện trực quan thông qua các bảng thống kê mô tả sai lệch chuẩn (SD), hệ số biến động (S%) và hệ thống đồ thị tương quan hồi quy tuyến tính trong phần mềm R. Biểu đồ phân tán (scatter plot) cho thấy mối tương quan đồng biến chặt chẽ giữa tuổi rừng với sinh khối và hàm lượng cacbon đất, đồng thời biểu thị tương quan nghịch giữa sinh khối với dung trọng đất.

Nguyên nhân của hiện tượng chua hóa đất bắt nguồn từ việc Keo tai tượng sinh trưởng nhanh đã hấp thụ một lượng rất lớn các cation kiềm (Ca²⁺, Mg²⁺, K⁺) tích lũy vào thân, cành và lá. Khi lượng rơi rụng chưa kịp khoáng hóa hoàn toàn hoặc bị rửa trôi do lượng mưa lớn (lên tới 3.056,7 mm/năm tại Thừa Thiên Huế), đất mặt bị thiếu hụt cation bazơ trầm trọng, làm giải phóng các ion Al³⁺ và H⁺ gây chua đất.

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

Nhằm nâng cao năng suất rừng trồng và bảo vệ độ phì nhiêu đất hướng tới phát triển lâm nghiệp bền vững, nghiên cứu đề xuất 4 nhóm giải pháp kỹ thuật cụ thể:

  1. Áp dụng kỹ thuật khai thác bảo tồn vật liệu hữu cơ: Tuyệt đối nghiêm cấm biện pháp đốt dọn thực bì sau khai thác trắng. Chủ rừng cần để lại toàn bộ cành ngọn, lá mục và vỏ cây rải đều trên mặt đất nhằm hoàn trả tối thiểu 30% lượng cation bazơ và nitơ hữu cơ cho đất. Biện pháp này do các chủ hộ trồng rừng và công ty lâm nghiệp triển khai ngay trong giai đoạn 2026 - 2030.

  2. Cải tạo độ chua và bổ sung dinh dưỡng khoáng: Tiến hành bón lót định kỳ từ 350 - 500 kg vôi bột (CaCO3)/ha kết hợp bón 150 - 200 kg lân nung chảy (P2O5)/ha ở giai đoạn rừng 1 - 3 tuổi. Mục tiêu nâng độ pH đất tầng mặt lên trên ngưỡng 4,5 nhằm hạn chế độc tính nhôm di động và tăng cường khả năng hoạt động của vi khuẩn nốt sần. Chủ thể thực hiện là cán bộ kỹ thuật lâm sinh và ban quản lý rừng sản xuất.

  3. Chuyển hóa rừng gỗ nhỏ sang kinh doanh rừng gỗ lớn: Kéo dài luân kỳ khai thác từ chu kỳ ngắn 5 - 7 năm sang luân kỳ 10 - 15 năm. Mục tiêu nâng tổng sinh khối rừng vượt mức 130 - 150 tấn/ha, giúp hệ rễ phát triển sâu, giảm dung trọng đất xuống dưới 1,05 g/cm³ và tối ưu hóa tích lũy cacbon hữu cơ trong đất đạt trên 2,8%. Sở Nông nghiệp và Môi trường các địa phương cần định hướng chính sách hỗ trợ tín dụng dài hạn cho các hộ lâm nghiệp.

  4. Xây dựng mô hình rừng trồng hỗn giao cây bản địa: Giảm tỷ lệ trồng keo thuần loài bằng cách thiết lập các đai rừng hỗn giao, đưa 25 - 35% cây họ Đậu bản địa hoặc cây gỗ lớn (như Lim xẹt, Mỡ, Re gừng) xen kẽ với Keo tai tượng. Giải pháp này giúp đa dạng hóa quần xã vi sinh vật đất, giảm áp lực suy thoái dinh dưỡng khoáng và hạn chế sâu bệnh hại trên diện rộng giai đoạn 2026 - 2035.

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

Nội dung và bộ số liệu thực nghiệm trong luận văn mang lại giá trị tham khảo thiết thực cho 4 nhóm đối tượng chính:

  • Cơ quan quản lý nhà nước và nhà hoạch định chính sách lâm nghiệp: Cung cấp cơ sở dữ liệu định lượng về tăng trưởng sinh khối và biến động lập địa để xây dựng quy hoạch phân vùng trồng rừng kinh tế, tiêu chuẩn kỹ thuật phục hồi đất thoái hóa tại khu vực miền Bắc và Bắc Trung Bộ.
  • Kỹ sư lâm sinh và các chủ rừng, doanh nghiệp trồng rừng: Nắm vững phương pháp tính toán sinh khối cá lẻ, quy trình bón phân cân đối, kỹ thuật quản lý vật rơi rụng và kiểm soát độ chua đất nhằm tối ưu hóa chu kỳ kinh doanh gỗ nguyên liệu.
  • Giảng viên, nghiên cứu sinh và sinh viên ngành Lâm nghiệp - Môi trường: Tiếp cận phương pháp luận nghiên cứu hiện đại, nắm vững kỹ thuật phân tích mô hình ảnh hưởng hỗn hợp tuyến tính (LME) trên phần mềm R để ứng dụng vào xử lý dữ liệu đo đếm lặp trong sinh thái học rừng.
  • Các tổ chức phi chính phủ và chuyên gia dự án tín chỉ cacbon (REDD+/FSC): Sử dụng phương trình sinh khối và hệ số tích lũy SOC làm tài liệu kiểm chứng khoa học khi đánh giá trữ lượng cacbon rừng và chứng chỉ quản lý rừng bền vững.

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

1. Trồng Keo tai tượng thuần loài liên tục có làm suy thoái đất lâm nghiệp không?
Rừng Keo tai tượng giúp cải thiện rõ rệt dung trọng đất (đạt 1,02 - 1,21 g/cm³) và tăng lượng đạm, cacbon hữu cơ tầng mặt nhờ vi khuẩn nốt sần cố định đạm. Tuy nhiên, nếu khai thác trắng chu kỳ ngắn kèm theo đốt thực bì, đất sẽ bị suy giảm cation bazơ nghiêm trọng và tăng độ chua nhanh chóng.

2. Tại sao độ pH đất dưới tán rừng Keo tai tượng lại giảm mạnh ở các tuổi đầu?
Giai đoạn 5 - 7 tuổi là thời kỳ keo sinh trưởng sinh khối nhanh nhất (tăng trưởng trên 8 tấn/ha/năm). Cây hấp thu ồ ạt canxi, magie và kali từ dung dịch đất để tạo thân gỗ, đồng thời rễ tiết ra axit hữu cơ, khiến đất mặt trở nên rất chua với độ pH giảm xuống dưới 4,0.

3. Tốc độ tăng trưởng sinh khối Keo tai tượng giữa Tuyên Quang và Thừa Thiên Huế có sự khác biệt lớn không?
Tốc độ tăng trưởng sinh khối bình quân hàng năm giữa hai vùng tương đương nhau, đạt 8,48 tấn/ha/năm tại Tuyên Quang và 8,24 tấn/ha/năm tại Thừa Thiên Huế. Tuy nhiên, tổng sinh khối tích lũy ở Tuyên Quang cao hơn nhờ hệ thống rừng điều tra có nhiều cấp tuổi lớn (lên đến 16 tuổi).

4. Ưu điểm nổi bật của mô hình ảnh hưởng hỗn hợp (LME) so với ANCOVA trong phân tích đất rừng là gì?
Mô hình LME phân tách rạch ròi biến cố định (kiểu sử dụng đất) và biến ngẫu nhiên (lập địa - Site), xử lý triệt để hiện tượng sai số do giả lặp trên các ô tiêu chuẩn cùng lâm phần. Điều này giúp kết luận thống kê chính xác tuyệt đối mà phân tích ANCOVA truyền thống không làm được.

5. Cần lưu ý điều gì nhất về mặt phân bón khi tái trồng Keo tai tượng luân kỳ 2?
Cần đặc biệt chú trọng bón lót bổ sung vôi bột (350 - 500 kg/ha) để khử chua tầng mặt và bổ sung phân lân nung chảy (150 - 200 kg/ha). Điều này giúp hóa giải hiện tượng lân bị cố định bởi ion nhôm di động trong điều kiện đất chua pH < 4,5.

Kết luận

  • Luận văn đã xác định chính xác quy luật sinh trưởng sinh khối của Keo tai tượng thuần loài với mức tăng trưởng trung bình đạt 8,24 - 8,48 tấn/ha/năm tại hai vùng sinh thái Tuyên Quang và Thừa Thiên Huế.
  • Quá trình gây trồng Keo tai tượng mang lại tác động tích cực trong việc làm giảm dung trọng đất (xuống mức 1,02 - 1,21 g/cm³) và nâng cao hàm lượng cacbon hữu cơ (SOC đạt 1,85 - 2,52%) cùng đạm tổng số.
  • Nghiên cứu chỉ ra nguy cơ hiện hữu về tình trạng chua hóa đất tầng mặt (pH dao động 3,73 - 4,97) do cây hấp thụ mạnh các cation bazơ kiềm trong suốt chu kỳ sinh trưởng.
  • Đóng góp khoa học then chốt là việc tiên phong ứng dụng thành công mô hình ảnh hưởng hỗn hợp tuyến tính (LME) bằng ngôn ngữ R trong việc giải bài toán tương tác đất - cây rừng tại Việt Nam.
  • Trong giai đoạn 2026 - 2030, các địa phương cần nhanh chóng chuyển đổi cơ cấu lâm sinh sang mô hình kinh doanh gỗ lớn kết hợp trồng hỗn giao bản địa để bảo vệ vẹn toàn độ phì đất rừng. Hãy áp dụng ngay các giải pháp kỹ thuật quản lý lập địa bền vững để tối đa hóa hiệu quả kinh tế và sinh thái cho các vùng trồng keo!