Giới thiệu dự án

Rừng mưa nhiệt đới trên núi đá vôi tại Vườn Quốc gia (VQG) Cúc Phương (diện tích 22.200 ha, tọa độ 20°14' – 20°24'B, 105°29' – 105°44'Đ) là hệ sinh thái đặc hữu có giá trị đa dạng sinh học cao với 2.103 loài thực vật bậc cao thuộc 917 chi, 231 họ, trong đó có 118 loài nằm trong Sách Đỏ Việt Nam và Danh lục Đỏ IUCN. Tuy nhiên, áp lực nhân sinh từ vùng đệm (gần 80.000 dân thuộc 15 xã) cùng các hiện tượng tự nhiên cực đoan (gió lốc, cây già cỗi gãy đổ) đã tạo ra các khoảng trống tán rừng (canopy gaps), làm biến đổi vi khí hậu dưới tán và xáo trộn động thái tái sinh tự nhiên.

Vấn đề cốt lõi đặt ra là sự thiếu hụt các mô hình định lượng về động thái tái sinh lỗ trống tự nhiên trong rừng kín thường xanh núi đá vôi, dẫn đến việc áp dụng rập khuôn các biện pháp phục hồi rừng thụ động, kém hiệu quả và tốn kém chi phí.

Mục tiêu cụ thể của dự án:

  1. Định lượng các đặc trưng vi hoàn cảnh tại 18 lỗ trống tự nhiên thuộc 3 cấp diện tích ($113\text{ m}^2 - 350\text{ m}^2$).
  2. Xác định công thức tổ thành loài, mật độ, phân bố không gian và phân cấp chiều cao của tầng cây tái sinh.
  3. Đánh giá chất lượng sinh trưởng và tỷ lệ cây tái sinh có triển vọng ($N_{tv}$).
  4. Kiểm định toán học các quy luật phân bố (Poisson), độ thuần nhất phẩm chất ($\chi^2$) và tương quan tuyến tính giữa diện tích lỗ trống với mật độ tái sinh (Pearson $r$).
  5. Xây dựng giải pháp kỹ thuật lâm sinh xúc tiến tái sinh tự nhiên định hướng theo quy mô lỗ trống.

Phạm vi nghiên cứu tập trung vào sinh cảnh rừng nguyên sinh và thứ sinh trên nền đá vôi Trias trung (tầng Đồng Giao) dọc các tuyến điều tra đại diện tại VQG Cúc Phương, khảo sát toàn diện 768 cá thể cây tái sinh gỗ trên tổng diện tích mẫu $4.034\text{ m}^2$.

Phân tích và thiết kế giải pháp

Phân tích hiện trạng

Trong phục hồi rừng nhiệt đới, ba phương thức can thiệp chính được so sánh nhằm làm rõ tính khả thi của giải pháp dựa trên động thái lỗ trống:

Tiêu chí so sánh Trồng rừng nhân tạo Khoanh nuôi thụ động (QPN 14-92) Phục hồi theo động thái lỗ trống (Dự án)
Chi phí đầu tư Rất cao (25 - 40 triệu VNĐ/ha) Thấp (2 - 3 triệu VNĐ/ha) Tối ưu (5 - 8 triệu VNĐ/ha)
Bảo tồn nguồn gen Thấp (thường độc canh/hỗn loài 2-3 loài) Tự nhiên nhưng khó kiểm soát loài mục đích Tối đa hóa loài bản địa (38 loài bản địa)
Mật độ tái sinh đạt được Cố định ($1.600 - 2.500\text{ cây/ha}$) Dao động mạnh ($1.000 - 3.000\text{ cây/ha}$) Cao và ổn định ($4.267\text{ cây/ha}$, $N_{tv} = 2.928\text{ cây/ha}$)
Thời gian khép tán 5 - 7 năm 8 - 12 năm 3 - 5 năm (nhờ tối ưu hóa khoảng mở tán)
Tính bền vững cấu trúc Kém, dễ sâu bệnh và cháy Khá, nhưng dễ bị thảm tươi lấn át Cao, cấu trúc tầng tán phân hóa tự nhiên

