Giới thiệu dự án

Suy thoái đa dạng sinh học (ĐDSH) toàn cầu đang diễn ra với tốc độ báo động, đe dọa trực tiếp đến cấu trúc hệ sinh thái (HST) và an ninh sinh thái rừng nhiệt đới. Theo thống kê của Tổ chức Bảo tồn Thiên nhiên Quốc tế (IUCN), hơn 8.321 loài thực vật toàn cầu đang đối mặt với nguy cơ tuyệt chủng. Tại Việt Nam – một trong 16 quốc gia có tính đa dạng sinh học cao nhất thế giới – các loài cây gỗ quý đặc hữu đang suy giảm nghiêm trọng do khai thác cạn kiệt và biến đổi sinh cảnh.

Cinnamomum balansae H. Lecomte, 1913 (tên bản địa: Gù hương, Vù hương, Kháo gừng) là loài cây gỗ lớn thuộc họ Re (Lauraceae), được phân hạng Sẽ nguy cấp (VU A1c) trong Sách đỏ Việt Nam (2007), Nguy cấp (EN A1cd, B1+2c ver 2.3) theo Danh lục đỏ IUCN và thuộc Nhóm IIA theo Nghị định 32/2006/NĐ-CP của Chính phủ.

Vấn đề thực tế (Problem Statement)

Tại khu vực huyện Võ Nhai, tỉnh Thái Nguyên (đặc biệt trong đệm và lõi Khu bảo tồn thiên nhiên Thần Sa – Phượng Hoàng), quần thể Cinnamomum balansae đang đứng trước nguy cơ tuyệt chủng cục bộ với các thách thức:

  • Tình trạng khai thác kiệt quệ: 60% người dân điều tra xác nhận buôn bán thương phẩm gỗ với đơn giá lên tới 20.000.000 – 24.000.000 VNĐ/m³, phục vụ chạm khắc lũa, làm nhà và chiết xuất tinh dầu.
  • Áp lực tái sinh tự nhiên: 0/50 hộ dân điều tra từng gây trồng loài; khả năng tái sinh từ hạt trong tự nhiên suy giảm mạnh, phụ thuộc chủ yếu vào tái sinh chồi kém bền vững.
  • Thiếu hụt dữ liệu định lượng: Chưa có bộ dữ liệu chuẩn hóa về tổ thành sinh thái ($IV_i$), mật độ tầng cao, động thái tái sinh và đặc tính thổ nhưỡng làm cơ sở khoa học để thiết lập quy hoạch bảo tồn nội vi (in-situ) và ngoại vi (ex-situ).

Mục tiêu dự án

  1. Khảo sát và đánh giá 100% tri thức bản địa cùng mức độ sử dụng lâm sinh của 50 người dân địa phương tại các xã trọng điểm huyện Võ Nhai.
  2. Định lượng các đặc điểm hình thái thân, lá ($n=100$), hoa và quả ($n=100$) của Cinnamomum balansae so với cơ sở dữ liệu chuẩn quốc gia.
  3. Xác định đặc tính sinh thái tầng cây gỗ (chỉ số quan trọng $IV_i$), độ tàn che ánh sáng, cấu trúc thảm tươi - cây bụi, đặc tính lý - hóa học đất và động thái tái sinh tầng dưới tán.
  4. Đề xuất quy chuẩn kỹ thuật lâm sinh cho bảo tồn in-situ tại Thần Sa – Phượng Hoàng và thiết lập khung nhân giống phục hồi sinh thái ex-situ.

Phương pháp tiếp cận & Phạm vi

  • Phương pháp: Ứng dụng điều tra theo tuyến (57 điểm đo trên toàn huyện, khoảng cách tuyến 3 km), lập 8 ô tiêu chuẩn (OTC) điển hình diện tích $1.000\text{ m}^2$ ($R=18\text{ m}$) tích hợp 5 ô dạng bản (ODB) diện tích $25\text{ m}^2$ ($R=2.82\text{ m}$) phân bố 4 hướng và trung tâm.
  • Phạm vi: Rừng tự nhiên trên núi đất và núi đá vôi tại huyện Võ Nhai (tọa độ $105^\circ 51'05'' - 106^\circ 08'38''\text{E}$, $21^\circ 45'12'' - 21^\circ 56'30''\text{N}$, độ cao $100 - 800\text{ m}$).

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

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

