Tổng quan nghiên cứu

Theo số liệu thống kê từ Bộ Nông nghiệp và Phát triển nông thôn, diện tích rừng tự nhiên tại Việt Nam đạt khoảng 10,4 triệu ha vào cuối năm 2006. Tuy nhiên, trước sức ép gia tăng dân số, tập quán canh tác nương rẫy và áp lực khai thác lâm sản, diện tích và chất lượng rừng trên cả nước đối mặt với nguy cơ suy giảm đáng kể. Riêng tại tỉnh Bắc Giang, diện tích rừng bị thiệt hại được ghi nhận ở mức 29,67 ha trong giai đoạn tính đến tháng 5/2009. Thực trạng này đặt ra yêu cầu cấp thiết phải làm sáng tỏ các quy luật cấu trúc nội tại của quần xã thực vật rừng nhằm thiết lập cơ sở khoa học vững chắc cho công tác quản lý và phục hồi sinh thái.

Nghiên cứu tập trung giải quyết vấn đề thiếu hụt các dẫn liệu định lượng có tính hệ thống về đặc điểm lâm học của rừng tự nhiên tại vùng Đông Bắc Việt Nam. Mục tiêu cụ thể là xác định các đặc trưng cấu trúc cơ bản của tầng cây cao (tổ thành loài, mật độ, hình thái phân bố không gian, quy luật phân bố theo đường kính $D_{1.3}$ và chiều cao $H_{vn}$), đánh giá động thái tầng cây tái sinh, từ đó đề xuất hệ thống giải pháp kỹ thuật lâm sinh quản lý bền vững.

Phạm vi nghiên cứu được thực hiện tại Khu bảo tồn thiên nhiên Tây Yên Tử, trải rộng trên địa bàn các xã An Lạc, Thanh Sơn, Thanh Luận (huyện Sơn Động) và xã Lục Sơn (huyện Lục Nam), tỉnh Bắc Giang. Khu bảo tồn có tổng diện tích tự nhiên 16.462 ha với độ che phủ đạt 73%, bao gồm phân khu bảo vệ nghiêm ngặt 6.716 ha và phân khu phục hồi sinh thái 9.724 ha. Công trình hoàn thành vào tháng 6/2010 mang ý nghĩa thực tiễn to lớn, hỗ trợ định hướng bảo tồn 728 loài thực vật, 226 loài động vật rừng quý hiếm và tối ưu hóa hiệu quả khoanh nuôi trên 200 ha rừng tái sinh.

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 xây dựng trên nền tảng học thuyết kiểu rừng của G.F. Morozov (1912), nhìn nhận rừng là một thể thống nhất hữu cơ giữa quần lạc sinh vật và hoàn cảnh lập địa. Đồng thời, đề tài kế thừa hệ thống phân loại trạng thái rừng của Loetschau (1966) được chuẩn hóa theo quy phạm kỹ thuật QPN 6-84 của Viện Điều tra Quy hoạch Rừng. Ngoài ra, lý thuyết sinh thái không gian của Clark và Evans (1954) cùng bộ tiêu chí và chỉ số quản lý rừng bền vững của Tổ chức Gỗ Nhiệt đới Quốc tế (ITTO) và ASEAN C&I được tích hợp làm khung tham chiếu.

Năm khái niệm then chốt được vận dụng xuyên suốt bao gồm: Cấu trúc tổ thành loài (xác định qua tỷ lệ số cây $N%$ và chỉ số giá trị quan trọng $IV%$), Tiết diện ngang lâm phần ($\Sigma G$ tính bằng $\text{m}^2/\text{ha}$), Quy luật phân bố số cây theo cấp đường kính ($N/D_{1.3}$) và cấp chiều cao ($N/H_{vn}$), Hình thái phân bố không gian (kiểu cụm, ngẫu nhiên, cách đều), và Động thái tái sinh tự nhiên dưới tán rừng.

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

