Giới thiệu dự án

Rừng nhiệt đới tại Việt Nam là hệ sinh thái có tính đa dạng sinh học cao nhưng đang đối mặt với nguy cơ suy thoái nghiêm trọng. Theo thống kê lịch sử lâm nghiệp, trong hơn 50 năm qua, Việt Nam đã mất hơn 5 triệu hecta rừng tự nhiên (từ 14,3 triệu ha năm 1943 xuống còn 9,5 triệu ha năm 1993, trung bình mất 100.000 ha/năm). Hiện nay, diện tích rừng gỗ tự nhiên giàu và trung bình chỉ còn khoảng 1,4 triệu ha (chiếm xấp xỉ 13% tổng diện tích rừng), tập trung chủ yếu ở các vùng núi cao, hiểm trở.

Vườn Quốc gia (VQG) Bến En (tỉnh Thanh Hóa) có diện tích vùng lõi 16.634 ha và vùng đệm 12.000 ha, là nơi lưu trữ 1.389 loài thực vật bậc cao có mạch thuộc 650 chi, 173 họ. Nơi đây sở hữu 29 loài nằm trong Danh lục đỏ IUCN (2006) và 42 loài trong Sách đỏ Việt Nam (2007). Tuy nhiên, áp lực khai thác gỗ trái phép phục vụ thương mại đã làm suy giảm nghiêm trọng các quần thể gỗ quý nhóm "tứ thiết" và gỗ quý bản địa. Điển hình, từ năm 2007 đến 2011, đã có ít nhất 112 cây Lim xanh bị chặt hạ trái phép với khối lượng thất thoát hơn $110,813\text{ m}^3$, bên cạnh các loài Chò chỉ và Sao hải nam đang đứng trước nguy cơ tuyệt chủng cục bộ.

Vấn đề nghiên cứu và Pain Points

  • Thiếu hụt dữ liệu định lượng sinh thái: Các chương trình phục hồi rừng trước đây thiếu tham số cấu trúc vi mô như không gian dinh dưỡng ($S_{dd}$), cự ly sinh thái tối ưu và tổ hợp loài đi kèm tự nhiên.
  • Tái sinh tự nhiên bị nghẽn: Cây mạ và cây con dưới tán rừng tự nhiên gặp tỷ lệ tử vong cao do chu kỳ sai quả cách quãng (2–5 năm/lần) và thiếu hụt tán sáng phù hợp theo từng giai đoạn tuổi.
  • Xung đột bảo tồn và sinh kế: Vùng đệm VQG Bến En gồm 16 xã, 2 thị trấn với 8.023 hộ dân ($41.500$ nhân khẩu, mật độ 57 người/$\text{km}^2$), dẫn đến áp lực lớn về đất đai và lâm sản.

Mục tiêu dự án

  1. Xác định chính xác tọa độ, phân bố không gian và mật độ quần thể của 3 loài cây quý: Lim xanh (Erythrophloeum fordii Oliv.), Chò chỉ (Parashorea chinensis Wang Hsie), và Sao hải nam (Hopea hainanensis Merr. et Chun) trên tổng chiều dài 10,6 km tuyến điều tra tại tiểu khu 616 và 619 (khu vực Sông Chàng).
  2. Định lượng các chỉ tiêu sinh thái thân cây (đường kính ngang ngực $D_{1.3}$, đường kính tán $D_t$, chiều cao vút ngọn $H_{vn}$, chiều cao dưới cành $H_{dc}$).
  3. Thiết lập công thức tổ thành loài (CTTT) tầng cây gỗ và xác định diện tích dinh dưỡng ($S_{dd}$) bằng phương pháp ô tiêu chuẩn 7 cây Thomasius.
  4. Đánh giá động thái tái sinh tự nhiên dưới tán rừng và xung quanh gốc cây mẹ ($1 \times D_t$).
  5. Xây dựng gói giải pháp lâm sinh học tích hợp (Bảo tồn nguyên vị In-situ, Bảo tồn chuyển vị Ex-situ, kỹ thuật gieo ươm).

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

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

