Tổng quan nghiên cứu

Hệ sinh thái rừng tự nhiên tại vùng chuyển tiếp giữa Tây Nguyên và Đông Nam Bộ giữ vai trò xung yếu trong việc duy trì cân bằng sinh thái, bảo vệ nguồn nước và bảo tồn nguồn gen sinh học quý hiếm. Vườn quốc gia Bù Gia Mập sở hữu tổng diện tích tự nhiên 26.032 ha, trong đó diện tích đất có rừng đạt 25.536 ha, chiếm tỷ lệ che phủ lên tới 98,1%. Nơi đây ghi nhận sự hiện diện phong phú của 724 loài thực vật bậc cao có mạch thuộc 362 chi, 109 họ và 359 loài động vật có xương sống. Tuy nhiên, các tác động lịch sử từ khai thác chọn và áp lực nương rẫy canh tác du canh trước đây đã làm suy thoái đáng kể cấu trúc lâm phần, làm xáo trộn quy luật tái sinh tự nhiên và làm giảm sút trữ lượng các loài cây gỗ bản địa có giá trị kinh tế cao.

Nghiên cứu tập trung giải quyết vấn đề cấp bách về việc thiếu hụt cơ sở dữ liệu định lượng phục vụ công tác phục hồi rừng tại phân khu phục hồi sinh thái. Mục tiêu cụ thể là làm sáng tỏ quy luật cấu trúc tầng tán, phân bố số cây theo đường kính, phân bố theo chiều cao và đặc điểm tái sinh tự nhiên tại hai trạng thái rừng chiếm diện tích lớn là trạng thái phục hồi nghèo (IIB) và trạng thái trung bình (IIIA2). Nghiên cứu được triển khai thực nghiệm trên phạm vi 22 ô tiêu chuẩn định vị tại tỉnh Bình Phước. Kết quả nghiên cứu cung cấp các chỉ số định lượng then chốt, hỗ trợ ban quản lý thiết lập các giải pháp lâm sinh chính xác, nâng cao độ tàn che từ mức 0,43 lên trên 0,70 và thúc đẩy tiến trình diễn thế rừng hướng đến trạng thái ổn định lâu dài.

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 vận dụng hệ sinh thái học rừng nhiệt đới kết hợp các nguyên lý lâm học thực nghiệm của nhiều nhà khoa học trong và ngoài nước. Khung lý thuyết cấu trúc tầng tán kế thừa quan điểm phân chia cấu trúc thẳng đứng của Richards và Thái Văn Trừng, phân định lâm phần thành các tầng vượt tán (A1), tầng ưu thế sinh thái (A2), tầng dưới tán (A3), tầng cây bụi (B) và tầng thảm tươi (C). Để lượng hóa cấu trúc quần xã thực vật, nghiên cứu tích hợp các mô hình toán học giải tích như hàm phân bố giảm Mayer, hàm phân bố Weibull và hàm phân bố khoảng cách nhằm mô phỏng quy luật phân bố số cây theo đường kính ngang ngực và chiều cao vút ngọn.

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

  • Tiết diện ngang lâm phần: Tổng diện tích mặt cắt thân cây tại vị trí 1,3 m tính trên một đơn vị diện tích rừng.
  • Chỉ số giá trị quan trọng: Chỉ số tổng hợp giữa tỷ lệ phần trăm mật độ và tỷ lệ phần trăm tiết diện ngang nhằm xác định vị thế sinh thái của từng loài cây trong quần thể.
  • Chỉ số đa dạng sinh học Simpson và Shannon-Wiener: Các đại lượng toán học phản ánh mức độ phong phú và tính đồng đều của các loài trong quần xã.
  • Hình thái phân bố không gian: Quy luật sắp xếp vị trí tương đối giữa các cá thể cây rừng trên mặt phẳng nằm ngang theo chỉ số Clark-Evans.

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

Nguồn dữ liệu sơ cấp được thu thập thông qua điều tra thực địa theo phương pháp lập ô tiêu chuẩn điển hình trên tuyến đại diện cho hai trạng thái rừng IIB và IIIA2. Cỡ mẫu nghiên cứu bao gồm 22 ô tiêu chuẩn với diện tích mỗi ô là 2.000 m² (kích thước 40m x 50m), tương đương tổng diện tích lấy mẫu tầng cây cao là 4,4 ha. Trong mỗi ô tiêu chuẩn, tiến hành bố trí 5 ô dạng bản diện tích 16 m² (kích thước 4m x 4m) tại 4 góc và 1 ô ở tâm, tạo thành hệ thống 110 ô dạng bản để theo dõi cây tái sinh, cây bụi và thảm tươi. Phương pháp chọn mẫu điển hình được lựa chọn nhằm phản ánh chính xác các đặc trưng lâm học đặc thù của từng trạng thái diễn thế mà không bị phân tán bởi các yếu tố nhiễu địa hình cực đoan.