Cỡ mẫu nghiên cứu gồm 9 ô tiêu chuẩn (ÔTC) với quy mô $5.000\text{ m}^2$/ô ($0,5\text{ ha}$/ô), tổng diện tích điều tra đạt $4,5\text{ ha}$. Tầng cây tái sinh được thu thập trên 2 dải đo vuông góc xuất phát từ tâm các cạnh ô tiêu chuẩn, chiều rộng $2,5\text{ m}$, đạt diện tích $368\text{ m}^2$/ô.

Phương pháp chọn mẫu điển hình đại diện được áp dụng nhằm bao quát các trạng thái rừng chủ yếu tại các đai độ cao từ dưới 200m đến trên 500m. Lý do lựa chọn các phương pháp phân tích chuyên sâu:

  • Sử dụng chỉ số $IV%$ của Daniel Marmillod kết hợp $N%$ nhằm phản ánh chính xác cả số lượng cá thể lẫn ưu thế sinh khối của từng loài trong lâm phần.
  • Áp dụng hàm phân bố xác suất Weibull và hàm Khoảng cách kết hợp kiểm định Chi bình phương ($\chi^2$ Pearson) nhằm mô hình hóa khách quan các quy luật kết cấu lâm phần thực nghiệm.
  • Sử dụng tiêu chuẩn $U$ của Clark và Evans nhằm lượng hóa chuẩn xác kiểu phân bố không gian của cây rừng trên mặt đất.

Nguồn dữ liệu thực địa được thu thập qua đo đếm trực tiếp đường kính ngang ngực $D_{1.3}$ bằng thước kẹp kính cho cây có $D_{1.3} \ge 6\text{ cm}$, đo chiều cao vút ngọn $H_{vn}$ bằng thước Blumeleiss và đo đường kính tán $D_t$ theo 2 chiều trực giao. Dữ liệu khí tượng thủy văn được tổng hợp từ trạm quan trắc Bắc Giang với lượng mưa trung bình $1.483,3\text{ mm}$/năm và độ ẩm bình quân 82%. Toàn bộ dữ liệu được xử lý bằng phần mềm SPSS 13.0 và Excel trong khung thời gian từ năm 2008 đến tháng 6/2010.

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

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

Kết quả điều tra thực nghiệm tại 9 ô tiêu chuẩn mang lại 3 phát hiện lâm học nổi bật:

Thứ nhất, có sự phân hóa rõ rệt về mật độ và trữ lượng giữa 3 trạng thái rừng theo phân loại Loetschau:

  • Trạng thái IIB (rừng non phục hồi đai cao dưới 200m) có mật độ dao động từ 684 đến 940 cây/ha, tổng tiết diện ngang $\Sigma G$ đạt $12,190 - 13,911\text{ m}^2/\text{ha}$.
  • Trạng thái IIIA2 (rừng phục hồi sau khai thác kiệt, đai cao 200 - 500m) ghi nhận mật độ 666 - 758 cây/ha, $\Sigma G$ đạt $15,100 - 15,626\text{ m}^2/\text{ha}$.
  • Trạng thái IIIA3 (rừng tác động nhẹ đai cao trên 500m) đạt mật độ 666 - 846 cây/ha, $\Sigma G$ đạt giá trị cao nhất từ $16,037$ đến $19,682\text{ m}^2/\text{ha}$ (tăng hơn 41% so với trạng thái IIB).

