Tổng quan nghiên cứu

Khu bảo tồn thiên nhiên Pù Hu nằm ở phía Tây tỉnh Thanh Hóa, trải rộng trên địa giới hành chính của 11 xã thuộc hai huyện Quan Hóa và Mường Lát. Khu bảo tồn sở hữu tổng diện tích tự nhiên 27.502,89 ha, trong đó phân khu bảo vệ nghiêm ngặt chiếm 16.246,74 ha và phân khu phục hồi sinh thái là 11.238,15 ha. Địa hình khu vực có sự chia cắt sâu sắc với độ cao biến thiên từ dưới 50 m lên đến đỉnh cao nhất là đỉnh Pù Hu ở mức 1.470 m. Trong cấu trúc cảnh quan, hệ sinh thái rừng núi đất đóng vai trò chủ đạo khi chiếm tới 98,55% tổng diện tích đất lâm nghiệp. Khu hệ thực vật tại đây ghi nhận 508 loài thuộc 323 chi và 102 họ, bao gồm 28 loài thực vật quý hiếm nằm trong Sách Đỏ Việt Nam và Nghị định 48/2002/NĐ-CP như Sến mật, Lát hoa, Kim giao.

Tuy nhiên, áp lực dân sinh với 3.894 hộ gia đình và 19.360 nhân khẩu sinh sống tại vùng lõi và vùng đệm (mật độ dân số 36 người/km², tỷ lệ tăng dân số dân tộc H’Mông đạt 3,3% và người Thái đạt 2,9%) đã khiến tài nguyên rừng bị tác động nghiêm trọng qua tập quán canh tác nương rẫy và khai thác lâm sản tự do. Vấn đề nghiên cứu trọng tâm là làm rõ quy luật biến đổi cấu trúc lâm phần và đa dạng thực vật theo đai độ cao, từ đó khắc phục tình trạng suy thoái sinh thái.

Mục tiêu cụ thể của luận văn là định lượng hóa cấu trúc tầng cây cao, tầng cây tái sinh, phân loại trạng thái rừng núi đất từ đai 600 m đến 1.200 m, và đề xuất các giải pháp kỹ thuật lâm sinh quản lý bền vững. Nghiên cứu mang ý nghĩa then chốt khi cung cấp cơ sở khoa học chuẩn xác giúp tăng độ che phủ rừng lên trên 75%, bảo tồn nguồn gen bản địa và định hình mô hình phục hồi rừng tiếp cận tự nhiên cho vùng Tây Bắc Thanh Hóa.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Luận văn vận dụng học thuyết hệ sinh thái rừng của Tansley (1935), Sukachev (1957) và quan điểm sinh thái phát sinh thảm thực vật rừng Việt Nam của Thái Văn Trừng (1978). Cấu trúc rừng được tiếp cận toàn diện theo quan điểm lâm sinh của Phùng Ngọc Lan (1986) và quan điểm sản lượng của Husch (1982), coi cấu trúc là sự sắp xếp tổ chức nội bộ của thực vật theo không gian và thời gian.

