Tổng quan nghiên cứu

Hệ thống rừng tự nhiên Việt Nam đã trải qua nhiều giai đoạn biến động mạnh mẽ, từ độ che phủ 43% vào năm 1943 giảm xuống mức đáy 28,4% vào năm 1990 trước khi phục hồi lên mức 39,1% với 13,2 triệu ha vào năm 2009. Tại vùng cực Bắc tổ quốc, bốn huyện vùng cao Quản Bạ, Yên Minh, Đồng Văn và Mèo Vạc thuộc tỉnh Hà Giang sở hữu tổng diện tích tự nhiên 234.763,0 ha, trong đó đất lâm nghiệp chiếm tới 140.392,3 ha. Mặc dù số liệu thống kê ghi nhận diện tích đất có rừng đạt khoảng 77.302 ha (tương đương 55% diện tích tự nhiên), song phần lớn diện tích lại là rừng thứ sinh nghèo kiệt, phân bố chủ yếu trên địa hình núi đá vôi Karst và núi đất có tỷ lệ đá lẫn từ 30% đến 70%. Tình trạng khai thác lạm dụng và tập quán đốt nương làm rẫy kéo dài đã làm xáo trộn nghiêm trọng cấu trúc sinh thái, làm suy giảm khả năng giữ nước ngọt phục vụ đời sống của hơn 175.944 người dân địa phương.

Vấn đề cốt lõi đặt ra là sự thiếu hụt các luận cứ khoa học định lượng về quy luật cấu trúc lâm phần, dẫn đến việc các biện pháp can thiệp kỹ thuật trước đây đạt hiệu quả rất thấp hoặc chỉ dừng lại ở việc khoanh nuôi tự nhiên thụ động. Luận văn thạc sĩ khoa học lâm nghiệp của tác giả Tô Đức Hiện, dưới sự hướng dẫn của TS. Nguyễn Trọng Bình tại Trường Đại học Lâm nghiệp, đã giải quyết căn cơ bài toán này. Nghiên cứu tập trung giải mã các đặc điểm cấu trúc tổ thành, cấu trúc không gian tầng tán, mật độ và quy luật tái sinh tự nhiên của các trạng thái rừng tự nhiên (trạng thái IIA, IIB, IIIA2 trên núi đất lẫn đá và trạng thái IIA trên núi đá vôi) cùng các dạng rừng trồng hỗn giao và thuần loài. Ý nghĩa của công trình thể hiện ở việc thiết lập hệ thống cơ sở dữ liệu lâm sinh chuẩn xác, làm nền tảng đưa ra các giải pháp kỹ thuật thúc đẩy phục hồi thảm thực vật, nâng cao độ che phủ và bảo tồn tính đa dạng sinh học tại khu vực Công viên Địa chất Toàn cầu UNESCO Cao nguyên đá Đồng Văn.

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 mưa nhiệt đới của các học giả kinh điển trên thế giới kết hợp với các công trình nghiên cứu hệ sinh thái rừng nhiệt đới ẩm Việt Nam. Về mặt sinh thái, cấu trúc rừng phản ánh quy luật tổ hợp, sắp xếp giữa các thành phần thực vật trong không gian và thời gian dưới áp lực chọn lọc tự nhiên và tương tác môi trường sống. Hệ thống phân loại trạng thái rừng của Loestchau (năm 1961) bổ sung bởi Viện Điều tra Quy hoạch Rừng được áp dụng để phân định các trạng thái rừng thứ sinh từ giai đoạn phục hồi non (IIA, IIB) đến giai đoạn rừng giàu sau khai thác (IIIA2).