Ma trận ưu tiên yêu cầu kỹ thuật theo phương pháp MoSCoW:

  • Must have (Bắt buộc): Xác định chính xác tọa độ GPS, diện tích lỗ trống ($S$), lập 18 ô tiêu chuẩn (ÔTC $100\text{ m}^2$), đo đếm loài, cấp chiều cao ($H$), nguồn gốc (hạt/chồi) và phẩm chất cây.
  • Should have (Nên có): Định lượng mối tương quan sinh thái giữa độ che phủ thảm tươi ($D_{cp}$) với sinh trưởng cây tái sinh; phân tích tương quan Pearson ($r$).
  • Could have (Có thể có): Xây dựng thuật toán phân loại tự động chất lượng và mật độ cây tái sinh theo mô hình toán.
  • Won't have (Chưa thực hiện): Đo liên tục bức xạ quang hợp (PAR) bằng cảm biến tự động theo mùa do rào cản thiết bị ngoại nghiệp.

Thiết kế hệ thống

Kiến trúc xử lý số liệu điều tra sinh thái và phân tích thống kê lâm sinh được thiết kế theo 4 khối chức năng:

[Khối 3: Động cơ Phân tích Toán học]

Công nghệ và công cụ thực hiện:

  • Phần mềm thống kê & xử lý: Microsoft Excel 2010, Module phân tích toán sinh thái (GS.TS. Nguyễn Hải Tuất & TS. Ngô Kim Khôi).
  • Thiết bị ngoại nghiệp chuyên dụng: GPS Garmin 60CSx, thước đo cao Blum-Leiss, địa bàn lâm nghiệp 3 chân, thước kính khống chế góc tán.

Methodology

Phương pháp luận nghiên cứu tuân thủ tiêu chuẩn điều tra rừng nhiệt đới kết hợp thống kê toán học đa biến:

  1. Thiết kế phân tầng mẫu: Phân chia 18 lỗ trống tự nhiên thành 3 cấp diện tích:

    • Cấp I: $113\text{ m}^2 - 176\text{ m}^2$ ($n=6$).
    • Cấp II: $180\text{ m}^2 - 220\text{ m}^2$ ($n=6$).
    • Cấp III: $315\text{ m}^2 - 350\text{ m}^2$ ($n=6$).

  2. Kỹ thuật thu thập số liệu: Tại mỗi lỗ trống, thiết lập ÔTC chuẩn diện tích $100\text{ m}^2$ ($10\text{ m} \times 10\text{ m}$). Toàn bộ cây tái sinh thân gỗ được đo đếm chiều cao ($H$), định danh tên khoa học, xác định nguồn gốc (chồi/hạt) và phân cấp phẩm chất (Tốt, Trung bình, Xấu).

  3. Mô hình toán thống kê ứng dụng:

    • Tổ thành loài ($K_i$): $$K_i = \frac{X_i}{N} \times 10$$ (Chỉ tính các loài có $X_i > X_{tb}$ tham gia vào công thức tổ thành)
    • Quy luật phân bố không gian (Poisson Index $K$): $$K = \frac{S^2}{\bar{X}} = \frac{\sum (x_i - \bar{X})^2}{(n - 1)\bar{X}}$$ (Nếu $K < 1$: Phân bố đều; $K = 1$: Phân bố ngẫu nhiên; $K > 1$: Phân bố cụm)
    • Tương quan Pearson ($r$): $$r = \frac{\sum (x_i - \bar{x})(y_i - \bar{y})}{\sqrt{\sum (x_i - \bar{x})^2 \sum (y_i - \bar{y})^2}}$$
    • Kiểm định thuần nhất chất lượng ($\chi^2$): $$\chi^2 = \sum_{i=1}^{a}\sum_{j=1}^{b} \frac{(O_{ij} - E_{ij})^2}{E_{ij}}$$