Công tác bảo tồn nguồn gen lâm nghiệp quý hiếm tại miền Bắc Việt Nam hiện đang áp dụng hai mô hình chính:

Tiêu chí Bảo tồn In-situ truyền thống (Rừng đặc dụng) Bảo tồn Ex-situ đơn lẻ (Vườn thực vật) Giải pháp Lâm sinh tích hợp (Đề tài đề xuất)
Bảo toàn cấu trúc gen Cao, giữ nguyên quần xã sinh thái Trung bình, dễ lai tạp hoặc thoái hóa Tối ưu hóa: Kết hợp bảo tồn cấu trúc $IV_i$ tại chỗ & tuyển chọn nguồn giống
Khả năng can thiệp phục hồi Thấp, chủ yếu bảo vệ nghiêm ngặt thụ động Cao, can thiệp nhân giống vô tính/hữu tính Cao: Xúc tiến tái sinh tự nhiên kết hợp tra dặm bằng cây ươm vườn
Sự tham gia của cộng đồng Hạn chế, xung đột lợi ích kinh tế rừng Không có sự tham gia của cư dân bản địa Tối đa hóa: Tận dụng tri thức bản địa, hỗ trợ kỹ thuật nhân giống hộ gia đình
Chi phí duy trì Rất cao cho bộ máy tuần tra bảo vệ Tốn kém kinh phí xây dựng trạm thí nghiệm Chi phí thấp, hiệu quả phân tán theo mô hình nông - lâm kết hợp

Yêu cầu kỹ thuật theo mô hình MoSCoW

  • Must have: Thuật toán xử lý chỉ số sinh thái Daniel Marmillod ($IV_i \ge 5%$), bảng chuẩn đoán lý tính - hóa tính mẫu đất (pH, N, P, K, Mùn), đo đạc hình thái lá/quả ($n=100$).
  • Should have: Phân tầng tái sinh hạt/chồi, đánh giá cây triển vọng ($H \ge 1\text{ m}$), phân tích ma trận tương quan độ tàn che.
  • Could have: Bản đồ số hóa GIS phân bố không gian quần thể Cinnamomum balansae.
  • Won't have: Giải trình tự chuỗi DNA phân tử toàn bộ hệ gen (chuyển tiếp sang pha nghiên cứu tiếp theo).

Thiết kế hệ thống thu thập và xử lý dữ liệu

Stack công nghệ và tiêu chuẩn khoa học

  • Hệ thống phân loại thực vật: Hệ thống Takhtajan (1997), Sách đỏ Việt Nam (2007), IUCN Red List ver 2.3 (1994).
  • Thiết bị đo đạc lâm sinh: Máy định vị GPS Garmin eTrex 30x (sai số $<3\text{ m}$), thước đo cao Blume-Leiss, thước kẹp Caliper Mitutoyo, kính đo độ tàn che lõm.
  • Hệ thống xử lý tính toán sinh thái: Microsoft Excel 2016 / Python Pandas Engine tích hợp công thức phân tích lâm học chuẩn của Daniel Marmillod và Thái Văn Trừng.

Triển khai và kết quả thực nghiệm

Thuật toán và mô hình toán học sinh thái

Chỉ số mức độ quan trọng ($IV_i%$) của loài cây gỗ tầng cao được tính toán dựa trên thuật toán Daniel Marmillod:

$$IV_i (%) = \frac{A_i + D_i + RF_i}{3}$$

Trong đó:

  • Độ phong phú tương đối ($A_i%$): $$A_i (%) = \frac{N_i}{\sum_{j=1}^s N_j} \times 100$$
  • Độ ưu thế tương đối ($D_i%$): Dựa trên tổng tiết diện thân ngang tầng ngực ($G_i$): $$D_i (%) = \frac{G_i}{\sum_{j=1}^s G_j} \times 100 \quad \text{với} \quad G_i = \sum \frac{\pi \cdot D_{1.3}^2}{4}$$
  • Tần suất xuất hiện tương đối ($RF_i%$): $$RF_i (%) = \frac{F_i}{\sum_{j=1}^s F_j} \times 100 \quad \text{với} \quad F_i = \frac{\text{Số ô xuất hiện loài } i}{\text{Tổng số ô điều tra}} \times 100$$

