Tổng quan nghiên cứu

Tài nguyên rừng tự nhiên tại Việt Nam đóng vai trò then chốt trong bảo vệ môi trường sinh thái, an ninh nguồn nước và phát triển kinh tế bền vững. Theo số liệu thống kê của Bộ Nông nghiệp và Phát triển nông thôn, diện tích rừng tự nhiên toàn quốc vào năm 2006 đạt khoảng 10,4 triệu ha, trong khi tỷ lệ che phủ rừng từng suy giảm nghiêm trọng từ mức 43% năm 1943 xuống còn 28% vào năm 1993. Trước sức ép gia tăng dân số, tình trạng khai thác quá mức và biến đổi khí hậu, việc bảo tồn các hệ sinh thái rừng tự nhiên trở thành nhiệm vụ cấp bách hàng đầu.

Khu bảo tồn thiên nhiên Tây Yên Tử thuộc tỉnh Bắc Giang được thành lập chính thức theo Quyết định số 117/QĐ-UB ngày 22/7/2002 của Ủy ban nhân dân tỉnh Bắc Giang. Khu bảo tồn trải rộng trên địa bàn hai huyện Lục Nam và Sơn Động với tổng diện tích rừng tự nhiên là 16.462 ha, đạt độ che phủ lên tới 73%. Trong đó, phân khu bảo vệ nghiêm ngặt chiếm 6.716 ha, phân khu phục hồi sinh thái chiếm 9.724 ha và phân khu dịch vụ du lịch sinh thái là 22 ha. Khu vực này sở hữu tính đa dạng sinh học vượt trội với 728 loài thực vật bậc cao thuộc 189 chi, 86 họ cùng 226 loài động vật có xương sống.

Tuy nhiên, các công trình nghiên cứu trước đây tại khu vực Tây Yên Tử còn phân tán và thiếu tính định lượng hệ thống về quy luật cấu trúc lâm phần. Vấn đề nghiên cứu trọng tâm của luận văn là làm sáng tỏ các quy luật cấu trúc cơ bản của tầng cây cao và tầng cây tái sinh. Mục tiêu cụ thể là xác định tổ thành loài, quy luật phân bố không gian, mô hình hóa toán học cấu trúc đường kính và chiều cao ở các trạng thái rừng chính, từ đó cung cấp luận cứ khoa học chuẩn xác để xây dựng các giải pháp kỹ thuật quản lý 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

Nghiên cứu được xây dựng trên nền tảng học thuyết về kiểu lâm phần của G. Morozov (1912) và lý thuyết phân loại rừng ba cấp của P. Pogrebnyak, trong đó xem xét rừng là một thể thống nhất hữu cơ giữa quần xã thực vật và hoàn cảnh lập địa. Đồng thời, nghiên cứu vận dụng hệ thống phân loại trạng thái rừng của Loetschau (1966) đã được Viện Điều tra Quy hoạch rừng cải tiến và chuẩn hóa theo quy phạm QPN 6-84. Khung lý thuyết này phân định rõ các giai đoạn diễn thế phục hồi của rừng tự nhiên dựa trên trữ lượng gỗ và tổng tiết diện ngang.

Bên cạnh đó, đề tài ứng dụng lý thuyết phân tích cấu trúc tổ thành loài thông qua chỉ số giá trị quan trọng IV% của Daniel Marmillod và quan điểm ưu thế sinh thái của Thái Văn Trừng (1978). Các khái niệm cốt lõi được vận dụng xuyên suốt bao gồm: cấu trúc hình thái không gian, động thái tái sinh tự nhiên, và khung tiêu chuẩn quản lý rừng bền vững của Tổ chức Gỗ Nhiệt đới Quốc tế (ITTO) cùng Bộ tiêu chí và chỉ số quản lý rừng vùng ASEAN gồm 7 tiêu chí nền tảng.

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

