Giới thiệu dự án

Bối cảnh và thực trạng ngành lâm nghiệp

Rừng tự nhiên nhiệt đới tại Việt Nam đóng vai trò then chốt trong việc bảo tồn nguồn gen, điều hòa nguồn nước, chống xói mòn và cân bằng hệ sinh thái. Tuy nhiên, dữ liệu thống kê lâm nghiệp giai đoạn lịch sử cho thấy tỷ lệ che phủ rừng toàn quốc đã suy giảm nghiêm trọng từ 43% (năm 1943) xuống còn 26% (năm 1993) do sức ép từ chiến tranh, khai thác chọn kiệt quệ và tập quán du canh nương rẫy. Tại huyện Văn Chấn, tỉnh Yên Bái – khu vực cửa ngõ vùng Tây Bắc – diện tích rừng tự nhiên bị chia cắt và thoái hóa mạnh. Xã Sùng Đô có tổng diện tích tự nhiên 4.081,50 ha, trong đó đất quy hoạch lâm nghiệp chiếm 2.262,42 ha (55,43%). Trạng thái rừng phục hồi sau nương rẫy (trạng thái IIa) chiếm diện tích lớn với 833,30 ha, đóng vai trò xung yếu trong phòng hộ đầu nguồn và sinh kế của đồng bào dân tộc bản địa.

       CƠ CẤU ĐẤT LÂM NGHIỆP XÃ SÙNG ĐÔ (TỔNG: 2.262,42 HA)

Vấn đề nghiên cứu (Problem Statement)

Rừng phục hồi trạng thái IIa tại vùng cao Tây Bắc đang đối mặt với các rào cản sinh thái nghiêm trọng:

  • Cấu trúc tầng tán bị xáo trộn, trữ lượng cây gỗ thương phẩm thấp, chủ yếu là các loài cây tiên phong ưa sáng có giá trị kinh tế không cao.
  • Quá trình tái sinh tự nhiên diễn ra không đồng đều do sự cạnh tranh khốc liệt về ánh sáng và không gian dinh dưỡng từ tầng cây bụi, thảm tươi dày đặc.
  • Hoạt động chăn thả gia súc tự do (trâu, bò) và áp lực thu hái gỗ củi làm suy giảm tỷ lệ cây con chất lượng cao mọc từ hạt.
  • Thiếu các dẫn liệu định lượng chính xác về cấu trúc tổ thành ($IV%$), mật độ phân bố và động thái phát triển chiều cao làm cơ sở khoa học để thiết kế các biện pháp kỹ thuật lâm sinh (silvicultural treatments) thích ứng theo từng đai địa hình.

Mục tiêu đề tài

  1. Định lượng cấu trúc tầng cây gỗ cao: Xác định mật độ ($N/ha$), đường kính ngang ngực ($D_{1.3}$), chiều cao vút ngọn ($H_{vn}$), độ tàn che và chỉ số quan trọng ($IV%$) của các loài cây gỗ tại các vị trí chân, sườn, đỉnh dông.
  2. Đánh giá toàn diện đặc tính tái sinh tự nhiên: Phân tích mật độ tầng cây tái sinh, tỷ lệ cây triển vọng, chất lượng (Tốt - Trung bình - Xấu), nguồn gốc phát sinh (Hạt vs. Chồi) và cấu trúc phân bố không gian.
  3. Xác lập giải pháp lâm sinh ứng dụng: Đề xuất quy trình tác động cơ giới - sinh học (khoanh nuôi xúc tiến tái sinh, tra dặm cây bản địa mục đích, điều tiết tán che) nhằm rút ngắn chu kỳ diễn thế phục hồi rừng đỉnh cực.

