Giới thiệu dự án

Sến mật (Madhuca pasquieri (Dubard) H.J.Lam), thuộc họ Hồng xiêm (Sapotaceae), là loài cây gỗ lớn bản địa quý hiếm thuộc nhóm "Tứ thiết" (Đinh - Lim - Sến - Táu) của lâm nghiệp Việt Nam. Với đặc tính cơ lý vượt trội (gỗ cứng, tỷ trọng cao, vân mịn, chống mối mọt) cùng hàm lượng lipid trong hạt đạt 30% - 55% phục vụ y dược và công nghiệp, Sến mật là đối tượng chịu áp lực khai thác cạn kiệt trong nhiều thập kỷ. Loài đã được xếp vào cấp Nguy cấp (EN A1a,c,d) trong Sách Đỏ Việt Nam và Danh lục Đỏ IUCN.

Tại Vườn Quốc gia (VQG) Pù Mát - khu dự trữ sinh quyển thế giới diện tích vùng lõi 94.275 ha thuộc tỉnh Nghệ An, các quần thể Sến mật tự nhiên đang đối mặt với sự suy thoái nghiêm trọng về mật độ và cấu trúc tái sinh. Nghiên cứu "Nghiên cứu một số đặc điểm lâm học của loài cây Sến mật (Madhuca pasquieri) tại Vườn Quốc gia Pù Mát - Tỉnh Nghệ An" do tác giả Nguyễn Văn Anh thực hiện (dưới sự hướng dẫn của TS. Hồ Ngọc Sơn, Đại học Nông Lâm Thái Nguyên) giải quyết bài toán thiếu hụt cơ sở dữ liệu định lượng về sinh thái, cấu trúc quần xã và động thái tái sinh của loài.

Mục tiêu nghiên cứu cụ thể

  1. Xác định tọa độ và đai cao phân bố: Định vị chính xác sự hiện diện của Sến mật theo không gian địa lý và đai cao từ 900m đến 1.500m trên 3 tuyến trọng điểm: Cao Vều, Tam Đình - Tam Hợp, Thác Kèm.
  2. Khảo sát đặc điểm hình thái và vật hậu: Giải mã chỉ số sinh học thân, vỏ, tán, lá, hoa, quả và nhịp điệu sinh học theo mùa.
  3. Định lượng cấu trúc quần xã và hoàn cảnh rừng: Xác định hệ số tổ thành tầng cây cao ($N_a$), cấu trúc tầng thứ, mật độ cây bụi thảm tươi và mối tương quan với các loài cây đi kèm qua ô tiêu chuẩn 6 cây.
  4. Đánh giá hiện trạng tái sinh tự nhiên: Đo lường mật độ cây tái sinh ($N/\text{ha}$), tỷ lệ phân cấp chất lượng cây con và khả năng kế cận tầng tán.
  5. Xây dựng giải pháp kỹ thuật lâm sinh: Đề xuất khung can thiệp bảo tồn nguyên vị (in-situ) và kỹ thuật gieo ươm phát triển loài.

Phạm vi và giới hạn

  • Không gian: 9 Ô tiêu chuẩn (OTC) điển hình (diện tích $1.500,\text{m}^2/\text{OTC}$, tổng diện tích điều tra $13.500,\text{m}^2$) và 6 phẫu diện đất chuyên sâu tại VQG Pù Mát.
  • Thời gian thực hiện: Giai đoạn điều tra thực địa tập trung từ 31/12/2014 đến 20/04/2015 kết hợp phương pháp kế thừa dữ liệu vật hậu chu kỳ năm.

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

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

Nghiên cứu lâm học loài bản địa quý hiếm tại Việt Nam đã có những bước tiến đối với các loài như Lát hoa (Chukrasia tabularis), Trai lý (Garcinia fragraeoides), hay Phỉ ba mũi (Cephalotaxus fortunei). Tuy nhiên, dữ liệu riêng về Sến mật tại dải Trường Sơn Bắc còn phân tán và thiếu tính hệ thống.