Nguồn dữ liệu sơ cấp được thu thập thông qua hệ thống 9 ô tiêu chuẩn định hình với diện tích 5.000 m² mỗi ô (tương đương 0,5 ha/ô), đạt tổng diện tích điều tra trực tiếp là 45.000 m². Phương pháp chọn mẫu điển hình đại diện được áp dụng nghiêm ngặt trên 3 trạng thái rừng chủ yếu gồm rừng non phục hồi (trạng thái IIB), rừng trung bình phục hồi sau khai thác chọn (trạng thái IIIA2) và rừng giàu trên núi cao (trạng thái IIIA3). Lý do lựa chọn phương pháp chọn mẫu điển hình là nhằm phản ánh chính xác các đặc trưng lâm học đặc thù của từng đai cao và địa hình phân cắt phức tạp tại vùng núi Yên Tử.

Trên mỗi ô tiêu chuẩn, tiến hành đo đếm toàn bộ các cá thể cây gỗ có đường kính ngang ngực từ 6 cm trở lên bằng thước kẹp kính chuyên dụng, xác định chiều cao vút ngọn bằng thước đo cao Blumeleiss và đo đường kính tán cây theo hai chiều vuông góc. Tầng cây tái sinh được khảo sát chi tiết trên 2 dải vuông góc rộng 2,5 m với tổng diện tích 368 m² mỗi ô tiêu chuẩn, phân chia thành 5 cấp chiều cao sinh trưởng từ dưới 1,0 m đến trên 4,0 m.

Về phương pháp phân tích, đề tài sử dụng phần mềm thống kê sinh học chuyên dụng kết hợp với phần mềm SPSS để mô phỏng phân bố số cây theo cỡ đường kính và chiều cao thông qua hàm mật độ xác suất Weibull và hàm phân bố Khoảng cách. Việc lựa chọn hai hàm toán học này xuất phát từ khả năng mô tả linh hoạt các dạng phân bố lệch trái và giảm liên tục của rừng nhiệt đới hỗn loài khác tuổi. Mức độ phù hợp giữa phân bố lý thuyết và phân bố thực nghiệm được kiểm định bằng tiêu chuẩn Chi bình phương của Pearson ở mức ý nghĩa 5%. Cấu trúc phân bố không gian cây rừng trên mặt đất được lượng hóa chính xác qua chỉ số khoảng cách trung bình đến cây cận kề và tiêu chuẩn U của Clark - Evans. Toàn bộ timeline thu thập và xử lý dữ liệu được thực hiện đồng bộ trong giai đoạn hoàn thành công trình năm 2010.

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

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

Thứ nhất, kết quả phân loại trạng thái rừng phản ánh rõ nét sự biến động về mật độ và sinh khối giữa các đai lập địa. Tại trạng thái rừng non phục hồi IIB (phân bố ở độ cao dưới 200m), mật độ lâm phần dao động từ 684 đến 940 cây/ha với tổng tiết diện ngang từ 12,190 đến 13,911 m²/ha. Trạng thái rừng phục hồi sau khai thác IIIA2 (độ cao 200 - 500m) có mật độ từ 666 đến 758 cây/ha, tổng tiết diện ngang đạt 15,100 đến 15,626 m²/ha. Tại trạng thái rừng giàu IIIA3 trên núi cao hiểm trở (độ cao trên 500m, độ dốc trên 35 độ), mật độ cây biến động từ 666 đến 846 cây/ha nhưng tổng tiết diện ngang đạt mức cao nhất, dao động từ 16,037 đến 19,682 m²/ha.

Thứ hai, tổ thành loài cây gỗ biến đổi sâu sắc theo mức độ phục hồi của hệ sinh thái. Tầng cây cao ghi nhận sự hiện diện phong phú từ 27 đến 48 loài thực vật trong mỗi ô tiêu chuẩn. Ở trạng thái IIB và IIIA2, loài Dẻ đỏ chiếm ưu thế sinh thái vượt trội với tỷ lệ số cây N% đạt từ 7,45% đến 10,90% và chỉ số giá trị quan trọng IV% đạt từ 7,18% đến 12,94%. Ngược lại, tại trạng thái IIIA3, cấu trúc ưu thế chuyển dịch rõ rệt sang các loài cây gỗ quý vùng núi cao như Bứa (IV% đạt 10,06%), Nhung nhựa (IV% đạt 9,96%), Trâm tía (IV% đạt 8,95%) và Thị rừng (IV% đạt 7,40%).

