Tổng quan nghiên cứu

Tài nguyên rừng nhiệt đới Việt Nam đang chịu áp lực suy thoái kéo dài, kéo theo mất cân bằng sinh thái và giảm đa dạng sinh học. Trong bối cảnh đó, phục hồi rừng theo hướng "tiếp cận tự nhiên" đòi hỏi phải hiểu rõ cấu trúc rừng, đặc biệt là quy luật phân bố không gian của cây trên mặt phẳng nằm ngang — mảng nghiên cứu gần như còn bỏ ngỏ tại Việt Nam trước năm 2010.

Luận văn thạc sĩ Lâm học của tác giả Hoàng Thị Thu Hòa (Trường Đại học Lâm nghiệp, bảo vệ ngày 28/08/2010, hướng dẫn khoa học: TS. Trần Quang Bảo) tập trung nghiên cứu quy luật phân bố không gian của tầng cây cao trên kiểu rừng kín lá rộng thường xanh tại hai Vườn quốc gia Ba Vì và Cúc Phương. Đây là hai khu rừng đặc dụng được bảo vệ nghiêm ngặt, đại diện cho trạng thái rừng ít bị tác động nhất ở miền Bắc: VQG Ba Vì có tổng diện tích khoảng 7.377 ha với ba đỉnh cao nhất là đỉnh Vua (1.270 m), Tản Viên (1.227 m) và Ngọc Hoa (1.131 m); VQG Cúc Phương có tổng diện tích tự nhiên 52.852 ha, gồm 22.200 ha vùng lõi và hơn 30.000 ha vùng đệm.

Mục tiêu nghiên cứu gồm hai nội dung cụ thể: (1) xác định đặc điểm cấu trúc của kiểu rừng kín lá rộng thường xanh tại hai VQG; (2) tính toán các chỉ số phản ánh quy luật phân bố không gian theo chiều ngang của kiểu rừng này. Đối tượng là tầng cây cao thuộc trạng thái IIIA3 (Ba Vì) và trạng thái IV (Cúc Phương), số liệu kế thừa từ chu kỳ II chương trình theo dõi diễn biến tài nguyên rừng toàn quốc của Viện Điều tra Quy hoạch rừng năm 2007.

Ý nghĩa của đề tài nằm ở chỗ lượng hóa được "mẫu chuẩn" phân bố không gian của rừng tự nhiên, làm cơ sở để các nhà lâm học thiết kế biện pháp tác động đưa các lâm phần suy thoái tiệm cận mô hình chuẩ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 dựa trên hai nền tảng lý thuyết chính. Thứ nhất là học thuyết hệ sinh thái rừng của Tansley (1935) kết hợp học thuyết sinh địa quần học của Sucasev (1957), được Thái Văn Trừng (1978) phát triển thành hệ thống phân loại sinh thái phát sinh chia rừng Việt Nam thành 14 kiểu thảm thực vật — hệ thống mà luận văn lựa chọn vì tính chặt chẽ và khả năng áp dụng cho cả rừng nguyên sinh, thứ sinh lẫn rừng nhân tạo. Thứ hai là lý thuyết thống kê không gian điểm (spatial point pattern analysis), xuất phát từ phạm trù ba kiểu phân bố cá thể trên mặt đất rừng: cách đều (regular), ngẫu nhiên (Poisson/random) và cụm — lan truyền (contagious/aggregated).

Các khái niệm chính được vận hành hóa gồm: cấu trúc rừng (theo Phùng Ngọc Lan, 1986 — quy luật sắp xếp tổ hợp thành phần quần thể theo không gian và thời gian); tổ thành loài; mức độ quan trọng IV% theo Daniel Marmillod (loài có IV% ≥ 5% mới có ý nghĩa sinh thái, nhóm dưới 10 loài có tổng IV% ≥ 40% được coi là nhóm loài ưu thế); phân hóa giá trị cao — thấp; và tự tương quan không gian.

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

