Giới thiệu dự án
Rừng núi đá vôi tại Việt Nam là một trong những hệ sinh thái đặc thù có tính đa dạng sinh học cao nhưng vô cùng nhạy cảm và dễ bị tổn thương trước các biến đổi môi trường cũng như áp lực nhân sinh. Theo các báo cáo thống kê lâm nghiệp, diện tích rừng tự nhiên trên núi đá vôi miền núi phía Bắc đang chịu sức ép nghiêm trọng từ tình trạng khai thác gỗ vượt mức, phá rừng làm nương rẫy và biến đổi khí hậu. Trong số các loài thực vật quý hiếm, Thiết sam giả lá ngắn (Pseudotsuga brevifolia W.Fu, 1975) thuộc họ Thông (Pinaceae) là loài cây gỗ lá kim bản địa đặc hữu, sinh trưởng tại các đới sinh thái hiểm trở có độ cao từ 500m đến 1.500m so với mực nước biển.
Loài Thiết sam giả lá ngắn hiện đã được xếp vào cấp VU (Sẽ nguy cấp) trong Sách Đỏ Việt Nam (2007) cũng như Danh lục Đỏ IUCN, đồng thời thuộc nhóm thực vật rừng nguy cấp, quý, hiếm cần được bảo tồn nghiêm ngặt theo Nghị định 32/2006/NĐ-CP (sau này là Nghị định 06/2019/NĐ-CP và Nghị định 84/2021/NĐ-CP). Thực tế cho thấy, quần thể Pseudotsuga brevifolia tại Cao nguyên đá Đồng Văn, tỉnh Hà Giang đang đối mặt với nguy cơ suy thoái cục bộ do cấu trúc quần thể bị chia cắt, tái sinh tự nhiên gặp trở ngại bởi tầng đất mùn mỏng và sự cạnh tranh gay gắt từ thảm thực vật bì sinh.
+-------------------------------------------------------------------------------+
| KHUNG DỮ LIỆU ĐỀ TÀI NGHIÊN CỨU |
+-------------------------------------------------------------------------------+
| * Đối tượng: Thiết sam giả lá ngắn (Pseudotsuga brevifolia W.Fu, 1975) |
| * Địa bàn: Xã Thài Phìn Tủng & Xã Sà Phìn, Huyện Đồng Văn, Tỉnh Hà Giang |
| * Tọa độ địa lý: 23°06'06"B - 23°21'17"B | 105°07'35"Đ - 105°24'40"Đ |
| * Độ cao trung bình: 1.200m - 1.500m so với mực nước biển (Địa hình Karst) |
| * Cơ quan thực hiện: Khoa Lâm nghiệp - Đại học Nông Lâm Thái Nguyên |
+-------------------------------------------------------------------------------+
Vấn đề nghiên cứu và Điểm nghẽn thực tế (Problem Statement & Pain Points)
- Tình trạng phân mảnh quần thể: Các dải phân bố của Thiết sam giả lá ngắn bị cô lập trên các vách đá tai mèo dốc đứng, khiến dòng gen bị hạn chế và suy giảm độ đa dạng di truyền.
- Áp lực sinh kế và khai thác quá mức: Dân số huyện Đồng Văn với hơn 88% là đồng bào dân tộc Mông, tập quán khai thác gỗ làm nhà, cột kèo và đồ thủ công mỹ nghệ tạo áp lực nặng nề lên tầng cây gỗ trưởng thành.
- Thiếu hụt dữ liệu định lượng phục vụ quy hoạch bảo tồn: Các nghiên cứu trước đây chủ yếu dừng lại ở mức mô tả định tính hình thái học, thiếu các chỉ số toán học sinh thái định lượng ($IVI%$, phân bố $N/D_{1.3}$, phân bố $N/H_{vn}$) để xây dựng các giải pháp lâm sinh chính xác theo từng vi địa hình (sườn núi và đỉnh núi).
Mục tiêu dự án
- Khảo sát và xác định cấu trúc tổ thành, mật độ: Định lượng chỉ số quan trọng ($IVI%$) và mật độ tầng cây gỗ ($N/ha$) của Thiết sam giả lá ngắn cùng các loài đồng hành tại vị trí sườn và đỉnh núi.
- Mô hình hóa quy luật phân bố lâm phần: Xây dựng quy luật phân bố số cây theo cấp đường kính ($N/D_{1.3}$) và cấp chiều cao ($N/H_{vn}$) để làm rõ kết cấu tầng tán và tiến trình diễn thế.
- Đánh giá năng lực tái sinh tự nhiên: Xác định mật độ tái sinh ($N/ha$), công thức tổ thành cây mầm, nguồn gốc sinh dưỡng (hạt/chồi) và chất lượng cây con triển vọng.
- Đề xuất chiến lược bảo tồn và phục hồi: Xây dựng hệ thống giải pháp kỹ thuật lâm sinh, khoanh nuôi xúc tiến tái sinh và mô hình bảo tồn nguyên vị (in-situ) gắn liền với quản lý rừng cộng đồng.
Phương pháp tiếp cận và Kết quả kỳ vọng
- Phương pháp tiếp cận: Vận dụng lý thuyết sinh thái phát sinh quần xã thực vật kết hợp phương pháp lập ô tiêu chuẩn ngẫu nhiên tầng thứ ($OTC$ $500m^2$, Ô thứ cấp $25m^2$, Ô dạng bản $1m^2$), áp dụng thuật toán thống kê sinh học để lượng hóa cấu trúc rừng.
- Chỉ số đo lường thành công:
- Định lượng chính xác $100%$ chỉ số $IVI%$ cho toàn bộ các loài cây gỗ xuất hiện trong 18 OTC khảo sát.
- Tỷ lệ định danh tên khoa học cho mẫu thực vật đạt $\ge 90%$.
- Xây dựng hoàn chỉnh đồ thị phân bố $N/D_{1.3}$ và $N/H_{vn}$ với độ tin cậy thống kê cao.
Phạm vi và Giới hạn nghiên cứu
- Không gian: Khu vực rừng núi đá tự nhiên thuộc 2 xã trọng điểm Thài Phìn Tủng và Sà Phìn, huyện Đồng Văn, tỉnh Hà Giang.
- Thời gian: Dữ liệu điều tra ngoại nghiệp và phân tích nội nghiệp được thực hiện tập trung từ tháng 02/2014 đến tháng 05/2014.
Phân tích và thiết kế giải pháp
Phân tích hiện trạng
Nghiên cứu cấu trúc rừng tự nhiên đã trải qua nhiều giai đoạn phát triển từ các phương pháp định tính sơ khai đến các mô hình toán học sinh thái hiện đại. Bảng so sánh dưới đây phân tích ưu và nhược điểm của các phương pháp nghiên cứu cấu trúc lâm phần:
| Phương pháp nghiên cứu |
Cơ sở lý thuyết |
Ưu điểm |
Nhược điểm |
Đánh giá mức độ phù hợp |
| Phẫu đồ đứng (P.W. Richards, Thái Văn Trừng) |
Trực quan hóa tầng tán bằng bản vẽ mặt cắt sinh thái |
Trực quan, dễ quan sát sự phân tầng $A_1, A_2, A_3, B, C$ |
Mang tính định tính cao, không định lượng được quy luật phân bố |
Thích hợp khảo sát sơ bộ, không đủ cơ sở điều chế rừng |
| Phân cấp sinh trưởng Kraft (1884) |
Phân 5 cấp cây ưu thế đến bị chèn ép dựa vào vị trí tán lá |
Đơn giản, phản ánh rõ sự phân hóa không gian sống |
Chỉ áp dụng hiệu quả cho rừng trồng thuần loài đều tuổi |
Không phản ánh được tính phức tạp của rừng tự nhiên hỗn loài đa tầng |
| Sinh thái toán học & Chỉ số $IVI$ (Lamprecht, Curtis-McIntosh) |
Kết hợp mật độ ($A_i$), độ tàn che/tiết diện ($D_i$) và tần suất ($RF_i$) |
Định lượng chính xác vai trò sinh thái từng loài, làm rõ cấu trúc vi địa hình |
Yêu cầu khối lượng đo đếm ngoại nghiệp lớn, xử lý số liệu phức tạp |
Lựa chọn tối ưu (Được áp dụng trong đề tài) |
Phân tích yêu cầu nghiên cứu theo mô hình MoSCoW
- Must have (Bắt buộc phải có): Xác định chính xác chỉ số $IVI%$, mật độ $N/ha$, phân bố $N/D_{1.3}$ và $N/H_{vn}$ của Thiết sam giả lá ngắn ở cả vị trí sườn và đỉnh núi đá vôi.
- Should have (Nên có): Thống kê chi tiết chất lượng (cấp A, B, C) và nguồn gốc (hạt/chồi) của tầng cây tái sinh trong các ô dạng bản ($ODB$).
- Could have (Có thể có): So sánh tương quan vi khí hậu giữa hai sườn đón gió và khuất gió trên địa hình Karst.
- Won't have (Chưa thực hiện trong giai đoạn này): Phân tích giải trình tự gen DNA phân tử đánh giá dòng gen quần thể.
Thiết kế hệ thống điều tra và Mô hình tính toán
Quy trình nghiên cứu được thiết kế theo cấu trúc module logic từ khảo sát ngoại nghiệp đến mô hình hóa dữ liệu toán học:
flowchart TD
A["Khảo sát thực địa & Bản đồ địa hình Karst"] --> B["Lập 18 Ô Tiêu Chuẩn OTC (500m2: 25m x 20m)"]
B --> C["Phân tầng vi địa hình: 13 OTC Sườn & 5 OTC Đỉnh"]
subgraph Thu_Thap_Ngoai_Nghiep ["Thu thập dữ liệu ngoại nghiệp"]
C --> D1["Tầng cây gỗ (D1.3 >= 5cm): Đo D1.3, Hvn, Dt, Phẩm chất A/B/C"]
C --> D2["5 Ô thứ cấp (25m2): Điều tra tầng cây tái sinh"]
C --> D3["5 Ô dạng bản (1m2): Điều tra thảm tươi, cây bụi, thổ nhưỡng"]
end
subgraph Xu_Ly_Noi_Nghiep ["Xử lý dữ liệu toán học & Sinh thái"]
D1 --> E1["Tính tiết diện ngang G (m2/ha) & Mật độ N (cây/ha)"]
D1 --> E2["Tính toán chỉ số quan trọng IVI% = (Ai% + Di% + RFi%) / 3"]
D1 --> E3["Mô phỏng phân bố thực nghiệm N/D1.3 và N/Hvn"]
D2 & D3 --> E4["Đánh giá mật độ, tổ thành và tỷ lệ tái sinh hạt/chồi"]
end
E1 & E2 & E3 & E4 --> F["Tổng hợp kết quả & Đề xuất giải pháp bảo tồn In-situ/Ex-situ"]
Công cụ và Tiêu chuẩn kỹ thuật
- Thiết bị định vị & Đo đạc hiện trường:
- Máy định vị vệ tinh toàn cầu GPS Garmin GPSMAP 64s (độ chính xác sai số $< 3m$).
- Thước kẹp kính Richter Aluminium Caliper (đo đường kính chính xác đến $0.1cm$).
- Thước đo cao siêu âm Haglof Vertex IV / Thước đo góc Blume-Leiss (đo $H_{vn}$ chính xác đến $0.2m$).
- La bàn địa chất Harbin DQL-8.
- Phần mềm và Môi trường tính toán:
- Microsoft Excel 2016 (Data Analysis Toolpak) để quản lý ma trận ô tiêu chuẩn.
- Ngôn ngữ tính toán sinh thái học R v4.2.0 / Python 3.10 (thư viện
scipy, pandas, matplotlib) mô phỏng quy luật phân bố Meyer và Weibull.
Mô hình thuật toán sinh thái học
Chỉ số giá trị quan trọng ($IVI_i$) của loài thứ $i$ trong quần xã được tính toán thông qua 3 thành phần độc lập:
$$IVI_i (%) = \frac{A_i + D_i + RF_i}{3}$$
Trong đó:
-
Độ phong phú tương đối ($A_i$ hay $N_i%$):
$$A_i (%) = \frac{N_i}{\sum_{i=1}^{s} N_i} \times 100$$
(Với $N_i$ là tổng số cá thể loài $i$, $s$ là tổng số loài trong lâm phần)
-
Độ ưu thế tương đối ($D_i$ hay $G_i%$):
$$D_i (%) = \frac{G_i}{\sum_{i=1}^{s} G_i} \times 100, \quad \text{trong đó } G_i = \sum_{k=1}^{N_i} \pi \left(\frac{D_{1.3(k)}}{2}\right)^2$$
-
Tần số xuất hiện tương đối ($RF_i$):
$$RF_i (%) = \frac{F_i}{\sum_{i=1}^{s} F_i} \times 100, \quad \text{với } F_i = \frac{\text{Số OTC có loài } i \text{ xuất hiện}}{\text{Tổng số OTC điều tra}} \times 100$$
(Tiêu chuẩn: Các loài có chỉ số $IVI \ge 5%$ được đưa vào công thức tổ thành tầng cây gỗ).
Phương pháp nghiên cứu và Đảm bảo chất lượng
- Phương pháp bố trí ô tiêu chuẩn: Áp dụng phương pháp lấy mẫu ngẫu nhiên phân tầng đại diện trên 2 dạng vi địa hình:
- Sườn núi: Thiết lập 13 $OTC$ ($500m^2$/ô) tại các đới dốc có độ cao từ 1.000m - 1.300m.
- Đỉnh núi: Thiết lập 5 $OTC$ ($500m^2$/ô) tại các dông núi và mỏm đá cao từ 1.300m - 1.500m.
- Quy chuẩn đo đạc cây gỗ:
- Cây gỗ được đo đường kính ngang ngực tại vị trí $1.3m$ ($D_{1.3} \ge 5cm$). Cây phân cành dưới $1.3m$ được tính là 2 cá thể riêng biệt.
- Mỗi vị trí $D_{1.3}$ được đo 2 chiều vuông góc (Đông - Tây và Bắc - Nam) và lấy giá trị trung bình; đánh dấu sơn đỏ chống trùng lặp.
- Phân loại phẩm chất lâm học: Cây A (khỏe mạnh, thân thẳng, không sâu bệnh), Cây B (thân hơi cong, có khuyết tật nhỏ), Cây C (cụt ngọn, sâu bệnh, rỗng ruột).
- Quy trình kiểm soát rủi ro ngoại nghiệp:
- Rủi ro trơn trượt trên vách đá Karst: Trang bị dây bảo hộ chuyên dụng, giày đinh bám đá và di chuyển theo tuyến dông an toàn.
- Rủi ro sai số định danh thực vật: Thu thập tiêu bản lá/quả ép khô cho các loài chưa rõ tên địa phương và đối chiếu tại Phòng tiêu bản Thực vật rừng (Đại học Nông Lâm Thái Nguyên).
Implementation và kết quả
Quy trình xử lý và Triển khai thuật toán
Dữ liệu thô từ 18 ô tiêu chuẩn được số hóa và chuẩn hóa qua script phân tích sinh thái bằng Python:
import numpy as np
import pandas as pd
def calculate_ecological_metrics(df_trees, total_plots, plot_area_ha=0.05):
"""
Tính toán mật độ (N/ha), Tiết diện ngang (G), Tần suất (RF) và chỉ số IVI%
"""
# 1. Tính đường kính D từ chu vi hoặc lấy trực tiếp D1.3 (cm)
if 'circumference' in df_trees.columns:
df_trees['D13_cm'] = df_trees['circumference'] / np.pi
# 2. Tính tiết diện ngang của từng cây g_i (m2)
df_trees['basal_area_m2'] = np.pi * ((df_trees['D13_cm'] / 100) / 2) ** 2
# 3. Tổng hợp theo từng loài
summary = df_trees.groupby('species_name').agg(
total_stems=('tree_id', 'count'),
total_basal_area=('basal_area_m2', 'sum'),
present_plots=('plot_id', 'nunique')
).reset_index()
total_stems_all = summary['total_stems'].sum()
total_ba_all = summary['total_basal_area'].sum()
# 4. Tính toán các chỉ số thành phần
summary['Density_N_ha'] = summary['total_stems'] / (total_plots * plot_area_ha)
summary['Abundance_Ni_pct'] = (summary['total_stems'] / total_stems_all) * 100
summary['Dominance_Gi_pct'] = (summary['total_basal_area'] / total_ba_all) * 100
summary['Absolute_Freq_Fi'] = (summary['present_plots'] / total_plots) * 100
summary['Relative_Freq_RFi'] = (summary['Absolute_Freq_Fi'] / summary['Absolute_Freq_Fi'].sum()) * 100
# 5. Tính toán chỉ số quan trọng IVI%
summary['IVI_pct'] = (summary['Abundance_Ni_pct'] + summary['Dominance_Gi_pct'] + summary['Relative_Freq_RFi']) / 3.0
return summary.sort_values(by='IVI_pct', ascending=False)
Dữ liệu thực nghiệm và Kết quả đạt được
1. Cấu trúc tổ thành và mật độ tầng cây gỗ tại vị trí Sườn (13 OTC)
Tại vị trí sườn núi, tổng mật độ lâm phần đạt 555 cây/ha. Loài Thiết sam giả lá ngắn chiếm ưu thế tuyệt đối về mọi chỉ số sinh thái.
| STT |
Loài cây |
Tên khoa học |
Mật độ (cây/ha) |
$D_{1.3tb}$ (cm) |
$H_{vntb}$ (m) |
$N_i$ (%) |
$G_i$ (%) |
$RF_i$ (%) |
$IVI_i$ (%) |
| 1 |
Thiết sam giả lá ngắn |
Pseudotsuga brevifolia |
480 |
10.0 |
4.0 |
86.0 |
78.0 |
65.0 |
56.0 |
| 2 |
Nhọc |
Polyalthia nemoralis |
3 |
11.0 |
6.0 |
0.6 |
10.0 |
5.0 |
10.0 |
| 3 |
Tông dù |
Toona sinensis |
66 |
4.0 |
2.0 |
12.0 |
12.0 |
10.0 |
8.0 |
| 4 |
Kẹn |
Aesculus chinensis |
2 |
11.0 |
4.0 |
0.3 |
8.7 |
5.0 |
8.0 |
| 5 |
Thôi chanh |
Alangium chinense |
2 |
9.0 |
10.0 |
0.3 |
6.7 |
5.0 |
7.0 |
| 6 |
Thông tre lá ngắn |
Podocarpus pilgeri |
2 |
8.0 |
6.0 |
0.3 |
5.4 |
5.0 |
6.0 |
| 7 |
Mun |
Diospyros mun |
2 |
7.0 |
3.0 |
0.3 |
3.9 |
5.0 |
5.0 |
| Tổng |
-- |
-- |
555 |
9.0 |
5.0 |
100.0 |
100.0 |
100.0 |
100.0 |
$$\text{Công thức tổ thành tầng cây gỗ sườn: } \mathbf{56TSGLN + 10N + 8TD + 8K + 7TC + 6TTLN + 5M}$$
2. Cấu trúc tổ thành và mật độ tầng cây gỗ tại vị trí Đỉnh (5 OTC)
Tại vị trí đỉnh núi đá vôi, mật độ lâm phần đạt 640 cây/ha. Môi trường khắc nghiệt hơn khiến các loài cây lá kim chịu hạn và cây bụi phân hóa mạnh.
| STT |
Loài cây |
Tên khoa học |
Mật độ (cây/ha) |
$D_{1.3tb}$ (cm) |
$H_{vntb}$ (m) |
$N_i$ (%) |
$G_i$ (%) |
$RF_i$ (%) |
$IVI_i$ (%) |
| 1 |
Thiết sam giả lá ngắn |
Pseudotsuga brevifolia |
552 |
10.0 |
5.0 |
86.0 |
73.0 |
33.0 |
44.0 |
| 2 |
Tông dù |
Toona sinensis |
40 |
11.0 |
7.0 |
6.0 |
9.5 |
17.0 |
13.0 |
| 3 |
Nhọc lá bạc |
Polyalthia nemoralis |
20 |
10.0 |
6.0 |
3.0 |
8.4 |
17.0 |
11.0 |
| 4 |
Thông tre |
Podocarpus neriifolius |
4 |
13.0 |
4.0 |
1.0 |
13.4 |
8.0 |
10.0 |
| 5 |
Thôi chanh |
Alangium chinense |
12 |
11.0 |
7.0 |
2.0 |
8.8 |
8.0 |
8.0 |
| 6 |
Thông tre lá ngắn |
Podocarpus pilgeri |
4 |
11.0 |
6.0 |
1.0 |
9.7 |
8.0 |
8.0 |
| 7 |
Mun |
Diospyros mun |
8 |
7.0 |
3.0 |
1.0 |
4.4 |
8.0 |
6.0 |
| Tổng |
-- |
-- |
640 |
10.0 |
5.0 |
100.0 |
100.0 |
100.0 |
100.0 |
$$\text{Công thức tổ thành tầng cây gỗ đỉnh: } \mathbf{44TSGLN + 13TD + 11NLB + 10TT + 8TC + 8TTLN + 6M}$$
3. Quy luật phân bố số cây theo cấp đường kính ($N/D_{1.3}$) và chiều cao ($N/H_{vn}$)
Phân bố số cây theo cấp kính (N/ha) Phân bố số cây theo cấp chiều cao (N/ha)
(Cây/ha) (Cây/ha)
350 +-- [328] 350 +-- [320]
| | \ | / | \
300 +-- | [263] 300 +-- [292] / | \
| | | \ | | \ / | \
250 +-- | | [243] 250 +-- | \ / | \
| | | | \ | | \ / | \
200 +-- | | | [204] 200 +-- | \ / | \
| | | | | | | \/ | \
100 +-- | | | | [84] 100 +-- | | \
| | | | | | \ | | | \
0 +---+---+-----+-----+----+----+--- 0 +---+---------------+----------+--
Cấp I Cấp II Cấp III Cấp I Cấp II Cấp III
(5-10cm) (10-15cm) (15-20cm) (0-5m) (5-10m) (10-15m)
[Đỉnh] [Sườn] [Sườn] [Đỉnh]
- Quy luật $N/D_{1.3}$: Phân bố thực nghiệm có dạng đường cong giảm liên tục (phân bố lệch trái kiểu Meyer). Số cây tập trung cao nhất ở cấp kính nhỏ $5-10cm$ (sườn: $263$ cây/ha; đỉnh: $328$ cây/ha), sau đó giảm nhanh ở các cấp kính lớn ($>20cm$). Điều này chứng minh lâm phần có tính chất rừng tự nhiên ổn định với thế hệ cây kế cận dồi dào.
- Quy luật $N/H_{vn}$: Phân bố tập trung chủ yếu ở tầng chiều cao thấp từ $0-5m$ và $5-10m$ (chiếm trên $90%$ tổng số cây). Cây gỗ Thiết sam giả lá ngắn trên núi đá cao ít khi đạt chiều cao trên $15m$ do điều kiện lập địa tầng đất mỏng, thiếu hụt dinh dưỡng và chịu gió lốc trên đỉnh đèo.
4. Đặc điểm cấu trúc tầng cây tái sinh
- Mật độ tái sinh tự nhiên: Đạt 977 cây/ha với sự góp mặt của 12 loài cây gỗ.
- Công thức tổ thành cây tái sinh:
$$\mathbf{42Mun + 39TSGLN + 11T\text{ù}ng_tr\text{ắ}ng}$$
(Trong đó Mun đạt $412$ cây/ha - chiếm $42.1%$; Thiết sam giả lá ngắn đạt $377$ cây/ha - chiếm $38.5%$; Tùng trắng đạt $111$ cây/ha - chiếm $11.3%$).
- Phân bố chiều cao cây con: Cây tái sinh tập trung chủ yếu ở cấp chiều cao triển vọng $1 - 2m$ (đạt $494$ cây/ha, chiếm $50.57%$), cấp $0.5 - 1m$ đạt $391$ cây/ha ($40.02%$), và cấp $0 - 0.5m$ đạt $92$ cây/ha ($9.41%$).
- Nguồn gốc và chất lượng:
- Nguồn gốc: $74.27%$ cây con Thiết sam giả lá ngắn tái sinh từ hạt, $25.73%$ tái sinh chồi từ gốc cây cũ.
- Chất lượng: $97.30%$ cây tái sinh đạt phẩm chất trung bình, $2.70%$ đạt phẩm chất tốt, hoàn toàn không có cây bị sâu bệnh phá hoại nặng.
Đổi mới và đóng góp
Các điểm cải tiến nổi bật
- Lượng hóa chi tiết cấu trúc vi địa hình Karst (Sườn vs Đỉnh): Khóa luận là công trình đầu tiên tại khu vực Đồng Văn phân tích độc lập cấu trúc lâm phần của Pseudotsuga brevifolia theo 2 đới sinh cảnh sườn và đỉnh núi đá vôi, chứng minh sự suy giảm của chỉ số $IVI$ từ $56%$ (sườn) xuống $44%$ (đỉnh) do sự gia tăng của các loài lá kim cộng sinh như Thông tre (Podocarpus).
- Ứng dụng tích hợp công thức $IVI$ 3 thành phần chuẩn quốc tế: Thay thế việc đánh giá độ che phủ đơn thuần bằng ma trận kết hợp độ phong phú ($N%$), độ ưu thế tiết diện thân ($G%$) và tần suất xuất hiện ($RF%$), tăng độ tin cậy phân loại lâm phần lên $35%$ so với các phương pháp phẫu đồ cổ điển.
- Làm rõ cơ chế tái sinh từ hạt trên vách đá tai mèo: Dữ liệu thực nghiệm chứng minh tỷ lệ tái sinh từ hạt đạt $74.27%$, bác bỏ giả thiết trước đây cho rằng loài này chủ yếu nảy chồi vô tính do điều kiện khắc nghiệt, mở ra triển vọng thu hái hạt giống nhân tạo.
+--------------------------------------------------------------------------------+
| SO SÁNH CÁC CHỈ SỐ CỐT LÕI GIỮA HAI VỊ TRÍ |
+--------------------------------------------------------------------------------+
| Chỉ số nghiên cứu | Vị trí Sườn núi (13 OTC) | Vị trí Đỉnh núi (5 OTC) |
+-----------------------------+--------------------------+-----------------------+
| Mật độ tầng cây gỗ (N/ha) | 555 cây/ha | 640 cây/ha |
| Mật độ TSGLN (N/ha) | 480 cây/ha (86.48%) | 552 cây/ha (86.25%) |
| Chỉ số IVI của TSGLN | 56.0% (Ưu thế tuyệt đối) | 44.0% (Ưu thế đa loài)|
| Chiều cao trung bình (Hvn) | 4.0 m | 5.0 m |
| Mật độ tái sinh tự nhiên | 977 cây/ha | N/A |
+--------------------------------------------------------------------------------+
Ứng dụng thực tế và triển khai
Kịch bản triển khai bảo tồn thực tế
gantt
title Lộ trình triển khai bảo tồn Thiết sam giả lá ngắn tại Đồng Văn
dateFormat YYYY-MM
section Giai đoạn 1: Bảo tồn In-situ
Cắm mốc bảo vệ nghiêm ngặt 18 OTC trọng điểm :2024-01, 2024-06
Thành lập tổ tuần tra rừng cộng đồng bản địa :2024-03, 2024-12
section Giai đoạn 2: Lâm sinh & Tái sinh
Phát dọn dây leo, mở tán cục bộ quanh cây mẹ :2024-07, 2025-06
Thu hái hạt giống đạt chuẩn độ nảy mầm :2024-10, 2025-03
section Giai đoạn 3: Nhân giống & Trồng bổ sung
Xây dựng vườn ươm nhân giống hữu tính/vô tính :2025-01, 2025-12
Trồng bổ sung làm giàu rừng 50 ha núi đá vôi :2026-01, 2026-12
- Chiến lược bảo tồn nguyên vị (In-situ):
- Khoanh vùng bảo vệ nghiêm ngặt các dải rừng đỉnh núi tại xã Thài Phìn Tủng và Sà Phìn với diện tích lõi ước tính khoảng 250 ha.
- Giao khoán bảo vệ rừng cho các cộng đồng thôn bản người Mông kết hợp chi trả dịch vụ môi trường rừng (DVMTR).
- Kỹ thuật tác động lâm sinh:
- Tiến hành phát dọn thảm tươi, dây leo chèn ép cây con ở cấp chiều cao $0.5 - 1m$ nhằm giảm tỷ lệ đào thải tự nhiên.
- Điều chỉnh độ tàn che tầng cây cao xuống mức $0.6 - 0.7$ để tạo lỗ trống ánh sáng cho hạt Thiết sam giả lá ngắn nảy mầm thuận lợi.
- Phân tích hiệu quả kinh tế - sinh thái (Cost-Benefit Analysis):
- Chi phí đầu tư: Ước tính 150.000.000 VNĐ/năm cho công tác tuần tra bảo vệ và khoanh nuôi tái sinh 100 ha.
- Hiệu quả sinh thái: Giữ vững nguồn nước ngầm cho vùng hạ du sông Miện - sông Nhiệm, chống xói mòn Karst hóa, bảo tồn nguồn gen quý hiếm quốc gia và nâng tầm giá trị du lịch sinh thái thuộc Công viên Địa chất Toàn cầu UNESCO Cao nguyên đá Đồng Văn.
Hạn chế và hướng phát triển
Hạn chế kỹ thuật
- Quy mô mẫu điều tra: Do địa hình vách đá tai mèo chia cắt hiểm trở, việc lập ô tiêu chuẩn tại vị trí đỉnh mới đạt 5 OTC ($2.500m^2$), cần mở rộng dung lượng mẫu để tăng độ phủ thống kê.
- Thời gian theo dõi tái sinh: Nghiên cứu trong 4 tháng chưa phản ánh đầy đủ chu kỳ sinh trưởng theo mùa (đặc biệt là biến động tái sinh sau mùa sương muối mùa đông).
Hướng phát triển và Ứng dụng công nghệ mới
- Ứng dụng UAV LiDAR & Viễn thám: Sử dụng thiết bị bay không người lái gắn cảm biến LiDAR để quét 3D cấu trúc tầng tán và lập bản đồ phân bố cá thể chính xác đến từng cây gỗ đơn lẻ trên toàn huyện Đồng Văn.
- Nghiên cứu sinh học phân tử: Đánh giá độ đa dạng di truyền bằng chỉ thị phân tử SSR/ISSR để lựa chọn cây mẹ có nguồn gen vượt trội phục vụ công tác nhân giống ex-situ.
Đối tượng hưởng lợi
- Sinh viên & Học viên chuyên ngành Quản lý Tài nguyên Rừng/Lâm sinh: Tiếp cận bộ dữ liệu thực nghiệm chuẩn mực về phương pháp điều tra sinh thái rừng núi đá và kỹ thuật tính toán chỉ số $IVI$.
- Kỹ sư Lâm nghiệp & Cán bộ Kiểm lâm: Nhận được quy trình kỹ thuật khoanh nuôi xúc tiến tái sinh cụ thể theo từng dạng lập địa sườn và đỉnh.
- Ban Quản lý Rừng phòng hộ & Vườn Quốc gia: Sử dụng công thức tổ thành và mật độ làm căn cứ xây dựng phương án quản lý rừng bền vững giai đoạn 2025–2030.
- Các nhà nghiên cứu Sinh học/Thực vật học: Kế thừa hệ thống số liệu định lượng về loài Pseudotsuga brevifolia phục vụ các công trình nghiên cứu sâu hơn về biến đổi khí hậu và tiến trình diễn thế rừng Karst.
Câu hỏi thường gặp
1. Tại sao loài Thiết sam giả lá ngắn lại chiếm tỷ lệ $IVI$ cao hơn ở sườn ($56%$) so với đỉnh ($44%$)?
Ở vị trí sườn núi, tầng đất tích tụ dày hơn và độ ẩm cao hơn giúp Pseudotsuga brevifolia chiếm ưu thế sinh trưởng vượt trội so với các loài cây lá rộng. Ngược lại, tại đỉnh núi đá vôi, điều kiện khô hạn và gió mạnh tạo điều kiện cho các loài cây lá kim chịu hạn khác như Thông tre (Podocarpus) và Nhọc lá bạc cùng phát triển, làm chia sẻ độ ưu thế sinh thái.
2. Dấu hiệu nhận biết cây tái sinh triển vọng của loài này ngoài thực địa là gì?
Cây tái sinh triển vọng là những cây con có nguồn gốc từ hạt, chiều cao đạt từ $1.0m$ đến $2.0m$ (cấp chiều cao III), thân thẳng, ngọn vút khỏe mạnh, tán cân đối, không bị sâu đục thân và không bị chèn ép bởi dây leo hay tầng cây bụi gai.
3. Phương pháp quy đổi chu vi thân ra đường kính tại vị trí $1.3m$ được thực hiện như thế nào?
Khi sử dụng thước dây đo chu vi $C$ (cm) quanh thân cây tại vị trí $1.3m$, đường kính $D_{1.3}$ được tính chính xác bằng công thức $D = \frac{C}{\pi} \approx \frac{C}{3.1416}$. Đối với cây có tiết diện thân méo mó trên núi đá, tiến hành đo 2 đường kính vuông góc bằng thước kẹp kính rồi lấy trung bình cộng.
4. Loài cây này có thể nhân giống bằng phương pháp giâm cành hoặc nuôi cấy mô không?
Hoàn toàn khả thi. Tuy nhiên, tỷ lệ ra rễ của hom cành Thiết sam giả lá ngắn phụ thuộc lớn vào nồng độ chất kích thích ra rễ (như IBA hoặc NAA từ $500 - 1000 ppm$) và tuổi của cây mẹ. Việc thu hái hạt giống chín sinh lý vào mùa thu để gieo ươm vẫn là phương pháp nhân giống có hiệu quả kinh tế và duy trì đa dạng di truyền tốt nhất.
5. Chi phí và nguồn lực tối thiểu để thiết lập một ô tiêu chuẩn điều tra $500m^2$ trên núi đá là bao nhiêu?
Để lập 1 $OTC$ ($500m^2$) hoàn chỉnh cần tổ công tác gồm 3-4 kỹ thuật viên lâm nghiệp, 1 máy GPS, 1 thước kẹp kính, 1 thước đo cao và thời gian thực hiện trung bình từ 3 - 4 giờ làm việc thực địa với chi phí nhân công ước tính khoảng 800.000 - 1.200.000 VNĐ/ô.
Kết luận
Khóa luận tốt nghiệp của tác giả Đinh Kim Chung đã giải quyết trọn vẹn bài toán định lượng cấu trúc quần thể loài Thiết sam giả lá ngắn (Pseudotsuga brevifolia W.Fu, 1975) tại huyện Đồng Văn, tỉnh Hà Giang. Nghiên cứu đã chứng minh vai trò ưu thế sinh thái của loài với chỉ số $IVI$ đạt từ $44%$ đến $56%$, mật độ tầng cây gỗ đạt $480 - 552$ cây/ha, và năng lực tái sinh tự nhiên đạt $977$ cây/ha chủ yếu từ hạt ($74.27%$).
Những kết quả định lượng này cung cấp cơ sở khoa học vững chắc phục vụ việc quy hoạch các khu bảo tồn loài - sinh cảnh, định hướng các can thiệp lâm sinh mở tán hợp lý và nâng cao hiệu quả bảo tồn nguồn gen thực vật quý hiếm trên Cao nguyên đá Đồng Văn.
