Tổng quan nghiên cứu

Thực trạng tài nguyên rừng Việt Nam qua hơn 60 năm lịch sử đã chứng kiến sự suy giảm nghiêm trọng về diện tích và chất lượng do sức ép gia tăng dân số và phương thức canh tác nương rẫy. Nếu như trước năm 1945 độ che phủ rừng cả nước đạt khoảng 43% diện tích đất tự nhiên, thì đến năm 2005 diện tích rừng chỉ còn khoảng 12.700 nghìn ha, độ che phủ giảm xuống mức 37% theo công bố tại Quyết định số 1970/QĐ/BNN-KL của Bộ Nông nghiệp và Phát triển nông thôn. Đặc biệt, tại một số khu vực Tây Bắc như Sơn La và Lai Châu, độ che phủ có thời điểm chạm đáy ở mức 8% đến 9%. Hiện tượng phân mảnh sinh cảnh và khai thác quá mức đã làm phá vỡ cấu trúc lâm phần, đe dọa trực tiếp đến tính đa dạng sinh học và chức năng phòng hộ đầu nguồn.

Trước yêu cầu cấp bách về bảo tồn nguồn gen quý hiếm và quản lý bền vững các hệ sinh thái vùng núi cao, luận văn thạc sĩ chuyên ngành Lâm học của tác giả Nguyễn Duy Thịnh, dưới sự hướng dẫn khoa học của Tiến sĩ Phạm Ngọc Giao tại Trường Đại học Lâm nghiệp (năm 2007), đã triển khai đề tài nghiên cứu đặc điểm cấu trúc lâm phần tại Khu bảo tồn thiên nhiên Hoàng Liên - Văn Bàn, tỉnh Lào Cai. Khu bảo tồn được thành lập theo Quyết định số 702/QĐ-UBND ngày 27/3/2007 của Ủy ban nhân dân tỉnh Lào Cai, sở hữu diện tích tự nhiên rộng lớn với 23.476,89 ha đất có rừng, chiếm 91,5% tổng diện tích khu bảo tồn, trong đó rừng tự nhiên đạt 23.321,99 ha (chiếm 90,9%).

Mục tiêu cụ thể của luận văn là định lượng hóa các quy luật cấu trúc tầng cây gỗ và động thái tái sinh tự nhiên trên phạm vi 21.629,40 ha thuộc Phân khu Bảo tồn nguyên vẹn. Nghiên cứu tập trung tại các tiểu khu trọng điểm như 447, 462, 500, 506, 510, 518, 528 (xã Nậm Xé) và 522, 526, 527, 528, 532, 533, 534, 535, 537 (xã Nậm Xây). Kết quả nghiên cứu cung cấp luận cứ khoa học chuẩn xác để bảo vệ quần thể 891 loài thực vật bậc cao có mạch cùng 60 loài thú, 310 loài chim bản địa, đồng thời hỗ trợ Ban quản lý xây dựng phương án điều chế rừng tối ưu.

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 sinh thái học rừng mưa nhiệt đới của Odum (1971) và nguyên lý quản lý rừng bền vững của Baur (1964). Luận văn kế thừa quan điểm hệ thống của Bertram Husch và Đồng Sĩ Hiền (1974), định nghĩa lâm phần như một đơn vị sinh vật học hoàn chỉnh mà ở đó các cây gỗ cùng sinh trưởng, tạo lập tàn che và vận động theo những quy luật kết cấu xác định về không gian thẳng đứng và mặt phẳng nằm ngang.

Bên cạnh đó, đề tài vận dụng hệ thống phân loại thảm thực vật nhiệt đới gió mùa của Thái Văn Trừng (1978) kết hợp với bảng phân chia trạng thái rừng theo Loeschau cải tiến của Viện Điều tra Quy hoạch Rừng. Các khái niệm học thuật cốt lõi được lượng hóa gồm có:

  1. Tiết diện ngang lâm phần (G, tính bằng mét vuông trên hecta).
  2. Trữ lượng gỗ lâm phần (M, tính bằng mét khối trên hecta).
  3. Cấu trúc tổ thành loài theo công thức tiết diện ngang.
  4. Quy luật phân bố mật độ cá thể theo đường kính ngang ngực (N-D1,3) và chiều cao vút ngọn (N-Hvn).
  5. Phân bố không gian tái sinh theo chỉ số khoảng cách Clark & Evans (1954) kết hợp mô hình xác suất Poisson.

