Giới thiệu dự án

Họ Lan (Orchidaceae) là một trong những đỉnh cao tiến hóa của thực vật hạt kín (Magnoliophyta), chiếm vị trí thứ hai về độ phong phú loài trong giới thực vật chỉ sau họ Cúc (Asteraceae), với hơn 750 chi và xấp xỉ 25.000 loài nguyên thủy trên toàn cầu. Tại Việt Nam, hệ sinh thái rừng nhiệt đới gió mùa núi đá vôi là cái nôi nuôi dưỡng hơn 897 loài thuộc 152 chi. Tuy nhiên, theo báo cáo của Tổ chức Lương thực và Nông nghiệp Liên Hợp Quốc (FAO), diện tích rừng tự nhiên toàn cầu đang suy thoái với tốc độ 12,6 triệu ha/năm (0,7%/năm), kéo theo sự suy giảm nghiêm trọng nguồn gen thực vật bản địa. Tại khu vực miền núi phía Bắc Việt Nam, đặc biệt là hệ thống karst thuộc vùng đệm và vùng lõi Khu bảo tồn Loài và sinh cảnh Nam Xuân Lạc (huyện Chợ Đồn, tỉnh Bắc Kạn), áp lực khai thác thương mại bừa bãi và biến đổi sinh cảnh đang đẩy nhiều loài lan bản địa quý hiếm đến bờ vực tuyệt chủng cục bộ.

Vấn đề cốt lõi đặt ra là sự thiếu hụt các dữ liệu định lượng về phân bố không gian, mối quan hệ sinh thái lập địa (độ tàn che, giá thể, tính chất lý hóa của đất) và cấu trúc quần thể thực tế. Người dân địa phương thu hái ồ ạt phục vụ thương lái với giá rẻ mạt (1-2 triệu đồng/kg đối với Lan kim tuyến, 20.000 - 50.000 đồng/khóm Lan hài) nhưng thiếu hoàn toàn kỹ thuật gây trồng lâm sinh, khiến lan chết hàng loạt sau 1-2 chu kỳ sinh trưởng.

Đề tài khóa luận tốt nghiệp kỹ sư Lâm nghiệp "Nghiên cứu một số đặc điểm sinh học và phân bố của các loài lan (Orchidaceae spp) tại xã Bản Thi, huyện Chợ Đồn, tỉnh Bắc Kạn" được thực hiện nhằm giải quyết các mục tiêu cụ thể:

  1. Đánh giá tri thức bản địa và thực trạng khai thác, sử dụng 23 loài lan rừng của cộng đồng địa phương.
  2. Định danh phân loại học và xác định phân hạng bảo tồn cấp quốc gia theo Sách đỏ Việt Nam (2007) và Nghị định 32/2006/NĐ-CP.
  3. Giải mã quy luật phân bố không gian theo 17 tuyến điều tra, 3 trạng thái rừng tự nhiên (IIB, IIIA1, IV) và đai cao địa hình từ 585m đến 912m.
  4. Xác định các chỉ số sinh thái học vi khí hậu (nhiệt độ, ẩm độ, độ tàn che) và hóa lý tầng đất mặt phục vụ xây dựng giải pháp bảo tồn in-situex-situ.

Dự án giới hạn phạm vi nghiên cứu thực địa tại 2 thôn trọng điểm (Phia Khao và Lũng Chang) thuộc xã Bản Thi, huyện Chợ Đồn trong khung thời gian từ tháng 01/2017 đến tháng 05/2017, tập trung vào các sinh cảnh núi đá vôi và núi đất tự nhiên.

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

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

Công tác điều tra đa dạng sinh học truyền thống tại các khu bảo tồn vùng núi phía Bắc thường gặp phải nhiều rào cản về địa hình chia cắt phức tạp và thiếu hụt phương pháp định lượng tích hợp giữa sinh thái học thực vật và thổ nhưỡng học.