Phạm vi và giới hạn nghiên cứu

  • Địa bàn: Xã Sùng Đô, huyện Văn Chấn, tỉnh Yên Bái (độ cao trung bình 1.015 m, độ dốc trung bình $52,5^\circ$).
  • Đối tượng: Trạng thái rừng phục hồi sau nương rẫy (trạng thái IIa) phân bố tại 3 điều kiện lập địa/địa hình: Chân đồi, Sườn dốc và Đỉnh dông.
  • Quy mô mẫu: 9 ô tiêu chuẩn (OTC) điển hình diện tích $2.000\text{ m}^2$ ($40\text{ m} \times 50\text{ m}$) và 45 ô tiêu chuẩn dạng bản (ODB) diện tích $25\text{ m}^2$ ($5\text{ m} \times 5\text{ m}$).

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

Phân tích hiện trạng và giải pháp lâm sinh

Phương thức phục hồi Ưu điểm Nhược điểm Đánh giá tính khả thi kinh tế - sinh thái
Trồng rừng mới thuần loại Tăng trưởng sinh khối nhanh, dễ cơ giới hóa chăm sóc. Chi phí đầu tư cao ($25 - 35\text{ triệu VNĐ/ha}$), cấu trúc đơn tầng kém bền vững, nguy cơ cháy rừng và xói mòn đất dốc cao. Thấp đối với vùng phòng hộ đầu nguồn có độ dốc $>50^\circ$.
Khoanh nuôi đóng cửa rừng thuần túy Chi phí đầu tư thấp ($0,5 - 1\text{ triệu VNĐ/ha/năm}$). Diễn thế chậm ($>20 - 30\text{ năm}$), tầng thảm tươi lấn át cây tái sinh mục đích, thành phần loài dễ lệch về cây gỗ tạp. Trung bình; chỉ phù hợp nơi có nguồn giống mẹ dồi dào và mật độ cây tái sinh tốt.
Khoanh nuôi xúc tiến tái sinh có làm giàu (Đề xuất) Tận dụng tối đa cây tái sinh tự nhiên ($53,08%$ cây triển vọng), bổ sung cây bản địa có giá trị, rút ngắn thời gian thành rừng $40 - 50%$. Đòi hỏi điều tra thực địa chi tiết và kỹ thuật điều tiết tán che, phát luỗng chính xác. Rất cao; tối ưu hóa chi phí đầu tư, đảm bảo tính đa dạng sinh học và phòng hộ bền vững.
                       MA TRẬN ƯU TIÊN YÊU CẦU LÂM SINH (MoSCoW)

Phương pháp luận và thiết kế hệ sinh thái điều tra

Phương pháp tiếp cận dựa trên lý thuyết sinh thái phát sinh quần thể của Thái Văn Trừng (1978) kết hợp phương pháp dãy phát triển tự nhiên (dùng không gian thay thế thời gian).

                      MÔ HÌNH THU THẬP VÀ XỬ LÝ DỮ LIỆU ĐIỀU TRA

Công cụ và thiết bị kỹ thuật sử dụng:

  • Đo trắc địa và lâm học: Thước kẹp kính Haglöf Mantax Precision (độ chính xác $1\text{ mm}$), thước đo cao Laser Vertex IV/Hypsometer (độ chính xác $0,1\text{ m}$), thước dây bọc sợi thủy tinh Yamayo $50\text{ m}$, máy định vị GPS Garmin GPSMAP 64s.
  • Phần mềm xử lý thống kê sinh thái: R v4.3.2 (packages: vegan, spatstat), Python v3.10 (thư viện: pandas, numpy, scipy).

Implementation và kết quả

Các công thức và thuật toán định lượng sinh thái

1. Chỉ số giá trị quan trọng (Important Value Index - $IV%$ theo Daniel Marmillod):