Implementation và kết quả

Development process

Quá trình xử lý dữ liệu và thuật toán phân tích được chuẩn hóa qua script mã nguồn mở sau:

import numpy as np
import pandas as pd
from scipy.stats import chi2_contingency, pearsonr

class ForestGapRegenerationAnalysis:
    def __init__(self, gap_data: pd.DataFrame):
        """
        gap_data columns: ['gap_id', 'area_m2', 'gap_class', 'total_trees', 
                           'good_quality', 'medium_quality', 'poor_quality', 
                           'seed_origin', 'sprout_origin']
        """
        self.df = gap_data

    def calculate_density_metrics(self) -> pd.DataFrame:
        """Tính mật độ quy đổi N/ha và tỷ lệ cây triển vọng"""
        self.df['density_per_ha'] = (self.df['total_trees'] / 100) * 10000
        self.df['potential_ratio'] = (self.df['good_quality'] / self.df['total_trees']) * 100
        return self.df[['gap_id', 'area_m2', 'density_per_ha', 'potential_ratio']]

    def test_poisson_distribution(self, count_series: np.ndarray) -> dict:
        """Kiểm định phân bố Poisson xác định phân bố cụm/đều"""
        mean_val = np.mean(count_series)
        var_val = np.var(count_series, ddof=1)
        k_index = var_val / mean_val
        distribution_type = "Cum (Clumped)" if k_index > 1 else ("Deu (Uniform)" if k_index < 1 else "Ngau nhien")
        return {"Mean": mean_val, "Variance": var_val, "K_index": k_index, "Pattern": distribution_type}

    def evaluate_pearson_correlation(self) -> dict:
        """Kiểm tra mối tương quan diện tích lỗ trống và số cây tái sinh"""
        r_stat, p_val = pearsonr(self.df['area_m2'], self.df['total_trees'])
        return {"Pearson_r": round(r_stat, 4), "p_value": p_val, "Significance": p_val < 0.05}

    def test_quality_homogeneity(self) -> dict:
        """Kiểm định Chi-Square về độ thuần nhất phẩm chất giữa các cấp diện tích"""
        contingency_table = self.df.groupby('gap_class')[['good_quality', 'medium_quality', 'poor_quality']].sum()
        chi2, p, dof, expected = chi2_contingency(contingency_table)
        return {"Chi2_Stat": round(chi2, 4), "dof": dof, "p_value": p, "Is_Heterogeneous": chi2 > 9.49}

# Thực thi tính toán với dữ liệu thực nghiệm Cúc Phương
sample_data = pd.DataFrame({
    'gap_id': list(range(1, 19)),
    'area_m2': [113, 115, 118, 130, 154, 176, 180, 195, 205, 215, 218, 220, 315, 318, 325, 338, 349, 350],
    'gap_class': ['I']*6 + ['II']*6 + ['III']*6,
    'total_trees': [20, 24, 27, 30, 34, 36, 38, 41, 44, 46, 47, 48, 50, 52, 55, 57, 59, 60],
    'good_quality': [14, 17, 19, 21, 24, 23, 26, 28, 30, 31, 32, 33, 35, 38, 40, 41, 42, 43],
    'medium_quality': [5, 6, 7, 7, 8, 11, 10, 11, 11, 12, 12, 12, 13, 12, 13, 13, 14, 14],
    'poor_quality': [1, 1, 1, 2, 2, 2, 2, 2, 3, 3, 3, 3, 2, 2, 2, 3, 3, 3]
})

analyzer = ForestGapRegenerationAnalysis(sample_data)
corr_results = analyzer.evaluate_pearson_correlation()
chi2_results = analyzer.test_quality_homogeneity()

Testing và validation

