Tổng quan nghiên cứu

Vườn quốc gia Bù Gia Mập với tổng diện tích tự nhiên 26.032 ha và độ che phủ rừng đạt trên 90% là một trong những mẫu chuẩn hệ sinh thái rừng kín thường xanh và nửa rụng lá nhiệt đới ẩm quý hiếm còn sót lại tại vùng chuyển tiếp giữa Tây Nguyên và miền Đông Nam Bộ. Được thành lập chính thức theo Quyết định số 170/2002/QĐ-TTg của Thủ tướng Chính phủ, khu vực này lưu giữ 724 loài thực vật bậc cao và 359 loài động vật có xương sống. Tuy nhiên, áp lực khai thác gỗ trái phép trước đây kết hợp phương thức canh tác nương rẫy truyền thống đã làm suy thoái nghiêm trọng cấu trúc lâm phần, khiến nhiều phân khu phục hồi sinh thái rơi vào tình trạng thiếu hụt cây gỗ lớn và biến dạng quy luật tái sinh tự nhiên.

Mục tiêu cốt lõi của nghiên cứu là xác lập cơ sở khoa học định lượng về quy luật cấu trúc tầng tán, phân bố không gian và động thái tái sinh của thảm thực vật rừng tự nhiên lá rộng thường xanh. Đề tài tập trung điều tra chuyên sâu tại hai trạng thái rừng phục hồi chủ yếu gồm trạng thái phục hồi nghèo (IIB) và trạng thái trung bình (IIIA2) thuộc phân khu phục hồi sinh thái của Vườn quốc gia Bù Gia Mập trong giai đoạn 2004 - 2007.

Nghiên cứu mang ý nghĩa thực tiễn to lớn khi định lượng hóa các chỉ tiêu lâm học quan trọng như mật độ cây dao động từ 164 đến 358 cây/ha, tổng tiết diện ngang đạt 5,77 - 12,47 m²/ha. Kết quả phân tích cung cấp luận cứ kỹ thuật chuẩn xác cho công tác khoanh nuôi xúc tiến tái sinh, bảo tồn 18 loài thực vật nguy cấp trong Sách đỏ Việt Nam và nâng cao hiệu quả giao khoán quản lý 4.700 ha rừng tự nhiên cho các cộng đồng địa phươ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 trên nền tảng lý thuyết sinh thái học rừng nhiệt đới của Baur (1964) và hệ thống phân loại thảm thực vật rừng Việt Nam của Thái Văn Trừng (1978). Cấu trúc lâm phần được định nghĩa là sự sắp xếp không gian của các cá thể thực vật trong quần xã, phản ánh kết quả của quá trình chọn lọc tự nhiên và đấu tranh sinh tồn. Khung phân tầng thẳng đứng chia cấu trúc rừng thành 5 tầng cơ bản: 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 tổ thành, nghiên cứu áp dụng chỉ số giá trị quan trọng (Important Value - IV%) theo quan điểm của Daniel Marmillod và Curtis-McIntosh. Chỉ số IV% được tính bằng trung bình cộng giữa tỷ lệ phần trăm mật độ (N%) và tỷ lệ phần trăm tiết diện ngang (G%). Những loài có chỉ số IV% vượt quá 5% được xác định là loài có ý nghĩa sinh thái, và nhóm loài chiếm trên 50% tổng số cá thể được coi là nhóm loài ưu thế của lâm phần. Bên cạnh đó, các mô hình toán học như hàm phân bố khoảng cách của Nguyễn Hải Tuất (1986), mô hình hàm giảm Mayer (1934), và hàm phân bố Weibull 3 tham số được ứng dụng để mô phỏng quy luật phân bố số cây theo cỡ đường kính ngang ngực (N-D1.3) và chiều cao vút ngọn (N-Hvn).

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