$$IV_i% = \frac{N_i% + G_i%}{2}$$ Trong đó:

  • $N_i% = \frac{n_i}{N} \times 100$: Tỷ lệ phần trăm số cá thể của loài $i$ so với tổng số cá thể lâm phần.
  • $G_i% = \frac{g_i}{G} \times 100$: Tỷ lệ phần trăm tiết diện ngang thân cây ($g_i = \sum \frac{\pi D_{1.3}^2}{4}$) của loài $i$ so với tổng tiết diện ngang lâm phần.
  • Ngưỡng sinh thái: Loài có $IV% \ge 5%$ được xác định là loài tham gia vào công thức tổ thành tầng cây cao.

2. Kiểm định mô hình phân bố không gian mặt phẳng ngang (Clark & Evans Index):

$$R = \frac{\bar{r}}{E(r)} = \frac{\bar{r}}{\frac{1}{2\sqrt{\lambda}}}; \quad U = \frac{\bar{r} - E(r)}{\sigma_r} = \frac{\bar{r} - \frac{1}{2\sqrt{\lambda}}}{\frac{0,26136}{\sqrt{n \cdot \lambda}}}$$ Trong đó:

  • $\bar{r}$: Khoảng cách trung bình thực tế giữa cá thể tái sinh đến cá thể lân cận gần nhất ($m$).
  • $\lambda = \frac{N}{S}$: Mật độ cá thể trên một đơn vị diện tích ($\text{cây/m}^2$).
  • $n$: Tổng số lần đo khoảng cách.
  • Quy tắc đánh giá:
    • $-1,96 \le U \le 1,96$: Phân bố ngẫu nhiên (Poisson distribution, $p > 0,05$).
    • $U > 1,96$: Phân bố đều (Uniform distribution).
    • $U < -1,96$: Phân bố cụm/đám (Aggregated/Clumped distribution).

Script Python xử lý tự động hóa chỉ số $IV%$ và kiểm định Clark-Evans:

import numpy as np
import pandas as pd

def calculate_importance_value(df_trees: pd.DataFrame) -> pd.DataFrame:
    """
    Tính chỉ số quan trọng IV% cho từng loài cây gỗ theo Daniel Marmillod.
    Yêu cầu dataframe có các cột: ['Species', 'DBH_cm']
    """
    df = df_trees.copy()
    # Tính tiết diện ngang g = pi * (D/200)^2 (đơn vị: m2)
    df['Basal_Area_m2'] = np.pi * (df['DBH_cm'] / 200.0) ** 2
    
    total_trees = len(df)
    total_ba = df['Basal_Area_m2'].sum()
    
    grouped = df.groupby('Species').agg(
        n_trees=('DBH_cm', 'count'),
        sum_ba=('Basal_Area_m2', 'sum')
    ).reset_index()
    
    grouped['N_pct'] = (grouped['n_trees'] / total_trees) * 100.0
    grouped['G_pct'] = (grouped['sum_ba'] / total_ba) * 100.0
    grouped['IV_pct'] = (grouped['N_pct'] + grouped['G_pct']) / 2.0
    
    return grouped.sort_values(by='IV_pct', ascending=False)

def clark_evans_test(distances: np.ndarray, density_lambda: float) -> dict:
    """
    Thực hiện kiểm định phân bố không gian Clark & Evans (1954).
    """
    n = len(distances)
    r_bar = np.mean(distances)
    e_r = 1.0 / (2.0 * np.sqrt(density_lambda))
    sigma_r = 0.26136 / np.sqrt(n * density_lambda)
    u_stat = (r_bar - e_r) / sigma_r
    
    if u_stat < -1.96:
        pattern = "Phân bố cụm (Clustered)"
    elif u_stat > 1.96:
        pattern = "Phân bố đều (Regular)"
    else:
        pattern = "Phân bố ngẫu nhiên (Random)"
        
    return {
        "Mean_Distance_m": float(r_bar),
        "Expected_Distance_m": float(e_r),
        "R_ratio": float(r_bar / e_r),
        "U_statistic": float(u_stat),
        "Spatial_Pattern": pattern
    }

