Tổng quan nghiên cứu

Trong hệ sinh thái rừng nhiệt đới, chu trình sinh địa hóa vận hành liên tục thông qua hai quá trình đối lập nhưng thống nhất: hấp thu dinh dưỡng để kiến tạo sinh khối và hoàn trả vật chất hữu cơ về lại đất thông qua thảm vật rơi rụng. Tại Vườn Quốc gia Cúc Phương – khu bảo tồn thiên nhiên rộng hơn 22.200 ha với hơn 2.000 loài thực vật bậc cao – lớp thảm mục đóng vai trò như một cơ chế tự bón phân vĩnh cửu. Tuy nhiên, trước áp lực phát triển kinh tế và nhu cầu gỗ nguyên liệu ngày càng tăng, nhiều diện tích rừng tự nhiên giàu trữ lượng đang bị thay thế bởi các mô hình rừng trồng thuần loài có chu kỳ khai thác ngắn (khoảng 5 đến 10 năm) và thiếu biện pháp thâm canh cải tạo đất. Thực trạng này dẫn đến nguy cơ thoái hóa đất nghiêm trọng, suy giảm 20% đến 35% độ phì nhiêu và làm giảm sút năng suất cây trồng ở các luân kỳ tiếp theo.

Đề tài nghiên cứu được thực hiện nhằm giải quyết bài toán khoa học và thực tiễn cấp bách: định lượng chính xác lượng vật rơi rụng mới, lượng tồn đọng thảm khô tích lũy, tốc độ phân giải chất hữu cơ và động thái hoàn trả các nguyên tố khoáng đa lượng (Nitơ, Phốt pho, Kali) giữa các trạng thái rừng. Nghiên cứu tiến hành so sánh đối chứng giữa ba kiểu rừng đặc trưng tại Cúc Phương gồm: Rừng tự nhiên trạng thái phục hồi IIIA3 tại Tiểu khu 20, Rừng trồng thuần loài Xoan đào 14 tuổi và Rừng trồng hỗn giao Phay vi - Xoan đào - Côm bàng 14 đến 15 tuổi tại Vườn thực vật. Quá trình theo dõi thực địa được triển khai trong giai đoạn chuyển tiếp khí hậu Xuân - Hạ (từ tháng 3 đến tháng 7). Kết quả nghiên cứu cung cấp cơ sở dữ liệu định lượng tin cậy giúp các nhà quản lý rừng thiết lập giải pháp kỹ thuật lâm sinh, nâng cao 30% đến 40% tỷ lệ tái sinh tự nhiên từ hạt và bảo tồn độ phì của đất rừng một cách 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

Nghiên cứu được xây dựng dựa trên nền tảng lý thuyết tuần hoàn vật chất và năng lượng trong hệ sinh thái rừng do các học giả kinh điển như Bray, Rodin và Basilevic thiết lập. Theo lý thuyết này, vật rơi rụng (litterfall) là tập hợp toàn bộ sinh khối chết của thực vật rụng xuống bề mặt đất gồm lá, cành, hoa, quả, vỏ cây và xác bã sinh vật, đóng vai trò là mắt xích khởi đầu cho chu trình khoáng hóa. Tốc độ chuyển hóa vật chất được xác định thông qua hệ số phân giải thảm mục K của Rodin (tỷ số giữa lượng thảm mục tồn dư trung bình trên mặt đất với lượng thảm mục rơi rụng hàng năm), phản ánh đặc trưng sinh thái của từng đới rừng.

Bên cạnh đó, nghiên cứu kế thừa quan điểm của các nhà khoa học lâm nghiệp Việt Nam như Hoàng Xuân Tý, Nguyễn Ngọc Lung và Vũ Đại Hải về chức năng bảo vệ đất của thảm mục hữu cơ. Lớp thảm khô được định nghĩa như một ngân hàng tiết kiệm dinh dưỡng khoáng, đồng thời là tấm khiên thủy văn có tác dụng triệt tiêu 80% đến 90% động năng của hạt mưa, ngăn ngừa hiện tượng nén chặt đất mặt và giảm thiểu xói mòn rửa trôi. Khung phân tích tập trung vào ba khái niệm then chốt: động thái rơi rụng mới (lượng sinh khối bổ sung định kỳ), lượng tồn đọng thảm khô (lượng vật chất tích lũy chưa phân hủy) và khả năng tự bón phân (lượng nguyên tố dinh dưỡng đa lượng N, P, K được giải phóng sau quá trình khoáng hóa).

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

