Giới thiệu dự án

Họ Ngọc lan (Magnoliaceae) là một trong những nhánh thực vật hạt kín (Angiospermae) nguyên thủy bậc nhất thế giới với 17 chi và khoảng 300 loài có giá trị sinh thái, khoa học và kinh tế lâm sản đặc biệt cao. Tại Việt Nam, loài Giổi na (Magnolia grandis (Hu & W.Kumar)) được xếp vào cấp Cực kỳ nguy cấp (Critically Endangered - CR) theo Danh lục Đỏ IUCN (2011) do suy giảm sinh cảnh nghiêm trọng và áp lực khai thác quá mức. Tại huyện Quản Bạ, tỉnh Hà Giang – khu vực ranh giới sinh thái núi đá vôi Đông Bắc – các ghi nhận thực địa ban đầu của tổ chức Quốc tế FFI (Fauna & Flora International) mới chỉ dừng lại ở phạm vi hẹp, gắn biển định danh một số cây mẹ để thu hái giống mà chưa có cơ sở dữ liệu định lượng về cấu trúc quần thể, phân bố không gian và các yếu tố suy thoái lâm sinh.

Vấn đề nghiên cứu và Điểm nghẽn thực tế

  1. Thiếu hụt dữ liệu quần thể cấp cơ sở: Chưa có tọa độ không gian chính xác, mật độ cây mẹ và tỷ lệ tái sinh tự nhiên tại 3 xã trọng điểm: Tùng Vài, Tả Ván, Cao Mã Pờ.
  2. Xung đột sinh kế nông lâm: Hoạt động canh tác Thảo quả (Amomum tsao-ko) dưới tán rừng tự nhiên làm phát quang tầng cây tái sinh, đe dọa khả năng phục hồi tự nhiên của loài.
  3. Tổn thương sinh học và kỹ thuật bảo tồn sai cách: Đóng đinh biển sắt cố định vào thân cây gây hoại tử tầng phát sinh giác - libe; hiện tượng rụng quả do động vật gặm nhấm ăn hạt non làm giảm tỷ lệ nảy mầm tự nhiên.

Mục tiêu dự án

  1. Mục tiêu 1: Điều tra, định vị không gian và lập bản đồ phân bố 100% các cá thể Giổi na phát hiện tại 3 xã Tùng Vài, Tả Ván và Cao Mã Pờ thuộc huyện Quản Bạ.
  2. Mục tiêu 2: Định lượng cấu trúc lâm phần, độ tàn che, tổ thành loài đi cùng thông qua hệ thống Ô tiêu chuẩn ($500\text{ m}^2$) và Ô dạng bản ($25\text{ m}^2$).
  3. Mục tiêu 3: Đánh giá toàn diện các yếu tố tác động tiêu cực sinh học - nhân vi sinh và đề xuất quy trình kỹ thuật bảo tồn ngoại vi (ex-situ) lẫn nội vi (in-situ).

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

Nghiên cứu tập trung tại hệ sinh thái rừng tự nhiên núi đất và núi đá vôi thuộc 3 xã Tùng Vài, Tả Ván, Cao Mã Pờ (độ cao từ $1.262\text{ m}$ đến $1.623\text{ m}$ so với mực nước biển), diện tích rừng phòng hộ $28.185,06\text{ ha}$ của huyện Quản Bạ.

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 Nhược điểm Khả năng mở rộng
Khảo sát truyền thống (FFI) Nhanh, xác định được cây mẹ có quả Chỉ tập trung tuyến mòn, bỏ sót tầng tái sinh, thiếu số hóa bản đồ Thấp, dễ sai số tọa độ
Mục trắc cảm quan lâm sinh Không yêu cầu thiết bị phức tạp Sai số độ tàn che $> 25%$, không đo được phân tầng vi khí hậu Kém, không có khả năng kiểm chứng số
Giải pháp tích hợp Digital Forestry (Đề tài) Tích hợp GPS Garmin, GLAMA v1.2, Canopeo v2.0, MapInfo v17.0 Yêu cầu nhân lực có kỹ năng phân tích GIS và lâm sinh thực nghiệm Cao, dễ dàng đồng bộ cơ sở dữ liệu quốc gia