Kết quả nghiên cứu thực nghiệm

                      MẬT ĐỘ TẦNG CÂY GỖ VÀ CÂY TÁI SINH THEO ĐỊA HÌNH
  Cây/ha
         [Vị trí Chân đồi]           [Vị trí Sườn dốc]        [Vị trí Đỉnh dông]
         

1. Cấu trúc tầng cây gỗ cao

  • Mật độ tầng cao: Dao động từ $111 - 117\text{ cây/ha}$ (Chân: $117\text{ cây/ha}$; Sườn: $113\text{ cây/ha}$; Đỉnh: $111\text{ cây/ha}$). Xu hướng mật độ giảm dần từ chân lên đỉnh do tầng đất mỏng và khả năng giữ ẩm kém ở vị trí đỉnh.
  • Tổ thành tầng cây gỗ: Ghi nhận từ 10 đến 13 loài/OTC. 11 loài ưu thế tham gia vào công thức tổ thành gồm: Giẻ (Castanopsis sieboldii), Sâng (Pometia pinnata), Re (Cinnamomum sp.), Kháo (Machilus sp.), Sồi (Lithocarpus sp.), Mản (Helicia cochinchinensis), Vạng (Endospermum chinense), Xoan nhừ (Choerospondias axillaris), Ba gạc (Rauvolfia verticillata), Kẹn, Bồ đề tía.
  • Công thức tổ thành tổng hợp:
    • Chân: $40,62\text{LK} + 11,19\text{G} + 9,42\text{S} + 8,23\text{Sg} + 6,29\text{R} + 5,49\text{Khv}$
    • Sườn: $64,34\text{LK} + 18,02\text{G} + 9,24\text{Try} + 7,28\text{M} + 6,41\text{Sg} + 6,31\text{Kh} + 6,11\text{Bg} + 5,65\text{Sđ} + 5,32\text{V}$
    • Đỉnh: $45,32\text{LK} + 8,50\text{G} + 7,39\text{Xn} + 6,66\text{Sg} + 5,95\text{M} + 5,65\text{Ng} + 5,75\text{R} + 5,43\text{V}$

2. Đặc điểm định lượng tầng cây tái sinh tự nhiên

Vị trí lập địa Mật độ tái sinh ($N/\text{ha}$) Cây triển vọng ($N/\text{ha}$) Tỷ lệ triển vọng (%) Tỷ lệ Tốt (%) Tỷ lệ TB (%) Tỷ lệ Xấu (%) Nguồn gốc Hạt (%) Nguồn gốc Chồi (%)
Chân đồi 2.613 1.333 50,88% 48,28% 25,53% 26,20% 48,58% 51,42%
Sườn dốc 2.453 1.307 53,28% 44,42% 35,73% 19,85% 57,88% 42,12%
Đỉnh dông 2.907 1.600 55,07% 46,17% 35,12% 16,02% 62,84% 37,76%
Trung bình 2.658 1.413 53,08% 46,29% 32,13% 20,69% 56,43% 43,77%

3. Cấu trúc phân bố tầng tái sinh theo cấp chiều cao

                 PHÂN BỐ SỐ CÂY TÁI SINH THEO CẤP CHIỀU CAO (N = 2.658 CÂY/HA)
  • Phân tích động thái: Đồ thị phân bố chiều cao có dạng đường cong giảm cấp mạnh ($J$-shaped curve). Số lượng cây con tập trung áp đảo ở cấp $<0,5\text{ m}$ ($51,32%$), sau đó giảm nhanh khi đạt các cấp chiều cao $>1,5\text{ m}$. Hiện tượng này minh chứng cho sự chọn lọc đào thải tự nhiên khốc liệt khi cây tái sinh bắt đầu vươn lên cạnh tranh ánh sáng với tầng cây bụi, thảm tươi.

