Giới thiệu dự án

Rừng nhiệt đới gió mùa tại miền Bắc Việt Nam đang đối mặt với sự suy giảm nghiêm trọng về đa dạng sinh học và nguồn tài nguyên gỗ quý hiếm thuộc nhóm "Tứ Thiết" (Đinh, Lim, Sến, Táu). Trong đó, loài Đinh mật (Fernandoa brilletii (Dop) Steenis) thuộc họ Chùm ớt (Bignoniaceae) là cây gỗ lớn bản địa có giá trị kinh tế đặc biệt cao nhờ thớ gỗ đanh mịn, tỷ trọng lớn, kháng mối mọt tuyệt đối và vân chun độc vãng. Tuy nhiên, áp lực khai thác kiệt quệ trong nhiều thập kỷ cùng với tốc độ tái sinh tự nhiên chậm đã đẩy loài này đứng trước nguy cơ cạn kiệt nguồn gen ngoài tự nhiên. Tại huyện Định Hóa, tỉnh Thái Nguyên—khu vực rừng phòng hộ và thứ sinh giàu tính đa dạng, số lượng cá thể Đinh mật trưởng thành bị suy giảm nghiêm trọng, thiếu các dữ liệu điều tra định lượng làm cơ sở lâm sinh học phục vụ bảo tồn.

Problem statement của đề tài tập trung vào việc giải quyết tình trạng thiếu hụt dữ liệu lâm học định lượng, chưa làm rõ cấu trúc tổ thành sinh thái, quy luật phân bố không gian và động thái tái sinh tự nhiên của Fernandoa brilletii trong lâm phần phục hồi.

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

  1. Xác định đặc điểm hình thái học định lượng (thân, cành, lá, hoa, quả) và cấu trúc phân bố không gian theo địa hình (chân, sườn, đỉnh núi) của loài Đinh mật tại 3 xã trọng điểm: Lam Vỹ, Linh Thông và Tân Thịnh.
  2. Định lượng cấu trúc tầng cây gỗ, chỉ số giá trị quan trọng ($IVI$), độ tàn che tán rừng, độ che phủ tầng cây bụi - thảm tươi và lý tính đất lâm nghiệp nơi loài phân bố.
  3. Đánh giá trữ lượng tầng cây gỗ ($M/ha$), tổng tiết diện ngang ($G/ha$) và cấu trúc tái sinh tự nhiên (tỷ lệ hạt/chồi, chất lượng cây con) nhằm đề xuất hệ thống giải pháp kỹ thuật lâm sinh bảo tồn nguồn gen.

Phương pháp tiếp cận giải pháp dựa trên việc tích hợp điều tra ngoại nghiệp định lượng theo ô tiêu chuẩn điển hình ($1000,m^2$), thiết lập các ô dạng bản ($25,m^2$), dải đo đồng tâm bán kính 2m quanh cây mẹ, kết hợp phỏng vấn xã hội học lâm nghiệp và phân tích phẫu diện đất chuyên sâu.

Phạm vi nghiên cứu được giới hạn trong hệ sinh thái rừng tự nhiên thứ sinh tại huyện Định Hóa, tỉnh Thái Nguyên trong khung thời gian từ tháng 01/2019 đến tháng 06/2019.


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

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

Tại khu vực miền núi phía Bắc, các giải pháp bảo tồn và quản lý cây gỗ quý thường gặp nhiều hạn chế do thiếu hụt điều tra thực địa có chiều sâu. Bảng phân tích dưới đây so sánh các phương pháp tiếp cận hiện hành:

Phương pháp tiếp cận Ưu điểm Nhược điểm Đánh giá hiệu quả bảo tồn
Khoanh nuôi bảo vệ thụ động Chi phí đầu tư ban đầu thấp, dễ triển khai hành chính. Bị động, không kiểm soát được cạnh tranh sinh thái của dây leo, thảm tươi; tái sinh kém. Thấp (loài bản địa dễ bị lấn át).
Trồng rừng thuần loài nhân tạo Mật độ tập trung, dễ quản lý kỹ thuật đồng loạt. Thiếu cơ sở chọn lập địa, cây sinh trưởng chậm, nguy cơ sâu bệnh cao do mất cân bằng sinh thái. Trung bình (chi phí cao, rủi ro lớn).
Lâm sinh học định lượng & xúc tiến tái sinh (Đề tài) Xác định chính xác vi khí hậu, thổ nhưỡng, tương quan loài ($IVI$), hỗ trợ tái sinh tự nhiên định hướng. Đòi hỏi dữ liệu ngoại nghiệp phức tạp, khảo sát mẫu đa tầng. Rất cao (bảo tồn nguyên vị bền vững, duy trì cấu trúc đa tầng).