Thứ ba, nghiên cứu đã định lượng hóa thành công hình thái phân bố không gian của các cá thể cây rừng. Kiểm định tiêu chuẩn U của Clark - Evans trên 9 ô tiêu chuẩn cho thấy có 5 ô biểu hiện dạng phân bố ngẫu nhiên (trị số U nằm trong khoảng từ -0,419 đến 1,582, tập trung chủ yếu ở trạng thái IIIA3) và 4 ô biểu hiện dạng phân bố cụm rõ nét (trị số U nhỏ hơn -1,96, điển hình ở trạng thái IIB với U đạt -2,089).

Thứ tư, quy luật cấu trúc đường kính ngang ngực và chiều cao vút ngọn đều tuân theo các quy luật toán học chặt chẽ. Cấu trúc số cây theo cỡ đường kính biểu hiện đường cong giảm liên tục hoặc dạng một đỉnh lệch, mô phỏng xuất sắc qua hàm Weibull và hàm Khoảng cách với giá trị Chi bình phương tính toán luôn nhỏ hơn giá trị Chi bình phương tra bảng ở độ tin cậy 95%.

Thảo luận kết quả

Sự phân hóa rõ rệt về tổ thành và cấu trúc tầng tán bắt nguồn từ sự tương tác giữa điều kiện lập địa khắc nghiệt và lịch sử tác động của con người. Tại các đai cao dưới 200m, hoạt động canh tác nương rẫy và khai thác trong quá khứ đã làm suy giảm các loài cây mục đích lớn, tạo điều kiện cho các loài cây tiên phong ưa sáng như Vàng trứng, Phân mã chiếm tỷ lệ cao. Trong khi đó, ở đai cao trên 500m với địa hình dốc đứng trên 35 độ và nhiều vách đá vôi hiểm trở, rừng ít bị con người tác động nên vẫn bảo tồn được các quần xã thứ sinh ổn định với sự xuất hiện của các loài quý hiếm như Vàng kiếng, Nanh vàng, Chò chỉ.

Về mặt biểu diễn dữ liệu, toàn bộ cấu trúc tổ thành loài được trực quan hóa thông qua biểu đồ đường kết hợp giữa hệ số N% và IV%, cho phép các nhà lâm học nhận diện trực tiếp sự chênh lệch giữa số lượng cá thể và mức độ chiếm lĩnh không gian dinh dưỡng của từng loài cây. Đồng thời, cấu trúc đường kính và chiều cao được mô tả sinh động qua các biểu đồ phân bố thực nghiệm dạng cột kết hợp đường cong mô phỏng lý thuyết Weibull, thể hiện sự tập trung mật độ cao ở cấp kính nhỏ từ 6 cm đến 18 cm và giảm dần ở các cấp kính lớn trên 40 cm. Kết quả chuyển dịch phân bố không gian từ dạng cụm sang ngẫu nhiên chứng minh rừng Tây Yên Tử đang trải qua quá trình tự tỉa thưa tự nhiên mạnh mẽ để mở rộng không gian sống.

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