Nghiên cứu ứng dụng các khái niệm và mô hình định lượng nền tảng:

  • Tổ thành loài theo chỉ số giá trị quan trọng (Important Value - IV%) của Daniel Marmillod và Vũ Đình Huề (1984), xác định loài ưu thế sinh thái khi IV% vượt ngưỡng 5% và nhóm ưu hợp thực vật chiếm trên 50% tổng cá thể.
  • Phân loại trạng thái rừng theo hệ thống Loetschau (1960) được Viện Điều tra Quy hoạch rừng cải tiến, phân định rõ các nhóm kiểu rừng non phục hồi (IIa, IIb), rừng thứ sinh bị tác động (IIIa1, IIIa2, IIIa3, IIIb) và rừng nguyên sinh thành thục (IVa, IVb).
  • Đánh giá đa dạng sinh học qua chỉ số phong phú loài $R$ của Kjayarama (2000), chỉ số đa dạng Shannon - Wiener ($H'$) và chỉ số ưu thế Simpson.
  • Mô hình hóa phân bố số cây theo đường kính ($N/D_{1.3}$) và chiều cao ($N/H_{vn}$) bằng các hàm phân bố xác suất lý thuyết như Weibull, hàm Khoảng cách (Distance Distribution) và hàm giảm Meyer.

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 20 ô tiêu chuẩn (OTC) điển hình tạm thời, mỗi ô có diện tích 1.000 m² (tổng diện tích điều tra đạt 20.000 m² tương đương 2,0 ha). Các OTC được bố trí đại diện trên 4 đai độ cao chính gồm: 600 m, 800 m, 1.000 m và 1.200 m so với mực nước biển (mỗi đai lập 5 OTC độc lập). Phương pháp chọn mẫu điển hình được lựa chọn nhằm phản ánh đầy đủ các trạng thái lập địa đặc trưng của sườn núi đất và hạn chế tính áp đặt chủ quan.

Quy trình đo đếm thực địa được tiến hành trên toàn bộ cây gỗ tầng cao có đường kính ngang ngực $D_{1.3} \ge 6$ cm bằng thước kẹp kính theo 2 chiều trực giao (Đông - Tây và Nam - Bắc), chiều cao vút ngọn $H_{vn}$ đo bằng thước đo cao Blumeleiss, đường kính tán $D_t$ đo qua hình chiếu mặt đất, cùng việc thống kê toàn diện đường kính gốc $D_0$ và $H_{vn}$ của tầng cây tái sinh.

Dữ liệu được xử lý bằng phần mềm SPSS 13.0 và Excel. Đề tài sử dụng kiểm định $\chi^2$ (Chi-square) của Pearson với bậc tự do $k = m - 3$ để đánh giá mức độ tương thích giữa phân bố lý thuyết và thực nghiệm. Tiêu chuẩn Kolmogorov - Smirnov ($K_{\alpha} = 1,36$ tại mức ý nghĩa $\alpha = 0,05$) được áp dụng để so sánh sự khác biệt phân bố giữa các đai cao, kết hợp tiêu chuẩn $t$-Student để so sánh phương sai chỉ số đa dạng loài. Việc sử dụng toán tin giúp chuyển đổi các mô tả định tính truyền thống sang mô hình toán học chuẩn xác. Timeline nghiên cứu thực địa và phân tích dữ liệu được hoàn thành trọn vẹn trong chu kỳ 12 tháng.

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

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

Hệ sinh thái rừng núi đất Pù Hu phân bố chủ yếu từ độ cao 600 m đến 1.400 m, chiếm 98,55% diện tích tự nhiên với 3 kiểu rừng chính: kiểu rừng kín thường xanh mưa ẩm nhiệt đới ít bị tác động ở đai cao (>600 m), kiểu rừng thường xanh bị tác động trên đai cao, và kiểu rừng thứ sinh hỗn giao cây lá rộng xen tre nứa phục hồi sau nương rẫy ở đai thấp (<600 m).

Cấu trúc tổ thành loài thay đổi sâu sắc theo đai cao. Ở đai 600 m - 800 m, ưu thế thuộc về các loài nhiệt đới ưa sáng mọc nhanh như Ràng ràng (Ormosia balansae), Gáo (Anthocephalus chinensis), Vàng anh (Saraca dives), Sổ bà (Dillenia indica). Lên đai 1.000 m - 1.200 m, tổ thành chuyển dịch rõ nét sang các loài cây gỗ lớn á nhiệt đới thuộc các họ Dẻ (Fagaceae), Re (Lauraceae), Mộc lan (Magnoliaceae) và Dâu tằm (Moraceae), tiêu biểu là Sến mật (Madhuca pasquieri), Vàng tâm (Manglietia fordiana), Dổi xanh (Michelia mediocris) với tổng giá trị quan trọng IV% nhóm ưu thế đạt trên 55%.

Mức độ đa dạng loài có xu hướng giảm dần theo đai độ cao. Chỉ số phong phú loài $R$ đạt giá trị cao nhất tại đai 600 m ($R = 3,85 - 4,20$) và giảm đều xuống mức thấp nhất tại đai 1.200 m ($R = 2,90 - 3,15$). Tương tự, chỉ số Shannon - Wiener $H'$ tại đai 600 m đạt trên 3,60 nát/cá thể, trong khi ở đai 1.200 m giảm xuống dưới 2,85 nát/cá thể.

Quy luật phân bố số cây theo cỡ kính ($N/D_{1.3}$) và cỡ chiều cao ($N/H_{vn}$) trên cả 4 đai cao đều thích hợp nhất với hàm phân bố Weibull và hàm Khoảng cách (giá trị $\chi^2$ thực nghiệm đều nhỏ hơn $\chi^2_{0.05}$). Đường cong phân bố có dạng một đỉnh lệch trái rõ rệt với tham số hình dạng $\alpha < 3$, trong đó số lượng cá thể tập trung chủ yếu ở cỡ kính nhỏ từ 8 cm đến 20 cm (chiếm từ 65% đến 72% tổng số cây điều tra), số cây có $D_{1.3} > 40$ cm chỉ chiếm khoảng 4% đến 6%.

Thảo luận kết quả

Nguyên nhân căn bản dẫn đến sự biến đổi cấu trúc tổ thành và đa dạng loài là quy luật gradient nhiệt ẩm. Càng lên cao, nhiệt độ không khí giảm trung bình 0,5°C/100 m, lượng mưa dao động 1.400 - 1.500 mm phân bổ không đều và tốc độ gió sườn núi tăng cao, tạo nên điều kiện sống khắc nghiệt hơn, chọn lọc tự nhiên chỉ giữ lại các loài thực vật á nhiệt đới có biên độ sinh thái thích hợp.

Kết quả này hoàn toàn tương đồng với các phát hiện về cấu trúc rừng kín thường xanh tại Kon Hà Nừng của Lê Sáu (1996) và Hương Sơn của Đào Công Khanh (1996), nơi hàm Weibull chứng minh tính ưu việt vượt trội trong mô phỏng kết cấu lâm phần tự nhiên hỗn loài dị tuổi.

Dữ liệu nghiên cứu được biểu diễn trực quan thông qua biểu đồ phân bố tần số $N/D_{1.3}$ dạng đường cong giảm liên tục, kết hợp đồ thị hàm Weibull thể hiện đỉnh tái sinh mạnh ở cỡ kính 12 - 16 cm. Bảng ma trận so sánh kiểm định Kolmogorov - Smirnov giữa các cặp đai cao cho thấy giá trị sai lệch lớn nhất $\hat{D}$ giữa đai 600 m và đai 1.200 m đều vượt qua ngưỡng giới hạn $D_{(\alpha)} = 0,185$, khẳng định sự khác biệt có ý nghĩa thống kê về cấu trúc mật độ giữa vùng núi thấp và đai cao.

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

  1. Khoanh nuôi phục hồi và làm giàu rừng tại đai thấp dưới 800 m:

    • Hành động: Áp dụng biện pháp khoanh nuôi xúc tiến tái sinh tự nhiên có trồng bổ sung các loài cây bản địa bản địa mục đích như Dổi xanh, Trám đen, Lát hoa trên diện tích rừng trạng thái IIIa1 và IIIa2.
    • Target metric: Nâng độ tàn che lâm phần từ mức 0,3 - 0,5 lên $\ge 0,7$, mật độ cây tái sinh mục đích đạt trên 2.500 cây/ha sau 3 năm.
    • Timeline: Giai đoạn 2026 - 2028.
    • Chủ thể thực hiện: Ban quản lý Khu bảo tồn thiên nhiên Pù Hu chủ trì, phối hợp với cộng đồng dân cư các xã vùng đệm.
  2. Bảo tồn nghiêm ngặt đai cao từ 1.000 m đến 1.470 m:

    • Hành động: Thiết lập các trạm chốt kiểm lâm cố định và tăng cường tuần tra ứng dụng công nghệ định vị GPS nhằm bảo vệ nguyên vẹn các quần xã Sến mật, Vàng tâm, Kim giao.
    • Target metric: Giảm 100% các vụ vi phạm khai thác gỗ quý trái phép, bảo toàn 16.246,74 ha vùng bảo vệ nghiêm ngặt.
    • Timeline: Thực hiện liên tục hàng năm từ năm 2026.
    • Chủ thể thực hiện: Hạt Kiểm lâm Khu bảo tồn thiên nhiên Pù Hu phối hợp cùng lực lượng kiểm lâm liên huyện Quan Hóa - Mường Lát.
  3. Phát triển sinh kế và giao khoán rừng gắn với chi trả dịch vụ môi trường rừng:

    • Hành động: Giao khoán bảo vệ rừng ổn định lâu dài cho 3.894 hộ dân gắn với chính sách chi trả dịch vụ môi trường rừng lưu vực sông Mã và sông Luồng; hỗ trợ kỹ thuật thâm canh 9.000 ha rừng Luồng kinh tế.
    • Target metric: Tăng thu nhập từ lâm nghiệp của người dân lên 25 - 30%, giảm tỷ lệ thiếu lương thực giáp hạt từ 2 - 5 tháng xuống dưới 1 tháng/năm.
    • Timeline: Giai đoạn 2026 - 2030.
    • Chủ thể thực hiện: Ủy ban nhân dân hai huyện Quan Hóa, Mường Lát và Chi cục Kiểm lâm tỉnh Thanh Hóa.
  4. Xây dựng mạng lưới ô định vị giám sát diễn thế sinh thái dài hạn:

    • Hành động: Định vị cố định 20 ô tiêu chuẩn điều tra bằng mốc bê tông kiên cố trên 4 đai cao để theo dõi định kỳ động thái tăng trưởng đường kính, trữ lượng và biến động đa dạng sinh học.
    • Target metric: Số hóa 100% cơ sở dữ liệu giám sát cấu trúc lâm phần và cập nhật danh mục thực vật định kỳ.
    • Timeline: Chu kỳ đo lường 3 năm/lần bắt đầu từ năm 2027.
    • Chủ thể thực hiện: Phòng Khoa học và Hợp tác quốc tế thuộc Khu bảo tồn thiên nhiên Pù Hu liên kết với Trường Đại học Lâm nghiệp.

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

  • Ban quản lý các Khu bảo tồn thiên nhiên và Vườn quốc gia: Tiếp nhận tài liệu khoa học và phương pháp luận định lượng để xây dựng phương án quản lý rừng bền vững giai đoạn 2026 - 2035, tối ưu hóa phân khu chức năng theo đai độ cao cho các khu bảo tồn vùng núi phía Bắc và Bắc Trung Bộ.
  • Nhà nghiên cứu lâm học và chuyên gia sinh thái: Sử dụng nguồn dữ liệu thực chứng và bộ tham số hàm Weibull, hàm Khoảng cách để làm đối chứng cho các công trình nghiên cứu cấu trúc, mô hình hóa sinh khối và động thái thảm thực vật nhiệt đới.
  • Giảng viên, học viên cao học và sinh viên chuyên ngành Quản lý tài nguyên rừng: Khai thác làm tài liệu tham khảo chuẩn mực về kỹ thuật điều tra lâm học, ứng dụng phần mềm SPSS trong phân tích toán sinh học, và phương pháp phân loại trạng thái rừng tự nhiên.
  • Cơ quan quản lý nhà nước và các tổ chức bảo tồn quốc tế (IUCN, WWF): Làm căn cứ thực tiễn để hoạch định chính sách bảo vệ nguồn gen nguy cấp, xây dựng dự án chi trả dịch vụ hệ sinh thái rừng và phân bổ ngân sách hỗ trợ sinh kế giảm nghèo cho đồng bào dân tộc thiểu số vùng cao.

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

Tại sao luận văn lựa chọn hàm phân bố Weibull để mô phỏng cấu trúc lâm phần tại Pù Hu? Hàm Weibull có tính linh hoạt rất cao nhờ tham số hình dạng $\alpha$ và tham số tỷ lệ $\gamma$. Tại Pù Hu, với tham số $\alpha < 3$, hàm Weibull mô tả chính xác dạng phân bố lệch trái của rừng tự nhiên phục hồi, nơi mật độ cây nhỏ chiếm ưu thế trên 65%, vượt trội hoàn toàn so với phân bố chuẩn Gauss.

Cấu trúc tổ thành loài thực vật thay đổi như thế nào giữa đai 600 m và đai 1.200 m? Tại đai 600 m, lâm phần ưu thế bởi các loài cây nhiệt đới ưa sáng mọc nhanh như Ràng ràng, Gáo, Vàng anh. Khi lên đai 1.200 m, các loài cây á nhiệt đới chịu lạnh và ưa ẩm như Sến mật, Vàng tâm, Dổi xanh thuộc họ Dẻ và Re chiếm ưu thế với chỉ số IV% trên 55%.

Áp lực dân sinh tác động cụ thể ra sao đến thảm thực vật khu bảo tồn? Với 3.894 hộ dân sống xen kẽ trong vùng lõi và vùng đệm, việc du canh nương rẫy và khai thác lâm sản phụ tự do trước năm 1998 đã làm suy thoái thảm thực vật nguyên sinh, biến nhiều diện tích thành rừng thứ sinh nghèo kiệt (IIIa1, IIIa2) hoặc thảm tre nứa đơn ưu.

Phương pháp điều tra 20 ô tiêu chuẩn trên 4 đai cao có đảm bảo độ tin cậy thống kê? Có, việc thiết lập 20 OTC với tổng diện tích 20.000 m² (2,0 ha) trải đều trên 4 đai cao (600 m, 800 m, 1.000 m, 1.200 m) đã thu thập dữ liệu hàng nghìn cá thể cây gỗ, đáp ứng đầy đủ dung lượng mẫu theo lý thuyết điều tra lâm sinh và đạt độ tin cậy thống kê trên 95%.

Giải pháp nào giúp giải quyết hài hòa giữa bảo tồn đa dạng sinh học và phát triển kinh tế địa phương? Giải pháp đột phá là kết hợp giao khoán bảo vệ 21.523,44 ha đất có rừng gắn với chi trả dịch vụ môi trường rừng lưu vực sông Mã, kết hợp thâm canh 9.000 ha rừng Luồng thương phẩm, giúp tạo nguồn thu nhập ổn định cho 8.928 lao động địa phương.

Kết luận

  • Xác lập được đặc điểm phân bố của hệ sinh thái rừng núi đất chiếm 98,55% diện tích tự nhiên Khu bảo tồn thiên nhiên Pù Hu, phân định rõ 3 kiểu rừng chính chịu sự chi phối chặt chẽ của đai độ cao.
  • Ứng dụng thành công hàm Weibull và hàm Khoảng cách để lượng hóa chính xác quy luật phân bố số cây theo đường kính và chiều cao, phản ánh lâm phần đang trong giai đoạn phục hồi mạnh với trên 65% số cây thuộc cấp kính nhỏ.
  • Chứng minh quy luật suy giảm độ phong phú loài ($R$) và chỉ số đa dạng Shannon - Wiener ($H'$) từ đai thấp 600 m lên đai cao 1.200 m, cùng sự dịch chuyển tổ thành từ nhóm loài nhiệt đới sang nhóm ưu hợp á nhiệt đới (Fagaceae, Lauraceae, Magnoliaceae).
  • Đánh giá toàn diện hiện trạng suy thoái cấu trúc rừng dưới áp lực sinh kế của 3.894 hộ dân, làm cơ sở khoa học phân loại các trạng thái rừng thứ sinh phục vụ kinh doanh lâm sinh.
  • Đề xuất hệ thống 4 nhóm giải pháp kỹ thuật lâm sinh và chính sách sinh kế khả thi, tạo tiền đề nâng cao chất lượng rừng và bảo vệ 28 loài thực vật quý hiếm trong Sách Đỏ.

Đóng góp lớn nhất của luận văn là hoàn thiện cơ sở dữ liệu định lượng về sinh thái rừng núi đất Tây Bắc Thanh Hóa, chuyển hóa các quan sát định tính thành mô hình toán học chuẩn xác. Kế hoạch tiếp theo giai đoạn 2026 - 2030 là triển khai ứng dụng các mô hình kỹ thuật phục hồi rừng vào thực tiễn quản lý tại 11 xã của Khu bảo tồn thiên nhiên Pù Hu. Các nhà quản lý, nhà nghiên cứu và cơ quan lâm nghiệp hãy áp dụng ngay khung phương pháp luận và các giải pháp lâm sinh của đề tài nhằm thúc đẩy công tác bảo tồn đa dạng sinh học bền vững.