Dữ liệu được thu thập theo quy chuẩn ô định vị nghiên cứu sinh thái (OĐVNCST) của Viện Điều tra Quy hoạch rừng. Mỗi OĐVNCST rộng 100 ha (cạnh 1 km); lấy 1/4 về phía Đông Bắc (25 ha) làm ô điều tra cơ bản, trên đó thiết lập lưới ô vuông 50×50 m để phân chia lô trạng thái. Trong mỗi ô điều tra cơ bản mở 3 ô đo đếm đại diện, mỗi ô 1 ha (100×100 m), chia tiếp thành 25 phân ô 20×20 m. Sai số khép kín cho phép dưới 1/200. Hai ô định vị được sử dụng là số hiệu 335–40B (xã Ba Vì, huyện Ba Vì) và 343–43B (xã Cúc Phương, huyện Nho Quan, Ninh Bình).

Chọn mẫu theo phương pháp điển hình có định hướng: ô được đặt tại trạng thái rừng ít bị tác động nhất nhằm bảo đảm mẫu phản ánh trạng thái gần tự nhiên, phù hợp với mục tiêu tìm mẫu chuẩn thay vì suy rộng cho toàn bộ diện tích. Mỗi cây gỗ tầng cao được đo D1.3, Hvn, Hdc, đường kính tán theo bốn hướng Đông – Tây – Nam – Bắc, phẩm chất, tên loài (yêu cầu định danh ≥ 90% số cây) và đặc biệt là tọa độ gốc cây — dữ liệu then chốt cho phân tích không gian.

Xử lý số liệu gồm ba nhóm: (i) sàng lọc thô, loại bỏ mọi quan sát nằm ngoài khoảng X̄ ± 3S; (ii) mô phỏng quy luật tần số bằng hàm Weibull và hàm Khoảng cách, kiểm tra bằng tiêu chuẩn χ²; tương quan Hvn/D1.3 thử nghiệm 4 dạng phương trình (Logarit, Power, và các dạng đa thức); (iii) phân tích không gian bằng bốn chỉ số ANN, Ripley's K-function, General G (Getis-Ord) và Moran's I, cả ở mức toàn cục lẫn cục bộ (Gi*, Local I). Lý do chọn bộ chỉ số này là chúng bổ trợ nhau: ANN và K-function mô tả mạng hình vị trí, còn General G và Moran's I mô tả sự phân hóa giá trị thuộc tính (D1.3, Hvn) trong không gian.

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ấu trúc tổ thành và trạng thái. Kết quả phân loại xác định rừng Ba Vì thuộc trạng thái IIIA3 — rừng đã bị tác động vừa phải, đang phục hồi tốt, xu hướng diễn thế thứ sinh; rừng Cúc Phương thuộc trạng thái IV — rừng giàu, trữ lượng cao, diễn thế ổn định. Tổ thành tầng cây cao được xác định song song theo hai cách (tỷ lệ % số cây và IV%), cho thấy ưu thế rõ rệt của các họ nhiệt đới điển hình: họ Dẻ (Fagaceae) và họ Re (Lauraceae). Tại Cúc Phương, tầng vượt tán A1 gồm Chò xanh, Chò nhai, Chò chỉ, Chò đãi, Sấu; tầng ưu thế sinh thái A2–A3 gồm Vàng anh, Gội, Mạy tèo, Sồi, Mun núi đá, Dẻ gai.

Thứ hai, về quy luật tần số. Phân bố N/D1.3 được mô phỏng tốt nhất bằng hàm Khoảng cách, còn phân bố N/Hvn phù hợp với hàm Weibull. Kết quả này trùng khớp với kết luận của Trần Cẩm Tú (1999) cho rừng tự nhiên phục hồi sau khai thác tại Hương Sơn — Hà Tĩnh, củng cố nhận định rằng rừng tự nhiên hỗn giao khác tuổi có đường cong N/D dạng giảm với đỉnh ở cỡ kính nhỏ (thường quanh 12 cm), khác biệt căn bản với rừng trồng thuần loài đều tuổi.