Dựa trên kết quả nghiên cứu cấu trúc lâm phần, 4 nhóm giải pháp kỹ thuật lâm sinh và quản trị cần được triển khai đồng bộ:

  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 cây bản địa: Tiến hành xử lý thực bì cục bộ, chặt bỏ dây leo và cây sâu bệnh, trồng dặm bổ sung từ 300 đến 500 cây bản địa quý hiếm mỗi hecta (như Lim xanh, Pơ mu, Gụ, Dinh, Sến mật) trên 200 ha rừng thuộc trạng thái IIB và IIIA2 trong giai đoạn 2026 - 2028 do Ban Quản lý Khu bảo tồn thiên nhiên Tây Yên Tử chủ trì.
  2. Hoàn thiện hệ thống ranh giới và cắm mốc bảo vệ nghiêm ngặt: Tiến hành số hóa tọa độ GIS và cắm 150 cột mốc ranh giới kiên cố bao quanh 6.716 ha phân khu bảo vệ nghiêm ngặt, mục tiêu ngăn chặn 100% tình trạng chặt phá rừng và săn bắt động vật hoang dã trái phép trong lộ trình 2026 - 2027 do Hạt Kiểm lâm Sơn Động phối hợp chính quyền địa phương thực hiện.
  3. Mở rộng mô hình đồng quản lý và giao khoán bảo vệ rừng bền vững: Ký kết hợp đồng giao khoán quản lý bảo vệ rừng đối với 9.724 ha phân khu phục hồi sinh thái cho 5 tổ đội chuyên trách và các cộng đồng thôn bản người Tày, Nùng, Sán Dìu, bảo đảm mục tiêu nâng cao 15% đến 20% thu nhập từ lâm sản ngoài gỗ và dịch vụ môi trường rừng trong giai đoạn 2026 - 2030 dưới sự điều phối của Ủy ban nhân dân huyện Sơn Động.
  4. Xây dựng mạng lưới định vị ô tiêu chuẩn giám sát động thái đa dạng sinh học: Thiết lập cố định 15 ô tiêu chuẩn giám sát định kỳ với chu kỳ đo đạc 2 năm một lần nhằm theo dõi tốc độ tăng trưởng đường kính, sinh khối và động thái tầng cây tái sinh mục đích từ năm 2026 đến năm 2030 do các viện nghiên cứu lâm nghiệp phối hợp lực lượng kỹ thuật tại chỗ thực hiện.

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

  1. Ban quản lý các khu rừng đặc dụng và lực lượng kiểm lâm: Cung cấp cơ sở dữ liệu định lượng về chỉ số IV%, mật độ lâm phần và cấu trúc tầng tán để xây dựng phương án quản lý rừng bền vững và thiết kế các tuyến tuần tra bảo vệ cho 16.462 ha rừng tự nhiên.
  2. Các nhà nghiên cứu, nghiên cứu sinh và học viên cao học ngành Lâm học và Sinh thái học: Kế thừa quy trình ứng dụng hàm Weibull, hàm Khoảng cách và tiêu chuẩn Clark - Evans trong mô hình hóa cấu trúc quần xã thực vật nhiệt đới vùng Đông Bắc.
  3. Cơ quan quản lý nhà nước và tổ chức quy hoạch lâm nghiệp: Sử dụng các căn cứ khoa học từ luận văn để phê duyệt quy hoạch bảo tồn đa dạng sinh học cấp tỉnh, phân bổ ngân sách chi trả dịch vụ môi trường rừng và định hình chính sách phát triển lâm nghiệp xã hội tại tỉnh Bắc Giang.
  4. Giảng viên và sinh viên các trường đại học khối Nông - Lâm - Môi trường: Sử dụng làm tài liệu tham khảo giảng dạy chuyên đề điều tra quy hoạch rừng, phương pháp thiết lập ô tiêu chuẩn 5.000 m² và kỹ thuật phân tích số liệu lâm sinh học bằng phần mềm thống kê chuyên ngành.

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

  1. Hàm phân bố Weibull có ưu thế gì trong mô tả cấu trúc đường kính cây rừng tại Tây Yên Tử? Hàm Weibull thể hiện tính linh hoạt vượt trội nhờ khả năng mô tả chính xác cả dạng phân bố giảm liên tục lẫn dạng phân bố một đỉnh lệch trái. Kết quả kiểm định Chi bình phương với độ tin cậy 95% khẳng định hàm Weibull phản ánh chân thực kết cấu phân tầng và mật độ tập trung ở các cấp đường kính từ 6 cm đến 20 cm của lâm phần tự nhiên.

  2. Chỉ số giá trị quan trọng IV% phản ánh vai trò sinh thái của loài cây khác biệt như thế nào so với tỷ lệ số cây N%? Tỷ lệ số cây N% đơn thuần phản ánh số lượng cá thể, trong khi chỉ số IV% là tổng hợp giữa tỷ lệ cá thể và tỷ lệ tiết diện ngang. Nhờ đó, chỉ số IV% giúp xác định chuẩn xác các loài cây gỗ lớn chiếm ưu thế sinh thái thực sự như Dẻ đỏ (IV% đạt 12,94%) hay Bứa (IV% đạt 10,06%), loại trừ sai lệch từ các loài cây bụi nhỏ có mật độ cá thể cao.

  3. Nguyên nhân khiến cây rừng ở trạng thái IIB phân bố theo dạng cụm trên mặt đất là gì? Ở trạng thái rừng non phục hồi IIB, các loài cây tiên phong ưa sáng thường tái sinh tập trung tại các khoảng trống tán rừng hoặc ven nương rẫy cũ. Sự cạnh tranh không gian dinh dưỡng ở giai đoạn đầu chưa gay gắt, kết hợp đặc tính phát tán hạt cục bộ dẫn tới trị số tiêu chuẩn U nhỏ hơn -1,96, thể hiện cấu trúc phân bố cụm điển hình.

  4. Tầng cây cao tại trạng thái rừng IIIA3 ở độ cao trên 500m có đặc điểm cấu trúc nổi bật nào? Tại đai cao trên 500m với độ dốc trên 35 độ, rừng IIIA3 sở hữu tổng tiết diện ngang cao nhất đạt 19,682 m²/ha. Cấu trúc tổ thành xuất hiện nhiều loài cây gỗ quý chịu hạn trên núi đá như Nhung nhựa, Trâm tía, Bứa, Chò chỉ với hình thái phân bố ngẫu nhiên chiếm ưu thế do quá trình chọn lọc tự nhiên lâu dài.

  5. Việc điều tra cây tái sinh trên diện tích 368 m² mỗi ô tiêu chuẩn mang lại ý nghĩa gì cho công tác lâm sinh? Phương pháp điều tra cây tái sinh trên hai dải trực giao phân theo 5 cấp chiều cao giúp đánh giá chính xác mật độ và khả năng thay thế thế hệ cây già cỗi. Tỷ lệ cây tái sinh mục đích phát triển tốt ở cấp chiều cao trên 1,0 m là chỉ dấu sinh học then chốt để quyết định áp dụng biện pháp khoanh nuôi đơn thuần hay xúc tiến tái sinh có trồng bổ sung.