Phương pháp tiếp cận Ưu điểm Hạn chế Khả năng ứng dụng tại Bến En
Bảo vệ nghiêm ngặt thuần túy Giữ nguyên hiện trạng rừng nguyên sinh, không can thiệp nhân tạo. Bị động trước áp lực lâm tặc; cây tái sinh thiếu ánh sáng để vượt tán, tỷ lệ chết tự nhiên cao. Kém hiệu quả nếu không kết hợp hỗ trợ kỹ thuật lâm sinh.
Trồng rừng thuần loài nhân tạo Mật độ cao, dễ chăm sóc, tăng nhanh độ che phủ ngắn hạn. Dễ bùng phát sâu bệnh hại, triệt tiêu đa dạng sinh học, không phù hợp tập tính loài tầng vượt tán. Không phù hợp cho Lim xanh và Sao hải nam do nhu cầu che bóng giai đoạn đầu.
Lâm sinh sinh thái đa tầng (Đề tài ứng dụng) Tái lập cấu trúc rừng tự nhiên dựa trên công thức tổ thành ($CTTT$) và không gian dinh dưỡng ($S_{dd}$). Yêu cầu điều tra thực địa công phu, dữ liệu xử lý thống kê toán học phức tạp. Tối ưu nhất: Đảm bảo tương thích sinh thái và khả năng sống sót $>85%$.

Yêu cầu nghiên cứu theo mô hình MoSCoW

  • Must have: Mật độ cây cá thể (cây/km, cây/ha), công thức tổ thành tầng cây cao ($CTTT$), không gian dinh dưỡng cây mẹ ($S_{dd}$), đánh giá tầng tái sinh có triển vọng ($H > 1\text{ m}$).
  • Should have: Tương quan địa hình - thổ nhưỡng (núi đất feralit đỏ vàng, núi đá vôi, thung lũng), vật hậu học chu kỳ ra hoa kết quả.
  • Could have: Kỹ thuật giâm hom vô tính kết hợp chất điều hòa sinh trưởng (IBA/NAA) cho nguồn gen khó gieo hạt.
  • Won't have: Thử nghiệm can thiệp chỉnh sửa gen phân tử hoặc giải trình tự toàn bộ genome quy mô lớn.

Thiết kế hệ thống điều tra và thu thập dữ liệu

Nghiên cứu ứng dụng cấu trúc phân tầng thu thập dữ liệu thực địa kết hợp xử lý thống kê sinh trắc học rừng:

[Hệ thống thực địa]

Thiết bị và công cụ tiêu chuẩn

  • Định vị vệ tinh: Máy định vị GPS Garmin chuyên dụng (sai số $< 5\text{ m}$).
  • Đo trắc lượng hình thái: Thước đo cao Blume-Leiss / Thước ngắm laser, Thước kẹp kính cơ học chính xác (Caliper), Thước dây bọc sợi thủy tinh $50\text{ m}$.
  • Thiết bị đo phương vị: Địa bàn địa chất Harbin DQL-8.
  • Phần mềm xử lý số liệu: Microsoft Excel kết hợp các module phân tích phương sai và hồi quy sinh thái.

Methodology

Quy trình nghiên cứu áp dụng phương pháp luận Lâm sinh học định lượng (Quantitative Silviculture):

  1. Phương pháp kế thừa: Tổng hợp bản đồ phân chia tiểu khu 616, 619; số liệu khí tượng - thủy văn trạm Như Xuân (nhiệt độ trung bình $23,3^\circ\text{C}$, lượng mưa $1.790\text{ mm/năm}$, độ ẩm $85%$).
  2. Phương pháp điều tra tuyến sơ thám & chi tiết: Ghi nhận toàn bộ cá thể bắt gặp, xác định tọa độ, độ cao tuyệt đối ($m$), trạng thái lập địa (núi đất, núi đá, thung lũng).
  3. Phương pháp ô tiêu chuẩn 7 cây của Thomasius: Xác định quan hệ cạnh tranh và không gian chiếm lĩnh giữa cây mục tiêu ($O$) và 6 cây lân cận gần nhất ($i = 1 \dots 6$).
  4. Phương pháp bình quân gia quyền: Tính toán trị số trung bình mẫu sinh trắc: $$\bar{X} = \frac{\sum f_i \cdot x_i}{n} \quad (\text{với } n > 30)$$