Thứ ba, về tương quan sinh trưởng. Trong bốn dạng phương trình thử nghiệm, dạng Logarit Hvn = a + b·log(D1.3) được chọn cho phần lớn ô đo đếm, khẳng định quan hệ chiều cao — đường kính là quan hệ đường cong chặt, độ dốc giảm dần khi đường kính tăng.

Thứ tư, về mạng hình phân bố không gian. Đây là đóng góp cốt lõi. Bốn chỉ số được tính cho ba cấp đối tượng: tổng thể (TT), nhóm loài cây ưu thế (NLCUT) và loài cây ưu thế (LCUT). Với ANN, ngưỡng phán quyết là ANN < 1 và Z-score < −1,96 → phân bố cụm; ANN > 1 và Z-score > 1,96 → phân bố phân tán; khoảng giữa là ngẫu nhiên. K-function bổ sung chiều thông tin quan trọng: cùng một lâm phần có thể cụm ở khoảng cách ngắn nhưng ngẫu nhiên hoặc phân tán ở khoảng cách lớn hơn — điều mà chỉ số ANN đơn lẻ không phát hiện được.

Thảo luận kết quả

Xu hướng chung là mức độ tụ cụm giảm dần khi rừng tiến tới trạng thái thành thục. Điều này phù hợp với quan niệm của Nguyễn Hải Tuất và Ngô Kim Khôi (1990–1994): quá trình cạnh tranh không gian dinh dưỡng làm hình thái phân bố chuyển dịch từ cụm → ngẫu nhiên → cách đều. Trạng thái IV tại Cúc Phương, với diễn thế đã ổn định, do đó có xu hướng tiệm cận phân bố ngẫu nhiên rõ hơn trạng thái IIIA3 tại Ba Vì.

Kết quả cũng lý giải vì sao các nghiên cứu trước cho kết luận không thống nhất: Lê Sáu (1996) và Trần Cẩm Tú (1999) ghi nhận chủ yếu phân bố ngẫu nhiên và cách đều ở rừng thành thục, trong khi Đinh Quang Diệp (1994) thấy 10/12 ô nghiên cứu rừng khộp tái sinh có phân bố cụm. Sự khác biệt không phải mâu thuẫn mà phản ánh các giai đoạn diễn thế khác nhau.

Dữ liệu loại này được trình bày hiệu quả nhất qua ba dạng hình: bản đồ điểm vị trí cây trong từng ô đo đếm (thể hiện trực quan mức độ tụ cụm), đồ thị K-function với ba đường — đường kỳ vọng, đường quan sát và đường bao giới hạn (đọc trực tiếp ngưỡng khoảng cách chuyển đổi kiểu phân bố), và bản đồ nhiệt Gi* / Local I đánh dấu các điểm nóng cây lớn, điểm lạnh cây nhỏ. Bảng tổng hợp so sánh chéo ba cấp đối tượng × ba chỉ số cho phép kiểm chứng tính nhất quán giữa các phương pháp.

Ý nghĩa thực tiễn: khi mức độ cụm ở loài ưu thế cao hơn ở tổng thể, đó là dấu hiệu tái sinh theo lỗ trống hoặc theo cơ chế phát tán hạt gần cây mẹ — thông tin trực tiếp phục vụ thiết kế mật độ và cự ly trồng bổ sung.

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

Áp dụng bộ chỉ số không gian vào quy trình điều tra rừng quốc gia. Kiến nghị Viện Điều tra Quy hoạch rừng bổ sung bắt buộc trường tọa độ gốc cây (X, Y) vào biểu điều tra thống kê tầng cây gỗ cho toàn bộ hệ thống ô định vị nghiên cứu sinh thái, mục tiêu phủ 100% ô định vị trong chu kỳ điều tra kế tiếp. Đây là chi phí biên rất thấp nhưng mở khóa toàn bộ lớp phân tích không gian.