Khung lý thuyết này giúp chuyển đổi các mô tả định tính truyền thống thành hệ thống mô hình toán học giải tích chặt chẽ, phù hợp với điều kiện địa hình chia cắt mạnh trên đai cao từ 600 m đến 2.913 m của đỉnh Lang Cung.

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 ô tiêu chuẩn điển hình tạm thời diện tích 1.000 m2 (kích thước 25 m x 40 m) được thiết lập song song với đường đồng mức. Phương pháp chọn mẫu điển hình phân tầng được áp dụng trên 4 trạng thái rừng chủ yếu (rừng giàu IIIa3, rừng trung bình IIIa2, rừng nghèo IIIa1 và rừng phục hồi IIb) phân bố tại 3 cấp đai cao: dưới 700 m, từ 700 m đến 1.700 m và trên 1.700 m. Lý do lựa chọn phương pháp này là nhằm phản ánh toàn diện tính phân hóa của thảm thực vật theo độ cao và kiểm soát sai số do độ dốc bình quân từ 30 đến 35 độ.

Trên mỗi ô tiêu chuẩn 1.000 m2, tác giả tiến hành đo đếm toàn diện tất cả các cây gỗ có đường kính D1,3 từ 6 cm trở lên bằng thước kẹp kính theo hai chiều vuông góc, đo chiều cao vút ngọn (Hvn) và chiều cao dưới cành (Hdc) bằng thước Blume-Leiss, đo đường kính tán (Dt) bằng thước dây và vẽ phẫu đồ rừng trên diện tích 400 m2. Để điều tra lớp cây tái sinh, mỗi ô tiêu chuẩn bố trí 5 ô dạng bản diện tích 25 m2 (5 m x 5 m) tại 4 góc và trung tâm (tổng diện tích 125 m2 mỗi ô tiêu chuẩn), phân loại cây tái sinh theo 2 cấp chiều cao (dưới 1 m và trên 1 m) cùng 3 cấp chất lượng (tốt, trung bình, xấu), đồng thời đo 8 khoảng cách đến cây lân cận gần nhất trong từng ô dạng bản.

Dữ liệu ngoại nghiệp thu thập trong năm 2007 được xử lý loại bỏ sai số bất thường bằng phương pháp Rumski. Tính thuần nhất giữa các mẫu điều tra được kiểm định bằng tiêu chuẩn chuẩn tắc U (đối với cỡ mẫu n > 30) và tiêu chuẩn phi tham số Kruskal-Wallis. Các hàm phân bố toán học gồm hàm giảm Meyer, phân bố khoảng cách và phân bố Weibull được sử dụng để nắn khớp thực nghiệm, sau đó kiểm tra độ phù hợp thông qua tiêu chuẩn Khi bình phương (Chi-square) ở các mức ý nghĩa alpha = 0,05 và 0,01. Mối tương quan Hvn-D1,3 và Dt-D1,3 được mô hình hóa qua các phương trình hồi quy phi tuyến và tuyến tính một lớp.

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ơ cấu diện tích và hiện trạng thảm thực vật, phân khu bảo tồn nguyên vẹn có tổng diện tích 21.629,40 ha (chiếm 84,3% diện tích khu bảo tồn). Trong đó, diện tích đất có rừng tự nhiên đạt 19.943,20 ha (chiếm 92,2% diện tích phân khu). Trạng thái rừng giàu IIIa3 chiếm ưu thế áp đảo với 15.520,10 ha (tương đương 71,7% diện tích phân khu bảo tồn nguyên vẹn), rừng trung bình IIIa2 chiếm 3.523,40 ha (16,3%), rừng nghèo IIIa1 chiếm 264,70 ha (1,2%), rừng phục hồi IIb đạt 410,30 ha (1,9%) và rừng hỗn giao gỗ vầu chiếm 224,70 ha (1,0%). Thảm thực vật được phân chia thành 4 kiểu chính, điển hình là rừng kín thường xanh mưa ẩm á nhiệt đới núi thấp với tổ thành ưu thế gồm các loài Dẻ gai, Sồi hồng, Kháo đá, Vối thuốc và Mạ sa.

Thứ hai, về quy luật phân bố số cây theo cỡ đường kính (N-D1,3), kết quả kiểm định cho thấy các lâm phần tự nhiên đều mang đặc trưng của rừng hỗn giao khác tuổi. Đường cong phân bố thực nghiệm có dạng giảm liên tục hoặc 1 đỉnh lệch trái rõ nét. Tỷ lệ số cây ở các cấp kính nhỏ từ 6 cm đến 14 cm chiếm tới 65% đến 72% tổng số cây điều tra. Hàm phân bố Meyer và hàm Weibull đều mô phỏng xuất sắc quy luật này với giá trị Khi bình phương tính toán luôn nhỏ hơn giá trị giới hạn tra bảng ở mức ý nghĩa alpha = 0,05, chứng minh quá trình đào thải tự nhiên diễn ra liên tục và ổn định.

