Tổng quan nghiên cứu

Vườn Quốc gia Cúc Phương là khu rừng đặc dụng có giá trị bảo tồn sinh học hàng đầu tại Việt Nam với tổng diện tích 22.200 ha, trong đó rừng tự nhiên chiếm 20.473 ha (tương đương 92,2% diện tích) mang đặc trưng của thảm thực vật rừng mưa nhiệt đới thường xanh trên núi đá vôi. Hệ thực vật tại đây vô cùng phong phú với 2.103 loài thuộc 917 chi và 231 họ thực vật bậc cao. Dù được bảo vệ nghiêm ngặt, các quá trình diễn thế tự nhiên cùng tác động xáo trộn trong lịch sử đã tạo nên sự phân hóa cấu trúc rõ rệt giữa các trạng thái rừng. Vấn đề cốt lõi được đặt ra là cần làm sáng tỏ quy luật cấu trúc lâm phần và động thái không gian của các loài cây gỗ dưới ảnh hưởng của xáo trộn sinh thái.

Mục tiêu cụ thể của công trình là xác định mật độ, tổ thành loài, phân bố số cây theo cấp đường kính, đánh giá tính đa dạng sinh học và phân tích mối quan hệ không gian cùng loài cũng như khác loài của các loài cây gỗ ưu thế. Nghiên cứu được triển khai tại vùng lõi Vườn Quốc gia Cúc Phương (thuộc địa giới ba tỉnh Ninh Bình, Hòa Bình và Thanh Hóa) với dữ liệu ngoại nghiệp thu thập từ ngày 11/11/2016 đến ngày 31/3/2017. Kết quả nghiên cứu mang lại giá trị khoa học vượt trội khi định lượng hóa các tương tác sinh thái với độ tin cậy 95%, tạo cơ sở dữ liệu chuẩn xác phục vụ công tác quy hoạch, phục hồi và phát triển bền vững hệ sinh thái rừng nhiệt đới.

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 của các lý thuyết sinh thái học rừng nhiệt đới kinh điển:

  • Giả thuyết Janzen-Connell: Giải thích cơ chế điều hòa mật độ khi các cá thể tái sinh ở xa cây mẹ có tỷ lệ sống sót cao hơn nhờ tránh được các sinh vật gây hại chuyên biệt tập trung quanh gốc cây mẹ.
  • Giả thuyết bảo vệ nhóm loài (Species Herd Protection): Mở rộng từ giả thuyết Janzen-Connell, khẳng định sự chung sống của nhiều loài khác nhau giúp phân tán và giảm thiểu rủi ro bị tấn công bởi sâu bệnh hại.
  • Lý thuyết ổ sinh thái (Niche Theory) và Lý thuyết trung lập (Neutral Theory): Lý thuyết ổ sinh thái nhấn mạnh sự phân hóa không gian và môi trường sống định hình đặc tính của loài, trong khi lý thuyết trung lập xem xét sự bình đẳng tương đối giữa các cá thể trong quá trình sinh trưởng và đào thải tự nhiên.

Các khái niệm then chốt bao gồm: Chỉ số giá trị quan trọng (IVI%) phản ánh vai trò sinh thái thông qua mật độ tương đối và tiết diện ngang tương đối; Phân bố hoàn toàn ngẫu nhiên (CSR) đóng vai trò làm mô hình không chuẩn tắc; và Hàm tương quan theo cặp g(r) một biến số g11(r) cùng hai biến số g12(r) nhằm đo lường mật độ điểm theo khoảng cách.

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

Nguồn dữ liệu được thu thập thông qua phương pháp lập ô tiêu chuẩn điển hình tạm thời trên 2 trạng thái rừng tự nhiên. Tổng cỡ mẫu nghiên cứu gồm 3 ô tiêu chuẩn (OTC) với diện tích 1 ha mỗi ô (kích thước 100 m × 100 m, tổng diện tích khảo sát 3 ha), gồm 1 ô rừng nguyên sinh (OTC 1) và 2 ô rừng thứ sinh (OTC 2 và OTC 3). Mỗi ô tiêu chuẩn 1 ha được chia thành lưới 100 ô đo đếm nhỏ diện tích 100 m2 (10 m × 10 m).