Implementation và kết quả

Development Process & Key Algorithms

Hệ thống xử lý tính toán sinh trắc học rừng được chuẩn hóa theo mã giải thuật toán học dưới đây nhằm xác định công thức tổ thành ($CTTT$) và diện tích không gian dinh dưỡng ($S_{dd}$):

import math
from typing import List, Dict, Tuple

class ForestBiometricsAnalyzer:
    @staticmethod
    def calculate_trophic_area(r_max: float) -> float:
        """
        Tính diện tích không gian dinh dưỡng của cây tâm theo Thomasius:
        S_dd = (pi * R_0^2) / 6
        """
        return (math.pi * (r_max ** 2)) / 6.0

    @staticmethod
    def compute_composition_formula(species_counts: Dict[str, int]) -> Tuple[Dict[str, float], str]:
        """
        Xác định hệ số tổ thành (k_i) và công thức tổ thành (CTTT)
        k_i = (n_i / N) * 10
        """
        total_trees = sum(species_counts.values())
        num_species = len(species_counts)
        mean_trees_per_species = total_trees / num_species
        
        dominant_species = {}
        other_trees = 0
        
        for sp, count in species_counts.items():
            if count >= mean_trees_per_species:
                k_i = round((count / total_trees) * 10, 2)
                dominant_species[sp] = k_i
            else:
                other_trees += count
                
        k_other = round((other_trees / total_trees) * 10, 2)
        
        # Sort dominant species descending
        sorted_species = sorted(dominant_species.items(), key=lambda x: x[1], reverse=True)
        cttt_str = " + ".join([f"{k}{sp}" for sp, k in sorted_species])
        if k_other > 0:
            cttt_str += f" + {k_other}Lkh"
            
        return dominant_species, cttt_str

# Test data từ lâm phần Lim xanh Bến En
analyzer = ForestBiometricsAnalyzer()
s_dd_lim = analyzer.calculate_trophic_area(r_max=7.0)
print(f"Diện tích dinh dưỡng Lim xanh: {s_dd_lim:.2f} m² (Thực tế mẫu: 35.74 m²)")

Testing và validation

Nghiên cứu đã hoàn thành $100%$ khối lượng ngoại nghiệp trên $10,6\text{ km}$ tuyến điều tra, ghi nhận $67$ cá thể trưởng thành của 3 loài mục tiêu, tiến hành giải đoán cấu trúc trên $6.000\text{ m}^2$ ô tiêu chuẩn và $480\text{ m}^2$ ô dạng bản tái sinh.

                      TỔNG HỢP MẬT ĐỘ PHÂN BỐ THEO ĐỊA HÌNH
        Lim (Đá)       Lim (Đất)     Chò chỉ        Sao hải nam

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

1. Mật độ và phân bố không gian

  • Lim xanh (Erythrophloeum fordii): Tổng số $48$ cây ($4,55\text{ cây/km}$). Phân bố áp đảo trên núi đất ($45$ cây, chiếm $93,75%$), chỉ có $3$ cây trên núi đá ($6,25%$). Tập trung nhiều nhất tại khu vực Rừng giống thuộc tiểu khu 619.
  • Chò chỉ (Parashorea chinensis): Tổng số $11$ cây ($1,04\text{ cây/km}$). Thích nghi đa dạng trên cả 3 dạng sinh cảnh: núi đất ($7$ cây, $63,64%$), thung lũng ($3$ cây, $27,27%$), núi đá ($1$ cây, $9,09%$).
  • Sao hải nam (Hopea hainanensis): Tổng số $8$ cây ($0,75\text{ cây/km}$). Phân bố đặc hữu hẹp, $100%$ bắt gặp tại sinh cảnh núi đất rừng giống tiểu khu 619 ($1,64\text{ cây/km}$ trên tuyến 01).