Thứ ba, về tương quan sinh trưởng lâm học, mối quan hệ giữa chiều cao vút ngọn (Hvn) và đường kính ngang ngực (D1,3) biểu hiện dưới dạng hàm logarit và đường cong parabol với hệ số tương quan đạt từ 0,75 đến 0,88. Mối quan hệ giữa đường kính tán (Dt) và D1,3 phát triển theo dạng tuyến tính bậc nhất với hệ số tương quan R > 0,70. Tầng tích tụ tán cây tập trung dày đặc nhất ở khoảng độ cao từ 10 m đến 18 m, các tán cây cao trên 20 m thưa dần và rất hiếm gặp tán cây vượt trên 30 m.

Thứ tư, về động thái tái sinh tự nhiên, mật độ cây tái sinh trung bình đạt từ 4.500 đến 8.200 cây/ha. Cơ cấu chất lượng cây tái sinh rất khả quan với tỷ lệ cây tốt và trung bình chiếm trên 78%, trong đó cây có nguồn gốc hạt chiếm hơn 85%. Chỉ số Clark & Evans (Q) dao động trong khoảng 0,65 đến 0,82 (giá trị Q < 1), kết hợp với mô hình xác suất Poisson chứng minh cây tái sinh phân bố theo kiểu cụm tập trung quanh các lỗ trống tán rừng và vệt trượt tự nhiên.

Thảo luận kết quả

Điều kiện khí hậu cận nhiệt đới núi cao với lượng mưa bình quân lên tới 3.500 mm/năm, độ ẩm tương đối trung bình 87% và 160,8 ngày sương mù trong năm đã tạo lập môi trường ẩm mát lý tưởng, thúc đẩy quá trình tái sinh và phân tầng liên tục của thảm thực vật. Tuy nhiên, địa hình dốc lớn từ 30 đến 35 độ kết hợp hoạt động kinh tế vùng đệm của 391 hộ dân (2.884 nhân khẩu với mật độ 8 người/km2), đặc biệt là tập quán mở rộng diện tích trồng cây thảo quả dưới tán rừng, đang làm suy thoái nghiêm trọng tầng cây tái sinh và chia cắt tầng tán.

Khi so sánh với các nghiên cứu kinh điển của Thái Văn Trừng (1978) tại miền Bắc và Bảo Huy (1993) tại Tây Nguyên, mô hình cấu trúc rừng tại Hoàng Liên - Văn Bàn phản ánh tính nguyên sinh cao hơn nhờ tỷ lệ trạng thái IIIa3 vượt trội (71,7%). Hàm Weibull thể hiện tính mềm dẻo vượt bậc so với các hàm mũ đơn giản khi mô tả các lâm phần có cấu trúc phân tầng phức tạp từ 3 đến 5 tầng tán.