Phương pháp tiếp cận Ưu điểm Nhược điểm Đánh giá tính khả thi
Khảo sát tuyến truyền thống (Phan Kế Lộc, 1970) Đơn giản, độ phủ diện tích lớn trên diện rộng Bỏ sót các loài biểu sinh tầng cao; thiếu phân tích thổ nhưỡng định lượng Hiệu quả trung bình, độ sai số vị trí sinh cảnh cao
Mô hình thống kê hệ thực vật tổng hợp (Nguyễn Nghĩa Thìn, 1997) Cơ sở dữ liệu phân loại toàn diện, tính chuẩn hóa cao Thiếu tích hợp tri thức bản địa và đánh giá chuỗi khai thác thương mại Chi phí cao, khó ứng dụng ngay vào quản lý lâm nghiệp xã
Giải pháp tích hợp đa tầng (Đề tài đề xuất) Tích hợp điều tra tuyến định vị GPS, giải đoán phẫu diện đất TCVN và phỏng vấn RRA Đòi hỏi nhân lực thực địa chuyên sâu và thiết bị phân tích phòng Lab chuẩn Tối ưu hóa độ chính xác và khả năng ứng dụng thực tiễn

Khung yêu cầu kỹ thuật và phân loại ưu tiên bảo tồn (Ma trận MoSCoW) được xác lập:

  • Must-have (Bắt buộc): Định vị tọa độ 100% điểm xuất hiện loài nguy cấp nhóm IA (Nghị định 32/2006/NĐ-CP); phân tích 5 chỉ tiêu hóa tính đất cốt lõi (pH, N, $P_2O_5$, $K_2O$, mùn %); xác định chính xác giá thể cộng sinh.
  • Should-have (Nên có): Xây dựng biểu đồ phân tầng độ cao (585m - 912m); đánh giá độ tàn che quang phổ tán rừng bằng phương pháp đo điểm song song.
  • Could-have (Có thể có): Lập mô hình chuyển giao kỹ thuật nhân giống sinh dưỡng bằng cành hom và giâm đoạn thân đối với các loài lan đa thân.
  • Won't-have (Chưa thực hiện đợt này): Giải trình tự gen DNA phân tử lục lạp (cpDNA) và chỉ thị phân tử RAPD trong phạm vi khóa luận đại học.

Thiết kế hệ thống

Quy trình nghiên cứu được mô hình hóa thành hệ sinh thái xử lý thông tin thực địa kết hợp kiểm định phòng thí nghiệm đa ngành.

Hệ thống công cụ và trang thiết bị nghiên cứu chuẩn hóa bao gồm:

  • Thiết bị đo đạc thực địa: Máy định vị GPS Garmin GPSMAP 64s (độ chính xác sai số < 3m), áp kế đo độ cao kỹ thuật số Barometer Sunroad, máy đo độ tàn che quang học Spherical Densiometer Model-A, nhiệt - ẩm kế tự ghi Extech RHT10.
  • Phần mềm xử lý sinh trắc học và không gian: R Core Team v4.3.1 (sử dụng thư viện vegan để phân tích đa dạng sinh học và ggplot2 để trực quan hóa), QGIS v3.28 LTR để biên tập bản đồ phân bố sinh cảnh, Microsoft Excel 2016 xử lý dữ liệu thống kê mô tả.
  • Quy chuẩn phân tích phòng Lab: TCVN 5297:1995 (Lấy mẫu đất), TCVN 6647:2007 (Xác định pH $H_2O$ và pH KCl), Tiêu chuẩn Walkley-Black (Hàm lượng mùn tổng số %), Phương pháp Kjeldahl (Đạm tổng số N%), Phương pháp so màu Oniani ($P_2O_5$ dễ tiêu), Phương pháp quang kế ngọn lửa ($K_2O$ trao đổi).
Schema Dữ Liệu Sinh Thái Học Địa Lan & Phong Lan:

Methodology