Tiêu chí phân tích Phương pháp điều tra truyền thống Giải pháp nghiên cứu tích hợp của đề tài
Độ bao phủ không gian Quan sát tuyến định tính, ghi nhận cá thể rời rạc Lập 9 OTC cố định ($30\text{m} \times 50\text{m}$) bám sát đường đồng mức tại 3 phân khu
Đánh giá vi hoàn cảnh Mô tả thổ nhưỡng ngoại quan đại trà Đào 6 phẫu diện đất sâu >85cm, phân tích lý hóa tính định lượng (pH, mùn, $\text{N}, \text{P}_2\text{O}_5, \text{K}_2\text{O}$)
Cấu trúc không gian Đo đếm tầng tán đơn lẻ Thiết lập ô đồng tâm 6 cây lân cận để xác định hệ số cây bạn đi kèm ($P_o, P_c$)
Đánh giá nhân tác Ước lượng cảm tính Phân ô $400,\text{m}^2$ định kỳ mỗi $500\text{m}$ tuyến, chấm điểm định lượng áp lực (thang 0–3)

Phân loại ưu tiên nghiên cứu và bảo tồn theo khung MoSCoW

  • Must have (Bắt buộc phải có): Tọa độ GPS chính xác 9 OTC; mật độ tầng cây cao và cây tái sinh; phân tích lý hóa tính tầng đất $0 - 90,\text{cm}$; công thức tổ thành tầng cây cao.
  • Should have (Nên có): Phân hạng mức độ gắn kết cây bạn đi kèm qua chỉ số $P_o \ge 30%$ và $P_c > 70%$; phân cấp chất lượng sống cây con (Tốt, Trung bình, Xấu).
  • Could have (Có thể mở rộng): Bản đồ GIS số hóa mật độ phân bố chi tiết toàn lâm phần; phân tích tương quan sinh trưởng chiều cao ($H_{vn}$) và đường kính ($D_{1.3}$).
  • Won't have (Tạm thời chưa thực hiện): Khảo nghiệm di truyền học phân tử phân biệt các dòng vô tính và định lượng chu kỳ năng suất hạt thương mại.

Thiết kế hệ thống và phương pháp luận

Hệ thống điều tra và xử lý dữ liệu lâm học được chuẩn hóa theo quy chuẩn sinh thái rừng nhiệt đới:

[Hệ thống tuyến khảo sát]

Công thức toán học và tiêu chuẩn kỹ thuật

  1. Mật độ tái sinh quy đổi: $$N/\text{ha} = \frac{10.000 \times n}{S}$$ (Trong đó: $n$ là tổng số cây tái sinh xuất hiện, $S$ là diện tích ô khảo sát tính bằng $\text{m}^2$)

  2. Hệ số tổ thành loài ($N_a$): Loài chính tham gia cấu trúc có số cây $n \ge N_{TB} = \frac{N}{m}$ ($N$: Tổng số cây, $m$: Tổng số loài): $$N_a = \frac{n}{N} \times 10$$

  3. Mức độ xuất hiện loài cây bạn đi kèm (Tiêu chuẩn Triệu Văn Hùng):

    • Nhóm I (Rất hay gặp): $P_o \ge 30%$ và $P_c > 70%$
    • Nhóm II (Hay gặp): $15% \le P_o \le 30%$ và $30% \le P_c \le 70%$
    • Nhóm III (Ít gặp): $P_o \le 15%$ và $P_c \le 30%$

Implementation và kết quả

Quy trình thu thập và thuật toán xử lý dữ liệu

Dữ liệu điều tra thực địa được chuẩn hóa và xử lý thông qua pipeline phân tích lâm học bằng ngôn ngữ thống kê R (phiên bản 4.3.1) / Python (3.10) với thư viện tidyversevegan.

import pandas as pd
import numpy as np