Mật độ tái sinh tự nhiên quy chuẩn trên 1 hecta ($N/\text{ha}$):

$$N/\text{ha} = \frac{10.000 \times n}{S}$$

($n$: số lượng cây tái sinh đếm được trong các ODB; $S$: tổng diện tích các ODB điều tra tính theo $\text{m}^2$).

import numpy as np
import pandas as pd

def calculate_marmillod_ivi(df_trees, total_plots=8):
    """
    Tính toán chỉ số quan trọng IVI (%) theo thuật toán Daniel Marmillod
    Input: DataFrame gồm các cột ['plot_id', 'species', 'd13_cm']
    """
    # 1. Tính tiết diện ngang thân cây G (m2)
    df_trees['basal_area_m2'] = np.pi * ((df_trees['d13_cm'] / 100) / 2) ** 2
    
    # 2. Tổng hợp theo loài
    species_summary = df_trees.groupby('species').agg(
        N_i=('d13_cm', 'count'),
        G_i=('basal_area_m2', 'sum'),
        plots_present=('plot_id', 'nunique')
    ).reset_index()
    
    # 3. Tính toán các thành phần Ai, Di, RFi
    total_N = species_summary['N_i'].sum()
    total_G = species_summary['G_i'].sum()
    
    species_summary['Ai_pct'] = (species_summary['N_i'] / total_N) * 100
    species_summary['Di_pct'] = (species_summary['G_i'] / total_G) * 100
    species_summary['Fi_pct'] = (species_summary['plots_present'] / total_plots) * 100
    species_summary['RFi_pct'] = (species_summary['Fi_pct'] / species_summary['Fi_pct'].sum()) * 100
    
    # 4. Chỉ số IVI
    species_summary['IVI_pct'] = (species_summary['Ai_pct'] + species_summary['Di_pct'] + species_summary['RFi_pct']) / 3
    return species_summary.sort_values(by='IVI_pct', ascending=False)

Kết quả đo đạc thực nghiệm và chỉ số kỹ thuật

1. Đặc điểm hình thái học

Khảo sát thực địa trên $n=100$ lá và $n=100$ quả trưởng thành cho thấy sự sai khác sinh thái học so với dữ liệu mô tả văn bản chuẩn:

Cơ quan Chỉ số đo đạc Giá trị nhỏ nhất Giá trị lớn nhất Giá trị trung bình Dữ liệu Sách đỏ (2007)
Chiều dài cuống ($cm$) 0.60 2.20 1.40 2.00 – 3.00
Chiều dài phiến lá ($cm$) 4.25 9.80 6.45 9.00 – 11.00
Chiều rộng phiến lá ($cm$) 2.10 5.20 3.45 3.50 – 5.00
Quả Chiều dài quả ($cm$) 0.90 1.90 1.35 ~1.30 – 1.50
Đường kính quả ($cm$) 0.60 1.30 0.94 ~1.00

Ghi chú: Kích thước phiến lá trung bình tại Võ Nhai nhỏ hơn chuẩn Sách đỏ $35.5%$, chứng minh sự thích nghi hình thái học với điều kiện khô hạn trên nền karst núi đá vôi.

2. Cấu trúc tầng cây gỗ và vai trò sinh thái ($IV_i$)

  • Công thức tổ thành chung tầng cây cao: $$\text{CTTT} = 12.94\text{Gh} + 23.59\text{Mt} + 15.12\text{Ngh} + 14.88\text{Tl} + 12.35\text{Tt} + \dots + 12.59\text{LK}$$ (Gh: Gù hương; Mt: Mạy tèo; Ngh: Nghiến; Tl: Trai lý; Tt: Trám trắng; LK: Loài khác).
  • Gù hương giữ vai trò loài ưu thế sinh thái quan trọng với $IV_i = 12.94% > 5%$.
  • Đường kính ngang ngực ($D_{1.3}$) đạt $70 - 90\text{ cm}$, chiều cao vút ngọn ($H_{vn}$) đạt $20 - 25\text{ m}$. Độ tàn che trung bình của lâm phần phân bố là 0.32 (chứng minh tính chất cây ưa sáng mạnh).

3. Cấu trúc và động thái tái sinh tầng dưới tán