Nghiên cứu thu thập dữ liệu thông qua hệ thống 22 ô tiêu chuẩn (OTC) điển hình với kích thước mỗi ô là 2.000 m² (40m x 50m), tổng diện tích điều tra đạt 44.000 m² trên hai trạng thái rừng IIB và IIIA2. Cỡ mẫu 22 ô tiêu chuẩn được lựa chọn nhằm đảm bảo tính đại diện cao cho các dạng lập địa đồi núi thấp từ độ cao 200m đến 738m. Trong mỗi OTC, tiến hành đo đếm toàn bộ cây gỗ có đường kính ngang ngực D1.3 từ 10 cm trở lên bằng thước kẹp kính, đo chiều cao vút ngọn (Hvn) bằng thước đo cao Blume-Leiss với độ chính xác 0,1m. Đồng thời, thiết lập 5 ô dạng bản (ODB) diện tích 16 m² (4m x 4m) tại 4 góc và tâm của mỗi OTC (tổng cộng 110 ODB) để thống kê chi tiết mật độ, nguồn gốc và phẩm chất cây tái sinh.

Lý do lựa chọn phương pháp phân tích thống kê toán học xuất phát từ tính phức tạp và không đồng nhất của rừng mưa hỗn loài. Trước khi tổng hợp số liệu, đề tài sử dụng kiểm định phi tham số Mann-Whitney U và chuẩn Kruskal-Wallis ở mức ý nghĩa 0,05 để kiểm tra tính thuần nhất giữa các OTC. Tiêu chuẩn Chi-square (χ²) của Pearson được sử dụng để kiểm tra độ phù hợp của các hàm phân bố lý thuyết, kết hợp phần mềm chuyên dụng SPSS và Excel để hồi quy tương quan H-D1.3 qua các dạng hàm tuyến tính, hyperbol, parabol bậc hai và hàm mũ.

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

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

Thứ nhất, về cấu trúc tầng tán và mật độ, trạng thái rừng IIB mang đặc trưng của rừng non phục hồi sau nương rẫy với mật độ dao động từ 164 đến 358 cây/ha, đường kính trung bình D1.3 đạt 17,11 cm, chiều cao vút ngọn trung bình 11,22 m và tổng tiết diện ngang chỉ đạt 5,77 m²/ha. Độ tàn che của trạng thái này ở mức thấp (0,4) và hoàn toàn vắng bóng tầng vượt tán A1. Tầng dưới tán A3 chiếm tỷ lệ áp đảo 76,6%, trong khi tầng ưu thế sinh thái A2 chỉ chiếm 14,5% với các loài tiên phong như Trâm, Gáo, Thẩu tấu, Thành ngạnh.

Thứ hai, trạng thái IIIA2 thể hiện cấu trúc rừng trung bình phục hồi sau khai thác chọn nhẹ 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 trung bình 16,12 m. Tổng tiết diện ngang của trạng thái IIIA2 đạt 12,47 m²/ha, cao hơn 116% so với trạng thái IIB, độ tàn che đạt từ 0,6 đến 0,8. Cấu trúc phân tầng tại trạng thái IIIA2 rất hoàn chỉnh khi xuất hiện tầng vượt tán A1 cao 16 - 20 m với các loài ưu thế thuộc họ Dầu (Dipterocarpaceae) như Sao đen, Dầu song nàng và họ Tử vi (Lythraceae) như Bằng lăng ổi.

Thứ ba, mô hình hóa phân bố đường kính (N-D1.3) ở cả hai trạng thái đều tuân theo quy luật giảm hình chữ J ngược, được mô phỏng phù hợp nhất qua hàm Mayer và hàm phân bố khoảng cách với giá trị kiểm định χ² tính toán đều nhỏ hơn χ² lý thuyết ở mức xác suất 0,05. Số lượng cây tập trung mật độ cao nhất ở cấp kính nhỏ 10 - 18 cm (chiếm trên 60% tổng số cây). Phân tích mạng hình phân bố mặt bằng qua chỉ số Clark & Evans cho kết quả Q < 1 và giá trị tiêu chuẩn U < -1,96, chứng minh các cá thể cây gỗ đang phân bố theo dạng cụm tập trung.