Phương pháp phân tích số liệu kết hợp giữa thống kê sinh học và mô hình hóa toán học trên nền tảng phần mềm Excel và SPSS. Trước khi gộp số liệu lâm phần, nghiên cứu áp dụng tiêu chuẩn phi tham số Mann-Whitney U để so sánh từng cặp ô tiêu chuẩn và kiểm định Kruskal-Wallis để kiểm tra tính thuần nhất tổng thể về đường kính và chiều cao. Tiêu chuẩn Chi bình phương Pearson được dùng để kiểm định mức độ tương thích giữa tần số thực nghiệm và các hàm phân bố lý thuyết. Lý do lựa chọn các công cụ thống kê này là vì cấu trúc rừng tự nhiên nhiệt đới có độ biến động cao, dữ liệu thường không tuân theo quy luật phân bố chuẩn hoàn hảo, đòi hỏi các thuật toán kiểm định phi tham số để bảo đảm tính khách quan 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

Thứ nhất, cấu trúc tầng cây cao có sự phân hóa rõ rệt giữa hai trạng thái rừng nghiên cứu. Tại trạng thái IIB, mật độ cây dao động từ 164 đến 358 cây/ha, đường kính ngang ngực bình quân đạt 17,11 cm, chiều cao vút ngọn bình quân đạt 11,22 m, tổng tiết diện ngang đạt 5,77 m²/ha và độ tàn che bình quân ở mức thấp 0,43. Ngược lại, trạng thái IIIA2 thể hiện tính chất của lâm phần phục hồi tốt hơn với mật độ từ 238 đến 344 cây/ha, đường kính bình quân đạt 21,95 cm, chiều cao bình quân đạt 16,12 m, tổng tiết diện ngang đạt 12,47 m²/ha và độ tàn che đạt 0,61. So sánh tương quan cho thấy tổng tiết diện ngang của trạng thái IIIA2 cao hơn trạng thái IIB tới 116,1%, phản ánh trữ lượng sinh khối vượt trội.

Thứ hai, cấu trúc phân tầng và tổ thành loài phản ánh rõ giai đoạn diễn thế. Trạng thái IIB chỉ hình thành 4 tầng rõ rệt, thiếu vắng tầng vượt tán A1; trong đó tầng dưới tán A3 chiếm ưu thế áp đảo với 76,6% tổng số cá thể, tầng ưu thế sinh thái A2 chỉ chiếm 14,5% và tầng cây bụi B chiếm 8,9%. Trong khi đó, trạng thái IIIA2 đã định hình cấu trúc 5 tầng hoàn chỉnh với tầng vượt tán A1 cao 16 đến 20 m gồm các cây họ Dầu và họ Tử vi, tầng A2 phát triển mạnh với độ tàn che 0,60 gồm các loài Bằng lăng, Trâm, Trường và Chò chai.

Thứ ba, quy luật phân bố số cây theo cỡ kính và chiều cao tuân thủ các mô hình toán học xác định. Phân bố số cây theo đường kính ở cả hai trạng thái đều mang dạng đường cong giảm liên tục hình chữ J, mô phỏng thích hợp nhất qua hàm phân bố khoảng cách và phương trình Mayer. Phân bố cây theo chiều cao ở trạng thái IIIA2 có dạng một đỉnh lệch rõ nét, được nắn chỉnh tối ưu bằng hàm Weibull với độ tin cậy thống kê cao khi kiểm định qua tiêu chuẩn Chi bình phương Pearson.

Thứ tư, lớp cây tái sinh tự nhiên thể hiện tiềm năng phục hồi lớn nhưng phân bố không gian chưa đồng đều. Cây tái sinh chủ yếu có nguồn gốc từ hạt chiếm trên 80%, phân bố chủ yếu theo dạng cụm dưới các khoảng trống tán rừng với chỉ số Clark-Evans nhỏ hơn 1.

Thảo luận kết quả

Sự khác biệt về mật độ, độ tàn che và tiết diện ngang giữa trạng thái IIB và IIIA2 phản ánh tác động trực tiếp của cường độ can thiệp con người trong quá khứ. Trạng thái IIB hình thành chủ yếu sau nương rẫy bỏ hóa hoặc khai thác trắng, dẫn đến sự xâm chiếm ban đầu của các loài cây tiên phong ưa sáng, mọc nhanh như Thành ngạnh, Thẩu tấu, Cò ke với kích thước thân gỗ nhỏ tập trung ở cỡ kính 12 đến 18 cm. Trạng thái IIIA2 trải qua quá trình khai thác chọn nhẹ hơn nên vẫn giữ lại được các cá thể gỗ lớn thế hệ trước có đường kính trên 40 cm, đóng vai trò nguồn gieo giống quan trọng cho các loài cây chịu bóng và định cư dài ngày.