Khảo sát chi tiết tại 40 ô dạng bản ($S = 1.000\text{ m}^2$) thuộc 8 OTC thu được các thông số định lượng:

OTC Số Mật độ tái sinh ($N/\text{ha}$) Tái sinh từ Hạt ($%$) Tái sinh từ Chồi ($%$) Cây Tốt ($%$) Cây Trung bình ($%$) Cây Xấu ($%$)
OTC 1 720 25.0 75.0 22.2 44.4 33.4
OTC 2 80 0.0 100.0 0.0 100.0 0.0
OTC 3 270 0.0 100.0 0.0 66.7 33.3
OTC 4 0 0.0 0.0 0.0 0.0 0.0
OTC 5 480 33.3 66.7 33.3 50.0 16.7
OTC 6 320 50.0 50.0 25.0 50.0 25.0
OTC 7 0 0.0 0.0 0.0 0.0 0.0
OTC 8 560 66.7 33.3 50.0 33.3 16.7
Trung bình 303 37.5 62.5 24.0 44.4 31.6
  • Mật độ tái sinh bình quân đạt 303 cây/ha, nhưng tỷ trọng tham gia vào công thức tái sinh chung chỉ đạt $0.41%$.
  • Tỷ lệ cây phẩm chất tốt chỉ đạt 24.0%, trong khi cây trung bình và xấu chiếm 76.0%. Hiện tượng hạt bị động vật (chim, sóc) ăn và lớp thảm mục dày ngăn cản hạt tiếp xúc đất khoáng là rào cản tái sinh lớn nhất.

Đổi mới và đóng góp khoa học

  1. Khám phá đặc điểm sinh thái thích nghi karst: Phát hiện sự thu hẹp kích thước bản lá ($6.45\text{ cm}$ so với chuẩn $9 - 11\text{ cm}$) cùng hệ số phân bố tập trung ở đai cao $100 - 300\text{ m}$ trên đất xám feralit phát triển trên đá vôi có độ ẩm cục bộ cao.
  2. Lượng hóa ma trận quan hệ sinh thái tầng cây gỗ: Xác lập công thức liên kết sinh thái thực nghiệm giữa Cinnamomum balansae với các loài chỉ thị (Mạy tèo, Nghiến, Trai lý), mở ra giải pháp trồng hỗn giao phỏng sinh thái thay thế trồng thuần loài kém hiệu quả.
  3. Phát hiện nghịch lý tái sinh: Chứng minh tỷ lệ tái sinh chồi chiếm áp đảo ($62.5%$) phản ánh áp lực khai thác chặt hạ trong quá khứ, giải thích nguyên nhân suy thoái cấu trúc quần thể nếu chỉ áp dụng bảo tồn đóng cửa rừng thụ động.

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

Kịch bản triển khai lâm sinh ứng dụng

Phân tích Chi phí - Lợi ích (Cost-Benefit & ROI)

  • Chi phí đầu tư ban đầu: Khoảng 450.000.000 VNĐ cho 50 ha mô hình bảo tồn và chuyển giao giống.
  • Giá trị kinh tế thương phẩm:
    • Gỗ xẻ hộp: $20.000.000 - 24.000.000\text{ VNĐ/m}^3$.
    • Tinh dầu lá và rễ: Chiết xuất dược liệu có giá trị xuất khẩu cao.
    • Ván sàn mỹ nghệ ($2.2\text{ m} \times 0.7\text{ m} \times 0.05\text{ m}$): $1.500.000 - 2.000.000\text{ VNĐ/tấm}$.
  • Thời gian hoàn vốn (ROI): 12 - 15 năm đối với khai thác gỗ bền vững; 3 - 5 năm đối với khai thác tỉa cành lá thu hồi tinh dầu dược liệu.

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

Hạn chế kỹ thuật

  • Chu kỳ ra hoa không đồng đều: Hiện tượng sai hạt cách niên (1 - 2 năm mới có hoa rộ) làm gián đoạn nguồn hạt ngoại nghiệp.
  • Dữ liệu thời gian thực: Chưa tích hợp cảm biến IoT quan trắc tiểu khí hậu và độ ẩm đất tự động tại các ô tiêu chuẩn dài hạn.