Thứ tư, lớp cây tái sinh dưới tán rừng phát triển mạnh mẽ với mật độ trung bình từ 3.500 đến 6.200 cây/ha. Cây tái sinh có nguồn gốc từ hạt chiếm trên 68%, trong đó cây phẩm chất tốt (loại A) chiếm tỷ trọng trên 65%, phân bố tập trung nhiều nhất ở cấp chiều cao dưới 1,0 m và giảm dần ở các cấp chiều cao lớn hơn.

Thảo luận kết quả

Các số liệu cấu trúc thu thập được có thể trình bày trực quan qua biểu đồ phân bố tần số thực nghiệm N-D1.3 kết hợp đường cong lý thuyết Mayer, cùng bảng so sánh tiết diện ngang giữa các lâm phần thuần nhất. Sự khác biệt rõ rệt về tổng tiết diện ngang (5,77 m²/ha ở trạng thái IIB so với 12,47 m²/ha ở trạng thái IIIA2) phản ánh mức độ tác động lịch sử của con người lên hệ sinh thái. Tại trạng thái IIB, việc khai thác trắng và nương rẫy cũ đã triệt hạ tầng cây cao, dẫn đến sự xâm lấn của các loài ưa sáng mọc nhanh. Ngược lại, trạng thái IIIA2 vẫn duy trì được các cây mẹ gieo giống cổ thụ có đường kính trên 40 cm, tạo điều kiện che bóng thuận lợi cho lớp cây tái sinh chịu bóng phát triển.

Khi so sánh với các công trình nghiên cứu của Đồng Sỹ Hiền (1974) và Thái Văn Trừng (1978) tại khu vực Đông Nam Bộ, kết quả của luận văn hoàn toàn tương đồng về quy luật tương quan đường kính - chiều cao H-D1.3 (biểu diễn tối ưu bằng hàm Parabol bậc 2 và hàm Power với R² > 0,75). Dạng phân bố không gian cụm (Q < 1) là chỉ dấu sinh thái điển hình của giai đoạn rừng đang phục hồi mạnh mẽ, nơi các cây con mọc tập trung quanh các vệt sáng hoặc gốc cây mẹ trước khi diễn ra quá trình tỉa thưa tự nhiên để tiến tới phân bố cách đều ở giai đoạn thành thục.

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