Kết luận

  • Luận văn đã xác định toàn diện các đặc trưng cấu trúc lâm phần của 16.462 ha rừng tự nhiên tại Khu bảo tồn thiên nhiên Tây Yên Tử qua 3 trạng thái diễn thế chính gồm IIB, IIIA2 và IIIA3.
  • Định lượng hóa chính xác cấu trúc tổ thành loài với 27 đến 48 loài thực vật mỗi ô tiêu chuẩn, chứng minh vai trò ưu thế sinh thái của Dẻ đỏ ở đai thấp và các loài Bứa, Nhung nhựa, Trâm tía ở đai cao.
  • Ứng dụng thành công hàm phân bố Weibull và hàm Khoảng cách trong mô phỏng quy luật cấu trúc đường kính và chiều cao với độ tin cậy thống kê trên 95%.
  • Làm rõ quy luật chuyển dịch không gian từ phân bố cụm ở rừng non sang phân bố ngẫu nhiên ở rừng già thông qua kiểm định tiêu chuẩn U của Clark - Evans.
  • Đề xuất hệ thống 4 giải pháp lâm sinh và quản trị đồng bộ nhằm bảo tồn nguồn gen quý hiếm của 728 loài thực vật và 226 loài động vật đặc hữu.

Trong giai đoạn 2026 - 2030, Ban Quản lý Khu bảo tồn và các cơ quan lâm nghiệp cần khẩn trương số hóa toàn bộ hệ thống ô tiêu chuẩn định vị và triển khai chương trình làm giàu 200 ha rừng phục hồi sinh thái. Các nhà quản lý và nghiên cứu lâm nghiệp hãy ứng dụng ngay khung phân tích cấu trúc định lượng này vào thực tiễn quy hoạch để bảo tồn trọn vẹn giá trị đa dạng sinh học và phát triển bền vững tài nguyên rừng tự nhiên.