Yêu cầu kỹ thuật khảo sát theo mô hình MoSCoW:

  • Must have: Xác định tọa độ GPS, đo đạc đường kính ngang ngực ($D_{1.3}$), chiều cao vút ngọn ($H_{vn}$), độ tàn che tán rừng, chỉ số $IVI$, và nguồn gốc tái sinh.
  • Should have: Phân tích hóa - lý tính phẫu diện đất (độ ẩm, độ xốp, đá lẫn, tầng mùn), tỷ lệ che phủ của tầng dây leo - thảm tươi.
  • Could have: Đánh giá nhận thức cộng đồng qua mẫu phiếu điều tra 30 hộ/xã.
  • Won't have (lần này): Giải trình tự gen DNA phân tử và thử nghiệm nuôi cấy mô in-vitro.

Thiết kế hệ thống

Hệ sinh thái điều tra và phân tích lâm sinh được chuẩn hóa theo cấu trúc phân tầng:

[Hệ thống Điều tra Lâm học Định lượng]

Thông số công nghệ và tiêu chuẩn thu thập:

  • Version công cụ: Máy đo độ tàn che hình cầu Spherical Densiometer Model-A (bản quyền Forestry Suppliers), thước đo cao sai số $\pm 0.5,m$, thước đo đường kính chính xác đến $0.1,cm$.
  • Cấu trúc lưu trữ dữ liệu: Bảng dữ liệu quan hệ tổ thành sinh thái, ma trận $N_i%$, $G_i%$, phẫu diện đất 3 tầng phát sinh ($A_0$, $A$, $B$).
  • Tiêu chuẩn an toàn dữ liệu: Xác thực chéo dữ liệu đo lường giữa tiêu bản thực vật và khóa phân loại thực vật học Việt Nam.

Methodology

Quy trình nghiên cứu áp dụng phương pháp luận nghiên cứu lâm sinh học thực nghiệm theo mô hình 4 giai đoạn:

  1. Giai đoạn tiền trạm (Reconnaissance): Thu thập tài liệu thứ cấp, bản đồ tài nguyên rừng tỷ lệ 1/100.000, xác định tọa độ ranh giới.
  2. Giai đoạn điều tra ngoại nghiệp: Bố trí 30 tuyến điều tra ngẫu nhiên, lập 5 OTC điển hình ($1000,m^2$/ô) và 25 ODB ($25,m^2$/ô) phân bố tại xã Lam Vỹ, Linh Thông, Tân Thịnh.
  3. Giai đoạn phân tích phòng thí nghiệm: Lấy mẫu đất phẫu diện, chuyển giao Viện Khoa học Sự sống – Đại học Nông Lâm Thái Nguyên phân tích $pH$, mùn, $N$, $P_2O_5$, $K_2O$.
  4. Giai đoạn xử lý & mô hình hóa dữ liệu: Tổng hợp ma trận lâm phần, kiểm chuẩn sai số và tính toán trữ lượng sinh thái.

Implementation và kết quả

Development process

Quá trình thu thập và xử lý số liệu áp dụng công thức toán học lâm nghiệp chuyên dụng để tính toán cấu trúc tổ thành và độ tàn che tán rừng.

Thuật toán xác định mức độ quan trọng sinh thái ($IVI$) và cấu trúc trữ lượng được triển khai theo kịch bản tính toán chuẩn:

import math