def calculate_stand_composition(tree_df):
    """
    Tính toán hệ số tổ thành tầng cây cao theo tiêu chuẩn điều tra quy hoạch rừng
    Input: tree_df chứa cột 'species' và 'count'
    Output: Công thức tổ thành chuẩn dạng văn bản
    """
    total_trees = tree_df['count'].sum()
    num_species = len(tree_df)
    n_mean = total_trees / num_species
    
    formula_parts = []
    other_species_count = 0
    
    for _, row in tree_df.iterrows():
        species = row['species']
        count = row['count']
        coeff = round((count / total_trees) * 10, 1)
        
        if count >= n_mean and coeff >= 1.0:
            formula_parts.append(f"{int(round(coeff))}{species}")
        elif coeff >= 0.5:
            formula_parts.append(f"(+){species}")
        else:
            other_species_count += count
            
    if other_species_count > 0:
        other_coeff = round((other_species_count / total_trees) * 10, 1)
        formula_parts.append(f"{int(round(other_coeff))}LK")
        
    return " + ".join(formula_parts)

# Pipeline kiểm thử trên dữ liệu thực tế OTC 09 (Thác Kèm 3)
data_otc9 = pd.DataFrame({
    'species': ['SM', 'Re', 'De', 'Tram', 'Bo de'],
    'count': [11, 8, 7, 5, 4]
})
print("Công thức tổ thành tầng cây cao OTC 09:", calculate_stand_composition(data_otc9))

Dữ liệu thực nghiệm lâm phần và thổ nhưỡng

1. Phân bố không gian và mật độ cá thể Sến mật tại 9 OTC

STT Ô Tiêu Chuẩn (OTC) Tọa độ GPS (VN-2000) Độ cao (m) Cây trưởng thành (cây/OTC) Cây tái sinh (cây/OTC)
1 Cao Vều 1 0452975 - 2099280 980 5 0
2 Cao Vều 2 0455345 - 2104316 1.010 7 0
3 Cao Vều 3 0455126 - 2104264 1.080 8 3
4 Tam Đình - Tam Hợp 1 0452900 - 2099367 1.120 0 1
5 Tam Đình - Tam Hợp 2 0452867 - 2099324 1.200 5 3
6 Tam Đình - Tam Hợp 3 0454548 - 2099233 1.280 5 2
7 Thác Kèm 1 0454502 - 2099146 1.305 5 2
8 Thác Kèm 2 0454486 - 2099208 1.450 10 3
9 Thác Kèm 3 0454462 - 2099175 1.500 11 5
Tổng 9 OTC ($13.500,\text{m}^2$) -- 980 - 1.500 56 cây 19 cây

Nhận định định lượng: Mật độ cây trưởng thành đạt bình quân $6,22,\text{cây/OTC}$ (tương đương $41,48,\text{cây/ha}$), tăng dần theo đai cao và đạt cực đại tại Thác Kèm 3 ($11,\text{cây/OTC}$ tại độ cao $1.500\text{m}$). Ngược lại, mật độ tái sinh tự nhiên cực kỳ thấp: tổng 9 OTC chỉ ghi nhận 19 cá thể tái sinh ($14,07,\text{cây/ha}$), trong đó tại đai dưới $1.050\text{m}$ (Cao Vều 1, 2) hoàn toàn vắng bóng cây tái sinh ($0,\text{cây/OTC}$).

Mật độ cây Sến mật theo đai cao (m)

2. Đặc tính lý hóa thổ nhưỡng lập địa nơi Sến mật sinh trưởng