Toàn bộ các cây gỗ có đường kính ngang ngực (D1.3) từ 2,5 cm trở lên được đánh số, đo chu vi bằng thước dây chính xác 0,5 cm và xác định tọa độ không gian (x, y) bằng máy đo khoảng cách laser độ chính xác 0,1 cm kết hợp la bàn. Phương pháp phân tích mô hình điểm không gian sử dụng hàm tương quan theo cặp g(r) được lựa chọn thay vì hàm Ripley's K truyền thống nhằm loại bỏ hiệu ứng tích lũy khoảng cách, giúp nhận diện chính xác các tương tác trực tiếp ở cự ly hẹp từ 0 đến 50 m. Quy trình kiểm định thống kê thực hiện 199 lần mô phỏng Monte Carlo để xây dựng khoảng tin cậy 95% trên phần mềm chuyên dụng Programita trong giai đoạn từ tháng 11/2016 đến tháng 3/2017.

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

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

Quá trình điều tra và xử lý dữ liệu thực địa đã ghi nhận các kết quả định lượng nổi bật:

  • Mật độ và tiết diện ngang tầng cây gỗ: Rừng nguyên sinh (OTC 1) có mật độ trung bình 825 cây/ha với tổng tiết diện ngang đạt 36,00 m2/ha. Ở trạng thái rừng thứ sinh, mật độ biến động mạnh từ 601 cây/ha (OTC 2, tiết diện ngang 19,78 m2/ha) đến 1.006 cây/ha (OTC 3, tiết diện ngang đạt 36,34 m2/ha).
  • Cấu trúc tổ thành và đa dạng loài: Rừng nguyên sinh OTC 1 ghi nhận 106 loài cây gỗ với công thức tổ thành gồm Vàng anh (IV% = 9,32%), Nang trứng (7,49%), Phân mã (6,85%) và Thừng mực (5,06%). Tại OTC 2 có 82 loài, ưu thế thuộc về Vàng anh (12,58%) và Tèo nông (9,90%). Tại OTC 3 ghi nhận 89 loài, trong đó Nhò vàng chiếm ưu thế tuyệt đối với IV% đạt 25,72%, theo sau là Cà lồ (18,39%). Chỉ số Shannon-Wiener (H') đạt 0,860 tại OTC 1, tăng lên 0,964 tại OTC 2 và 0,941 tại OTC 3; chỉ số Simpson (D) lần lượt là 0,963; 0,953 và 0,825.
  • Cấu trúc đường kính (N/D): Phân bố số cây theo cấp đường kính ở cả 3 ô đều có dạng đường cong chữ J ngược giảm dần, tập trung cao ở các cấp kính nhỏ từ 5 đến 15 cm. Riêng OTC 1 có 135 cây đạt đường kính trên 25 cm và xuất hiện các cá thể cổ thụ đạt đường kính xấp xỉ 160 cm, trong khi rừng thứ sinh chủ yếu tập trung cây đường kính trung bình từ 20 đến 60 cm.
  • Quan hệ không gian: Phân bố cùng loài dạng cụm chiếm ưu thế ở hầu hết các loài ưu thế trong cự ly 1 đến 34 m (Phân mã cụm ở 3-8 m; Nang trứng cụm ở 4-8 m; Nhò vàng cụm từ 1-34 m). Đáng chú ý, 100% các cặp loài ưu thế khi phân tích hàm g12(r) đều thể hiện mối quan hệ không gian độc lập trong dải khoảng cách từ 0 đến 50 m.

Thảo luận kết quả

Các phát hiện nghiên cứu phản ánh trực quan quy luật phục hồi và diễn thế rừng nhiệt đới. Sự xuất hiện dày đặc của các loài tiên phong ưa sáng như Nhò vàng (mật độ 392 cây/ha tại OTC 3) và Vàng anh là minh chứng rõ ràng cho tác động xáo trộn tán rừng trong quá khứ, tạo điều kiện thuận lợi cho hạt giống nảy mầm. Dữ liệu khi trực quan hóa qua biểu đồ tọa độ không gian (Point pattern map) và biểu đồ phân bố N/D cho thấy cấu trúc rừng nguyên sinh có độ ổn định sinh thái cao với tầng cây vượt tán đường kính lớn, trong khi rừng thứ sinh thể hiện tính linh hoạt trong quá trình tự điều chỉnh mật độ.

Chỉ số đa dạng Shannon-Wiener ở rừng thứ sinh (H' = 0,964 tại OTC 2) cao hơn rừng nguyên sinh (H' = 0,860 tại OTC 1) do sự cộng hưởng giữa các loài bản địa sót lại và các loài tái sinh ưa sáng xâm nhập vào các khoảng trống. Kiểu phân bố cụm ở cự ly dưới 35 m là hệ quả tất yếu của cơ chế phát tán hạt hạn chế quanh gốc cây mẹ. Việc 100% các cặp loài ưu thế có quan hệ độc lập (nằm trọn trong khoảng tin cậy 95% của mô hình không) khẳng định sự cạnh tranh trực tiếp giữa các loài trưởng thành đã giảm thiểu tối đa nhờ môi trường sống đồng nhất và sự phân hóa ổ sinh thái hoàn hảo, tương đồng với các kết luận nghiên cứu lâm học tại khu vực Đông Nam Á.

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

Dựa trên các bằng chứng khoa học thực nghiệm, 4 nhóm giải pháp kỹ thuật và quản lý được đề xuất cụ thể:

  • Thiết lập mạng lưới 10 ô tiêu chuẩn định vị vĩnh viễn (PSP) diện tích 1 ha tại vùng lõi Vườn Quốc gia Cúc Phương trong giai đoạn 2026-2028, do Ban Quản lý Vườn chủ trì, nhằm theo dõi liên tục động thái diễn thế và tốc độ tăng trưởng đường kính theo chu kỳ 5 năm một lần.
  • Áp dụng biện pháp lâm sinh khoanh nuôi xúc tiến tái sinh tự nhiên kết hợp tỉa thưa cục bộ tại khoảng 1.500 ha rừng thứ sinh nghèo kiệt (mật độ dưới 600 cây/ha) trước năm 2028, ưu tiên bảo tồn các loài cây gỗ bản địa tái sinh có giá trị cao như Chò xanh, Chò chỉ, Đăng và Vàng anh.
  • Triển khai chương trình đồng quản lý rừng và tạo sinh kế bền vững cho 100% cộng đồng dân cư thuộc các xã vùng đệm (nơi có mật độ dân số từ 23 đến 559 người/km2), đặt mục tiêu giảm thiểu trên 85% các vụ việc xâm hại tài nguyên và khai thác lâm sản ngoài gỗ trái phép đến năm 2030.
  • Xây dựng hệ thống cơ sở dữ liệu số hóa WebGIS quản lý cấu trúc không gian rừng cho toàn bộ 20.473 ha rừng tự nhiên, hoàn thành tích hợp thuật toán phân tích điểm thống kê trong vòng 18 tháng nhằm hỗ trợ công tác dự báo cháy rừng và giám sát suy thoái sinh thái.

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

Công trình là tài liệu tham khảo chuyên sâu dành cho các nhóm đối tượng:

  • Ban quản lý các vườn quốc gia và khu bảo tồn thiên nhiên: Ứng dụng mô hình phân bố không gian và tổ thành loài để xây dựng phương án kỹ thuật phục hồi cho hơn 20.000 ha rừng đặc dụng có điều kiện lập địa karst tương tự.
  • Cán bộ nghiên cứu và giảng viên ngành Lâm học, Sinh thái học: Sử dụng phương pháp luận phân tích mô hình điểm không gian bằng hàm tương quan theo cặp g(r) và kiểm định Monte Carlo 199 lần lặp để phát triển các đề tài nghiên cứu sinh thái chuyên sâu.
  • Học viên sau đại học và sinh viên ngành Quản lý tài nguyên rừng: Tiếp cận bộ số liệu thực nghiệm chuẩn xác từ 3 ha điều tra chi tiết với hơn 2.400 cá thể cây gỗ, làm tài liệu học tập về phương pháp điều tra lâm học hiện đại.
  • Cơ quan hoạch định chính sách và tổ chức bảo tồn: Tham khảo các chỉ số định lượng về mật độ, độ tàn che và sinh khối để xây dựng khung tiêu chuẩn kỹ thuật phục hồi rừng tự nhiên cấp quốc gia.

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

  • Tại sao chỉ số đa dạng loài Shannon-Wiener ở rừng thứ sinh lại cao hơn rừng nguyên sinh? Hiện tượng này xuất hiện do xáo trộn rừng tạo ra các khoảng trống tán, kích thích các loài cây tiên phong ưa sáng như Nhò vàng, Vàng anh nảy mầm cùng sinh sống với các loài bản địa chịu bóng, khiến số lượng loài tạm thời tăng lên đạt H' = 0,964 tại ô thứ sinh so với 0,860 tại ô nguyên sinh.

  • Hàm tương quan theo cặp g(r) có ưu điểm gì vượt trội so với hàm Ripley's K? Hàm Ripley's K có tính chất tích lũy khoảng cách nên dễ che lấp các tương tác ở cự ly ngắn. Ngược lại, hàm g(r) là đạo hàm của hàm K, giúp phân tích mật độ chuẩn hóa tại từng khoảng cách cụ thể từ 0 đến 50 m một cách độc lập và chính xác.

  • Nguyên nhân nào khiến các loài cây gỗ ưu thế chủ yếu phân bố theo kiểu cụm ở khoảng cách dưới 35 m? Phân bố cụm phát sinh chủ yếu do đặc tính phát tán hạt hạn chế quanh gốc cây mẹ của các loài nhiệt đới. Khi điều kiện lập địa đồng nhất, cây con tập trung nảy mầm dày đặc xung quanh cây mẹ trong bán kính từ 1 đến 34 m.

  • Mối quan hệ không gian độc lập giữa các cặp loài ưu thế phản ánh điều gì? Kết quả kiểm định 100% các cặp loài có quan hệ độc lập (đường thực nghiệm nằm trong khoảng tin cậy 95%) phản ánh các loài cây trưởng thành đã đạt trạng thái cân bằng sinh thái, không xảy ra hiện tượng cạnh tranh loại trừ khốc liệt ở quy mô không gian nghiên cứu.

  • Đồ thị phân bố số cây theo cấp đường kính dạng chữ J ngược có ý nghĩa lâm học như thế nào? Đường cong chữ J ngược với số lượng lớn cây ở cấp kính nhỏ dưới 15 cm minh chứng lâm phần có khả năng tự tái sinh và duy trì cấu trúc quần xã ổn định, đảm bảo nguồn cây kế cận liên tục thay thế tầng cây già cỗi.

Kết luận

  • Luận văn đã định lượng hóa toàn diện cấu trúc mật độ (601 đến 1.006 cây/ha) và tiết diện ngang (19,78 đến 36,34 m2/ha) trên hai trạng thái rừng nguyên sinh và thứ sinh tại Cúc Phương.
  • Chứng minh tính đa dạng loài đạt mức cao với 106 loài tại rừng nguyên sinh và 82 đến 89 loài tại rừng thứ sinh, phản ánh tiềm năng phục hồi mạnh mẽ sau xáo trộn.
  • Khẳng định phân bố dạng cụm ở cự ly dưới 34 m là quy luật phổ biến do phát tán hạn chế, đồng thời các loài cây ưu thế duy trì mối quan hệ không gian độc lập.
  • Cung cấp cơ sở khoa học vững chắc để thực hiện lộ trình khoanh nuôi bảo vệ 1.500 ha rừng thứ sinh và số hóa dữ liệu quản lý rừng đặc dụng giai đoạn 2026-2030.
  • Kêu gọi các đơn vị quản lý tài nguyên và nhà khoa học đẩy mạnh ứng dụng mô hình toán sinh thái không gian vào thực tiễn quản lý rừng bền vững tại Việt Nam.