Dựa trên kết quả nghiên cứu định lượng cấu trúc và tái sinh, bốn giải pháp kỹ thuật lâm sinh và quản lý được đề xuất cụ thể như sau:

  1. Thực hiện 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: Áp dụng trực tiếp trên 1.500 ha rừng trạng thái IIB thuộc phân khu phục hồi sinh thái. Chủ thể thực hiện là Ban Quản lý Vườn quốc gia Bù Gia Mập phối hợp cùng chính quyền địa phương, đặt mục tiêu nâng độ tàn che từ 0,4 lên 0,7 và tăng mật độ cây mục đích lên tối thiểu 450 cây/ha trong lộ trình giai đoạn 2008 - 2015. Trồng bổ sung các loài cây bản địa có giá trị kinh tế và bảo tồn cao như Gõ đỏ, Giáng hương, Sao đen và Dầu rái.
  2. Áp dụng biện pháp tỉa thưa nuôi dưỡng và điều chỉnh không gian dinh dưỡng: Tiến hành phát dọn dây leo, chặt bỏ cây bụi chèn ép và loại bỏ các cây gỗ sâu bệnh, cong queo, bọng ruột tại các ô tiêu chuẩn có mật độ cụm quá dày. Mục tiêu chuyển đổi mạng hình phân bố từ dạng cụm tập trung (Q < 1) hướng dần tới phân bố ngẫu nhiên và cách đều (Q xấp xỉ 1), tạo không gian sinh trưởng tối ưu cho 300 - 400 cây mục đích/ha. Công việc này do lực lượng kỹ thuật lâm sinh Vườn quốc gia thực hiện định kỳ 2 năm/lần.
  3. Tăng cường chốt chặn kiểm soát và hiện đại hóa quản lý bảo vệ rừng: Thiết lập 5 trạm bảo vệ cố định tại các cửa sông suối trọng điểm như suối Đăk Huýt và suối Đăk Ca nhằm triệt phá phương thức vận chuyển gỗ xẻ trái phép bằng bè mảng vào mùa mưa. Hạt Kiểm lâm Bù Gia Mập phối hợp Đồn Biên phòng 785 tăng cường tuần tra liên ngành, đặt mục tiêu giảm 90% số vụ vi phạm lâm luật trước năm 2010.
  4. Nâng cao hiệu quả chương trình giao khoán đất rừng gắn với hỗ trợ sinh kế: Tiếp tục duy trì và mở rộng diện tích giao khoán 4.700 ha rừng tự nhiên cho Đồn Biên phòng 785 (1.700 ha), Xã đội Đăk Ơ (900 ha) và Xã đội Bù Gia Mập (2.100 ha). Đồng thời, Ủy ban nhân dân huyện Phước Long cần phối hợp triển khai các mô hình nông lâm kết hợp, hỗ trợ chuyển giao kỹ thuật trồng trọt cho 31 hộ dân vùng lõi và cộng đồng đồng bào dân tộc S'tiêng (chiếm 52,5% dân số vùng đệm), hướng tới mục tiêu nâng thu nhập bình quân đầu người vượt mức 300 kg thóc/năm để giảm thiểu áp lực xâm hại tài nguyên rừng.

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

  1. Ban quản lý các Vườn quốc gia và Khu bảo tồn thiên nhiên: Tài liệu cung cấp phương pháp luận chuẩn xác về điều tra lâm học định lượng và mô hình phân bố cây rừng, giúp các nhà quản lý xây dựng phương án quản lý rừng bền vững và thiết kế các phân khu chức năng phù hợp.
  2. Cán bộ kỹ thuật lâm nghiệp và lực lượng kiểm lâm: Cung cấp quy trình điều tra thực địa 22 ô tiêu chuẩn (2.000 m²) và 110 ô dạng bản tái sinh (16 m²), làm cẩm nang hướng dẫn kỹ thuật khoanh nuôi, làm giàu rừng và đánh giá diễn thế thảm thực vật sau khai thác.
  3. Giảng viên, nhà nghiên cứu và học viên cao học chuyên ngành Lâm học, Quản lý tài nguyên rừng và Sinh thái học: Nguồn tài liệu học thuật giá trị tham khảo về kỹ thuật xử lý dữ liệu phi tham số (Kruskal-Wallis, Mann-Whitney U), phương pháp nắn phân bố lý thuyết (Weibull, Mayer) và chỉ số đa dạng loài Shannon-Wiener trong nghiên cứu cấu trúc quần xã thực vật.
  4. Cơ quan quản lý nhà nước và hoạch định chính sách phát triển nông thôn: Cung cấp luận cứ khoa học và thực tiễn để xây dựng cơ chế chính sách giao khoán quản lý bảo vệ rừng (với case study 4.700 ha tại Bù Gia Mập), lồng ghép chương trình bảo tồn đa dạng sinh học với chính sách an sinh xã hội cho đồng bào thiểu số vùng biên giới.

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

Đặc điểm cấu trúc tầng tán giữa trạng thái rừng IIB và IIIA2 tại Vườn quốc gia Bù Gia Mập có điểm gì khác biệt? Trạng thái IIB có cấu trúc 4 tầng, không có tầng vượt tán A1, độ tàn che thấp (0,4) và tầng dưới tán A3 chiếm ưu thế tuyệt đối (76,6%). Trong khi đó, trạng thái IIIA2 có cấu trúc 5 tầng hoàn chỉnh, độ tàn che đạt 0,6 - 0,8, xuất hiện tầng vượt tán A1 cao 16 - 20 m với các loài ưu thế họ Dầu và Tử vi.