Thứ hai, cấu trúc tổ thành loài thể hiện tính đa dạng sinh học cao nhưng có sự chuyển dịch nhóm loài ưu thế:

  • Trạng thái IIB có 34 - 48 loài, loài Dẻ đỏ chiếm ưu thế tuyệt đối với $N% = 10,90%$ và $IV% = 12,94%$, cùng sự tham gia của Vạng trứng ($IV% = 9,93%$) và Phân mã ($IV% = 7,59%$).
  • Trạng thái IIIA2 gồm khoảng 44 loài, Dẻ đỏ tiếp tục ưu thế ($N% = 11,11%$, $IV% = 12,10%$), theo sau là Côm trâu ($IV% = 7,77%$) và Tô hợp đỏ ($IV% = 7,56%$).
  • Trạng thái IIIA3 xuất hiện 27 - 29 loài thân gỗ lớn, các loài ưu thế chuyển dịch sang Bứa ($IV% = 10,06%$), Nhung nhựa ($IV% = 9,96%$, $N% = 9,09%$) và Trâm tía ($IV% = 7,95%$).

Thứ ba, quy luật phân bố không gian và cấu trúc tầng tán phản ánh giai đoạn diễn thế phục hồi:

  • Tại trạng thái IIB, 75% số ô tiêu chuẩn (3/4 ô) có kiểu phân bố cụm với tiêu chuẩn $U$ âm rõ rệt ($U = -1,972$ đến $-2,089$). Ngược lại, tại trạng thái IIIA3, 100% số ô tiêu chuẩn (3/3 ô) đạt kiểu phân bố ngẫu nhiên với $U = 0,250$ đến $1,582$ (thỏa mãn điều kiện $|U| < 1,96$).
  • Cấu trúc $N/D_{1.3}$ và $N/H_{vn}$ có dạng đường cong giảm một đỉnh hoặc phân bố chữ J, trong đó hàm Weibull thể hiện sự tương thích lý thuyết xuất sắc với độ tin cậy thống kê cao.

Thảo luận kết quả

Sự chuyển dịch hình thái phân bố từ dạng cụm ở rừng non (IIB) sang dạng ngẫu nhiên ở rừng già ổn định (IIIA3) phản ánh quy luật tự điều tiết sinh thái. Trong giai đoạn đầu, cây con tái sinh tập trung quanh cây mẹ hoặc khoảng trống tán rừng. Qua thời gian, quá trình cạnh tranh không gian dinh dưỡng và ánh sáng dẫn đến sự tự tỉa thưa, phân tán cây rừng về trạng thái ngẫu nhiên. Kết quả này tương đồng với các nghiên cứu của Lê Sáu (1996) tại Tây Nguyên và Trần Cẩm Tú (1999) tại Hương Sơn.

Dữ liệu nghiên cứu được biểu diễn tối ưu thông qua hệ thống đồ thị và bảng ma trận chuyên dụng. Đường biểu diễn tương quan giữa $N%$ và $IV%$ giúp trực quan hóa vị trí ưu thế sinh thái của các loài mục đích. Các biểu đồ phân bố tần số thực nghiệm $N/D_{1.3}$ ghép cặp đường cong hàm Weibull và hàm Khoảng cách chứng minh tính chất phân tầng ổn định của rừng nhiệt đới. Việc phân tích bảng phân loại 5 cấp chiều cao tái sinh cho thấy mật độ cây con đạt chất lượng tốt chiếm tỷ lệ cao, đảm bảo năng lực kế thừa tự nhiên mà không cần can thiệp trồng mới nhân tạo trên diện rộng.

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

Dựa trên kết quả thực nghiệm, 4 nhóm giải pháp kỹ thuật lâm sinh và quản lý được đề xuất cụ thể:

Thực hiện khoanh nuôi xúc tiến tái sinh kết hợp làm giàu rừng trên diện tích 500 ha thuộc trạng thái IIB và IIIA2. Bổ sung các loài cây bản địa có giá trị kinh tế và bảo tồn cao như Lim xanh, Gụ, Sến mật, Pơ mu với mật độ 300 - 400 cây/ha. Thời gian thực hiện từ năm 2011 đến năm 2015 do Ban Quản lý Khu bảo tồn thiên nhiên Tây Yên Tử chủ trì.