Khi biểu diễn dữ liệu thực nghiệm qua biểu đồ phân bố tần số đường kính, đường cong thực nghiệm của trạng thái IIIA2 thể hiện dạng răng cưa nhẹ do sự thiếu hụt cục bộ ở các cỡ kính trung bình từ 24 đến 32 cm. Hiện tượng này hoàn toàn phù hợp với các kết luận của các công trình lâm học nhiệt đới khi giải thích hậu quả của việc chặt chọn tập trung vào các cây gỗ thương phẩm có giá trị. Bảng ma trận tương quan sinh thái loài cho thấy mối liên kết dương chặt chẽ giữa các loài cây họ Dầu và họ Bằng lăng ở tầng ưu thế, tạo tiền đề vững chắc cho việc định hướng tái cấu trúc rừng tự nhiên thông qua các can thiệp lâm sinh chọn lọc.

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

Thứ nhất, áp dụng biện pháp khoanh nuôi xúc tiến tái sinh tự nhiên có tác động làm giàu rừng trên toàn bộ 14.553 ha diện tích rừng hỗn giao và rừng thứ sinh nghèo kiệt. Mục tiêu cụ thể là nâng độ tàn che trung bình từ 0,43 lên 0,70 và gia tăng mật độ cây tái sinh mục đích lên trên 2.500 cây/ha. Thời gian triển khai thực hiện liên tục trong chu kỳ 5 năm từ 2008 đến 2013, do Ban Quản lý Vườn quốc gia Bù Gia Mập làm chủ thể điều hành chính.

Thứ hai, thực hiện kỹ thuật tỉa thưa nuôi dưỡng và vệ sinh rừng cục bộ tại các lâm phần trạng thái IIIA2 có mật độ vượt quá 350 cây/ha. Tiến hành loại bỏ các cây phẩm chất xấu, cây bị sâu bệnh hoặc rỗng ruột chiếm khoảng 15% đến 20% lâm phần nhằm giải phóng không gian dinh dưỡng cho các loài cây bản địa có giá trị kinh tế cao như Sao đen, Dầu song nàng, Giáng hương và Gõ đỏ. Hoạt động này cần được giám sát nghiêm ngặt bởi Hạt Kiểm lâm Vườn quốc gia trong giai đoạn 2008-2010.

Thứ ba, thiết lập 05 chốt trạm bảo vệ rừng cơ động tại các khu vực đầu nguồn suối Đăk Huýt và Đăk Ca nhằm kiểm soát triệt để tình trạng khai thác gỗ trái phép bằng cưa xăng và vận chuyển lâm sản đường thủy. Phấn đấu giảm thiểu 100% các vụ vi phạm lâm luật trong vùng lõi trong vòng 2 năm đầu triển khai.

Thứ tư, duy trì và mở rộng chính sách giao khoán quản lý bảo vệ rừng cho các đơn vị quân đội và cộng đồng địa phương, tiếp tục quản lý hiệu quả 4.700 ha rừng tự nhiên đã giao cho Đồn biên phòng 785 và lực lượng dân quân hai xã Đăk Ơ, Bù Gia Mập. Định kỳ hằng quý tiến hành kiểm tra nghiệm thu sinh khối và đa dạng sinh thái để chi trả dịch vụ môi trường rừng đầy đủ.

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

Thứ nhất, ban quản lý các vườn quốc gia và khu bảo tồn thiên nhiên khu vực Đông Nam Bộ và Tây Nguyên. Luận văn cung cấp phương pháp luận thực chứng và các chỉ số cấu trúc cụ thể để xây dựng đề án quản lý rừng bền vững, phân vùng phục hồi sinh thái khoa học cho hơn 26.000 ha diện tích tương đồng.

Thứ hai, các nhà khoa học, giảng viên và nghiên cứu sinh chuyên ngành Lâm học, Sinh thái học và Quản lý tài nguyên rừng. Tài liệu mang lại bộ số liệu thực nghiệm phong phú về 724 loài thực vật cùng quy trình ứng dụng các mô hình toán học giải tích như hàm Mayer, Weibull trong mô phỏng cấu trúc lâm phần nhiệt đới.

Thứ ba, học viên cao học và sinh viên các trường đại học lâm nghiệp. Luận văn là hình mẫu chuẩn mực về phương pháp điều tra thực địa 22 ô tiêu chuẩn kết hợp 110 ô dạng bản, cũng như phương pháp xử lý số liệu phi tham số qua các kiểm định Mann-Whitney và Kruskal-Wallis.

Thứ tư, các chuyên viên hoạch định chính sách thuộc Sở Nông nghiệp và Phát triển nông thôn, Chi cục Kiểm lâm các tỉnh. Nghiên cứu cung cấp luận cứ thực tiễn để xây dựng quy chuẩn kỹ thuật phục hồi rừng thứ sinh và chính sách giao khoán bảo vệ rừng gắn với sinh kế đồng bào dân tộc thiểu số tại chỗ.

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