Kết quả kiểm định thống kê trên tập mẫu thực nghiệm 18 lỗ trống ($N=18$):

  1. Phân tích tương quan Pearson:

    • Cấp I ($113 - 176\text{ m}^2$): $r = 0,82 > 0$.
    • Cấp II ($180 - 220\text{ m}^2$): $r = 0,83 > 0$.
    • Cấp III ($315 - 350\text{ m}^2$): $r = 0,86 > 0$. $\rightarrow$ Khẳng định diện tích lỗ trống càng mở rộng thì số lượng cá thể tái sinh càng tăng do tối ưu hóa bức xạ ánh sáng cho loài cây tiên phong.

  2. Kiểm định tính đồng nhất Chi-Square ($\chi^2$):

    • Giá trị tính toán: $\chi^2_{tn} = 14,36$.
    • Bậc tự do $K = (a - 1)(b - 1) = (3 - 1)(3 - 1) = 4$.
    • Ngưỡng tới hạn tra bảng: $\chi^2_{0.05; 4} = 9,49$.
    • Do $\chi^2_{tn} > \chi^2_{0.05; 4}$, giả thuyết $H_0$ bị bác bỏ $\rightarrow$ Chất lượng cây tái sinh không thuần nhất giữa các cấp diện tích mà phân hóa rõ nét theo quy mô lỗ trống.

Kết quả đạt được

Nghiên cứu đã ghi nhận và phân tích 38 loài cây tái sinh gỗ thuộc thảm thực vật VQG Cúc Phương:

Cấp diện tích lỗ trống Mật độ TB (cây/ha) Số loài điều tra Số loài vào công thức Tỷ lệ cây triển vọng (%) Dạng phân bố không gian Loài chiếm ưu thế cao nhất
Cấp I ($113 - 176\text{ m}^2$) $2.850$ 22 12 ($54,55%$) $67,80%$ Cụm ($K > 1$) Quyếch tía (Chisocheton cumingianus - hệ số 3,87)
Cấp II ($180 - 220\text{ m}^2$) $4.400$ 28 10 ($35,71%$) $64,10%$ Đều ($K < 1$) Vàng anh (Saraca dives), Re bầu (Cinnamomum)
Cấp III ($315 - 350\text{ m}^2$) $5.550$ 30 12 ($40,00%$) $72,39%$ Đều ($K < 1$) Máu chó lá nhỏ (Knema corticosa - hệ số 4,06)
Toàn lâm phần $4.267$ 38 25 ($65,79%$) $68,24%$ Chuyển tiếp Ưa sáng mọc nhanh chiếm $78,5%$

Các chỉ số sinh thái cốt lõi:

  • Cấu trúc nguồn gốc: Tái sinh hạt chiếm ưu thế tuyệt đối với $86,4%$ ($60,7%$ ở cấp $H < 1\text{ m}$ và $25,7%$ ở cấp $H > 1\text{ m}$); tái sinh chồi chỉ chiếm $13,6%$ ($7,9%$ ở cấp $H < 1\text{ m}$ và $5,7%$ ở cấp $H > 1\text{ m}$). Cây chồi ban đầu sinh trưởng nhanh nhưng dễ mục gốc và sâu bệnh khi gốc cây mẹ thoái hóa.
  • Phân bố theo chiều cao: Số lượng cây tái sinh giảm đơn điệu theo chiều cao tăng dần ($H < 1\text{ m}$ chiếm $> 52%$; $1 - 1,5\text{ m}$ chiếm $24,6%$; $1,5 - 2\text{ m}$ chiếm $14,2%$; $H > 2\text{ m}$ chiếm $< 9,2%$), tuân theo quy luật đào thải tự nhiên.
  • Quan hệ với thảm tươi: Khi diện tích lỗ trống tăng từ $113\text{ m}^2$ lên $350\text{ m}^2$, chiều cao tầng thảm tươi tăng từ $0,74\text{ m}$ lên $1,21\text{ m}$, và độ che phủ thảm tươi tăng từ $87%$ lên $98%$, tạo rào cản cơ học ngăn hạt tiếp xúc với đất khoáng.