Phương pháp nghiên cứu kết hợp chặt chẽ giữa điều tra sinh thái thực địa định lượng và phân tích hóa sinh trong phòng thí nghiệm:

  • Cỡ mẫu và phương pháp thu thập: Tại trạng thái rừng tự nhiên IIIA3, tác giả thiết lập 20 điểm thu mẫu cố định; tại mỗi mô hình rừng trồng (thuần loài và hỗn giao) bố trí 10 điểm thu mẫu. Tại mỗi điểm, đặt một bẫy lưới chuyên dụng diện tích 0,2 m2 (kích thước mắt lưới 1,2 x 1,2 mm), căng cách mặt đất 30 cm nhằm thu gom triệt để lá, hoa, quả rơi rụng với chu kỳ 15 ngày một lần. Đối với vật rơi rụng là cành nhỏ (đường kính dưới 2 cm), nghiên cứu bố trí 30 ô tiêu chuẩn dạng bản diện tích 4 m2 ở rừng tự nhiên và 15 ô ở mỗi kiểu rừng trồng. Lượng thảm mục tồn đọng được thu thập bằng khung tròn sắt tiêu chuẩn tại 20 điểm rừng tự nhiên và 10 điểm rừng trồng. Độ dày tầng thảm mục được đo đạc chính xác bằng thước khắc vạch milimet tại 150 điểm ở rừng tự nhiên và 50 điểm ở mỗi kiểu rừng trồng.
  • Phương pháp phân tích tốc độ phân giải: Sử dụng phương pháp túi lưới phân hủy (litter bag) kích thước 15 x 15 cm chứa 100 gam mẫu lá hoặc cành khô kiệt (đã sấy ở nhiệt độ 60°C đến khối lượng không đổi). Mẫu được đặt trực tiếp dưới tán rừng trong thời gian 45 ngày để xác định độ suy giảm khối lượng sinh học.
  • Phân tích hóa học: 18 mẫu thực vật đại diện (gồm 9 mẫu lá và 9 mẫu cành, mỗi mẫu 100 gam) được phân tích xác định hàm lượng Nitơ tổng số bằng phương pháp Kjeldahl, Phốt pho dễ tiêu bằng phương pháp so màu trắc quang và Kali trao đổi bằng phương pháp quang kế ngọn lửa phát xạ.
  • Thử nghiệm sinh học tái sinh: Bố trí các ô thí nghiệm diện tích 0,5 m2 với 100 hạt Xoan đào đã qua tuyển chọn kỹ lưỡng, gieo dưới các lớp phủ thảm mục có độ dày kiểm soát từ 0 cm (đối chứng) đến 6 cm.
  • Xử lý số liệu: Toàn bộ dữ liệu được tổng hợp và kiểm định sai khác thống kê bằng tiêu chuẩn phi tham số Kruskal-Wallis trên phần mềm Excel. Lý do lựa chọn kiểm định phi tham số là vì dữ liệu vật rơi rụng sinh thái trong điều kiện tự nhiên thường phân bố lệch chuẩn, cỡ mẫu giữa các trạng thái có sự biến thiên lớn, đòi hỏi kiểm định Kruskal-Wallis để đảm bảo tính chuẩn xác và độ tin cậy khoa học cao nhất.

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 đã đưa ra các phát hiện khoa học mang tính quy luật rõ rệt giữa cấu trúc rừng và động thái tuần hoàn dinh dưỡng:

  1. Động thái rơi rụng mới: Tổng lượng vật rơi rụng mới ước tính cho 1 ha rừng dao động từ 0,5151 tấn/ha đến 0,6275 tấn/ha. Trong đó, rừng tự nhiên đạt giá trị cao nhất với 0,6275 tấn/ha (bao gồm 0,5476 tấn/ha lá và 0,0799 tấn/ha cành nhỏ). Rừng trồng thuần loài Xoan đào đạt 0,5591 tấn/ha (0,5104 tấn/ha lá và 0,0487 tấn/ha cành). Rừng trồng hỗn giao đạt thấp nhất với 0,5151 tấn/ha (0,4739 tấn/ha lá và 0,0412 tấn/ha cành). Kết quả kiểm định Kruskal-Wallis cho thấy giá trị H tính toán đối với vật rơi rụng lá đạt 503,61 và đối với cành đạt trên 5,99, đều vượt trội so với giá trị tra bảng chuẩn ở mức ý nghĩa 0,05. Điều này khẳng định sự sai khác về lượng rơi rụng mới giữa 3 kiểu rừng có ý nghĩa thống kê rõ rệt.
  2. Trữ lượng thảm khô tồn đọng: Trái ngược với lượng rơi rụng mới, lượng thảm mục tồn đọng trên bề mặt đất đạt giá trị lớn nhất ở rừng trồng thuần loài Xoan đào với 4,886 tấn/ha, cao hơn rừng tự nhiên (4,644 tấn/ha) và vượt xa rừng trồng hỗn giao (4,063 tấn/ha). Kiểm định Kruskal-Wallis ghi nhận H = 547,79 > 5,99, chứng minh sự phân hóa rõ ràng về lượng tích lũy thảm mục giữa các mô hình rừng.
  3. Tốc độ phân giải sinh khối: Trong thời gian theo dõi 45 ngày dưới điều kiện nhiệt ẩm cao, tốc độ phân giải diễn ra rất mạnh mẽ. Tỷ lệ phân hủy của vật rơi rụng lá đạt từ 36,94% đến 61,08% (trong đó rừng trồng thuần loài phân giải nhanh nhất với 61,08%, rừng tự nhiên đạt 51,40% và rừng hỗn giao đạt 36,94%). Tỷ lệ phân hủy của cành nhỏ đạt từ 23,46% đến 41,40% (rừng tự nhiên phân giải cành nhanh nhất với 41,40%). Trung bình sau 45 ngày, có trên 30% lượng lá rụng và 25% lượng cành rụng được khoáng hóa hoàn toàn.
  4. Động thái dinh dưỡng khoáng: Hàm lượng Nitơ luôn chiếm tỷ lệ ưu thế tuyệt đối trong cả lá và cành ở cả 3 kiểu rừng. Cụ thể, hàm lượng Nitơ trong lá đạt từ 0,977% đến 1,248% (cao nhất ở rừng tự nhiên với 1,248%), trong khi ở cành đạt từ 0,318% đến 0,840%. Hàm lượng Phốt pho trong lá dao động từ 0,096% đến 0,105%, trong cành đạt 0,045% đến 0,090%. Hàm lượng Kali trong lá đạt 0,035% đến 0,076%, trong cành đạt 0,013% đến 0,036%. Lượng Nitơ bồi hoàn lý thuyết cho đất rừng tự nhiên ước tính lên tới 13,102 kg/ha/tháng.
  5. Ảnh hưởng đến tái sinh tự nhiên: Thí nghiệm nảy mầm hạt Xoan đào khẳng định lớp phủ thảm mục đóng vai trò điều tiết độ ẩm quyết định. Ở ô đối chứng trơ trụi (độ dày 0 cm), tỷ lệ hạt nảy mầm thấp nhất. Tỷ lệ nảy mầm đạt cực đại khi độ dày thảm mục duy trì ở mức 3 cm đối với rừng trồng và 4 cm đối với rừng tự nhiên. Khi độ dày thảm phủ vượt quá 5 cm đến 6 cm, tỷ lệ nảy mầm suy giảm mạnh do hạt bị ngăn cách tiếp xúc với đất khoáng và dễ bị nấm bệnh gây hại.