Đổi mới và đóng góp khoa học

  1. Khắc phục khoảng trống dữ liệu khu vực: Cung cấp bộ dữ liệu thực chứng hoàn chỉnh về rừng trạng thái IIa tại đai cao $>1.000\text{ m}$ huyện Văn Chấn, bổ sung cơ sở khoa học cho hệ thống phân loại thảm thực vật miền Bắc.
  2. Định lượng hóa năng lực tái sinh: Làm rõ tỷ lệ tái sinh từ hạt ($56,43%$) chiếm ưu thế so với chồi ($43,77%$), khẳng định quần xã có tiềm năng trẻ hóa tự nhiên lâu dài mà không bị hiện tượng thoái hóa gốc chồi như các khu vực khai thác chọn chu kỳ ngắn.
  3. So sánh tương quan với các công trình nghiên cứu tiền nhiệm:
                            SO SÁNH CÁC CHỈ SỐ SINH THÁI TÁI SINH

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

Quy trình kỹ thuật lâm sinh ứng dụng

                SƠ ĐỒ QUY TRÌNH KỸ THUẬT LÂM SINH TRẠNG THÁI IIA

Phân tích hiệu quả kinh tế - xã hội (Cost-Benefit & ROI)

  • Tiết kiệm ngân sách phục hồi rừng: Chi phí khoanh nuôi xúc tiến tái sinh có trồng làm giàu dao động từ $6 - 8\text{ triệu VNĐ/ha/5 năm}$, chỉ bằng $20 - 25%$ so với chi phí trồng rừng mới hoàn toàn ($30 - 35\text{ triệu VNĐ/ha}$).
  • Giá trị kinh tế lâu dài: Sau 15 - 20 năm, tầng cây gỗ mục đích đạt đường kính khai thác chọn ($D_{1.3} > 35\text{ cm}$), mang lại giá trị gỗ thương phẩm ước tính $120 - 150\text{ triệu VNĐ/ha}$ song song với nguồn thu từ lâm sản ngoài gỗ (vỏ quế, hạt giổi, quả trám).
  • Giá trị phòng hộ sinh thái: Giảm $65 - 80%$ lượng đất trôi xói mòn trên sườn dốc $>50^\circ$, bảo vệ an toàn các công trình thủy lợi kênh mương tại các thôn Nà Nọi, Ngã Hai, Làng Mảnh.

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

Hạn chế kỹ thuật

  • Thời gian nghiên cứu thực địa giới hạn trong 5 tháng (từ tháng 01/2015 đến tháng 05/2015), chưa theo dõi được biến động sinh khối và tỷ lệ sống của cây con qua trọn vẹn chu kỳ mùa khô khắc nghiệt nhiệt độ $<0^\circ\text{C}$.
  • Chưa thực hiện phân tích giải đoán ảnh viễn thám đa thời gian để xây dựng bản đồ số hóa không gian phân bố cho toàn bộ 833,30 ha trạng thái IIa.

Hướng nghiên cứu mở rộng

  • Ứng dụng GIS và ảnh vệ tinh đa phổ (Sentinel-2 / PlanetScope): Xây dựng thuật toán phân loại tự động trạng thái rừng và chỉ số thực vật ($NDVI, EVI$) nhằm giám sát tốc độ phục hồi rừng theo thời gian thực.
  • Thiết lập hệ thống ô định vị vĩnh cửu (Permanent Sampling Plots - PSP): Quan trắc liên tục sinh trưởng đường kính, chiều cao và động thái đào thải tự nhiên trong giai đoạn 5 - 10 năm.

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

                           CƠ CẤU ĐỐI TƯỢNG HƯỞNG LỢI

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

1. Yêu cầu kỹ thuật quan trọng nhất để triển khai khoanh nuôi xúc tiến tái sinh trạng thái IIa là gì?

Cần đánh giá mật độ cây triển vọng ($>1.000\text{ cây/ha}$) và kiểm soát độ tàn che tầng cao ở mức $0,6 - 0,7$. Nếu độ tàn che quá dày ($>0,8$), cây tái sinh thiếu sáng sẽ còi cọc; nếu mở tán quá mạnh ($<0,4$), tầng cây bụi thảm tươi sẽ bùng phát lấn át cây gỗ tái sinh.