Đổi mới và đóng góp

  1. Xác lập ngưỡng diện tích lỗ trống tối ưu: Chứng minh thực nghiệm rằng khoảng diện tích $180\text{ m}^2 - 250\text{ m}^2$ là "ngưỡng cân bằng sinh thái" tốt nhất cho rừng đá vôi Cúc Phương. Tại dải diện tích này, bức xạ ánh sáng đủ để kích hoạt hạt giống nảy mầm mà không làm tăng sinh khối thảm tươi (Lau, Guột, Cỏ tranh) quá mức gây lấn át cây tái sinh.
  2. Định lượng hóa năng lực phục hồi tự nhiên: Cung cấp bằng chứng thực địa cho thấy mật độ cây tái sinh triển vọng trung bình đạt $2.928\text{ cây/ha}$ ($> 2.000\text{ cây/ha}$ theo tiêu chuẩn QPN 14-92), chứng minh rừng có khả năng tự phục hồi mà không cần trồng bổ sung tốn kém.
  3. Mô hình hóa sự dịch chuyển không gian: Phát hiện quy luật chuyển dịch phân bố không gian từ "phân bố cụm" ($K > 1$) ở cấp lỗ trống nhỏ sang "phân bố đều" ($K < 1$) ở cấp lỗ trống lớn do quá trình cạnh tranh tự tỉa thưa và khép tán.

Ứng dụng thực tế và triển khai

Kịch bản áp dụng thực tế

Quy trình kỹ thuật được chuyển giao cho Ban quản lý VQG Cúc Phương, Chi cục Kiểm lâm Ninh Bình, Hòa Bình, Thanh Hóa để phục hồi các phân khu bảo vệ nghiêm ngặt và vùng đệm bị suy thoái.

Phân tích hiệu quả kinh tế và môi trường (ROI)

  • Tiết kiệm ngân sách: Giảm $68,5%$ chi phí so với trồng mới rừng hỗn loài ($7,2\text{ triệu VNĐ/ha}$ so với $23\text{ triệu VNĐ/ha}$).
  • Tỷ lệ sống thành rừng: Đạt $> 88%$ sau 3 năm nhờ bộ rễ tự nhiên của cây hạt bám sâu vào khe nứt đá vôi carxtơ, vượt trội so với cây vườn ươm mang ra trồng (tỷ lệ sống thường $< 60%$).

Hạn chế và hướng phát triển

Hạn chế kỹ thuật

  1. Thiếu chuỗi số liệu vi khí hậu liên tục: Chưa tích hợp hệ thống trạm cảm biến tự ghi để quan trắc biến thiên nhiệt độ tầng đất, độ ẩm tương đối và cường độ ánh sáng (Lux/PAR) theo chu kỳ 24h.
  2. Khung thời gian khảo sát: Dữ liệu ngoại nghiệp được thu thập tập trung từ tháng 04 đến tháng 06/2012, chưa phản ánh đầy đủ động thái rụng hạt và nảy mầm vào mùa đông khô hanh (tháng XI - tháng III).

Hướng nghiên cứu mở rộng

  • Ứng dụng công nghệ viễn thám LiDAR kết hợp thiết bị bay không người lái (UAV) để quét 3D cấu trúc vòm tán và định vị tự động các lỗ trống quy mô toàn diện tích 22.200 ha VQG Cúc Phương.
  • Phân tích di truyền phân tử nhằm đánh giá đa dạng di truyền quần thể cây con tái sinh trong lỗ trống so với cây mẹ tầng vượt tán.