Định hướng mở rộng

  • Ứng dụng công nghệ nuôi cấy mô tế bào (in-vitro) và giâm hom quy mô công nghiệp để vượt qua giới hạn khan hiếm hạt giống.
  • Thiết lập hệ thống GIS Web-map theo dõi diễn biến tài nguyên rừng Gù hương theo thời gian thực kết hợp giám sát viễn thám độ phân giải cao.

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

  • Sinh viên & Học viên chuyên ngành Lâm sinh / Quản lý Tài nguyên: Nắm vững phương pháp luận khảo sát cấu trúc lâm phần ($IV_i$, Marmillod) và quy trình điều tra thực vật học.
  • Kỹ sư Lâm nghiệp & Cán bộ Kiểm lâm: Sở hữu bộ số liệu định lượng về sinh thái đất và ánh sáng để thiết kế quy trình xúc tiến tái sinh tự nhiên chính xác.
  • Doanh nghiệp Lâm nghiệp & Cộng đồng Địa phương: Tiếp cận hướng đi mới về kinh tế rừng dưới tán, phát triển nguồn dược liệu tinh dầu giá trị gia tăng cao.
  • Nhà nghiên cứu Đa dạng Sinh học: Khung tham chiếu cho các công trình nghiên cứu sâu hơn về chi Cinnamomum tại vùng núi đá vôi Đông Bắc Việt Nam.

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

  1. Điều kiện lập địa tối ưu để trồng Gù hương là gì?
    Cây phát triển tối ưu trên đất feralit xám phát triển trên đá vôi hoặc sa thạch, tầng đất dày $>50\text{ cm}$, giàu mùn ($>4%$), pH chua vừa ($5.1 - 5.5$), độ tàn che ánh sáng giai đoạn cây non từ $0.3 - 0.4$.

  2. Tại sao tỷ lệ tái sinh tự nhiên từ hạt của Gù hương lại thấp (37.5%)?
    Do vỏ hạt mỏng có hương thơm hấp dẫn chim và gặm nhấm tiêu thụ; đồng thời lớp thảm mục rừng tự nhiên dày khiến hạt khó tiếp xúc tầng đất khoáng để nảy mầm.

  3. Gù hương có thể trồng thuần loài theo quy mô công nghiệp không?
    Không khuyến khích. Dữ liệu $IV_i$ chứng minh Gù hương có mối quan hệ phụ thuộc sinh thái mật thiết với các loài đi kèm (Mạy tèo, Trám, Nghiến). Cần trồng theo phương thức hỗn giao phỏng tự nhiên.

  4. Kỹ thuật lâm sinh nào quan trọng nhất để thúc đẩy tái sinh chồi?
    Chăm sóc gốc chồi, tỉa bỏ các chồi yếu, giữ lại 1 - 2 chồi thân thẳng, khỏe mạnh nhất mọc sát mặt đất, kết hợp xới xáo đất quanh gốc bán kính $1\text{ m}$.

  5. Chi phí sản xuất cây giống từ hạt và hom vô tính khác nhau như thế nào?
    Cây ươm từ hạt có chi phí thấp ($8.000 - 12.000\text{ VNĐ/cây}$) nhưng nguồn giống bấp bênh; giâm hom/nuôi cấy mô có chi phí cao hơn ($15.000 - 22.000\text{ VNĐ/cây}$) nhưng đảm bảo giữ nguyên phẩm chất di truyền ưu trội.


Kết luận

Đề tài nghiên cứu đặc điểm sinh học của loài Gù hương (Cinnamomum balansae) tại huyện Võ Nhai đã xác lập cơ sở khoa học định lượng vững chắc: khẳng định vai trò ưu thế tầng cây gỗ ($IV_i = 12.94%$), nhận diện điểm nghẽn tái sinh hạt ($37.5%$) và mật độ tầng dưới tán ($303\text{ cây/ha}$). Việc kết hợp đồng bộ giữa bảo tồn nguyên vị thông qua can thiệp lâm sinh (mở tán, phát dọn thảm mục) và nhân giống phục hồi ngoại vị gắn liền với sinh kế cộng đồng là giải pháp cấp thiết duy nhất nhằm bảo vệ nguồn gen quý hiếm này khỏi nguy cơ tuyệt chủng, đóng góp thiết thực cho chiến lược phát triển lâm nghiệp bền vững tại Việt Nam.