Ma trận ưu tiên yêu cầu kỹ thuật (MoSCoW)

  • Must have: Lập danh bạ tọa độ UTM/VN2000, đo đường kính $C_{1.3}$, chiều cao vút ngọn ($H_{vn}$), chiều cao dưới cành ($H_{dc}$), phân loại tái sinh hạt (TS) và chồi (TSC).
  • Should have: Phân tích độ che phủ tầng thảo mộc bằng thuật toán phân tích ảnh kỹ thuật số (Canopeo), đo độ tàn che tán rừng qua ảnh bán cầu (GLAMA).
  • Could have: Xác định công thức tổ thành tầng cây gỗ đi cùng theo công thức toán lâm sinh tiêu chuẩn.
  • Won't have: Giải trình tự DNA đa gen nhân (PHYA, LFY, GAI1) trong đợt khảo sát hiện trường này.

Thiết kế hệ thống và Công nghệ sử dụng

  • GIS & Mapping Engine: MapInfo Professional v17.0, MapSource v6.16, WinCmd32.
  • Phân tích hình ảnh lâm sinh di động:
    • GLAMA (Gap Light Analysis Mobile App) v1.2: Phân tích ảnh bán cầu vòm tán, định lượng độ tàn che canopy closure từ 0.00 đến 1.00.
    • Canopeo v2.0: Xử lý ảnh màu RGB dựa trên phân loại tỷ lệ xanh lá ($R/G$, $B/G$) xác định độ che phủ tầng cỏ bổi thảm tươi.
  • Thiết bị định vị và trắc địa: Garmin GPSMAP 64s (độ chính xác sai số $< 3\text{ m}$), địa bàn địa chất Harbin DQL-8, thước đo cao lâm sinh Blume-Leiss.

Methodology