Thảo luận kết quả

Sự phân hóa về lượng vật rơi rụng và tốc độ phân giải được chi phối trực tiếp bởi đặc điểm cấu trúc tổ thành loài và điều kiện vi khí hậu tại Vườn Quốc gia Cúc Phương. Trong thời gian nghiên cứu (tháng 3 đến tháng 7), nhiệt độ không khí trung bình duy trì từ 22°C đến 32°C kết hợp độ ẩm tương đối từ 90% đến 97%, tạo môi trường lý tưởng cho hoạt động của vi sinh vật đất và nấm phân giải xenlulozo.

Đặc biệt, yếu tố địa hình núi đá vôi karst có tác động sâu sắc đến sự sắp xếp không gian của lớp thảm khô. Tại các vị trí sườn dốc từ 8° đến 10° với tỷ lệ đá lộ đầu lên đến 80%, lớp thảm mục chỉ đạt độ dày mỏng từ 0,5 cm đến 1,2 cm (nhiều điểm trơ sỏi đá bằng 0 cm) do nước mưa cuốn trôi sinh khối xuống chân dốc. Ngược lại, tại các vùng lòng thung lũng tích tụ màu nâu đen, tầng thảm mục lắng đọng dày tới 3 cm đến 3,5 cm, cá biệt có điểm đạt 7 cm đến 8 cm. Sự phân bố không đều này giải thích vì sao sinh trưởng của tầng cây cao ở thung lũng (đường kính thân cây 30 cm đến 40 cm, chiều cao vút ngọn 20 m đến 25 m) luôn vượt trội so với cây trên sườn dốc (đường kính chỉ đạt 15 cm đến 20 cm, chiều cao 15 m đến 22 m).