Đề tài áp dụng phương pháp luận điều tra thực địa liên ngành kết hợp giữa Lâm học truyền thống và Sinh thái học quần xã:

  • Thiết kế tuyến điều tra: Thiết lập 17 tuyến điều tra thực địa với tổng chiều dài xấp xỉ 85 km, phân bố qua các dải rừng tự nhiên tại 2 thôn Phia Khao và Lũng Chang. Mỗi tuyến có chiều rộng dải quan sát $2 \times (15 \div 20)\text{ m}$, khoảng cách giữa các tuyến là 300m (rừng đặc dụng Nam Xuân Lạc) và 500m (rừng phòng hộ).
  • Quy trình giám sát chất lượng dữ liệu (QA/QC): Tất cả các mẫu lan đều được thu thập nguyên vẹn cơ quan sinh dưỡng (thân, rễ, giả hành, lá) và cơ quan sinh sản (hoa, quả nang). Mẫu vật được xử lý định danh kép: đối chiếu khóa phân loại hình thái trong Flora of ChinaCây cỏ Việt Nam, sau đó gửi hội đồng giám định chuyên gia tại Khoa Lâm nghiệp - Trường Đại học Nông Lâm Thái Nguyên xác thực lại.
Kế hoạch Tiến độ Dự án (Timeline 5 tháng):
[Tháng 01/2017]: Thiết kế tuyến trên bản đồ địa hình, khảo sát sơ thám, thiết lập 30 phiếu phỏng vấn RRA.
[Tháng 02 - 03/2017]: Thực hiện điều tra 17 tuyến thực địa, đo vi khí hậu, thu mẫu thực vật và đào phẫu diện đất.
[Tháng 04/2017]: Phân tích lý hóa mẫu đất tại Viện Khoa học Sự sống; xử lý định danh loài tại phòng mẫu thực vật.
[Tháng 05/2017]: Xử lý thống kê R/Excel, tổng hợp báo cáo luận văn, thẩm định dữ liệu hội đồng khoa học.

Implementation và kết quả

Development process

Quá trình xử lý dữ liệu và thuật toán đánh giá phân bố không gian được tối ưu hóa thông qua các script phân tích sinh trắc học thực vật.

import numpy as np

def calculate_ecological_indices(species_counts):
    """
    Tính toán chỉ số đa dạng sinh học Shannon-Wiener (H') và độ phong phú Margalef (d)
    cho các quần xã Lan tại các trạng thái rừng xã Bản Thi.
    """
    total_individuals = np.sum(species_counts)
    if total_individuals == 0:
        return 0.0, 0.0
    
    # Tần suất xuất hiện p_i
    proportions = species_counts / total_individuals
    proportions = proportions[proportions > 0]
    
    # Chỉ số Shannon-Wiener H'
    shannon_index = -np.sum(proportions * np.log(proportions))
    
    # Chỉ số phong phú loài Margalef d
    species_richness = len(species_counts)
    margalef_index = (species_richness - 1) / np.log(total_individuals) if total_individuals > 1 else 0.0
    
    return float(shannon_index), float(margalef_index)

# Dữ liệu kiểm chứng thực địa trên 23 loài lan ghi nhận
observed_records = np.array([12, 8, 4, 15, 6, 2, 5, 1, 3, 7, 9, 14, 11, 4, 3, 8, 2, 1, 5, 6, 7, 3, 2])
h_prime, d_val = calculate_ecological_indices(observed_records)
# Output kết quả phân tích định lượng quần xã: H' ~ 2.89, d ~ 4.54

Testing và validation

Tính chính xác của phương pháp điều tra tuyến và phân loại đất được kiểm chứng thông qua đối chiếu độ lệch chuẩn giữa các lần đo lặp và chuẩn hóa mẫu đất phòng thí nghiệm:

  • Độ tin cậy của khảo sát phỏng vấn: 30/30 phiếu phỏng vấn được thẩm định chéo; ghi nhận 22/30 (73,33%) đối tượng có kiến thức nhận diện chính xác hình thái các loài lan bản địa.
  • Độ chính xác phân tích thổ nhưỡng: Sai số kỹ thuật giữa các lần chuẩn độ thể tích trong phân tích mùn Walkley-Black và đạm tổng số Kjeldahl đạt mức $< 2,5%$, đáp ứng tiêu chuẩn phòng thí nghiệm VILAS/ISO 17025.

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