Áp dụng kỹ thuật tỉa thưa nuôi dưỡng và điều chỉnh cấu trúc không gian đối với các lâm phần trạng thái IIIA2 có mật độ vượt quá 800 cây/ha. Chặt bỏ cây sâu bệnh, dây leo bao phủ và cây phi mục đích nhằm giải phóng không gian dinh dưỡng cho các cây mục đích, duy trì tiết diện ngang ổn định ở mức $16 - 20\text{ m}^2/\text{ha}$. Hoạt động thực hiện định kỳ hàng năm từ năm 2011 do Hạt Kiểm lâm huyện Sơn Động phối hợp các tổ đội chuyên trách phụ trách.

Mở rộng cơ chế giao khoán quản lý bảo vệ rừng gắn với chi trả dịch vụ môi trường rừng cho 100% hộ gia đình đồng bào dân tộc thiểu số (Tày, Nùng, Sán Dìu) tại 4 xã vùng đệm (An Lạc, Thanh Sơn, Thanh Luận, Lục Sơn), hướng tới mục tiêu quản lý bền vững trên 10.000 ha trong giai đoạn 2011 - 2013 do UBND huyện Sơn Động và UBND huyện Lục Nam phê duyệt.

Thiết lập hệ thống 10 ô tiêu chuẩn định vị cố định ($1\text{ ha}$/ô) để giám sát dài hạn động thái diễn thế và sinh trưởng đường kính chu kỳ 3 năm/lần từ năm 2011 đến năm 2020 do Trường Đại học Lâm nghiệp phối hợp với Ban Quản lý Khu bảo tồn thực hiện.

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

Luận văn là tài liệu tham khảo giá trị cho 4 nhóm đối tượng chuyên môn:

Ban Quản lý các Khu bảo tồn thiên nhiên và Vườn quốc gia: Nắm bắt phương pháp luận điều tra cấu trúc định lượng; vận dụng trực tiếp vào việc xây dựng phương án quản lý rừng bền vững và bảo tồn quần thể động thực vật quý hiếm trên diện tích trên 16.000 ha.

Cán bộ kỹ thuật kiểm lâm và Chi cục Kiểm lâm cấp tỉnh: Tiếp cận quy trình phân loại trạng thái rừng chuẩn hóa theo QPN 6-84; ứng dụng thiết kế hồ sơ kỹ thuật khoanh nuôi phục hồi rừng, tuần tra giám sát bảo vệ sinh cảnh vùng giáp ranh Bắc Giang - Quảng Ninh.

Giảng viên, nghiên cứu sinh và học viên cao học chuyên ngành Lâm học: Tham khảo khung phân tích toán học ứng dụng (hàm Weibull, phân bố khoảng cách, chỉ số $IV%$, tiêu chuẩn $U$ Clark-Evans); sử dụng làm tài liệu giảng dạy và mở rộng nghiên cứu động thái rừng tự nhiên nhiệt đới.

Các tổ chức phi chính phủ và cơ quan hoạch định chính sách lâm nghiệp: Khai thác dữ liệu thực chứng về đa dạng sinh học và sinh thái nhân văn vùng đệm; vận dụng xây dựng các dự án đồng quản lý tài nguyên rừng và phát triển sinh kế cộng đồng bền vững.

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

Tại sao việc sử dụng chỉ số giá trị quan trọng $IV%$ lại phản ánh cấu trúc tổ thành chính xác hơn tỷ lệ số cây $N%$? Tỷ lệ $N%$ thuần túy đo đếm số lượng cá thể, dễ dẫn đến đánh giá sai lệch khi lâm phần có nhiều cây con đường kính nhỏ. Chỉ số $IV%$ tích hợp cả $N%$ và tỷ lệ tiết diện ngang $G%$, giúp lượng hóa chính xác sinh khối, độ che phủ và vai trò sinh thái thực sự của từng loài cây gỗ lớn.