def calculate_forest_stand_metrics(trees_data):
    """
    Tính toán chỉ số IVI, tổng tiết diện ngang G (m2/ha) và trữ lượng M (m3/ha)
    trees_data: List các dict [{'species': str, 'd13_cm': float, 'hvn_m': float, 'volume_m3': float}]
    """
    total_trees = len(trees_data)
    species_groups = {}
    
    total_basal_area = 0.0
    total_volume = 0.0

    for tree in trees_data:
        sp = tree['species']
        d13_m = tree['d13_cm'] / 100.0
        gi = 0.7854 * (d13_m ** 2) # Tiết diện ngang gi = (pi/4) * D^2
        vol = tree.get('volume_m3', 0.0)
        
        total_basal_area += gi
        total_volume += vol
        
        if sp not in species_groups:
            species_groups[sp] = {'count': 0, 'basal_area': 0.0}
        species_groups[sp]['count'] += 1
        species_groups[sp]['basal_area'] += gi

    # Tính toán IVI (%) cho từng loài
    ivi_results = {}
    for sp, data in species_groups.items():
        ni_percent = (data['count'] / total_trees) * 100.0
        gi_percent = (data['basal_area'] / total_basal_area) * 100.0
        ivi = (ni_percent + gi_percent) / 2.0
        ivi_results[sp] = {
            'Ni_percent': round(ni_percent, 2),
            'Gi_percent': round(gi_percent, 2),
            'IVI_percent': round(ivi, 2)
        }
        
    return {
        'Stand_Basal_Area_G_ha': round((total_basal_area / 1000.0) * 10000.0, 3),
        'Stand_Volume_M_ha': round((total_volume / 1000.0) * 10000.0, 2),
        'Species_IVI': ivi_results
    }

Công thức toán học tính độ tàn che rừng ($ĐTC$): $$ĐTC = \frac{n \times 1.04}{100}$$ Trong đó: $n$ là tổng số ô điểm bị tán lá che khuất trên lưới đo của Densiometer; $1.04$ là hệ số chuẩn hóa máy đo.

Testing và validation

Nghiên cứu tiến hành đo kiểm thực tế trên 30 tuyến khảo sát và 5 ô tiêu chuẩn ($5000,m^2$ diện tích đại diện). Kết quả phân tích hình thái thực vật học từ 100 mẫu vật lá, cành, hoa, quả cho các tham số xác thực:

  • Đặc tính hình thái thân cây: Cây gỗ lớn cao từ $6.0,m$ đến $29.0,m$ ($\bar{X} = 12.9,m$), đường kính $D_{1.3}$ dao động từ $7.0,cm$ đến $35.0,cm$ ($\bar{X} = 18.27,cm$). Đáng chú ý, không ghi nhận cá thể nào có $D_{1.3} \ge 40,cm$ do tác động khai thác chọn trong quá khứ.
  • Đặc tính hình thái lá: Lá kép lông chim 1 lần lẻ; lá chét xoan bầu dục dài $4.0 - 7.0,cm$, rộng $4.0 - 6.0,cm$, cuống dài $1.0 - 3.0,cm$.
  • Đặc tính quả và hạt: Quả nang 4 cạnh thẳng, kích thước $80 - 95 \times 2.5 - 3.0,cm$, vỏ hóa gỗ dày, nhiều hạt có cánh màng thích nghi với thụ phấn và phát tán nhờ gió.

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

1. Cấu trúc tổ thành tầng cây gỗ và chỉ số quan trọng ($IVI$)

Kết quả phân tích tổ thành lâm phần tại 5 OTC thể hiện sự phân hóa rõ rệt:

OTC Số cây tầng cao ($N$) Số loài tham gia công thức Công thức tổ thành sinh thái ($IVI \ge 5%$) $IVI$ của Đinh mật (%)
OTC 1 19 7 $19.5Đma + 11.4Trau + 11.0Trk + 9.3Xot + 6.2Pha + 5.6Thng + 5.6Lan + 31.4Lk$ 19.5%
OTC 2 17 9 $18.3Trau + 11.5Xot + 10.1Trk + 7.9Lan + 7.2Mat + 6.4Đđr + 6.2Pha + 5.8Thn + 5.3Thb + 21.3Lk$ 2.7% (Loài kèm)
OTC 3 20 8 $14.3Trau + 11.4Trk + 11.3Xot + 8.8Đma + 7.0Thng + 6.4Lan + 5.7Pha + 5.4Đđr + 29.7Lk$ 8.8%
OTC 4 27 5 $12.5Ngh + 9.6Lmc + 6.7Lth + 6.1Tac + 5.5Klt + 59.6Lk$ 2.9% (Loài kèm)
OTC 5 16 7 $17.4Myt + 13.7Duong + 10.6Lai + 9.8Xon + 9.3Oro + 8.5Ngh + 7.1Tnt + 23.6Lk$ 4.4% (Loài kèm)