2. Chỉ tiêu sinh trắc học thân cây trưởng thành ($n=10$ cây/loài)

Chỉ tiêu sinh trắc Lim xanh (E. fordii) Chò chỉ (P. chinensis) Sao hải nam (H. hainanensis)
$D_{1.3}$ trung bình (cm) $21,4 \pm 6,8$ (Min: $14$ - Max: $36$) $22,4 \pm 11,2$ (Min: $8$ - Max: $50$) $28,3 \pm 7,4$ (Min: $15$ - Max: $40$)
$D_t$ trung bình (m) $5,6 \pm 1,4$ $3,33 \pm 0,9$ $3,92 \pm 0,8$
$H_{vn}$ trung bình (m) $12,85 \pm 2,6$ (Max: $18,0$) $17,45 \pm 4,2$ (Max: $27,5$) $16,95 \pm 3,8$ (Max: $23,0$)
$H_{dc}$ trung bình (m) $6,20 \pm 1,8$ $8,38 \pm 2,7$ $11,45 \pm 3,1$
Không gian $S_{dd}$ ($m^2$) $35,74$ (Cự ly max: $7\text{ m}$) $42,54$ (Cự ly max: $8\text{ m}$) $32,59$ (Cự ly max: $8\text{ m}$)

3. Công thức tổ thành tầng cây gỗ và cây đi kèm

  • Lâm phần Lim xanh ($N=111$ cây, 46 loài): $$CTTT = 2,52\text{Lxa} + 1,08\text{Sau} + 0,54\text{Chc} + 0,54\text{Nhl} - 0,45\text{Got} - 0,45\text{Trc} - 0,36\text{Thmg} - 0,27\text{De} + 3,87\text{Lkh}$$ (Ưu hợp chính: Lim xanh - Sau sau Liquidambar formosana chiếm 36% tổ thành lâm phần).
  • Lâm phần Chò chỉ ($N=70$ cây, 31 loài): $$CTTT = 1,43\text{Dh} + 1,29\text{Choc} + 0,86\text{Cal} + 0,71\text{Dpt} + 0,57\text{Or} - 0,43\text{Lmc} + 3,86\text{Lkh}$$ (Đứng thứ hai sau Đỉnh hương Dysoxylum cauliflorum).
  • Lâm phần Sao hải nam ($N=67$ cây, 28 loài): $$CTTT = 1,49\text{Shn} + 1,34\text{Lxa} + 0,75\text{Chc} + 0,60\text{Lxe} - 0,45\text{Sl} - 0,45\text{Nga} - 0,45\text{Trc} + 4,48\text{Lkh}$$ (Sao hải nam đóng vai trò loài ưu thế chiếm 14,9% cá thể, cộng sinh mật thiết cùng Lim xanh và Chân chim).

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

  1. Chuẩn hóa định lượng diện tích dinh dưỡng ($S_{dd}$): Thay thế việc ước tính cảm quan cự ly trồng rừng bằng phương pháp toán học Thomasius. Xác định diện tích tối ưu cho Lim xanh ($35,74\text{ m}^2$), Chò chỉ ($42,54\text{ m}^2$), và Sao hải nam ($32,59\text{ m}^2$), tương ứng với cự ly trồng ban đầu từ $3,5\text{ m} \times 3,5\text{ m}$ đến $4,0\text{ m} \times 4,0\text{ m}$.
  2. Khám phá cấu trúc ưu hợp sinh thái cộng sinh: Xác định sự hiện diện đồng thời của Lim xanh và Sao hải nam trong các ô tiêu chuẩn hỗn giao (Lim xanh xuất hiện trong $40%$ ô điều tra của Chò chỉ và chiếm vị trí thứ 2 trong lâm phần Sao hải nam với $HSTT = 1,34$). Điều này bác bỏ quan điểm trồng đơn loài, cung cấp cơ sở cho mô hình phục hồi rừng hỗn giao đa loài bản địa.
  3. Cơ chế tái sinh theo giai đoạn tuổi: Chứng minh quy luật biến thiên sinh thái của Lim xanh và Chò chỉ: giai đoạn mạ non ($H < 50\text{ cm}$) chịu bóng cục bộ, cần độ tàn che $0,5 - 0,6$; giai đoạn trung niên ($H > 100\text{ cm}$) chuyển hóa hoàn toàn sang loài cây ưa sáng tầng vượt tán, đòi hỏi can thiệp tỉa thưa tầng che phủ trên.

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

                     QUY TRÌNH TRIỂN KHAI LÂM SINH ĐỒNG BỘ