Kết quả phân tích 6 phẫu diện đất xác định Sến mật thích hợp trên đất Feralit mùn phát triển trên đá biến chất/trầm tích:

  • Tầng dày đất: Trung bình $>85,\text{cm}$, tầng $A_0$ mỏng $0 - 1,5,\text{cm}$, tầng $A$ dày $1 - 25,\text{cm}$, đất tơi xốp, cấu trúc viên hạt mịn.
  • Thành phần cơ giới (TPCG): Thịt nhẹ đến thịt trung bình (Cát vật lý $\text{CVL} > 0,02,\text{mm}$ chiếm $56,17%$, Sét vật lý $\text{SVL} < 0,02,\text{mm}$ chiếm $43,83%$).
  • Độ chua thổ nhưỡng: $\text{pH}_{\text{KCl}} = 4,7$ (đất chua nhẹ đến trung bình).
  • Hàm lượng dinh dưỡng khoáng tiêu chuẩn:
    • Hàm lượng mùn tổng số: $2,28%$ (mức mùn trung bình khá).
    • Đạm dễ tiêu ($\text{NH}_4^+$): $6,44,\text{mg}/100\text{g}$ đất.
    • Lân dễ tiêu ($\text{P}_2\text{O}_5$): $3,80,\text{mg}/100\text{g}$ đất.
    • Kali dễ tiêu ($\text{K}_2\text{O}$): $6,20,\text{mg}/100\text{g}$ đất.

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

Nghiên cứu mang lại những đóng góp mới cho công tác bảo tồn gen thực vật rừng:

  1. Bộ cơ sở dữ liệu định lượng đầu tiên về Sến mật tại Pù Mát: Chuyển hóa các ghi nhận định tính trước đây thành hệ thống số liệu tọa độ, trắc lượng thân cây ($D_{1.3} = 0,5 - 0,7\text{m}$, $H_{vn} = 25 - 35\text{m}$) và hóa tính đất chuyên sâu.
  2. Phát hiện "nút thắt cổ chai" trong tái sinh tự nhiên: Chỉ ra tỷ lệ thiếu hụt nghiêm trọng của thế hệ cây con dưới tán rừng tự nhiên (chỉ $1,4$ cây tái sinh/$1.000,\text{m}^2$). Nguyên nhân do tán rừng quá dày rậm, hạt có dầu dễ bị nấm mốc hoặc động vật ăn hạt tiêu thụ, đòi hỏi can thiệp lâm sinh mở tán đón sáng.
Hạng mục so sánh Lát hoa (Chukrasia tabularis) (Nguyễn Bá Chất, 1996) Trai lý (Garcinia fragraeoides) (Lê Phương Triều, 2003) Sến mật (Madhuca pasquieri) (Nguyễn Văn Anh, 2015)
Đai cao phân bố chính $300 - 800,\text{m}$ $200 - 600,\text{m}$ (núi đá vôi) $900 - 1.500,\text{m}$ (sườn núi cao dốc $35^\circ - 50^\circ$)
Mật độ tái sinh tự nhiên Khá ($>1.500,\text{cây/ha}$) Trung bình ($500 - 800,\text{cây/ha}$) Cực kỳ thấp ($14,07,\text{cây/ha}$)
Phản ứng lập địa Ưa sáng mạnh từ nhỏ Ưa bóng giai đoạn non Ưa bóng nhẹ khi non, vươn tầng vượt tán ($H > 25\text{m}$) khi trưởng thành
Đặc tính đất thích hợp Feralit đỏ vàng phong hóa từ đá vôi Đất xen kẽ hốc đá vôi Feralit mùn núi cao, $\text{pH}_{\text{KCl}} = 4,7$, mùn $>2,2%$

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

Kế hoạch quản trị bảo tồn và giải pháp kỹ thuật lâm sinh

                          CHIẾN LƯỢC BẢO TỒN SẾN MẬT TẠI PÙ MÁT