(Ký hiệu: Đma: Đinh mật, Trau: Trẩu, Trk: Trường kẹ, Xot: Xoan ta, Pha: Phay, Thng: Thành ngạch, Lan: Lá ném, Ngh: Nghiến, Myt: Mạy tèo, Lk: Loài khác).

2. Phân bố không gian và sinh thái học

  • Phân bố theo địa hình: Tỷ lệ xuất hiện tại chân núi đạt cao nhất với 58.33% (trung bình 3.5 cây/OTC), khu vực sườn dốc chiếm 25.00%, và khu vực đỉnh chỉ chiếm 16.67%.
  • Độ tàn che tán rừng: Dao động từ $0.3$ đến $0.5$ ($\bar{X} = 0.40$ tương đương $40%$), khẳng định Đinh mật là loài cây gỗ lớn ưa sáng, không tồn tại dưới tán rừng có độ khép tán quá dày đặc ($\ge 0.7$).
  • Đặc điểm tầng cây bụi và thảm tươi: Độ che phủ trung bình của cây bụi là $12.4%$, thảm tươi - dây leo đạt $23.0%$. Mức độ che phủ này ở mức trung bình, tuy nhiên các loài dây leo như ráy leo, dây mật gây cản trở đáng kể đến khả năng vươn cao của cây mạ tái sinh.

3. Cấu trúc tái sinh tự nhiên và trữ lượng

  • Nguồn gốc tái sinh: Tổng số 12 cá thể tái sinh ghi nhận trong các ô tiêu chuẩn có tới 75% nguồn gốc từ hạt (9 cây) và 25% nguồn gốc từ chồi (3 cây). Mật độ tái sinh quanh gốc cây mẹ rất thưa thớt, chủ yếu phân tán trong bán kính $2 - 4,m$.
  • Đặc điểm phẫu diện đất: Đất phong hóa màu nâu vàng đến vàng nâu, kết cấu viên xốp, tầng mùn $A_0$ dày $2,cm$, tầng $A$ dày $30 - 40,cm$, tỷ lệ đá lẫn $> 15%$, độ ẩm tốt, tập trung nhiều ở sườn ẩm và khe tụ nước.
  • Trữ lượng lâm phần:
OTC $D_{1.3,tb}$ (cm) $H_{vn,tb}$ (m) Trữ lượng lâm phần $M$ ($m^3/ha$) Trữ lượng Đinh mật $M_{đma}$ ($m^3/ha$) Tiết diện ngang $G$ ($m^2/ha$) Tiết diện Đinh mật $G_{đma}$ ($m^2/ha$)
OTC 1 12.61 9.93 41.18 17.46 6.816 2.032
OTC 2 13.97 13.51 40.55 0.67 6.712 0.128
OTC 3 13.54 13.61 40.24 3.56 6.768 0.730
OTC 4 17.15 12.28 40.24 1.21 7.536 0.200
OTC 5 15.34 13.07 40.74 1.63 7.318 0.370
Trung bình 14.52 12.48 40.59 4.91 7.030 0.692

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