Các khái niệm chuyên ngành then chốt được vận dụng xuyên suốt gồm:

  • Chỉ số giá trị quan trọng (Important Value Index - IV% theo phương pháp Daniel Marmillod) nhằm lượng hóa mức độ ưu thế sinh thái của từng loài qua sự kết hợp giữa tỷ lệ số cây (N%) và tỷ lệ tổng tiết diện ngang (G%).
  • Cấu trúc không gian thẳng đứng và tầng phiến sinh thái mô tả sự phân tầng tán và vùng tích tụ tầng tán theo chiều cao.
  • Các mô hình hàm số toán học mô phỏng phân bố số cây theo đường kính ngang ngực (N/D1.3) và chiều cao vút ngọn (N/Hvn) gồm hàm giảm Meyer dạng số mũ, phân bố xác suất Weibull và phân bố khoảng cách.
  • Chỉ số đa dạng sinh học loài Shannon-Wiener (H) và chỉ số Simpson (D) nhằm đánh giá độ phong phú và tính đồng đều của quần xã thực vật.

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

Nguồn dữ liệu thực chứng được thu thập thông qua phương pháp điều tra điển hình tạm thời trên hệ thống 47 ô tiêu chuẩn (OTC). Cỡ mẫu được phân bổ gồm 33 OTC rừng tự nhiên (diện tích tiêu chuẩn 1.000 m2 kích thước 25m x 40m, hoặc 500 m2 tại các vách núi đá có độ dốc quá lớn) và 14 OTC rừng trồng (diện tích 500 m2 kích thước 20m x 25m). Các ô được phân tầng theo hai đai độ cao rõ rệt là đai dưới 700m và đai từ 700m đến 1.000m so với mực nước biển, đồng thời định vị tọa độ tuyệt đối bằng thiết bị định vị vệ tinh GPS. Trong mỗi OTC, thiết lập 5 ô dạng bản (ODB) diện tích 25 m2 (5m x 5m) tại 4 góc và 1 ô ở giao điểm hai đường chéo để thống kê toàn diện chất lượng và số lượng lớp cây tái sinh tự nhiên.

Phương pháp phân tích số liệu kết hợp chặt chẽ giữa xử lý thống kê sinh học và mô hình hóa toán học trên phần mềm chuyên dụng SPSS. Lý do lựa chọn hệ phương pháp này xuất phát từ yêu cầu loại bỏ tính chủ quan trong đánh giá cấu trúc rừng hỗn loài khác tuổi. Toàn bộ số liệu thô được lọc bỏ giá trị ngoại lai theo công thức Brooks và Carruthers. Kiểm định tính thuần nhất giữa các mẫu điều tra được thực hiện qua tiêu chuẩn phi tham số Mann-Whitney U đối với 2 mẫu và tiêu chuẩn Kruskal-Wallis đối với nhiều mẫu. Mức độ phù hợp giữa tần số thực nghiệm và tần số lý luận của các hàm phân bố được kiểm chứng nghiêm ngặt bằng tiêu chuẩn Khi bình phương (χ2) của Pearson. Các dạng phương trình tương quan sinh trưởng chiều cao - đường kính (Hvn/D1.3) được đánh giá qua hệ số tương quan R, tiêu chuẩn t-Student và giá trị kiểm định F-Fisher. Timeline thực hiện khảo sát ngoại nghiệp và tổng hợp dữ liệu hoàn tất trong giai đoạn 2007 - 2010.

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

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

Phân tích định lượng từ 47 ô tiêu chuẩn đã chỉ ra những quy luật cấu trúc mang tính đặc thù cao của thảm thực vật vùng cao núi đá Hà Giang:

Thứ nhất, tính đa dạng và cấu trúc tổ thành loài có sự phân hóa sâu sắc theo đai độ cao và điều kiện lập địa. Ở trạng thái rừng tự nhiên phục hồi IIA, số lượng loài tăng vọt từ 30 loài tại đai dưới 700m lên 50 loài tại đai từ 700m đến 1.000m. Nhóm loài ưu thế sinh thái của trạng thái IIA gồm 4 loài chiếm 53,68% tổng chỉ số IV% (Vối thuốc chiếm 33,13%, Cáng lò 9,21%, Sau sau 6,45%, Trâm lá dài 4,90%). Tại đai dưới 700m chịu áp lực khai thác mạnh, Vối thuốc chiếm tỷ trọng áp đảo với IV đạt 35,52%, trong khi ở đai cao trên 700m xuất hiện thêm nhiều loài mục đích giá trị cao như Re hương, Kháo xanh (IV đạt 5,48%).