Dữ liệu nghiên cứu về cấu trúc thảm mục và tần suất độ dày có thể được mô hình hóa trực quan thông qua các biểu đồ phân bố tần suất (Histogram) kết hợp bảng ma trận phân tích phương sai Kruskal-Wallis. Cách tiếp cận trực quan này giúp các nhà quản lý lâm nghiệp dễ dàng xác định ngưỡng che phủ an toàn, phân biệt rõ các khu vực tích tụ sinh khối hữu cơ cao với các khu vực xói mòn nghiêm trọng để đưa ra biện pháp can thiệp kỹ thuật chính xác.

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

Dựa trên các bằng chứng thực nghiệm về chu trình vật chất và động thái dinh dưỡng, luận văn đưa ra 4 nhóm giải pháp kỹ thuật có tính ứng dụng cao:

  1. Thiết lập quy chế bảo vệ nghiêm ngặt lớp thảm mục hữu cơ mặt đất: Ban Quản lý Vườn Quốc gia và chính quyền địa phương cần ban hành lệnh cấm triệt để hành vi quét dọn, thu gom lá khô cành rụng để làm chất đốt hoặc đốt dọn thực bì tùy tiện tại 15 xã vùng đệm. Biện pháp này nhằm bảo toàn trọn vẹn nguồn hoàn trả tự nhiên từ 10 kg đến 13 kg Nitơ/ha/tháng, giảm thiểu 40% đến 50% nguy cơ rửa trôi đất mặt trên sườn dốc đá vôi. Thời gian thực hiện: Triển khai liên tục hàng năm, lồng ghép vào các hương ước bảo vệ rừng cấp thôn bản.
  2. Ứng dụng kỹ thuật điều tiết độ dày thảm mục trong xúc tiến tái sinh rừng: Cán bộ kỹ thuật lâm sinh cần chủ động can thiệp vào tầng thảm khô trước mùa phát tán hạt giống (từ tháng 3 đến tháng 6). Duy trì độ dày lớp phủ thảm mục ở mức tối ưu từ 3 cm đến 4 cm nhằm giữ ẩm cho hạt; tiến hành xới vạt cục bộ ở những vị trí thảm mục tích tụ quá dày (trên 5 cm) để tạo điều kiện cho hạt giống tiếp xúc trực tiếp với đất khoáng, nâng cao 30% đến 40% tỷ lệ nảy mầm và định cư của cây tái sinh bản địa.
  3. Chuyển đổi căn bản cơ cấu rừng trồng từ thuần loài sang hỗn giao đa tầng tán: Các công ty lâm nghiệp và ban quản lý rừng phòng hộ cần chấm dứt việc trồng rừng độc canh chu kỳ ngắn. Khuyến nghị nhân rộng mô hình hỗn giao theo băng giữa các loài cây bản địa như Phay vi, Xoan đào và Côm bàng với mật độ 320 cây/ha. Mô hình này giúp duy trì độ tàn che trên 65%, phân bố lượng vật rơi rụng đồng đều quanh năm và giảm 50% nguy cơ bùng phát cháy rừng trong các đợt gió Tây Nam khô nóng.
  4. Xây dựng hệ thống quan trắc dinh dưỡng đất và vật liệu cháy định kỳ: Phòng Khoa học của Vườn Quốc gia cần thiết lập mạng lưới 30 đến 50 ô tiêu chuẩn định vị lâu dài để theo dõi động thái suy giảm khối lượng hữu cơ và độ ẩm thảm mục. Khi độ ẩm thảm khô xuống dưới mức 65% trong mùa khô, hệ thống quan trắc sẽ kích hoạt cảnh báo cháy rừng cấp 4 và cấp 5, hỗ trợ lực lượng kiểm lâm ứng phó kịp thời.

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

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

  • Các nhà nghiên cứu, nghiên cứu sinh chuyên ngành Lâm học, Sinh thái học rừng và Khoa học Đất: Tài liệu cung cấp quy trình phương pháp luận chuẩn mực về thiết lập bẫy thu vật rơi rụng, kỹ thuật túi lưới phân giải sinh khối và ứng dụng kiểm định thống kê phi tham số Kruskal-Wallis trong nghiên cứu sinh thái thực địa.
  • Ban Quản lý các Vườn Quốc gia và Khu bảo tồn thiên nhiên: Cán bộ quản lý có thể khai thác trực tiếp các số liệu định lượng về lượng bồi hoàn N, P, K và quy luật phân bố thảm mục theo địa hình karst để xây dựng phương án quản lý vật liệu cháy rừng và bảo tồn đa dạng sinh học tại hơn 30 khu bảo tồn có địa hình đá vôi tương đồng trên cả nước.
  • Kỹ sư lâm sinh và các đơn vị phát triển rừng trồng: Cung cấp giải pháp kỹ thuật lâm sinh cụ thể về mật độ trồng hỗn giao (320 cây/ha), công thức phối loài bản địa và phương pháp điều tiết thảm mục giúp tăng 25% đến 30% năng suất rừng trồng chu kỳ sau.
  • Giảng viên và sinh viên các trường đại học khối Nông Lâm nghiệp: Sử dụng như một tài liệu tham khảo điển hình cho các học phần Sinh thái học rừng nhiệt đới, Thổ nhưỡng rừng và Kỹ thuật xúc tiến tái sinh tự nhiên.

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