Phân tích hiệu quả kinh tế - xã hội

  • Bảo vệ vốn rừng đầu nguồn: Duy trì dung tích $250.000.000\text{ m}^3$ của hồ Sông Mực, đảm bảo nguồn nước tưới tiêu cho hơn $10.000\text{ ha}$ đất nông nghiệp tại hai huyện Nông Cống và Như Thanh.
  • Tối ưu chi phí bảo tồn: Tăng tỷ lệ sống của cây con trồng rừng từ $55%$ (phương pháp truyền thống) lên $>85%$ (ứng dụng cự ly $S_{dd}$ và cây che bóng bạn sinh thái), giảm $35%$ chi phí trồng dặm và chăm sóc trong 3 năm đầu.
  • Phát triển du lịch sinh thái: Hình thành các tuyến du lịch nghiên cứu khoa học tại khu vực hồ Sông Mực - Sông Chàng, tạo sinh kế thay thế cho $8.023$ hộ dân vùng đệm, giảm áp lực săn lùng lâm sản quý.

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

Hạn chế kỹ thuật

  • Thời gian nghiên cứu thực địa giới hạn theo chu kỳ đề tài, chưa theo dõi trọn vẹn $100%$ biên độ biến động vật hậu học (chu kỳ ra hoa sai quả của Chò chỉ kéo dài $4 - 5$ năm/lần).
  • Số lượng mẫu cây Sao hải nam phát hiện còn tương đối ít ($8$ cây trên $10,6\text{ km}$ tuyến), phản ánh nguy cơ đe dọa tuyệt chủng cao nhưng cũng gây khó khăn cho các phân tích hồi quy phi tuyến đa biến.

Hướng nghiên cứu tiếp theo

  1. Ứng dụng GIS và Viễn thám (Remote Sensing): Tích hợp ảnh vệ tinh độ phân giải cao Sentinel-2/Landsat-8 để giải đoán thảm phủ và tự động phân vùng lập địa thích hợp cho từng loài.
  2. Kỹ thuật nhân giống sinh học phân tử: Thử nghiệm vi nhân giống (In-vitro tissue culture) và nuôi cấy mô sẹo cho Sao hải nam và Lim xanh từ phôi non để giải quyết bài toán chu kỳ quả cách năm.
  3. Thiết lập ô định vị sinh thái lâu dài (Permanent Sample Plots - PSP): Theo dõi động thái tăng trưởng đường kính ($D_{1.3}$) và chiều cao ($H_{vn}$) liên tục trong chu kỳ 5 năm, 10 năm và 20 năm.

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

                 HỆ THỐNG GIÁ TRỊ PHÂN TẦNG CHO CÁC ĐỐI TƯỢNG

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

1. Yêu cầu lập địa kỹ thuật để gieo ươm và gây trồng Lim xanh là gì?

Lim xanh thích hợp trên đất feralit đỏ vàng phát triển từ đá phiến sét hoặc sa thạch, tầng đất dày $>50\text{ cm}$, ẩm nhưng thoát nước tốt, độ cao tuyệt đối $<300\text{ m}$. Cây con dưới 1 năm tuổi bắt buộc phải có dàn che bóng với độ che tàn $50 - 60%$, sau đó dỡ bỏ dần để cây đón ánh sáng hoàn toàn.

2. Không gian dinh dưỡng ($S_{dd}$) $35,74\text{ m}^2$ của Lim xanh có ý nghĩa thế nào khi thiết kế mật độ trồng?