[BẢO TỒN NGUYÊN VỊ - IN-SITU]                                       [GÂY ƯƠM & TRỒNG BỔ SUNG - EX-SITU]
  • Định vị cắm mốc 56 cây mẹ giống tại 3 phân khu                    • Thu hái hạt giống từ cây mẹ (Tháng 11 - 12)
  • Khoanh vùng bảo vệ nghiêm ngặt đai 1.200 - 1.500m                 • Xử lý hạt: Ngâm nước ấm 40-45°C trong 6h
  • Mở tán cục bộ (độ che phủ 0.5 - 0.6) quanh cây mẹ                • Gieo ươm bầu đất Feralit mùn + 10% phân chuồng
  • Tuần tra chống khai thác gỗ và hạt thương mại                     • Trồng dặm làm giàu rừng vào mùa mưa (Tháng 8 - 9)
  1. Khoanh vùng bảo vệ nghiêm ngặt:

    • Thiết lập 3 tiểu khu ưu tiên bảo vệ tại Thác Kèm (đai $1.300 - 1.500\text{m}$ - nơi mật độ cây mẹ cao nhất, $10 - 11,\text{cây/OTC}$) và Cao Vều (đai $1.000 - 1.100\text{m}$).
    • Giao khoán bảo vệ rừng gắn với trách nhiệm trực tiếp của các Trạm Quản lý Bảo vệ Rừng Tam Hợp, Thác Kèm và Phúc Sơn.
  2. Kỹ thuật hỗ trợ tái sinh tự nhiên (ANR):

    • Phát luỗng dây leo, cây bụi xâm lấn quanh gốc cây mẹ bán kính $15 - 20\text{m}$ vào tháng 9 - 10 trước mùa hạt rụng.
    • Xới nhẹ đất mặt để hạt rụng dễ tiếp xúc với tầng đất khoáng, khắc phục tình trạng lớp thảm mục dày làm hạt mất sức nảy mầm.
  3. Thu hái và nhân giống nhân tạo:

    • Thu hái quả chín vào tháng 11 - 12 (khi vỏ quả chuyển sang màu nâu vàng, hạt đạt kích thước $2 - 2,2,\text{cm} \times 1,5 - 1,8,\text{cm}$).
    • Tuyệt đối không phơi hạt trực tiếp dưới nắng gắt vì hàm lượng dầu cao ($30 - 55%$) dễ gây ôi hóa phôi; tiến hành gieo ngay sau khi tách hạt.

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

Hạn chế kỹ thuật

  • Thời gian quan trắc: Đề tài thực hiện trong 4 tháng (cuối mùa khô sang đầu mùa mưa), chưa theo dõi trọn vẹn chu kỳ ra hoa kết quả liên tục trong năm trên cùng một cây định vị.
  • Quy mô mẫu phân tích hóa tính: Chỉ mới thực hiện phân tích 6 mẫu đất đại diện do điều kiện thực địa vùng lõi hiểm trở, độ dốc lớn ($>45^\circ$).

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

  1. Thiết lập Ô định vị sinh thái lâu dài (PSP - Permanent Sample Plot): Giám sát sinh trưởng đường kính và chiều cao định kỳ 5 năm/lần.
  2. Ứng dụng chỉ thị phân tử ADN: Đánh giá độ đa dạng di truyền giữa các tiểu quần thể Sến mật tại Pù Mát so với các VQG khác (Cúc Phương, Bến En) nhằm tránh hiện tượng suy thoái cận huyết.
  3. Thử nghiệm nuôi cấy mô tế bào (Tissue Culture): Nhân giống vô tính quy mô lớn phục vụ các chương trình phục hồi rừng gỗ lớn bản địa quốc gia.

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

  • Cơ quan quản lý VQG Pù Mát và Kiểm lâm: Sở hữu bản đồ phân bố tọa độ và số liệu định lượng để bố trí lực lượng tuần tra, phòng chống lâm tặc xâm hại cây di sản.
  • Kỹ sư lâm nghiệp và Ban quản lý rừng phòng hộ: Tiếp cận quy trình kỹ thuật lập địa (độ cao $>900\text{m}$, đất Feralit $\text{pH} = 4,7$) để thiết kế các dự án trồng rừng bản địa bền vững.
  • Cộng đồng nghiên cứu khoa học: Dữ liệu chuẩn mực về sinh thái quần xã họ Sapotaceae tại vùng Trường Sơn Bắc, làm tài liệu tham khảo cho các luận văn thạc sĩ, tiến sĩ chuyên ngành Lâm học.
  • Cộng đồng địa phương vùng đệm (Con Cuông, Tương Dương, Anh Sơn): Định hướng khai thác bền vững hạt Sến mật lấy dầu thông qua mô hình lâm sản ngoài gỗ (NTFP) có kiểm soát, nâng cao sinh kế gắn với bảo vệ rừng.

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