Thứ hai, trạng thái rừng ẩm thường xanh núi đất xen đá IIIA2 thể hiện ưu thế sinh thái vượt trội của các loài cây gỗ lớn bản địa. Mật độ lâm phần trạng thái IIIA2 dao động từ 210 đến 300 cây/ha, đường kính thân bình quân đạt 25,4 - 28,6 cm và tổng tiết diện ngang đạt từ 17,46 m2/ha đến 31,70 m2/ha. Nhóm 5 loài ưu thế (Nghiến, Vối thuốc, Trai lý, Thị rừng, Cáng lò) chiếm tới 71,75% tổng giá trị IV%. Đáng chú ý, loài Nghiến giữ vai trò loài chỉ thị ưu thế tuyệt đối với IV đạt 37,91% tính chung, và đạt đỉnh 43,56% tại đai độ cao từ 700m đến 1.000m.

Thứ ba, rừng trên núi đá vôi trạng thái IIA phản ánh trạng thái suy thoái nghiêm trọng do điều kiện khắc nghiệt và tác động nhân sinh. Mật độ cây gỗ biến động rất mạnh từ 220 đến 420 cây/ha, tổng tiết diện ngang luôn ở mức thấp (dưới 10 m2/ha). Đường kính bình quân của cây gỗ sụt giảm từ 17,3 cm ở đai dưới 700m xuống chỉ còn 9,77 cm ở đai từ 700m đến 1.000m. Tuy nhiên, kiểu rừng này lại là nơi cư trú duy nhất của các loài hạt trần quý hiếm có nguy cơ tuyệt chủng như Bách vàng, Thông đỏ, Đỉnh tùng và Thông tre.

Thứ tư, trạng thái rừng phục hồi sau nương rẫy IIB tại đai 700m - 1.000m ghi nhận 42 loài tham gia tổ thành với mật độ 420 - 470 cây/ha, tổng tiết diện ngang dao động 11,00 - 14,95 m2/ha. Tổ thành được dẫn dắt bởi 4 loài ưu thế gồm Vối thuốc (16,55%), Dẻ cau (11,75%), Trám đen (8,51%) và Dẻ gai (6,23%), chiếm tổng cộng 43,03% chỉ số IV%.

Thảo luận kết quả

Nguyên nhân cốt lõi tạo nên sự phân hóa cấu trúc rừng tại bốn huyện vùng cao là sự giao thoa phức tạp giữa yếu tố địa chất Karst đặc thù và tập quán nhân sinh. Địa hình núi đá vôi có độ dốc lớn, khả năng giữ nước mặt cực kém kết hợp với tầng đất mỏng khiến tốc độ sinh trưởng của thảm thực vật diễn ra rất chậm. Đồng thời, áp lực từ cộng đồng dân cư với hơn 60% là đồng bào dân tộc Mông, tỷ lệ hộ nghèo lên tới 60,12% và tỷ lệ lao động nông nghiệp nhàn rỗi trên 40% đã đẩy mạnh hoạt động chặt củi, khai thác gỗ làm nhà và phát nương làm rẫy, khiến nhiều diện tích rừng tự nhiên bị đẩy lùi thành trảng cây bụi hoặc đất trống đồi núi trọc.

Khi so sánh với các nghiên cứu cấu trúc rừng thường xanh tại Tây Nguyên hoặc vùng đệm Bắc Trung Bộ, quy luật phân bố số cây theo cấp kính (N/D) và cấp chiều cao (N/H) tại Hà Giang vẫn tuân theo đường cong giảm của hàm Meyer và hàm Weibull, nhưng có độ dốc suy giảm mạnh hơn ở các cỡ kính lớn trên 40 cm. Điều này khẳng định rừng vùng cao phía Bắc Hà Giang đã bị khai thác kiệt các tầng cây vượt tán trong quá khứ.