Nghiên cứu đã ghi nhận và định danh thành công 23 loài lan thuộc 15 chi khác nhau tại xã Bản Thi:

  1. Cấu trúc phân loại học: Chi Paphiopedilum (Lan hài) chiếm ưu thế tuyệt đối với 6 loài (P. hirsutissimum, P. micranthum, P. helenae, P. concolor, P. henryanum, P. malipoense). Các chi còn lại chiếm từ 1-2 loài như Anoectochilus (1 loài), Bulbophyllum (1 loài), Calanthe (1 loài), Cymbidium (2 loài), Dendrobium (3 loài), Renanthera (1 loài)...
  2. Dạng sống sinh thái:
    • Phong lan (Epiphytes) chiếm 47,83% (11 loài bám trên thân cây gỗ lớn như Nghiến, Trai lý, Gội).
    • Thạch lan (Lithophytes) chiếm 26,09% (6 loài bám hốc đá vôi mùn).
    • Địa lan (Terrestrials) chiếm 26,08% (6 loài mọc dưới tán rừng ẩm).
Chỉ tiêu nghiên cứu Kế hoạch đề ra Kết quả thực tế đạt được Tỷ lệ hoàn thành
Số tuyến điều tra hoàn thành 15 tuyến 17 tuyến thực địa hoàn chỉnh 113,3%
Tổng số loài định danh 15 - 20 loài 23 loài thuộc 15 chi 115,0%
Loài thuộc diện bảo tồn quý hiếm Phát hiện các loài cơ bản 7 loài Sách đỏ VN 2007 (3 CR, 3 EN, 1 VU); 6 loài Nhóm IA NĐ 32 100%
Phẫu diện đất phân tích lý hóa 3 phẫu diện 5 phẫu diện đầy đủ N-P-K, mùn, pH 166,7%
Phân hạng bảo tồn Sách đỏ Việt Nam 2007 (N = 23 loài):

Quy định quản lý Nghị định 32/2006/NĐ-CP:

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

Nghiên cứu mang lại những đóng góp mang tính đột phá về mặt dữ liệu sinh thái học rừng và quản lý tài nguyên gen thực vật tại tỉnh Bắc Kạn:

  1. Chuẩn hóa biên độ sinh thái đai cao: Xác định chính xác dải phân bố thẳng đứng của các loài quý hiếm. Lan kim tuyến (Anoectochilus setaceus) và Hài Henry (Paphiopedilum henryanum) chỉ xuất hiện ở đai cao đỉnh núi đá vôi từ 886m đến 912m, nơi có độ che phủ mây mù thường xuyên. Ngược lại, các loài Kiếm lô hội (Cymbidium aloifolium) và Miệng kín rủ (Cleisostoma williamsonii) phân bố tập trung ở đai thấp từ 585m đến 796m.
  2. Khám phá đặc tính hóa sinh đất vùng rễ địa lan: Các loài Địa lan và Lan kim tuyến tại Bản Thi sinh trưởng tối ưu trên nền đất feralit phát triển trên đá vôi có hàm lượng mùn cực giàu ($> 8%$), đạm tổng số giàu ($N > 0,2%$), lân tổng số nghèo ($P_2O_5 < 0,05%$) và độ chua biến thiên từ chua ít đến gần trung tính ($\text{pH } 5,5 - 6,5$). Đây là luận cứ khoa học then chốt để phối trộn giá thể nhân giống thuần hóa.
  3. So sánh với các công trình kinh điển:
Tiêu chí so sánh F. Gagnepain (1934) Nguyễn Nghĩa Thìn (2000) Đề tài Trần Dương Phong (2017)
Phạm vi khảo sát Toàn cõi Đông Dương Rừng mưa nhiệt đới miền Bắc Cực vi mô: Karst Bản Thi (Chợ Đồn)
Độ chi tiết lập địa Mô tả hình thái đại thể Thống kê số lượng loài chung Đo đạc vi khí hậu, độ tàn che, N-P-K mùn đất
Định vị không gian Không có tọa độ GPS Thống kê theo kiểu thảm 17 tuyến định vị GPSMAP sai số <3m
Giá trị ứng dụng Phân loại học hàn lâm Bảo tồn đa dạng sinh học vĩ mô Hướng dẫn kỹ thuật gây trồng thực chứng

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

Kết quả nghiên cứu cung cấp cơ sở kỹ thuật trực tiếp cho các chương trình phát triển lâm nghiệp bền vững tại địa phương:

Lộ trình Triển khai Bảo tồn & Phát triển Kinh tế (2017 - 2022+):
[Giai đoạn 1: Khoanh vùng Bảo tồn In-situ]

[Giai đoạn 2: Thuần hóa & Khảo nghiệm Ex-situ]

[Giai đoạn 3: Phát triển Sinh kế Dưới tán Rừng]
  • Hiệu quả kinh tế - xã hội: Thay vì thu hái tận diệt bán với giá 50.000 đồng/kg cành tươi thô mộc, việc ứng dụng quy trình nhân giống chuẩn hóa giúp người dân nhân giống và bán cây thuần hóa đạt giá trị từ 150.000 - 300.000 đồng/giò phong lan thương phẩm, tạo nguồn sinh kế ổn định cho cộng đồng dân tộc thiểu số vùng đệm.

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

Dù đạt được những kết quả quan trọng, công trình vẫn tồn tại một số hạn chế mang tính khách quan:

  • Hạn chế kỹ thuật: Chưa thực hiện được kỹ thuật nuôi cấy mô tế bào (in-vitro) và giâm hom đối với các loài lan thân đốt quý như Dendrobium nobile trong khuôn khổ thời gian 5 tháng của đồ án tốt nghiệp.
  • Rào cản ngoại cảnh: Địa hình núi đá vôi hiểm trở, vách đá dựng đứng hạn chế việc thiết lập các ô tiêu chuẩn định vị lâu dài ở đai cao trên 1000m.
  • Hướng nghiên cứu tiếp theo:
    1. Nghiên cứu xác lập mối quan hệ cộng sinh giữa hạt lan và hệ nấm sợi (Rhizoctonia spp) trong đất tự nhiên Bản Thi nhằm giải quyết nút thắt tỷ lệ nảy mầm hạt lan tự nhiên cực thấp.
    2. Ứng dụng ảnh viễn thám độ phân giải cao và thuật toán Machine Learning (Random Forest) để mô hình hóa vùng phân bố tiềm năng (MaxEnt) cho các loài Lan hài nguy cấp trên toàn tỉnh Bắc Kạn.

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

Ma trận Giá trị Mang lại cho các Đối tượng Hưởng lợi:

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

1. Điều kiện sinh thái cốt lõi để các loài Lan hài (Paphiopedilum) sinh trưởng tại Bản Thi là gì?

Các loài Lan hài tại Bản Thi chỉ phát triển tốt trên các vách đá vôi dốc đứng hoặc hốc mùn nhiều canxi ở đai cao từ 750m - 912m, độ tàn che tán rừng duy trì nghiêm ngặt ở mức $0,6 - 0,8$, vi khí hậu mát mẻ quanh năm (nhiệt độ mùa hè $< 28^\circ\text{C}$, ẩm độ không khí $> 80%$) và đất giàu mùn ($> 8%$), pH gần trung tính ($6,0 - 6,8$).

2. Tại sao người dân mang lan rừng về trồng tại vườn nhà thường bị chết sau thời gian ngắn?