Xây dựng bộ mẫu chuẩn phân bố không gian theo kiểu rừng. Đề xuất Trường Đại học Lâm nghiệp phối hợp các VQG mở rộng nghiên cứu từ 2 lên tối thiểu 8–10 VQG đại diện cho các vùng sinh thái, trong 3–5 năm, nhằm hoàn thiện thư viện mẫu chuẩn cho ít nhất 5 kiểu thảm thực vật chính.

Thiết kế biện pháp lâm sinh dựa trên khoảng cách. Ban quản lý các VQG và công ty lâm nghiệp nên sử dụng ngưỡng khoảng cách xác định từ K-function để định cự ly tỉa thưa và trồng bổ sung, thay cho cự ly cơ giới đồng đều hiện hành. Với lâm phần đang phục hồi như trạng thái IIIA3, ưu tiên điều tiết tại các "điểm nóng" tụ cụm cây nhỏ do Gi* cục bộ chỉ ra.

Số hóa và tin học hóa dữ liệu điều tra. Khuyến nghị chuyển toàn bộ số liệu ô định vị sang định dạng GIS, tích hợp module tính ANN, K-function, General G, Moran's I, rút thời gian xử lý từ hàng tuần xuống dưới một ngày cho mỗi ô 1 ha.

Tăng cường tuyên truyền vùng đệm. Với VQG Cúc Phương, nơi vùng đệm rộng hơn 30.000 ha và đất canh tác nông nghiệp chiếm 836,09 ha trong vùng lõi mở rộng, cần chương trình vận động thường xuyên do chính quyền ba tỉnh Ninh Bình, Hòa Bình, Thanh Hóa chủ trì, nhằm giữ nguyên trạng diện tích rừng kín thường xanh 12.683,83 ha.

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

Học viên cao học và nghiên cứu sinh ngành Lâm học, Quản lý tài nguyên rừng. Luận văn cung cấp mẫu tham chiếu đầy đủ cho một đề tài định lượng: từ thiết kế ô định vị 100 ha, ô đo đếm 1 ha, tới quy trình sàng lọc X̄ ± 3S và kiểm định χ². Đặc biệt hữu ích khi xây dựng chương phương pháp cho đề tài có yếu tố phân tích không gian.

Cán bộ kỹ thuật Viện Điều tra Quy hoạch rừng và các chi cục kiểm lâm. Có thể dùng trực tiếp bộ biểu điều tra, quy cách phân ô 20×20 m và tiêu chuẩn phân loại trạng thái Loetschau đã hiệu chỉnh. Use case điển hình: chuẩn hóa quy trình thu thập tọa độ cây khi triển khai chu kỳ điều tra mới.

Ban quản lý vườn quốc gia và khu bảo tồn. Kết quả về tổ thành, mật độ, chỉ số đa dạng Simpson và Shannon–Wiener tại hai VQG là dữ liệu nền để xây dựng phương án quản lý rừng bền vững và đánh giá biến động sau mỗi chu kỳ 5 năm.

Giảng viên và người làm chính sách lâm nghiệp. Phần tổng quan hệ thống hóa hơn 50 công trình trong và ngoài nước về cấu trúc rừng, từ Richards (1952), Rollet (1971) đến Đồng Sĩ Hiền (1974), Bảo Huy (1993), là tài liệu giảng dạy sẵn dùng cho học phần Điều tra quy hoạch rừng, đồng thời làm luận cứ khoa học cho các quy định kỹ thuật về phục hồi rừng tự nhiên.

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