Để trực quan hóa các kết quả này trong công tác lâm sinh, dữ liệu cấu trúc không gian tầng tán được mô tả tối ưu thông qua phẫu diện đồ đứng và phẫu diện đồ ngang, làm lộ rõ tầng tích tụ tán lá ở độ cao 8 - 14m. Bên cạnh đó, quan hệ sinh trưởng giữa đường kính và chiều cao vút ngọn được thể hiện bằng các biểu đồ đường cong hồi quy phi tuyến, cho phép các kỹ sư lâm nghiệp tra cứu chiều cao chuẩn xác dựa trên đường kính đo đếm mà không cần đo cao toàn bộ lâm phần.

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

Dựa trên các phát hiện thực chứng về cấu trúc và quy luật tái sinh, bốn nhóm giải pháp kỹ thuật lâm sinh đồng bộ được đề xuất nhằm phát triển bền vững tài nguyên rừng:

  1. Thực thi khoanh nuôi xúc tiến tái sinh kết hợp làm giàu rừng tự nhiên: Áp dụng cho các trạng thái rừng phục hồi IIA và IIB. Chủ thể thực hiện là Ban Quản lý Rừng phòng hộ, Ban Quản lý 661 cùng các cộng đồng thôn bản. Mục tiêu cụ thể là nâng mật độ cây tái sinh mục đích lên mức 1.500 - 2.000 cây/ha, cải thiện độ tàn che của lâm phần đạt trên 0,6 trong khung thời gian 5 năm. Kỹ thuật can thiệp tập trung vào việc tỉa thưa cục bộ các cây sâu bệnh, phát dây leo chèn ép và trồng bổ sung các loài cây bản địa có giá trị cao như Vối thuốc, Cáng lò, Re hương, Trám đen.

  2. Quy chuẩn hóa mật độ và cơ cấu loài cho rừng trồng phòng hộ và sản xuất: Ứng dụng công thức mật độ tối ưu Keller để điều chỉnh mật độ rừng trồng thuần loài (Sa mộc, Thông mã vĩ) và rừng hỗn giao (Vối thuốc - Sa mộc - Thông) về ngưỡng 1.100 - 1.600 cây/ha tùy theo cấp tuổi. Cơ quan khuyến nông và lâm nghiệp địa phương cần chuyển giao quy trình gieo ươm cây con bản địa từ nguồn giống tại chỗ, đảm bảo tỷ lệ sống của cây trồng rừng đạt trên 85% sau 3 năm đầu thiết lập.

  3. Thiết lập đai phòng hộ nghiêm ngặt và hành lang bảo tồn loài nguy cấp trên núi đá vôi: UBND tỉnh Hà Giang cùng Chi cục Kiểm lâm cần ưu tiên bảo vệ tuyệt đối diện tích 6.426,0 ha rừng đặc dụng và 51.953,4 ha rừng phòng hộ tự nhiên hiện có. Triển khai ngay đề án bảo tồn gen đặc hữu đối với quần thể các loài Bách vàng, Thông đỏ, Đỉnh tùng và Nghiến tại huyện Đồng Văn và Mèo Vạc. Target metric là chấm dứt hoàn toàn tình trạng khai thác trộm gỗ quý, giảm thiểu 95% nguy cơ cháy rừng mùa khô trong giai đoạn 3 năm.

  4. Đổi mới chính sách giao khoán và phát triển sinh kế lâm nghiệp bền vững: Sở Nông nghiệp và PTNT cần phối hợp với chính quyền địa phương hoàn thiện giao đất lâm nghiệp gắn với cấp giấy chứng nhận quyền sử dụng đất cho 18.340 hộ gia đình. Tăng mức chi trả tiền giao khoán bảo vệ rừng từ các nguồn quỹ dịch vụ môi trường rừng, kết hợp phát triển mô hình nông lâm kết hợp trồng cây ăn quả, dược liệu dưới tán rừng (như Thảo quả, Hoàng tinh) nhằm nâng mức thu nhập kinh tế của người dân vùng đệm thêm 25 - 30% sau 4 năm thực hiện.

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