1. Điều kiện lập địa bắt buộc để gây trồng cây Sến mật đạt tỷ lệ sống cao là gì?

Sến mật yêu cầu đai cao từ $900\text{m}$ trở lên, đất Feralit mùn núi cao có tầng dày $>80,\text{cm}$, tơi xốp, thành phần cơ giới thịt nhẹ đến trung bình, thoát nước tốt, độ chua $\text{pH}_{\text{KCl}}$ dao động từ $4,5 - 5,5$, độ che phủ giai đoạn non đạt $0,5 - 0,6$.

2. Tại sao mật độ cây tái sinh của Sến mật tại Pù Mát lại ở mức báo động?

Có 3 nguyên nhân chính: (1) Tán rừng kín lá ($0,7 - 0,8$) khiến cây non thiếu ánh sáng để quang hợp; (2) Hạt Sến mật chứa $30 - 55%$ dầu béo, là thức ăn ưa thích của sóc, chuột và thú nhỏ rừng; (3) Hạt mất sức nảy mầm nhanh trong điều kiện ẩm ướt nếu không tiếp xúc trực tiếp với đất khoáng.

3. Sến mật có thể trồng phân tán ở vùng đồi thấp dưới 500m được không?

Không khuyến khích trồng đại trà ở đai thấp dưới $500\text{m}$ cho mục đích kinh doanh gỗ lớn chất lượng cao, vì ở đai thấp nhiệt độ mùa hè có thể vượt $42^\circ\text{C}$ kèm gió Lào khô nóng, làm suy giảm phẩm chất thân gỗ và tăng tỷ lệ sâu bệnh hại vỏ.

4. Giá trị kinh tế và công dụng dược liệu nổi bật của Sến mật là gì?

Gỗ Sến mật thuộc nhóm "Tứ thiết", siêu bền, chịu lực uốn nắn cao. Hạt cho dầu béo ăn được và ứng dụng chữa bệnh viêm loét dạ dày; lá chứa hoạt chất có thể nấu cao chữa bỏng và vết thương phần mềm hiệu quả.

5. Chi phí và lộ trình phục hồi 1 ha rừng có Sến mật phân bố?

Chi phí khoanh nuôi xúc tiến tái sinh tự nhiên có trồng bổ sung ước tính $15 - 20$ triệu VNĐ/ha trong 3 năm đầu (gồm phát luỗng, gieo dặm $300 - 500,\text{cây/ha}$, bón lót lân vi sinh). Rừng bắt đầu khép tán sau 6 - 8 năm và hình thành tầng gỗ lớn sau $25 - 30$ năm.


Kết luận

Công trình nghiên cứu của tác giả Nguyễn Văn Anh đã định lượng hóa toàn diện các đặc trưng lâm học của loài Sến mật (Madhuca pasquieri) tại VQG Pù Mát. Nghiên cứu khẳng định Sến mật là loài cây gỗ lớn ưu thế tầng tán tại đai cao $900 - 1.500\text{m}$, sinh trưởng trên đất Feralit mùn chua nhẹ ($\text{pH} = 4,7$) nhưng đang gặp đứt gãy nghiêm trọng ở khâu tái sinh tự nhiên ($14,07,\text{cây/ha}$).

Kết quả đề tài cung cấp luận cứ khoa học thực tiễn vững chắc để Ban quản lý VQG Pù Mát triển khai ngay các biện pháp khoanh vùng bảo vệ nghiêm ngặt 56 cây mẹ giống tại Thác Kèm, Cao Vều và Tam Hợp, kết hợp các giải pháp lâm sinh can thiệp mở tán đón sáng nhằm bảo tồn nguồn gen quý hiếm này cho hệ sinh thái rừng Việt Nam.