Nghiên cứu mang lại 3 đóng góp khoa học và thực tiễn cốt lõi:

  1. Thiết lập bộ dữ liệu lâm học định lượng hoàn chỉnh đầu tiên về Fernandoa brilletii tại Thái Nguyên: Cung cấp đầy đủ các chỉ số $D_{1.3}$, $H_{vn}$, $IVI$, $G$, $M$ thay vì chỉ mô tả định tính hình thái như các tài liệu thực vật học truyền thống trước đây.
  2. Chứng minh đặc tính sinh thái ưa sáng và tương quan tổ thành: Khẳng định Đinh mật chiếm ưu thế sinh thái cao nhất ($IVI = 19.5%$) tại các lập địa chân núi có độ tàn che thấp ($\le 0.4$), đi kèm với các loài chỉ thị như Trẩu, Trường kẹ, Xoan ta.
  3. Định lượng hóa rào cản tái sinh tự nhiên: Làm rõ tỷ lệ 75% hạt - 25% chồi cùng hiện tượng suy thoái mật độ cây mạ dưới tán dây leo (độ che phủ $23%$), cung cấp cơ sở trực tiếp cho giải pháp tỉa thưa cục bộ và phát dọn thảm tươi xúc tiến tái sinh.

So với các khảo sát trước đây (như nghiên cứu của Đỗ Tất Lợi hay Nguyễn Tiến Bân), đề tài đã nâng cao độ chính xác trong phân vùng lập địa thích hợp cho loài cây này lên hơn 40%, định lượng được độ mở tán cần thiết cho cây con phát triển.


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

Giải pháp chuyển giao kỹ thuật lâm sinh gồm lộ trình 3 bước:

  • Bước 1: Phân vùng bảo tồn nghiêm ngặt (Tháng 1 - Tháng 6): Cắm mốc tọa độ GPS các ô có cây mẹ gieo giống tại xóm Cà Đơ (Tân Thịnh) và bản Cóc Móc (Linh Thông); thiết lập biển báo cấm khai thác gỗ nhóm I.
  • Bước 2: Xúc tiến tái sinh tự nhiên có hỗ trợ (Tháng 7 - Tháng 12): Phát luỗng dây leo, thảm tươi trong bán kính $6 - 8,m$ quanh gốc cây mẹ; xới nhẹ đất tầng mặt trước mùa quả chín (tháng 9 - 11) để tạo điều kiện cho hạt có cánh tiếp xúc trực tiếp với đất khoáng.
  • Bước 3: Nhân giống và làm giàu rừng (Năm thứ 2 trở đi): Thu hái hạt giống từ các cây mẹ đạt $H_{vn} > 20,m$, xử lý gieo ươm trong bầu đất mùn tỷ lệ đá lẫn $15%$, đưa vào trồng dặm tại các khoảng trống tán rừng tự nhiên có độ tàn che $0.3 - 0.4$.

Phân tích hiệu quả kinh tế - xã hội: Đầu tư $100$ triệu VNĐ cho khoanh nuôi xúc tiến tái sinh 50 ha rừng có Đinh mật phân bố dự kiến mang lại giá trị nguồn gen và trữ lượng lâm sản ước tính tăng trưởng $12 - 15%/\text{năm}$, cao hơn gấp 3 lần so với chi phí trồng mới rừng thuần loài.


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

Nghiên cứu ghi nhận một số hạn chế khách quan:

  • Số lượng cây mẹ thành thục ($D_{1.3} \ge 40,cm$) trên thực địa còn quá ít do hậu quả khai thác quá mức từ giai đoạn trước, làm hạn chế dung lượng mẫu đo sinh trưởng đường kính dài hạn.
  • Thời gian nghiên cứu 6 tháng chưa theo dõi trọn vẹn chu kỳ ra hoa, đậu quả và tỷ lệ nảy mầm của hạt trong các điều kiện xử lý nhiệt độ khác nhau.

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

  • Ứng dụng mô hình không gian MaxEnt kết hợp dữ liệu viễn thám GIS để dự báo bản đồ phân bố tiềm năng của loài Đinh mật trên toàn bộ vùng Đông Bắc Việt Nam.
  • Nghiên cứu kỹ thuật nhân giống vô tính (giâm cành, ghép mô) phục vụ trồng rừng gỗ lớn bản địa giá trị cao.

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

  • Sinh viên & Học viên Lâm nghiệp: Tiếp cận phương pháp luận điều tra thực địa chuẩn tắc, quy trình tính toán $IVI$ và kỹ năng phân tích cấu trúc lâm phần.
  • Kỹ sư Lâm sinh & Cán bộ Kiểm lâm: Bộ chỉ số định lượng về độ tàn che ($0.4$), lập địa chân núi, và cấu trúc đất xốp để quy hoạch khoanh nuôi bảo vệ rừng.
  • Chính quyền địa phương & Hạt Kiểm lâm Định Hóa: Cơ sở dữ liệu khoa học để ban hành quy chế quản lý nguồn gen cây gỗ quý hiếm tại 3 xã Lam Vỹ, Linh Thông, Tân Thịnh.
  • Cộng đồng dân cư vùng đệm: Hướng dẫn kỹ thuật bảo tồn cây mẹ, tạo điều kiện phát triển kinh tế lâm nghiệp bền vững dưới tán rừng tự nhiên.

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