Nội dung và dữ liệu thực nghiệm của luận văn mang lại giá trị thực tiễn và học thuật cho 4 nhóm đối tượng cụ thể:

  1. Cơ quan quản lý nhà nước và lực lượng kiểm lâm địa phương: Sở Nông nghiệp và PTNT tỉnh Hà Giang, Chi cục Kiểm lâm và Ban Quản lý Rừng phòng hộ các huyện Quản Bạ, Yên Minh, Đồng Văn, Mèo Vạc sử dụng tài liệu làm căn cứ quy hoạch 3 loại rừng theo Chỉ thị 38/CT-TTg, xây dựng đề án phân bổ vốn đầu tư lâm nghiệp chính xác trên 140.392,3 ha đất rừng.

  2. Các nhà nghiên cứu, giảng viên và học viên chuyên ngành Lâm học: Kế thừa bộ số liệu chi tiết về 47 ô tiêu chuẩn, các công thức tổ thành loài theo chỉ số IV%, hệ số tương quan sinh trưởng Hvn/D1.3 và phương pháp ứng dụng thống kê toán học sinh học để thực hiện các nghiên cứu chuyên sâu về sinh thái rừng núi đá Karst.

  3. Ban quản lý Công viên Địa chất Toàn cầu Cao nguyên đá Đồng Văn và các tổ chức bảo tồn sinh thái: Khai thác dữ liệu phân bố và cấu trúc của các loài thực vật quý hiếm (Nghiến, Trai lý, Bách vàng, Thông tre) nhằm xây dựng các chương trình cứu hộ nguồn gen, đánh giá sức tải sinh thái và phát triển du lịch sinh thái bền vững.

  4. Các chủ rừng, doanh nghiệp lâm nghiệp và cộng đồng dân cư vùng cao: Ứng dụng quy trình chọn loài cây bản địa (Vối thuốc, Cáng lò, Sa mộc) và kỹ thuật điều chỉnh mật độ lâm phần tối ưu nhằm chuyển đổi mô hình rừng trồng năng suất thấp sang mô hình rừng gỗ lớn đa mục đích kết hợp bảo vệ nguồn nước ngầm.

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

1. Đặc điểm cấu trúc nổi bật nhất của rừng tự nhiên trên núi đá vôi tại vùng cao Hà Giang là gì?
Rừng núi đá vôi tại khu vực có cấu trúc phân tầng phức tạp nhưng mật độ cây gỗ biến động rất lớn, dao động từ 220 đến 420 cây/ha với tổng tiết diện ngang thấp dưới 10 m2/ha. Điểm đặc biệt là sự hiện diện của các loài cây bản địa quý hiếm như Nghiến chiếm tỷ lệ ưu thế sinh thái lên tới 43,56% ở đai 700 - 1.000m, cùng sự tồn tại độc duy của các loài hạt trần nguy cấp như Bách vàng và Thông đỏ.

2. Tại sao chỉ số giá trị quan trọng (IV%) lại chuẩn xác hơn tỷ lệ số cây đơn thuần khi đánh giá tổ thành rừng?
Chỉ số IV% kết hợp đồng thời tỷ lệ số cây (N%) đại diện cho mặt sinh thái và tỷ lệ tiết diện ngang (G%) đại diện cho mặt sinh khối, sản lượng. Trong thực tế rừng hỗn loài, một loài cây có số lượng cây non rất nhiều nhưng đường kính nhỏ sẽ không đóng vai trò chi phối bằng những loài cây gỗ lớn chiếm lĩnh không gian tầng tán như Vối thuốc hay Nghiến.