Phân bố không gian của cây rừng có mấy kiểu và phân biệt thế nào? Có ba kiểu: cách đều (regular), ngẫu nhiên (Poisson/random) và cụm (contagious/aggregated). Phân biệt bằng chỉ số ANN: giá trị nhỏ hơn 1 kèm Z-score dưới −1,96 là phân bố cụm; lớn hơn 1 kèm Z-score trên 1,96 là phân tán; nằm trong khoảng ±1,96 thì không đủ bằng chứng bác bỏ giả thuyết phân bố ngẫu nhiên.

Vì sao cần dùng nhiều chỉ số thay vì một? ANN chỉ xét khoảng cách tới cây gần nhất nên bỏ sót thông tin ở các cự ly xa hơn. K-function khắc phục bằng cách quét liên tục theo nhiều khoảng cách. Trong khi đó General G và Moran's I không xét vị trí mà xét giá trị thuộc tính — trả lời câu hỏi cây to có tụ lại một chỗ hay không. Bốn chỉ số cho bốn lát cắt bổ trợ.

Hàm nào mô phỏng tốt phân bố đường kính và chiều cao rừng tự nhiên? Với rừng tự nhiên hỗn giao khác tuổi, hàm Khoảng cách phù hợp cho N/D1.3 vì đường cong thường có đỉnh ngay sau cỡ kính nhỏ nhất, còn hàm Weibull mô phỏng tốt N/Hvn. Ngược lại, với rừng trồng thuần loài đều tuổi, hàm Weibull thường phù hợp cho cả hai phân bố.

Vì sao chọn Ba Vì và Cúc Phương làm địa bàn? Hai VQG này cùng thuộc kiểu rừng kín lá rộng thường xanh, đều được bảo vệ nghiêm ngặt, nhưng ở hai trạng thái diễn thế khác nhau: IIIA3 đang phục hồi và IV đã ổn định. Cặp so sánh này cho phép quan sát sự chuyển dịch kiểu phân bố theo tiến trình diễn thế trong cùng một kiểu rừng.

Kết quả có thể áp dụng cho rừng trồng không? Không trực tiếp. Rừng trồng có phân bố cách đều nhân tạo ngay từ đầu nên các chỉ số sẽ phản ánh thiết kế trồng chứ không phải quy luật tự nhiên. Giá trị ứng dụng nằm ở chiều ngược lại: dùng mẫu chuẩn tự nhiên làm đích để điều chỉnh rừng trồng và rừng phục hồi tiệm cận cấu trúc bền vững.

Kết luận

  • Luận văn là công trình tiên phong ở Việt Nam lượng hóa quy luật phân bố không gian theo chiều ngang của tầng cây cao, bằng bốn chỉ số thống kê không gian ANN, Ripley's K-function, General G và Moran's I, áp dụng đồng thời ở mức toàn cục và cục bộ.
  • Cấu trúc tần số rừng kín lá rộng thường xanh tại hai VQG tuân theo hàm Khoảng cách với N/D1.3 và hàm Weibull với N/Hvn; tương quan Hvn/D1.3 mô tả tốt nhất bằng dạng Logarit.
  • Mức độ tụ cụm giảm dần theo tiến trình diễn thế, từ trạng thái IIIA3 (Ba Vì) sang trạng thái IV (Cúc Phương), phù hợp cơ chế cạnh tranh không gian dinh dưỡng.
  • Đóng góp phương pháp luận quan trọng nhất là chứng minh giá trị của trường tọa độ gốc cây — dữ liệu chi phí thấp nhưng mở ra toàn bộ lớp phân tích không gian trước đây bỏ trống.
  • Hướng phát triển tiếp theo: mở rộng ra 8–10 VQG trong 3–5 năm, số hóa dữ liệu sang nền GIS, tiến tới bộ mẫu chuẩn cho các kiểu rừng chính.

Nếu bạn đang thực hiện đề tài về cấu trúc rừng, phân bố không gian hay điều tra quy hoạch rừng, luận văn này là tài liệu tham khảo nền tảng đáng đọc kỹ phần Chương 3 và Chương 4.