2. Vì sao tỷ lệ cây tái sinh ở vị trí đỉnh ($2.907\text{ cây/ha}$) lại cao hơn ở vị trí chân ($2.613\text{ cây/ha}$)?

Vị trí đỉnh có độ mở tán lớn hơn, tạo điều kiện ánh sáng dồi dào cho các loài cây tiên phong ưa sáng nảy mầm từ hạt. Ngược lại, vị trí chân đồi tuy đất ẩm giàu dinh dưỡng nhưng độ tàn che lớn kết hợp với tác động giẫm đạp, gặm chồi của gia súc tập trung quanh khu dân cư làm giảm mật độ cây con sống sót.

3. Cần tích hợp những loài cây bản địa nào khi trồng bổ sung làm giàu rừng?

Ưu tiên các loài cây bản địa có giá trị kinh tế và phòng hộ cao, phù hợp với điều kiện khí hậu á nhiệt đới núi cao Sùng Đô: Quế (Cinnamomum cassia), Giổi xanh (Michelia mediocris), Trám trắng (Canarium album), Kháo (Machilus odoratissima), Re (Cinnamomum sp.).

4. Tần suất bảo dưỡng, phát luỗng thảm tươi định kỳ như thế nào?

Trong 3 năm đầu sau tác động lâm sinh, cần thực hiện phát luỗng dây leo, cây bụi lấn át định kỳ 2 lần/năm (lần 1 vào đầu mùa mưa tháng 4 - 5; lần 2 vào cuối mùa sinh trưởng tháng 9 - 10) xung quanh gốc cây tái sinh với bán kính $0,8 - 1,0\text{ m}$.

5. Chi phí đầu tư và thời gian thu hồi vốn (ROI) của mô hình này?

Tổng suất đầu tư khoanh nuôi có trồng làm giàu khoảng $6 - 8\text{ triệu VNĐ/ha}$ phân bổ trong 3 - 5 năm. Doanh thu từ tỉa thưa củi gỗ nhỏ và lâm sản ngoài gỗ bắt đầu từ năm thứ 6 - 8; khai thác chọn cây gỗ lớn từ năm thứ 15 - 20, mang lại tỷ suất hoàn vốn nội bộ (IRR) ước tính đạt $14 - 18%/\text{năm}$.


Kết luận

Nghiên cứu về đặc điểm tái sinh rừng tự nhiên trạng thái IIa tại xã Sùng Đô, huyện Văn Chấn, tỉnh Yên Bái đã xác lập hệ thống dữ liệu sinh thái định lượng chuẩn xác:

  • Tầng cây gỗ cao đạt mật độ $111 - 117\text{ cây/ha}$ với 11 loài ưu thế ($IV% \ge 5%$).
  • Tầng cây tái sinh phát triển mạnh mẽ với mật độ trung bình $2.658\text{ cây/ha}$, trong đó tỷ lệ cây triển vọng đạt $53,08%$ ($1.413\text{ cây/ha}$) và nguồn gốc tái sinh từ hạt chiếm ưu thế ($56,43%$).
  • Quần xã rừng trạng thái IIa có năng lực tự phục hồi rất cao, hoàn toàn đáp ứng tiêu chí kỹ thuật để áp dụng phương thức Khoanh nuôi xúc tiến tái sinh kết hợp làm giàu rừng bằng cây bản địa, loại bỏ nhu cầu trồng mới tốn kém.

Các cơ quan quản lý lâm nghiệp địa phương, ban quản lý rừng phòng hộ và chủ rừng có thể ứng dụng trực tiếp các định mức kỹ thuật và quy trình phân vùng sinh thái từ nghiên cứu này để xây dựng phương án quản lý rừng bền vững giai đoạn 2026 – 2035.