Diện tích dinh dưỡng $S_{dd} = 35,74\text{ m}^2$ tương đương với bán kính nuôi dưỡng $R \approx 3,37\text{ m}$ xung quanh mỗi cây trưởng thành. Khi thiết kế trồng rừng phòng hộ hoặc làm giàu rừng, khoảng cách ban đầu nên bố trí $3,5\text{ m} \times 3,5\text{ m}$ (mật độ $\approx 800 - 1.000\text{ cây/ha}$ bao gồm cả cây che bóng đi kèm), sau đó tỉa thưa định kỳ để đạt mật độ ổn định $250 - 300\text{ cây/ha}$ ở tuổi thành thục.

3. Làm thế nào để giải quyết vấn đề chu kỳ sai quả cách năm của Sao hải nam và Chò chỉ?

Cần kết hợp song song hai giải pháp: (1) Xây dựng vườn giống lấy ghép và tuyển chọn cây mẹ vượt trội để thu hái hạt ngay trong các năm sai quả; (2) Nghiên cứu ứng dụng nhân giống sinh dưỡng vô tính bằng phương pháp giâm hom cành non có xử lý chất kích thích ra rễ IBA nồng độ $1,0%$, giúp sản xuất cây con liên tục không phụ thuộc vào mùa quả.

4. Chi phí đầu tư mô hình bảo tồn và chuyển vị nguồn gen tại VQG Bến En là bao nhiêu?

Chi phí phân bổ thành 3 hạng mục chính: (1) Khoanh nuôi xúc tiến tái sinh tự nhiên In-situ: khoảng $1,5 - 2,0\text{ triệu VNĐ/ha/năm}$; (2) Vườn ươm giống bán kiên cố quy mô $20.000\text{ cây/năm}$: khoảng $80 - 100\text{ triệu VNĐ}$; (3) Trồng làm giàu rừng hỗn giao: khoảng $15 - 20\text{ triệu VNĐ/ha}$ trong 3 năm đầu.

5. Khả năng tích hợp mô hình này cho các khu bảo tồn thiên nhiên khác tại miền Bắc?

Mô hình hoàn toàn tương thích với các khu bảo tồn có cùng điều kiện lập địa và sinh thái Bắc Trung Bộ và Đông Bắc như VQG Cúc Phương, VQG Ba Bể, Khu BTTN Kẻ Gỗ. Các kỹ sư chỉ cần thay đổi tham số công thức tổ thành ($CTTT$) dựa trên kết quả điều tra ô tiêu chuẩn tại địa bàn cụ thể.

Kết luận

Đề tài khóa luận tốt nghiệp "Nghiên cứu bảo tồn một số loài thực vật có nguồn gen quý tại Vườn Quốc Gia Bến En, Thanh Hóa" đã giải quyết thành công bài toán định lượng cấu trúc sinh thái và đề xuất khung giải pháp toàn diện cho 3 loài cây bản địa có giá trị bảo tồn đặc biệt: Lim xanh (Erythrophloeum fordii), Chò chỉ (Parashorea chinensis), và Sao hải nam (Hopea hainanensis).

Bằng việc ứng dụng phương pháp ô tiêu chuẩn 7 cây của Thomasius kết hợp điều tra tuyến chi tiết $10,6\text{ km}$, nghiên cứu đã làm sáng tỏ mối quan hệ cộng sinh giữa các loài mục tiêu với các cây bạn trong rừng qua công thức tổ thành ($CTTT$), đồng thời xác định chính xác diện tích dinh dưỡng tối ưu ($S_{dd}$ từ $32,59\text{ m}^2$ đến $42,54\text{ m}^2$). Đây là cơ sở khoa học then chốt để chuyển đổi từ phương thức quản lý bảo tồn thụ động sang chủ động can thiệp lâm sinh học, vừa bảo tồn nguyên vị nguồn gen quý, vừa mở rộng diện tích phục hồi rừng bền vững, đóng góp thiết thực vào chiến lược bảo tồn đa dạng sinh học quốc gia.