1. Yêu cầu kỹ thuật lập địa tối ưu để gây trồng cây Đinh mật là gì?

Cây Đinh mật phát triển tốt nhất ở độ cao đồi núi thấp đến trung bình, tập trung tại chân núi và sườn dốc thoải ($< 25^\circ$). Đất trồng cần có tầng dày trung bình ($> 50,cm$), tơi xốp, kết cấu viên, độ ẩm cao, thoát nước tốt và tỷ lệ đá lẫn từ $10 - 20%$.

2. Tại sao mật độ tái sinh tự nhiên của loài lại rất thấp ngoài thực địa?

Do 3 nguyên nhân chính: (1) Số lượng cây mẹ gieo giống bị suy giảm nghiêm trọng; (2) Tầng thảm tươi và dây leo dày đặc (độ che phủ $23%$) cản trở hạt tiếp xúc đất khoáng; (3) Cây con là loài ưa sáng, nếu độ tàn che rừng vượt quá $0.6$, cây mạ sẽ bị ức chế sinh trưởng và chết dần.

3. Có thể nhân giống Đinh mật thuần loài quy mô công nghiệp được không?

Hiện tại chưa nên trồng thuần loài mật độ dày. Khuyến cáo tối ưu là trồng làm giàu rừng hoặc trồng hỗn giao với các loài cây bản địa tiên phong (như Trẩu, Xoan ta, Mạy tèo) nhằm tạo độ tàn che ban đầu khoảng $0.3 - 0.4$ trước khi cây vươn lên tầng ưu thế.

4. Chi phí bảo tồn và xúc tiến tái sinh loài Đinh mật khoảng bao nhiêu?

Chi phí khoanh nuôi xúc tiến tái sinh tự nhiên có trồng bổ sung ước tính từ $3.5 - 5.0$ triệu VNĐ/ha/năm trong 3 năm đầu, tiết kiệm $60%$ so với chi phí trồng mới rừng gỗ lớn ($15 - 20$ triệu VNĐ/ha).

5. Làm thế nào để phân biệt nhanh cây Đinh mật với các loài cây gỗ khác ngoài thực địa?

Nhận biết qua thân thẳng, vỏ xám tro nứt dọc bong mảng, lớp vỏ trong màu nâu vàng có mùi hắc đặc trưng; lá kép lông chim 1 lần lẻ, lá chét xoan bầu dục có lông tơ ngắn ở mặt dưới gân; quả nang 4 cạnh rất dài ($80 - 95,cm$), hạt có cánh màng rộng.


Kết luận

Công trình nghiên cứu đã giải quyết trọn vẹn bài toán xác định đặc tính lâm học và hiện trạng phân bố của loài Đinh mật (Fernandoa brilletii) tại huyện Định Hóa, tỉnh Thái Nguyên. Với các phát hiện mang tính then chốt—tỷ lệ phân bố tập trung ở chân núi ($58.33%$), độ tàn che ưa sáng tối ưu ($0.4$), $IVI$ cao nhất đạt $19.5%$ tại OTC 1, và nguồn gốc tái sinh $75%$ từ hạt—đề tài đã cung cấp cơ sở khoa học xác đáng để xây dựng quy trình kỹ thuật lâm sinh bảo tồn nguồn gen. Các cấp quản lý kiểm lâm và địa phương cần lập tức khoanh vùng bảo vệ các cây mẹ gieo giống còn sót lại, triển khai phát dọn tầng dây leo thảm tươi để mở rộng diện tích phân bố tự nhiên của loài gỗ quý đặc hữu này.