Trong các báo cáo quản lý, toàn bộ dữ liệu cấu trúc này có thể được trực quan hóa thông qua phẫu đồ đứng và phẫu đồ ngang tỷ lệ 1:100 trên diện tích 400 m2, kết hợp biểu đồ phân bố tần số N-D1,3 dạng cột chồng và bảng ma trận tương quan đa biến giữa các chỉ tiêu đường kính, chiều cao và diện tích tán.

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

  1. Thiết lập hệ thống ranh giới bảo vệ nghiêm ngặt và số hóa dữ liệu lâm phần: Ban quản lý Khu bảo tồn thiên nhiên Hoàng Liên - Văn Bàn chủ trì phối hợp cùng Hạt Kiểm lâm huyện Văn Bàn triển khai rà soát, cắm mốc ranh giới và số hóa bản đồ chuyên đề cho toàn bộ 21.629,40 ha phân khu bảo tồn nguyên vẹn. Mục tiêu đạt tỷ lệ bảo vệ nguyên vẹn 100% diện tích rừng tự nhiên, chấm dứt hoàn toàn tình trạng khai thác trái phép các loài gỗ quý như Pơ mu, Bách tán Đài Loan và Giổi găng trong giai đoạn 2008-2010.

  2. Quy hoạch và kiểm soát triệt để diện tích canh tác thảo quả dưới tán rừng: Ủy ban nhân dân huyện Văn Bàn cùng Ủy ban nhân dân hai xã Nậm Xé và Nậm Xây phối hợp với Ban quản lý khu bảo tồn tiến hành đo đạc, lập hồ sơ quản lý các nương thảo quả hiện có. Nghiêm cấm hành vi phát dọn tầng cây tái sinh để trồng mới thảo quả, duy trì tỷ lệ cây non chất lượng tốt đạt trên 80% trên diện tích 19.943,20 ha đất có rừng, tiến hành kiểm tra định kỳ 6 tháng một lần từ năm 2008 đến 2012.

  3. Áp dụng các giải pháp kỹ thuật lâm sinh nuôi dưỡng và phục hồi rừng: Phòng Kỹ thuật thuộc Ban quản lý khu bảo tồn triển khai các biện pháp xúc tiến tái sinh tự nhiên kết hợp làm giàu rừng đối với 264,70 ha rừng nghèo IIIa1 và 410,30 ha rừng phục hồi IIb. Mục tiêu thúc đẩy tăng trưởng trữ lượng bình quân từ 3% đến 5%/năm, tạo điều kiện chuyển hóa các trạng thái rừng thứ sinh sang trạng thái rừng trung bình và giàu trong lộ trình 5 năm (2008-2013).

  4. Xây dựng mô hình đồng quản lý và hỗ trợ chuyển đổi sinh kế vùng đệm: Sở Nông nghiệp và Phát triển nông thôn tỉnh Lào Cai phối hợp cùng các tổ chức bảo tồn triển khai chương trình giao khoán bảo vệ rừng gắn với chi trả dịch vụ môi trường rừng cho 391 hộ dân (chủ yếu là đồng bào dân tộc Dao và Mông). Hỗ trợ phát triển các mô hình nông lâm kết hợp ngoài vùng lõi, giảm thiểu 90% áp lực xâm hại tài nguyên rừng giai đoạn 2008-2015.

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

  1. Ban quản lý các khu bảo tồn thiên nhiên và vườn quốc gia: Tài liệu cung cấp cơ sở dữ liệu định lượng chi tiết về 4 trạng thái rừng và hệ phương trình hồi quy tương quan sinh trưởng. Cán bộ quản lý có thể ứng dụng trực tiếp để xây dựng phương án quản lý rừng bền vững, quy chế tuần tra rừng đặc dụng và định vị các tiểu khu trọng điểm cần ưu tiên bảo vệ nghiêm ngặt.

  2. Giảng viên và nhà nghiên cứu chuyên ngành Lâm học, Sinh thái học: Luận văn là tài liệu tham khảo giá trị về phương pháp ứng dụng các mô hình toán học giải tích (Meyer, Weibull, Poisson, tiêu chuẩn Clark & Evans) trong mô phỏng cấu trúc lâm phần nhiệt đới. Giảng viên có thể sử dụng số liệu thực chứng để minh họa cho học phần Điều tra rừng và Lâm sinh học.

  3. Học viên cao học và sinh viên ngành Quản lý tài nguyên rừng: Cung cấp quy trình chuẩn mực từ thiết kế hệ thống ô tiêu chuẩn 1.000 m2, kỹ thuật lập 5 ô dạng bản 25 m2 điều tra tái sinh, đến phương pháp xử lý số liệu bằng phần mềm tin học lâm nghiệp, hỗ trợ định hình phương pháp luận cho các đề tài nghiên cứu thực địa.

  4. Các nhà hoạch định chính sách và chính quyền địa phương: Cung cấp bằng chứng thực tiễn về mối quan hệ giữa mật độ dân số vùng đệm (8 người/km2) với mức độ tác động tài nguyên rừng. Dữ liệu giúp các nhà quản lý xây dựng chính sách hỗ trợ phát triển sinh kế, quy hoạch vùng trồng cây dược liệu và phân bổ ngân sách bảo tồn hợp lý tại tỉnh Lào Cai.

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

1. Tại sao nghiên cứu cấu trúc lâm phần là cơ sở quyết định cho công tác bảo tồn rừng đặc dụng?
Nghiên cứu cấu trúc lâm phần giúp định lượng hóa quy luật phân bố số cây theo đường kính, chiều cao và không gian dinh dưỡng. Việc nắm vững các thông số này trên diện tích 21.629,40 ha phân khu bảo vệ nghiêm ngặt cho phép nhà quản lý dự báo chính xác xu hướng sinh trưởng của tầng cây gỗ và đề ra các biện pháp điều tiết lâm sinh kịp thời, bảo đảm sự ổn định lâu dài của hệ sinh thái.