Cấu trúc tầng tán của trạng thái rừng phục hồi IIB có điểm gì khác biệt so với rừng nguyên sinh? Trạng thái IIB là rừng phục hồi sau nương rẫy và khai thác trắng nên cấu trúc tầng tán bị đơn giản hóa, chỉ có 4 tầng và hoàn toàn thiếu vắng tầng vượt tán A1. Tầng dưới tán A3 chiếm ưu thế áp đảo tới 76,6% tổng số cây với chiều cao từ 8 đến 15 m, trong khi tầng ưu thế A2 chỉ chiếm 14,5% và độ tàn che đạt mức thấp 0,43.

Tại sao tổng tiết diện ngang của trạng thái IIIA2 lại cao gấp hơn hai lần trạng thái IIB? Tổng tiết diện ngang của trạng thái IIIA2 đạt 12,47 m²/ha so với 5,77 m²/ha ở trạng thái IIB là do trạng thái IIIA2 còn lưu giữ được các cây gỗ lớn thế hệ trước có đường kính trên 40 cm thuộc họ Dầu và họ Tử vi. Ngược lại, trạng thái IIB chủ yếu gồm cây tiên phong đường kính nhỏ tập trung ở cỡ kính 12 đến 18 cm.

Mô hình toán học nào mô tả chính xác nhất quy luật phân bố số cây theo đường kính tại khu vực nghiên cứu? Quy luật phân bố số cây theo cỡ kính ở cả hai trạng thái đều tuân theo dạng đường cong giảm liên tục hình chữ J. Mô hình hàm phân bố khoảng cách và phương trình Mayer biểu diễn tốt nhất cấu trúc này với độ tin cậy cao, đã được kiểm chứng không có sự sai khác qua tiêu chuẩn Chi bình phương Pearson ở mức xác suất 0,05.

Hiện trạng đa dạng thực vật bậc cao tại Vườn quốc gia Bù Gia Mập được ghi nhận cụ thể ra sao? Nghiên cứu ghi nhận 724 loài thực vật thuộc 362 chi, 109 họ thuộc 6 ngành thực vật bậc cao. Trong đó, ngành Ngọc lan chiếm ưu thế tuyệt đối với 686 loài (chiếm 94,8% tổng số loài), có mặt 18 loài thực vật quý hiếm nằm trong Sách đỏ Việt Nam như Cẩm lai Bà Rịa, Gõ đỏ, Trầm hương và Giáng hương quả to.

Kỹ thuật lâm sinh nào được ưu tiên hàng đầu để phục hồi các diện tích rừng nghèo kiệt? Biện pháp kỹ thuật then chốt là khoanh nuôi xúc tiến tái sinh tự nhiên kết hợp trồng bổ sung làm giàu rừng bằng các loài cây bản địa có giá trị. Đồng thời tiến hành điều tiết độ tàn che lên trên 0,70 và thiết lập mạng lưới tuần tra nghiêm ngặt nhằm ngăn chặn hoàn toàn việc sử dụng cưa xăng khai thác gỗ bất hợp pháp.

Kết luận

  • Luận văn đã định lượng hóa toàn diện các đặc trưng cấu trúc lâm phần và tái sinh tự nhiên trên diện tích 4,4 ha ô tiêu chuẩn tại Vườn quốc gia Bù Gia Mập, làm sáng tỏ sự phân hóa sinh thái giữa hai trạng thái IIB và IIIA2.
  • Chứng minh tính phù hợp vượt trội của các hàm toán học giải tích như Mayer, Weibull và hàm phân bố khoảng cách trong việc mô phỏng quy luật phân bố số cây theo đường kính và chiều cao của rừng tự nhiên nhiệt đới.
  • Xác lập cơ sở khoa học sinh thái cho thấy trạng thái IIIA2 có tổng tiết diện ngang đạt 12,47 m²/ha, cao hơn 116,1% so với trạng thái IIB (5,77 m²/ha), chỉ ra hướng đi phục hồi rừng thông qua nuôi dưỡng thế hệ cây kế cận.
  • Đề xuất hệ thống 4 nhóm giải pháp kỹ thuật lâm sinh và quản trị bảo tồn đồng bộ, bao gồm làm giàu 14.553 ha rừng hỗn giao, tỉa thưa nuôi dưỡng và giao khoán bảo vệ 4.700 ha rừng trọng điểm.
  • Khuyến nghị các nhà quản lý lâm nghiệp và cơ quan nghiên cứu tiếp tục triển khai quan trắc định kỳ mạng lưới ô tiêu chuẩn định vị trong giai đoạn 2008-2015 nhằm đánh giá động thái diễn thế rừng dài hạn.