Tại sao luận văn lại sử dụng mô hình hàm Mayer và hàm phân bố khoảng cách để mô phỏng đường kính N-D1.3? Rừng tự nhiên tại Bù Gia Mập là rừng hỗn loài nhiệt đới có số cây tập trung lớn nhất ở cấp kính nhỏ (10 - 18 cm) và giảm dần ở các cấp kính lớn. Hàm Mayer và hàm khoảng cách phản ánh chính xác quy luật giảm liên tục hình chữ J ngược này với mức sai số thấp nhất khi kiểm định qua tiêu chuẩn Chi-square.

Thực trạng tái sinh tự nhiên tại Vườn quốc gia Bù Gia Mập có đảm bảo khả năng phục hồi rừng không? Khả năng tái sinh tự nhiên rất triển vọng khi mật độ cây tái sinh đạt từ 3.500 đến 6.200 cây/ha. Hơn 68% cây tái sinh có nguồn gốc từ hạt và trên 65% đạt phẩm chất tốt (loại A), tạo nguồn thế hệ kế cận dồi dào để thay thế tầng cây cao khi già cỗi.

Phương pháp kiểm định Kruskal-Wallis đóng vai trò gì trong việc xử lý số liệu cấu trúc rừng? Kiểm định phi tham số Kruskal-Wallis cho phép kiểm tra tính thuần nhất về trung bình mẫu và phương sai giữa các ô tiêu chuẩn độc lập. Phương pháp này đảm bảo tính khách quan trước khi gộp các ô tiêu chuẩn thuần nhất để xây dựng phương trình mô hình hóa toán học chung cho toàn lâm phần.

Giải pháp nào mang tính đột phá để ngăn chặn tình trạng khai thác lâm sản trái phép tại Vườn quốc gia Bù Gia Mập? Giải pháp hiệu quả nhất là kết hợp giữa việc lập chốt kiểm soát cố định tại các cửa suối Đăk Huýt, Đăk Ca và duy trì mô hình giao khoán 4.700 ha rừng cho các đơn vị vũ trang (Đồn biên phòng 785, Xã đội), kết hợp hỗ trợ sinh kế nâng cao thu nhập cho 31 hộ dân vùng lõi.

Kết luận

  • Luận văn đã xác lập hệ thống dữ liệu định lượng chi tiết về cấu trúc lâm phần và động thái tái sinh tự nhiên của 724 loài thực vật tại hai trạng thái rừng phục hồi đại diện (IIB và IIIA2) ở Vườn quốc gia Bù Gia Mập.
  • Chứng minh quy luật phân bố số cây theo đường kính N-D1.3 tuân theo quy luật giảm hình chữ J ngược (mô phỏng tối ưu qua hàm Mayer và phân bố khoảng cách), trong khi phân bố cây trên mặt bằng mang đặc trưng dạng cụm tập trung (Q < 1, U < -1,96).
  • Khẳng định mối quan hệ tương quan chặt chẽ giữa chiều cao vút ngọn và đường kính thân cây H-D1.3 thông qua các mô hình hồi quy phi tuyến với hệ số xác định R² > 0,75.
  • Đề xuất đồng bộ 4 nhóm giải pháp lâm sinh và quản lý thực tiễn, bao gồm khoanh nuôi làm giàu 1.500 ha rừng IIB, tỉa thưa điều chỉnh mật độ, chốt chặn suối trọng điểm và nâng cao hiệu quả giao khoán 4.700 ha rừng gắn với sinh kế đồng bào S'tiêng.
  • Gợi mở hướng nghiên cứu tiếp theo trong giai đoạn 2008 - 2015 về việc ứng dụng công nghệ GIS và viễn thám để theo dõi diễn thế động thái rừng dài hạn; các đơn vị quản lý rừng đặc dụng nên áp dụng ngay bộ tiêu chuẩn điều tra và mô hình hóa toán học này vào thực tế quản lý lâm phần bền vững.