Quy trình nghiên cứu áp dụng phương pháp điều tra tuyến kết hợp lập ô tiêu chuẩn tròn diện tích $500\text{ m}^2$ ($R \approx 12.62\text{ m}$ trên mặt phẳng nằm ngang, hiệu chỉnh góc dốc $\alpha$ theo công thức $R' = \frac{R}{\sqrt{\cos\alpha}}$).

                  SƠ ĐỒ BỐ TRÍ Ô TIÊU CHUẨN TRÒN VÀ Ô DẠNG BẢN
                              (Diện tích 500m2)

                                  Hướng Bắc
                                [ÔDB 2 (ngoài)]
                                [ÔDB 1 (trong)]
                                [ÔDB 5 (trong)]
                                [ÔDB 6 (ngoài)]
             (Bán kính R = 12.62m, 4 ÔDB trong tán + 4 ÔDB ngoài tán 5x5m)

Thuật toán tính toán tổ thành lâm phần

Công thức tổ thành tầng cây gỗ ($K_i$) được xác định dựa trên số cây trung bình của các loài trong ô tiêu chuẩn:

$$\bar{X} = \frac{\sum_{i=1}^{N} n_i}{N}$$

Trong đó $n_i$ là số cá thể của loài $i$, $N$ là tổng số loài bắt gặp. Những loài có $n_i \ge \bar{X}$ sẽ trực tiếp tham gia công thức tổ thành:

$$K_i = \frac{n_i}{\sum n_{\text{tham gia}}} \times 10$$

Implementation và kết quả

Development process & Pipeline xử lý dữ liệu

Quy trình thu thập thực địa 12 tuyến điều tra (tổng chiều dài $36.85\text{ km}$) được chuẩn hóa bằng mã xử lý dữ liệu phân tích lâm phần:

import numpy as np
import pandas as pd

def calculate_stand_composition(species_counts: dict) -> dict:
    """
    Tính toán công thức tổ thành tầng cây gỗ theo quy chuẩn lâm sinh.
    """
    df = pd.DataFrame(list(species_counts.items()), columns=['Species', 'Count'])
    total_trees = df['Count'].sum()
    mean_trees_per_species = df['Count'].mean()
    
    # Lọc loài ưu thế tham gia công thức tổ thành
    dominant_species = df[df['Count'] >= mean_trees_per_species].copy()
    other_species_count = df[df['Count'] < mean_trees_per_species]['Count'].sum()
    
    # Tính hệ số tổ thành Ki trên thang điểm 10
    dominant_species['Ki'] = np.round((dominant_species['Count'] / total_trees) * 10, 1)
    other_ki = np.round((other_species_count / total_trees) * 10, 1)
    
    result = dict(zip(dominant_species['Species'], dominant_species['Ki']))
    if other_ki > 0:
        result['Other_Species'] = other_ki
    return result

# Dữ liệu thực nghiệm lâm phần ÔTC 01 (Tùng Vài)
sample_stand = {
    'Magnolia_grandis': 4,
    'Castanopsis_indica': 8,
    'Lithocarpus_fenestratus': 6,
    'Schima_wallichii': 5,
    'Machilus_odoratissima': 3,
    'Miscellaneous': 6
}
print("Công thức tổ thành:", calculate_stand_composition(sample_stand))

Testing và Validation

  • Độ tin cậy định vị: Kiểm tra chéo 12 phỏng vấn người dân địa phương (tổ tuần tra Voọc, tổ tuần rừng) trước khi thực địa, độ khớp tuyến đạt $91.6%$.
  • Độ chính xác độ che phủ vòm tán: So sánh giữa phương pháp mục trắc truyền thống và phần mềm Canopeo v2.0 / GLAMA v1.2 trên 6 ÔTC cho thấy sai số giảm từ $\pm 18.4%$ xuống còn $\pm 2.1%$.

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

Nghiên cứu đã ghi nhận và định lượng trực tiếp 345 cá thể Giổi na phân bố tự nhiên trên địa bàn 3 xã thuộc huyện Quản Bạ:

                            PHÂN BỐ CÁ THỂ GIỔI NA TẠI 3 XÃ
  Tùng Vài   ████████████████████████████████████ 127 cá thể (36.8%)
  Tả Ván     ███████████████████████████████████  123 cá thể (35.7%)
  Cao Mã Pờ  ███████████████████████████          95 cá thể (27.5%)
             0         30         60         90        120

Bảng tổng hợp cấu trúc quần thể Giổi na tại 3 xã nghiên cứu

Chỉ tiêu đánh giá lâm sinh Xã Tùng Vài ($n=127$) Xã Tả Ván ($n=123$) Xã Cao Mã Pờ ($n=95$) Toàn vùng ($N=345$)
Tầng cây gỗ trắc địa ($G$) 79 cây ($62.20%$) 58 cây ($47.15%$) 44 cây ($46.30%$) 181 cây ($52.46%$)
Tái sinh từ hạt ($TS$) 33 cây ($25.98%$) 47 cây ($38.21%$) 46 cây ($48.40%$) 126 cây ($36.52%$)
Tái sinh từ chồi ($TSC$) 15 cây ($11.81%$) 18 cây ($14.63%$) 5 cây ($5.30%$) 38 cây ($11.02%$)
Sinh trưởng Tốt ($T$) 113 cây ($88.98%$) 90 cây ($73.17%$) 70 cây ($73.68%$) 273 cây ($79.13%$)
Sinh trưởng Trung bình ($TB$) 11 cây ($8.66%$) 18 cây ($14.63%$) 20 cây ($21.05%$) 49 cây ($14.20%$)
Sinh trưởng Xấu/Sâu bệnh ($X$) 3 cây ($2.36%$) 15 cây ($12.20%$) 5 cây ($5.26%$) 23 cây ($6.67%$)
Sinh cảnh Nương Thảo quả 92 cây ($72.44%$) 62 cây ($50.41%$) 7 cây ($7.37%$) 161 cây ($46.67%$)
Sinh cảnh Rừng tự nhiên 34 cây ($26.77%$) 61 cây ($49.64%$) 88 cây ($92.63%$) 183 cây ($53.04%$)
Sinh cảnh Trảng cỏ/Khác 2 cây ($1.57%$) 0 cây ($0.00%$) 0 cây ($0.00%$) 1 cây ($0.29%$)
Dải độ cao phân bố (m) $1.262 - 1.411\text{ m}$ $1.337 - 1.415\text{ m}$ $1.269 - 1.623\text{ m}$ $1.262 - 1.623\text{ m}$

Đặc tính lâm phần tại 6 Ô tiêu chuẩn đại diện

  • Độ tàn che ($Crown\ Density$): Dao động từ $0.40$ (rừng phục hồi ven nương thảo quả) đến $0.85$ (rừng lá rộng thường xanh trung bình núi đá).
  • Độ che phủ thảm tươi ($Ground\ Cover$): Dao động từ $50%$ đến $85%$.
  • Hướng phơi chủ đạo: Đông Nam ($83.3%$) và Đông ($16.7%$), độ dốc sườn núi từ $18^\circ$ đến $42^\circ$.

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

  1. Chuẩn hóa bản đồ phân bố mật độ cao: Khắc phục tình trạng "điểm trắng số liệu" khi nâng tổng số cá thể Giổi na ghi nhận có tọa độ từ $< 30$ cây mẹ lên 345 cá thể đa tầng thế hệ (cây trưởng thành, tái sinh hạt, chồi).
  2. Phát hiện nghịch lý sinh thái nương Thảo quả: Chứng minh diện tích nương Thảo quả chiếm tới $72.44%$ cá thể Giổi na tại Tùng Vài do người dân giữ cây mẹ lấy bóng che, nhưng lại triệt tiêu $85%$ tiềm năng phát triển của cây con tái sinh do hoạt động phát dọn cơ học định kỳ.
  3. Phát hiện bất cập trong hoạt động quản lý bảo tồn thực địa: Chỉ rõ việc đóng đinh sắt định danh tấm biển bảo vệ vào vỏ cây mẹ gây vết thương hở, tạo điều kiện cho nấm hoại sinh xâm nhiễm mô mạch gỗ.
Tiêu chí so sánh Phương pháp cũ (FFI 2018-2020) Điều tra lâm sinh thông thường Phương pháp tiếp cận đề tài
Độ phủ không gian Tuyến chính, điểm lẻ Ô tiêu chuẩn cố định diện rộng Tuyến liên hoàn kết hợp 6 ÔTC tròn cải bằng dốc
Đo lường tán rừng Ước lượng mắt thường Thước đo tán cơ học Xử lý ảnh số GLAMA v1.2 & Canopeo v2.0
Đánh giá tái sinh Ghi nhận sự có mặt Đếm tổng số Phân loại nguồn gốc (Hạt vs. Chồi), cấp chiều cao
Hiệu quả phát hiện Cây mẹ $> 20\text{ cm}$ Ngẫu nhiên Tăng $280%$ số lượng cá thể kiểm kê

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

Quy trình triển khai giải pháp bảo tồn đa tầng

                  QUY TRÌNH BẢO TỒN VÀ PHỤC HỒI LOÀI GIỔI NA
  1. Bảo tồn nguyên vị (In-situ):
    • Tháo gỡ lập tức các đinh sắt đóng biển trên thân cây gỗ Giổi na, chuyển sang dùng đai luồn co giãn thân thiện với sinh trưởng đường kính tầng phát sinh.
    • Thiết lập hành lang bảo vệ quanh 181 cây mẹ tại 3 xã, ưu tiên khu vực núi đá Tả Ván và Cao Mã Pờ (độ cao $> 1.350\text{ m}$).
  2. Kỹ thuật thu hái và gieo ươm nhân giống (Ex-situ):
    • Thời gian thu hái hạt chuẩn xác: Tháng 9 đến tháng 10 hàng năm khi vỏ quả chuyển màu nâu tím nứt nẻ, áo hạt đỏ thẫm.
    • Xử lý hạt: Ngâm nước ấm $40^\circ\text{C}$ trong 6 giờ, đãi sạch lớp dầu áo hạt màu đỏ để tránh nấm mốc, gieo ươm trên giá thể đất tầng $A$ trộn $30%$ mùn núi đá.
  3. Mô hình Nông - Lâm kết hợp bền vững:
    • Tuyên truyền, ký cam kết với các hộ dân canh tác Thảo quả: Giữ lại toàn bộ cây con Giổi na tái sinh có chiều cao $> 50\text{ cm}$ trong quá trình làm cỏ nương.

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

Hạn chế kỹ thuật

  • Địa hình núi cao chia cắt mạnh hiểm trở tại dải biên giới Việt - Trung khiến một số vách đá dựng đứng chưa thể tiếp cận trực tiếp.
  • Chưa có điều kiện nghiên cứu chuyên sâu về nguyên nhân sinh lý học gây hiện tượng gãy/cụt ngọn chính ở cây trưởng thành.

Hướng phát triển tiếp theo

  • Ứng dụng thiết bị bay không người lái (UAV/Drone trang bị cảm biến LiDAR và Multispectral) để quét vòm tán phát hiện hoa Giổi na nở rộ vào tháng 5 - 6.
  • Nghiên cứu sinh học phân tử đánh giá mức độ đa dạng di truyền ($Heterozygosity$) giữa các quần thể Tùng Vài, Tả Ván và Cao Mã Pờ.

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

  • Cơ quan Quản lý Lâm nghiệp & Kiểm lâm: Hạt Kiểm lâm Quản Bạ, Chi cục Kiểm lâm Hà Giang sở hữu bộ dữ liệu GIS tọa độ chính xác của 345 cây Giổi na để quản lý, cấp mã số bảo vệ.
  • Cộng đồng địa phương tại 3 xã: Nâng cao nhận thức bảo tồn loài cây gỗ quý, kết hợp duy trì canh tác Thảo quả bền vững dưới tán rừng phòng hộ mà không vi phạm quy chế bảo vệ loài nguy cấp.
  • Giới nghiên cứu thực vật & Lâm sinh học: Cung cấp dữ liệu định lượng thực chứng về sinh thái học họ Magnoliaceae tại vùng núi cao đá vôi phía Bắc, phục vụ công tác đối sánh thực vật chí khu vực Đông Nam Á.

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

1. Loài Giổi na (Magnolia grandis) phân bố tối ưu ở điều kiện lập địa nào?

Loài phân bố tối ưu tại dải độ cao từ $1.262\text{ m}$ đến $1.623\text{ m}$, thích hợp nhất trên sườn núi đất và núi đá có hướng phơi Đông và Đông Nam, độ dốc từ $18^\circ - 42^\circ$, độ tàn che vòm tán rừng đạt từ $0.60 - 0.85$.

2. Tại sao tỷ lệ cây cụt ngọn ở Giổi na lại cao và xử lý như thế nào?

Hiện tượng cụt ngọn do gió lốc cực đoan vùng núi cao làm gãy chồi đỉnh hoặc sâu đục thân non. Khi ngọn chính bị gãy, cây kích hoạt các chồi nách thay thế. Cần bảo vệ tầng tán xung quanh để tạo đai cản gió tự nhiên cho cây con.

3. Tác động bất lợi lớn nhất đối với tái sinh tự nhiên của loài là gì?

Hai tác động chính gồm: (1) Hoạt động phát dọn thực bì canh tác Thảo quả triệt hạ cây con $< 50\text{ cm}$; (2) Động vật gặm nhấm (sóc) ăn hạt quả non vào tháng 8 - 9 làm suy giảm nguồn hạt giống nảy mầm tự nhiên.

4. Phương pháp nào xác định độ tàn che rừng chính xác nhất trong điều kiện thực địa?

Sử dụng ứng dụng ảnh bán cầu vòm tán GLAMA v1.2 chụp thẳng góc thiên đỉnh kết hợp phần mềm Canopeo v2.0 để xử lý ngưỡng màu tự động, loại bỏ hoàn toàn sai số chủ quan của phương pháp mục trắc bằng mắt.

5. Khả năng tái sinh chồi của loài Giổi na như thế nào?

Khả năng tái sinh chồi đạt $11.02%$ toàn quần thể (38 cá thể). Cây có gốc bị chặt hoặc thân đổ vẫn có khả năng đâm chồi khỏe mạnh, chứng minh sức sống sinh dưỡng của hệ rễ rất cao.

Kết luận

Khóa luận tốt nghiệp đã hoàn thành toàn diện việc điều tra, số hóa và lượng hóa cấu trúc sinh thái quần thể loài Giổi na (Magnolia grandis) Cực kỳ nguy cấp (CR) tại huyện Quản Bạ, tỉnh Hà Giang với 345 cá thể được định vị chính xác. Kết quả phân tích lâm phần tại 6 ÔTC và 12 tuyến điều tra khẳng định khả năng tái sinh tự nhiên của loài vẫn rất khả quan (chiếm $47.54%$ tổng cá thể), nhưng đang chịu sức ép nghiêm trọng từ kỹ thuật bảo tồn thủ công chưa đúng cách và áp lực sản xuất Thảo quả. Các cơ quan quản lý cần lập tức áp dụng hệ thống bảo tồn nguyên vị, tháo gỡ biển đóng đinh sắt, chuyển giao quy trình kỹ thuật gieo ươm hạt giống và xây dựng quy chế bảo vệ lâm phận nghiêm ngặt nhằm duy trì nguồn gen đặc hữu quý giá này cho hệ sinh thái rừng Việt Nam.