3. Mô hình toán học nào phản ánh tốt nhất quy luật phân bố số cây theo đường kính của lâm phần nghiên cứu?
Quy luật phân bố số cây theo cỡ kính (N/D1.3) của các trạng thái rừng tự nhiên hỗn loài tại Hà Giang được mô phỏng phù hợp nhất bằng hàm giảm Meyer (dạng số mũ) và phân bố xác suất Weibull. Khi kiểm định bằng tiêu chuẩn Khi bình phương (χ2) của Pearson, các hàm này đều cho giá trị sai số nhỏ nhất và mô tả chính xác mật độ tập trung cao ở các cấp kính nhỏ từ 8 - 16 cm.

4. Khí hậu và địa chất Karst tại 4 huyện phía Bắc Hà Giang tác động như thế nào đến tái sinh rừng?
Địa hình đá vôi chia cắt mạnh, tỷ lệ đá lẫn từ 30% đến 70% cùng lượng mưa phân bố không đều dẫn đến hiện tượng khô hạn sinh lý gay gắt trong mùa đông. Điều này khiến hạt giống khó nảy mầm và tỷ lệ sống của cây mạ rất thấp. Tái sinh tự nhiên chủ yếu diễn ra theo kiểu phân tán liên tục của cây chịu bóng dưới tán rừng khép tán hoặc tái sinh vệt ở các lỗ trống tán rừng nhỏ dưới 1.000 m2.

5. Loài cây bản địa nào được khuyến nghị ưu tiên trong các mô hình trồng rừng và làm giàu rừng tại địa phương?
Các loài cây bản địa có biên độ sinh thái rộng, sinh trưởng tốt và giá trị kinh tế - phòng hộ cao gồm Vối thuốc, Cáng lò, Dẻ cau, Trám đen và Sa mộc. Đối với khu vực sườn núi đất và thung lũng, Vối thuốc và Cáng lò là hai loài tiên phong cải tạo đất tốt nhất với hệ số tổ thành IV% luôn chiếm vị trí ưu thế hàng đầu trong các lâm phần phục hồi.

Kết luận

  • Luận văn đã định lượng toàn diện quy luật cấu trúc tổ thành, mật độ, tầng tán và tái sinh của 4 trạng thái rừng tự nhiên cùng hệ thống rừng trồng trên địa bàn 234.763,0 ha đất thuộc bốn huyện vùng cao Hà Giang.
  • Xác lập vai trò thống trị sinh thái của các loài cây bản địa mục đích, trong đó loài Nghiến đạt chỉ số IV% tới 37,91% ở trạng thái IIIA2 và Vối thuốc đạt IV% từ 16,55% đến 35,52% ở các trạng thái IIA, IIB.
  • Chứng minh tính ứng dụng cao của các mô hình toán học sinh học (hàm Meyer, Weibull, kiểm định Pearson χ2) trong việc mô phỏng chính xác cấu trúc lâm phần vùng núi đá Karst.
  • Cung cấp cơ sở khoa học lâm sinh chuẩn xác thay thế phương thức khoanh nuôi thụ động, định hình rõ giải pháp tỉa thưa, làm giàu rừng và quản lý mật độ tối ưu.
  • Đóng góp luận cứ kỹ thuật quan trọng phục vụ chiến lược quy hoạch bảo tồn hệ sinh thái rừng phòng hộ đầu nguồn và đa dạng sinh học tại Công viên Địa chất Toàn cầu UNESCO.

Trong giai đoạn 3 đến 5 năm tới, các cơ quan lâm nghiệp địa phương cần tập trung triển khai thực nghiệm các mô hình kỹ thuật lâm sinh được đề xuất vào các dự án trồng rừng thay thế và phục hồi rừng đầu nguồn. Các cấp quản lý, nhà nghiên cứu và chủ rừng hãy tích cực khai thác, ứng dụng các phát hiện từ luận văn để chung tay bảo tồn màu xanh bền vững trên miền cao nguyên đá cực Bắc tổ quốc.