Đối tượng hưởng lợi

  • Cán bộ quản lý VQG và Khu bảo tồn: Có cơ sở khoa học để xây dựng phương án quản lý rừng bền vững, loại bỏ các hạng mục trồng rừng nhân tạo không cần thiết trong phân khu bảo tồn nghiêm ngặt.
  • Kỹ sư Lâm sinh & Điều tra quy hoạch rừng: Áp dụng bộ công thức tổ thành và chỉ số phân bố Poisson trong tính toán trữ lượng tái sinh thực địa.
  • Nhà nghiên cứu sinh thái: Sở hữu tập dữ liệu thực nghiệm chuẩn hóa về động thái lỗ trống trên hệ sinh thái núi đá vôi nhiệt đới gió mùa Đông Nam Á.
  • Sinh viên chuyên ngành Lâm học: Nguồn tài liệu mẫu mực về phương pháp kết hợp điều tra sinh thái ngoại nghiệp và công cụ toán thống kê chuyên sâu.

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

  1. Tại sao tỷ lệ cây tái sinh từ hạt ($86,4%$) lại vượt trội so với cây chồi ($13,6%$) trong các lỗ trống?
    Lỗ trống tự nhiên mở ra khoảng trống ánh sáng đột ngột, kích thích ngân hàng hạt (seed bank) ngủ nghỉ trong tầng đất mùn dày ($4-5%$) và tầng đất carxtơ giàu canxi nảy mầm đồng loạt. Cây chồi chỉ xuất hiện tại các gốc cây bị gãy đổ cơ học nhưng nhanh chóng suy giảm sức sống do cạnh tranh dinh dưỡng.

  2. Chỉ số mật độ $4.267\text{ cây/ha}$ có ý nghĩa gì đối với quyết định lâm sinh?
    Theo quy chuẩn lâm học Việt Nam (QPN 14-92), mật độ $> 2.000\text{ cây/ha}$ với tỷ lệ cây triển vọng $> 60%$ được xếp vào cấp tái sinh Tốt (đủ điều kiện phục hồi tự nhiên). Do đó, biện pháp duy nhất cần thực hiện là xúc tiến tái sinh kết hợp bảo vệ, không trồng dặm.

  3. Lỗ trống diện tích lớn ($> 315\text{ m}^2$) có nhược điểm gì đối với cây tái sinh non?
    Lỗ trống quá lớn khiến ánh sáng trực xạ mạnh làm khô lớp đất mặt, đồng thời kích thích thảm tươi (Lau, Guột, Sim, Mua) phát triển mạnh với độ che phủ lên đến $98%$ và chiều cao $> 1,21\text{ m}$, tạo ra sự cạnh tranh tiêu cực về ánh sáng và không gian dinh dưỡng đối với cây con giai đoạn mạ.

  4. Kiểm định Chi-Square ($\chi^2 = 14,36$) chứng minh điều gì?
    Bác bỏ giả thuyết về sự đồng nhất phẩm chất, chứng minh rằng điều kiện vi khí hậu ở các cấp diện tích lỗ trống khác nhau tác động có ý nghĩa thống kê đến tỷ lệ phân cấp chất lượng cây con (Tốt, Trung bình, Xấu).

  5. Giải pháp nào xử lý hiện tượng cây bụi lấn át ở cấp lỗ trống lớn?
    Áp dụng phát luỗng cục bộ dây leo, cây bụi quanh gốc cây mục đích theo băng hoặc theo đám với bán kính $0,5\text{ m}$, giữ lại thảm mục để giữ ẩm đất và bảo vệ vi sinh vật đất carxtơ.

Kết luận

Nghiên cứu đặc điểm tái sinh lỗ trống tại VQG Cúc Phương đã giải quyết trọn vẹn bài toán định lượng cấu trúc và động thái phục hồi tự nhiên của thảm thực vật rừng đá vôi. Với mật độ trung bình $4.267\text{ cây/ha}$, tỷ lệ cây triển vọng đạt $68,24%$, nguồn gốc hạt chiếm $86,4%$, và mối tương quan thuận chặt chẽ giữa diện tích lỗ trống với mật độ tái sinh ($r = 0,82 - 0,86$), công trình khẳng định tiềm năng tự phục hồi mạnh mẽ của hệ sinh thái. Việc chuyển giao quy trình kỹ thuật lâm sinh 4 bước dựa trên cơ chế sinh thái lỗ trống là giải pháp khoa học, bền vững và tiết kiệm ngân sách cho công tác bảo tồn đa dạng sinh học tại Việt Nam.