Tại sao lượng vật rơi rụng mới ở rừng tự nhiên lại cao hơn các mô hình rừng trồng?

Rừng tự nhiên trạng thái IIIA3 sở hữu cấu trúc không gian phức tầng hoàn chỉnh với mật độ tầng cây cao đạt từ 110 đến 120 cây/ha và tầng cây nhỡ đạt 100 đến 110 cây/ha. Tổ thành loài đa dạng gồm Vàng anh, Cà lồ, Phân mã tạo ra sinh khối tán lá dày đặc và chu kỳ thay lá đan xen liên tục, mang lại tổng lượng rơi rụng mới đạt 0,6275 tấn/ha, vượt trội so với cấu trúc đơn tầng thưa thoáng của rừng trồng thuần loài.

Lớp thảm mục rơi rụng đóng vai trò gì trong việc duy trì độ phì đất rừng?

Dưới tác động của nhiệt ẩm cao tại Cúc Phương, hơn 30% lượng lá rụng và 25% lượng cành rụng được phân giải chỉ sau 45 ngày. Quá trình khoáng hóa này chuyển hóa sinh khối chết thành mùn và hoàn trả trực tiếp cho đất khoảng 13,102 kg Nitơ/ha/tháng cùng các nguyên tố Phốt pho, Kali dễ tiêu, giúp duy trì độ phì tự nhiên mà không cần bón phân hóa học.