2. Dạng hàm toán học nào mô phỏng chính xác nhất phân bố đường kính (N-D1,3) tại khu vực nghiên cứu?
Hàm giảm Meyer và hàm phân bố xác suất Weibull là hai mô hình mô phỏng phù hợp nhất cho phân bố N-D1,3 của rừng tự nhiên tại Hoàng Liên - Văn Bàn. Kết quả kiểm định Khi bình phương ở mức ý nghĩa alpha = 0,05 đều cho thấy sự tương thích cao giữa tần số lý thuyết và thực nghiệm, phản ánh cấu trúc rừng hỗn giao khác tuổi có mật độ cây nhỏ chiếm ưu thế trên 65%.

3. Hoạt động canh tác cây thảo quả dưới tán rừng gây ra những nguy cơ sinh thái nào?
Việc mở rộng diện tích trồng cây thảo quả dẫn đến hành vi phát dọn lớp cây bụi và cây tái sinh dưới tán rừng tự nhiên. Tình trạng này làm suy giảm mật độ cây con triển vọng, phá vỡ cấu trúc phân tầng tự nhiên, làm đứt gãy quá trình tái sinh kế thừa của các loài cây gỗ lớn và đe dọa trực tiếp môi trường sống của quần thể động vật hoang dã.

4. Phân khu bảo tồn nguyên vẹn Hoàng Liên - Văn Bàn gồm những trạng thái rừng chủ yếu nào?
Phân khu bảo tồn nguyên vẹn có tổng diện tích đất có rừng là 19.943,20 ha. Trong đó gồm có: trạng thái rừng giàu IIIa3 chiếm 15.520,10 ha; trạng thái rừng trung bình IIIa2 chiếm 3.523,40 ha; trạng thái rừng nghèo IIIa1 chiếm 264,70 ha; trạng thái rừng phục hồi IIb chiếm 410,30 ha và rừng hỗn giao gỗ vầu chiếm 224,70 ha.

5. Kiểu phân bố không gian của tầng cây tái sinh được xác định bằng phương pháp nào?
Đề tài sử dụng mô hình phân bố xác suất Poisson kết hợp chỉ số khoảng cách Clark & Evans (Q) thông qua việc đo cự ly giữa các cây tái sinh trong 5 ô dạng bản 25 m2. Kết quả chỉ số Q đạt dưới 1 chứng minh tầng cây tái sinh phân bố theo dạng cụm tập trung, phụ thuộc chặt chẽ vào độ mở tán rừng và điều kiện lập địa cục bộ.

Kết luận

  • Luận văn đã định lượng hóa toàn diện các quy luật cấu trúc lâm phần và động thái tái sinh trên quy mô 21.629,40 ha thuộc Phân khu Bảo tồn nguyên vẹn của Khu bảo tồn thiên nhiên Hoàng Liên - Văn Bàn.
  • Xác lập thành công các mô hình toán học giải tích gồm hàm Meyer, hàm Weibull và các phương trình hồi quy tương quan sinh trưởng Hvn-D1,3 có độ tin cậy thống kê cao.
  • Đánh giá chi tiết hiện trạng thảm thực vật với 891 loài thực vật bậc cao và ghi nhận ưu thế vượt trội của trạng thái rừng giàu IIIa3 với diện tích 15.520,10 ha (chiếm 71,7% phân khu).
  • Làm sáng tỏ quy luật phân bố cụm của tầng cây tái sinh (chỉ số Clark & Evans Q < 1) và cảnh báo nguy cơ suy thoái rừng do tập quán mở rộng canh tác thảo quả dưới tán.
  • Đề xuất hệ thống 4 nhóm giải pháp kỹ thuật lâm sinh, quản lý ranh giới và chính sách sinh kế vùng đệm đồng bộ cho giai đoạn 2008-2015.

Về định hướng tiếp theo, đơn vị quản lý cần tiếp tục duy trì mạng lưới ô tiêu chuẩn định vị để quan trắc diễn biến cấu trúc lâm phần theo chu kỳ 5 năm một lần. Độc giả, các cán bộ quản lý tài nguyên và nghiên cứu sinh quan tâm có thể khai thác trực tiếp cơ sở dữ liệu của luận văn để ứng dụng hiệu quả vào công tác bảo tồn đa dạng sinh học tại các khu rừng đặc dụng trên toàn quốc.