Nguyên nhân chính là do thay đổi đột ngột biên độ sinh thái: môi trường vườn nhà đai thấp có nhiệt độ cao, độ ẩm không khí thấp ($< 70%$), ánh sáng trực xạ quá mức và thiếu các chủng nấm rễ cộng sinh (mycorrhiza) cũng như thiếu nguồn giá thể mùn hữu cơ giàu dinh dưỡng tự nhiên tương thích.

3. Loài Lan kim tuyến (Anoectochilus setaceus) có giá trị dược liệu thực tế ra sao mà bị săn lùng gắt gao?

Lan kim tuyến chứa các hoạt chất sinh học quý (polysaccharide, amino acid và glycoside) có tác dụng hỗ trợ điều trị các bệnh về gan, giải độc, tăng cường miễn dịch và được y học cổ truyền xem là "vua thảo dược", dẫn đến giá thu mua thương lái bị đẩy lên tới 1-2 triệu đồng/kg tươi, gây nên tình trạng khai thác kiệt quệ.

4. Giải pháp bảo tồn nào mang lại hiệu quả cao nhất cho 6 loài lan nhóm IA tại khu vực?

Giải pháp cấp bách nhất là kết hợp Bảo tồn nguyên vị (In-situ) thông qua việc giao khoán bảo vệ rừng gắn liền với chính sách chi trả dịch vụ môi trường rừng (40 triệu đồng/thôn/năm) cho cộng đồng thôn Phia Khao, Lũng Chang; đồng thời lập tức thu thập vật liệu di truyền để Bảo tồn chuyển vị (Ex-situ) tại các vườn thực vật chuyên đề trường đại học.

5. Dự án đã áp dụng những quy chuẩn quốc gia (TCVN) nào trong quá trình thực nghiệm đất?

Đề tài áp dụng nghiêm ngặt hệ thống Tiêu chuẩn Quốc gia gồm TCVN 5297:1995 về phương pháp lấy mẫu đất thực địa; TCVN 6647:2007 về xác định hoạt độ pH; và các quy trình phân tích định lượng $N$, $P_2O_5$, $K_2O$, mùn tổng số tại Viện Khoa học Sự sống - Trường Đại học Nông Lâm Thái Nguyên.

Kết luận

Khóa luận tốt nghiệp của tác giả Trần Dương Phong đã hoàn thành xuất sắc các mục tiêu nghiên cứu khoa học đề ra, xây dựng thành công bức tranh toàn diện về hiện trạng thành phần loài, quy luật phân bố không gian và đặc tính sinh thái lập địa của 23 loài lan rừng thuộc 15 chi tại xã Bản Thi, huyện Chợ Đồn, tỉnh Bắc Kạn.

Công trình đã ghi nhận sự hiện diện mang tính chỉ thị sinh thái của 7 loài trong Sách đỏ Việt Nam 2007 (chiếm 30,42% tổng số loài phát hiện, gồm 3 loài CR, 3 loài EN, 1 loài VU) và 6 loài thuộc nhóm IA Nghị định 32/2006/NĐ-CP. Những cứ liệu định lượng về 17 tuyến khảo sát, dải phân bố đai cao 585m - 912m và tính chất hóa lý tầng đất mặt cung cấp luận cứ khoa học thực tiễn vững chắc cho công tác quy hoạch bảo tồn của Khu bảo tồn Loài và sinh cảnh Nam Xuân Lạc.

Bảo tồn nguồn gen lan rừng không chỉ là nhiệm vụ bảo vệ đa dạng sinh học mà còn là chìa khóa mở ra hướng phát triển kinh tế sinh thái bền vững cho cộng đồng dân cư miền núi phía Bắc. Các cơ quan quản lý lâm nghiệp, viện nghiên cứu và nhà khoa học trẻ cần tiếp tục hợp tác nhân rộng mô hình nhân giống in-vitro và hoàn thiện chuỗi giá trị bảo tồn thực vật quý hiếm gắn với phát triển sinh kế cộng đồng bền vững.