Vì sao độ dày thảm mục 3 cm đến 4 cm lại tối ưu nhất cho hạt Xoan đào nảy mầm?

Lớp thảm mục dày 3 cm ở rừng trồng và 4 cm ở rừng tự nhiên đóng vai trò như một lớp đệm giữ ẩm hoàn hảo, duy trì độ ẩm ổn định cho hạt hút nước trương nở và kích thích nấm vi sinh có ích hoạt động. Nếu thảm mục dày quá 5 cm, hạt sẽ bị treo lơ lửng, mất tiếp xúc với đất khoáng và dễ thối hỏng do nấm bệnh tích tụ.

Sự khác biệt về tốc độ phân hủy giữa lá và cành cây diễn ra như thế nào?

Sau 45 ngày theo dõi thực tế, tốc độ phân giải của lá đạt từ 36,94% đến 61,08%, cao hơn đáng kể so với cành nhỏ chỉ đạt từ 23,46% đến 41,40%. Nguyên nhân là do lá cây chứa hàm lượng Nitơ cao (0,977% đến 1,248%), màng tế bào mềm giúp vi sinh vật dễ dàng xâm nhập và bẻ gãy cấu trúc nhanh hơn nhiều so với mô gỗ giàu linhin của cành.

Việc ứng dụng kiểm định Kruskal-Wallis mang lại lợi ích khoa học gì cho đề tài?

Kiểm định phi tham số Kruskal-Wallis xử lý hiệu quả các tập dữ liệu sinh thái phân bố lệch chuẩn và không đồng nhất về cỡ mẫu (như so sánh giữa 20 điểm rừng tự nhiên và 10 điểm rừng trồng). Với giá trị H tính toán đạt 503,61 đối với lá và 547,79 đối với lượng tồn đọng, kiểm định khẳng định tính chính xác tuyệt đối của sự sai khác sinh thái giữa các kiểu rừng.

Kết luận

  • Luận văn đã định lượng hóa toàn diện động thái vật rơi rụng với tổng lượng mới đạt từ 0,5151 đến 0,6275 tấn/ha và lượng thảm mục tích lũy đạt từ 4,063 đến 4,886 tấn/ha tại ba kiểu rừng điển hình ở Vườn Quốc gia Cúc Phương.
  • Chứng minh tốc độ phân giải chất hữu cơ diễn ra nhanh chóng dưới điều kiện nhiệt đới ẩm, đạt tỷ lệ phân hủy trung bình trên 30% đối với lá và 25% đối với cành sau 45 ngày theo dõi.
  • Khẳng định vai trò quyết định của nguồn Nitơ hoàn trả tự nhiên (đạt 13,102 kg/ha/tháng ở rừng tự nhiên) trong việc duy trì chu trình sinh địa hóa và độ phì nhiêu của đất rừng.
  • Xác lập ngưỡng độ dày thảm mục tối ưu từ 3 cm đến 4 cm giúp tối đa hóa tỷ lệ nảy mầm và xúc tiến tái sinh tự nhiên cho loài cây bản địa Xoan đào.
  • Cung cấp cơ sở khoa học vững chắc để các nhà quản lý chuyển đổi phương thức lâm sinh từ trồng rừng độc canh sang mô hình hỗn giao đa tầng bản địa bền vững.

Để phát triển sâu hơn các kết quả của luận văn, trong giai đoạn 12 đến 24 tháng tới, các cơ quan nghiên cứu cần mở rộng phạm vi đánh giá chu trình dinh dưỡng khoáng trong toàn bộ 12 tháng của năm và phân tích chi tiết hoạt tính của hệ vi sinh vật đất chuyên trách phân giải xenlulozo. Hãy liên hệ với thư viện nhà trường hoặc cơ quan chủ quản để tiếp cận toàn văn công trình nghiên cứu và đồng hành ứng dụng các giải pháp kỹ thuật lâm sinh tiên tiến vào thực tiễn quản lý rừng bền vững ngay hôm nay!