Hàm phân bố Weibull có ưu thế gì trong mô phỏng cấu trúc đường kính rừng tự nhiên Tây Yên Tử? Hàm Weibull có tính linh hoạt cao nhờ các tham số hình dạng và tỷ lệ, mô tả chính xác cả dạng phân bố giảm liên tục lẫn phân bố lệch một đỉnh. Kết quả kiểm định Chi bình phương tại 9 ô tiêu chuẩn đều cho giá trị tính toán nhỏ hơn ngưỡng tới hạn, chứng minh sự tương thích lý thuyết chặt chẽ.

Sự chuyển dịch từ phân bố cụm sang phân bố ngẫu nhiên giữa các trạng thái rừng có ý nghĩa sinh thái học gì? Ở trạng thái non (IIB), cây rừng phân bố cụm do tái sinh theo vệt quanh cây mẹ ($U = -2,089$). Khi rừng tiến đến trạng thái ổn định (IIIA3), quá trình cạnh tranh dinh dưỡng thúc đẩy sự tự điều tiết mật độ, đưa quần xã về phân bố ngẫu nhiên ($U = 1,582$), khẳng định hệ sinh thái đang phục hồi cân bằng.

Điều kiện lập địa tại Khu bảo tồn Tây Yên Tử có đặc điểm gì chi phối cấu trúc rừng? Địa hình vùng nghiên cứu hiểm trở với độ dốc phổ biến trên $30^\circ - 40^\circ$, đỉnh Yên Tử cao 1.068m. Thổ nhưỡng chủ yếu là đất Feralit trên sa thạch và phiến thạch sét, tạo nên tính phân tầng rõ rệt từ rừng kín thường xanh mưa ẩm nhiệt đới vùng thấp đến đai á nhiệt đới núi cao.

Những loài cây bản địa nào cần được ưu tiên bảo tồn và phát triển tại Tây Yên Tử? Nghiên cứu chỉ ra các loài bản địa ưu thế có giá trị kinh tế và phòng hộ cao cần ưu tiên gồm: Dẻ đỏ ($IV% = 12,94%$), Nhung nhựa ($IV% = 9,96%$), Bứa ($IV% = 10,06%$), Lim xanh, Pơ mu, Thông tre và Sến mật nhằm duy trì tính đa dạng sinh học đặc hữu.

Kết luận

  • Luận văn đã phân loại và mô tả chi tiết 3 trạng thái rừng tự nhiên chính (IIB, IIIA2, IIIA3) trên tổng diện tích 16.462 ha tại Tây Yên Tử.
  • Xác định chính xác các loài ưu thế sinh thái thông qua chỉ số $IV%$ với Dẻ đỏ ($12,94%$), Bứa ($10,06%$) và Nhung nhựa ($9,96%$).
  • Chứng minh quy luật chuyển dịch không gian từ phân bố cụm ở trạng thái rừng non ($U = -2,089$) sang phân bố ngẫu nhiên ở trạng thái rừng phục hồi ổn định ($U = 1,582$).
  • Khẳng định độ tin cậy vượt trội của hàm Weibull trong việc mô hình hóa phân bố đường kính $N/D_{1.3}$ và chiều cao $N/H_{vn}$.
  • Đề xuất hệ thống 4 giải pháp kỹ thuật lâm sinh đồng bộ gắn kết quản lý kỹ thuật với bảo tồn nguồn gen và phát triển sinh kế cộng đồng.

Công trình đóng góp bộ cơ sở dữ liệu định lượng quan trọng cho công tác quản trị tài nguyên rừng tại tỉnh Bắc Giang. Kế hoạch tiếp theo cần tập trung mở rộng quan trắc định vị định kỳ 3 năm một lần trong giai đoạn 2011 - 2020. Các nhà quản lý và cơ quan chuyên môn hãy áp dụng ngay các mô hình cấu trúc và quy trình kỹ thuật lâm sinh này vào thực tiễn quy hoạch bảo vệ rừng bền vững.