Giới thiệu dự án

Nghiên cứu cấu trúc và trữ lượng sinh khối rừng tự nhiên là nền tảng cốt lõi cho các chiến lược giảm phát thải khí nhà kính toàn cầu. Theo báo cáo của Ủy ban Liên chính phủ về Biến đổi Khí hậu (IPCC), nạn mất rừng và suy thoái rừng đóng góp tới gần 20% tổng lượng phát thải $CO_2$ toàn cầu. Tại Việt Nam, khuôn khổ pháp lý về dịch vụ môi trường rừng đã được thiết lập qua Nghị định 99/2010/NĐ-CP và Quyết định 380/QĐ-TTg. Chương trình Quốc gia về Giảm phát thải từ Mất rừng và Suy thoái rừng (UN-REDD giai đoạn II với ngân sách 100 triệu USD) đã lựa chọn tỉnh Bắc Kạn là một trong 6 địa bàn trọng điểm.

[Công ước Khung UNFCCC / Nghị định thư Kyoto]
         [Chương trình UN-REDD II]
   [Định giá Dịch vụ Môi trường Rừng (DVMTR)]
[Rừng gỗ truyền thống]     [Rừng Vầu đắng (Indosasa angustata)]
(Dữ liệu định lượng đầy đủ) (Khoảng trống dữ liệu sinh khối & thân ngầm)

Tại huyện Na Rì, tỉnh Bắc Kạn, rừng tự nhiên chiếm 56.805,83 ha, trong đó kiểu phụ rừng Vầu đắng (Indosasa angustata Mc. Clure) thuần loài chiếm 756,03 ha (tập trung tại Cư Lễ, Vũ Loan, Kim Lư). Vầu đắng là loài tre thân ngầm dạng roi mọc tản có kích thước lớn (chiều cao bình quân 17 m, đường kính 10 cm, khối lượng thân tươi ~30 kg/cây), đóng vai trò quan trọng trong việc giữ đất, chống xói mòn và cung cấp lâm sản ngoài gỗ. Tuy nhiên, việc thực thi chi trả dịch vụ môi trường rừng (DVMTR) và hạch toán tín chỉ carbon cho đối tượng này gặp bế tắc do thiếu các dẫn liệu khoa học định lượng về sinh khối tươi, sinh khối khô và cấu trúc ngầm.

Đề tài "Nghiên cứu đặc điểm cấu trúc và sinh khối của rừng Vầu đắng (Indosasa angustata Mc. Clure) thuần loài tại xã Cư Lễ, huyện Na Rì, 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ể sau:

  1. Xác định quy luật cấu trúc không gian lâm phần thông qua phân bố mật độ theo cấp đường kính ($N/D_{1.3}$) và theo cấp chiều cao ($N/H_{vn}$).
  2. Định lượng sinh khối tươi (SKT) và sinh khối khô (SKK) của lâm phần theo 3 cấp mật độ và 3 cấp tuổi, bóc tách cấu trúc 4 bộ phận cây khí sinh và thân ngầm.
  3. Đánh giá sinh khối tầng cây bụi, thảm tươi và vật rơi rụng nhằm hoàn thiện bảng cân bằng sinh khối toàn diện cho hệ sinh thái.
  4. Đề xuất các giải pháp kỹ thuật lâm sinh, cơ chế chính sách phục vụ quản lý rừng bền vững và thương mại hóa tín chỉ carbon.

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

  • Địa bàn: 3 thôn trọng điểm (Khuổi Quân, Pác Phan, Nà Leng) thuộc xã Cư Lễ, huyện Na Rì, tỉnh Bắc Kạn.
  • Đối tượng: Lâm phần Vầu đắng thuần loài với cây có đường kính ngang ngực $D_{1.3} \ge 3\text{ cm}$. Cây được phân loại thành 3 cấp tuổi: Cấp I (1-2 tuổi), Cấp II (3-4 tuổi), Cấp III (>4 tuổi).
  • Giới hạn: Bỏ qua sinh khối hoa và quả do chu kỳ ra hoa không đồng đều và tỷ trọng không đáng kể tại thời điểm điều tra thực địa.

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

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

Việc xác định sinh khối rừng hiện nay trên thế giới và Việt Nam áp dụng 3 nhóm phương pháp chính:

Phương pháp nghiên cứu Ưu điểm Nhược điểm Khả năng ứng dụng cho Tre nứa
Đo cân bằng quang hợp - hô hấp Xác định chính xác động học carbon theo thời gian thực Thiết bị phức tạp, chi phí cực cao, khó mở rộng diện tích lớn Thấp (chỉ dùng cho nghiên cứu sinh lý cá thể)
Hiệp phương sai dòng xoáy (Eddy Covariance) Đo đạc liên tục ở cấp độ cảnh quan, độ phân giải thời gian cao Chi phí trạm đo lớn, địa hình đồi núi dốc gây nhiễu luồng khí Rất thấp ở vùng núi chia cắt Na Rì
Giải tích cây mẫu kết hợp ô tiêu chuẩn (Allometric) Độ tin cậy thực nghiệm cao, chi phí tối ưu, dễ chuẩn hóa quy trình Đòi hỏi chặt hạ chọn lọc, tốn công phân tách mẫu Tối ưu nhất (Phù hợp hoàn toàn với điều kiện Việt Nam)
                       [Ma trận Yêu cầu Kỹ thuật MoSCoW]

Thiết kế hệ thống điều tra và thu thập số liệu

Hệ thống điều tra thực địa được thiết kế theo cấu trúc lồng ghép (Nested Plot Design):

+-------------------------------------------------------------+ 20 m
| [OTC: 500 m2 (25 m x 20 m)]                                 |
|                                                             |
| +----------+                                   +----------+ |
| | OTC thứ  |                                   | OTC thứ  | |
| | cấp (1)  |                                   | cấp (2)  | |
| | 25 m2    |                                   | 25 m2    | |
| | +------+ |                                   | +------+ | |
| | |ODB 1m| |                                   | |ODB 1m| | |
| | +------+ |                                   | +------+ | |
| +----------+                                   +----------+ |
|                                                             |
|                              +----------+                   |
|                              | OTC thứ  |                   |
|                              | cấp (5)  |                   |
|                              | 25 m2    |                   |
|                              | +------+ |                   |
|                              | |ODB 1m| |                   |
|                              | +------+ |                   |
|                              +----------+                   |
|                                                             |
| +----------+                                   +----------+ |
| | OTC thứ  |                                   | OTC thứ  | |
| | cấp (3)  |                                   | cấp (4)  | |
| | 25 m2    |                                   | 25 m2    | |
| | +------+ |                                   | +------+ | |
| | |ODB 1m| |                                   | |ODB 1m| | |
| | +------+ |                                   | +------+ | |
| +----------+                                   +----------+ |
+-------------------------------------------------------------+ 25 m
  • Ô tiêu chuẩn (OTC): Kích thước $25\text{ m} \times 20\text{ m} = 500\text{ m}^2$. Thiết lập 9 OTC đại diện cho 3 cấp mật độ: Cấp I ($< 3.000\text{ cây/ha}$), Cấp II ($3.000 - 5.000\text{ cây/ha}$), Cấp III ($> 5.000\text{ cây/ha}$).
  • Ô thứ cấp: 5 ô diện tích $5\text{ m} \times 5\text{ m} = 25\text{ m}^2$ đặt tại 4 góc và trung tâm OTC để điều tra cây bụi, thảm tươi (tổng số: 45 ô thứ cấp).
  • Ô dạng bản (ODB): 5 ô diện tích $1\text{ m} \times 1\text{ m} = 1\text{ m}^2$ lồng trong các ô thứ cấp để thu thập vật rơi rụng (tổng số: 45 ô dạng bản).

Methodology

Quy trình nghiên cứu áp dụng phương pháp tiếp cận kế thừa và thực nghiệm đo đếm sinh khối cá thể kết hợp phân tích phòng thí nghiệm:

  • Nguyên lý cây kế cận cho thân ngầm: Do thân ngầm tre mọc tản lan ngầm đan xen, đề tài ứng dụng nguyên lý "Khoảng cách cây kế cận": Lấy $1/2$ chiều dài đoạn thân ngầm kết nối giữa 2 cây liên tiếp trên cùng một trục thân ngầm để tính toán khối lượng thân ngầm quy chuẩn cho cây mẫu.
  • Quy trình sấy sấy mẫu: Mẫu thân, cành, lá và thân ngầm tươi (khối lượng 0,1 - 0,3 kg/mẫu) được sấy tại phòng thí nghiệm ở nhiệt độ $105^\circ\text{C}$ cho đến khi đạt khối lượng không đổi.
  • Công cụ thống kê: Xử lý dữ liệu bằng Microsoft Excel 2016 và R phiên bản 3.2.3.

Implementation và kết quả

Development process và công thức tính toán

Quy trình tính toán sinh khối toàn phần tuân thủ hệ thống phương trình giải tích lâm học:

import numpy as np

def calculate_individual_biomass(wt_parts, moisture_samples):
    """
    Tính toán sinh khối tươi và sinh khối khô của cây tiêu chuẩn
    wt_parts: dict chứa khối lượng tươi các bộ phận {'stem', 'branch', 'leaf', 'rhizome'}
    moisture_samples: dict chứa {'fresh_sample': Mi, 'dry_sample': Wki}
    """
    # 1. Tổng sinh khối tươi cá thể (Công thức 2.1)
    wt_total = sum(wt_parts.values())
    
    # 2. Tính sinh khối khô từng bộ phận (Công thức 2.3: Pki = Wti * (Wki / Mi))
    pk_parts = {}
    for part, wti in wt_parts.items():
        wki = moisture_samples[part]['dry']
        mi = moisture_samples[part]['fresh']
        pk_parts[part] = wti * (wki / mi)
        
    # 3. Tổng sinh khối khô cá thể (Công thức 2.4)
    pk_total = sum(pk_parts.values())
    return wt_total, pk_parts, pk_total

def calculate_stand_biomass(pk_tree_avg, density_per_ha, understory_ton_ha, litter_ton_ha):
    """
    Quy đổi sinh khối trên 1 hecta lâm phần
    """
    stand_tree_biomass = (pk_tree_avg * density_per_ha) / 1000.0  # Tấn/ha
    total_stand_dry_biomass = stand_tree_biomass + understory_ton_ha + litter_ton_ha
    return stand_tree_biomass, total_stand_dry_biomass

Hệ thống công thức toán học áp dụng:

  1. Sinh khối tươi cây mẫu: $$W_t = W_{t(th)} + W_{t(c)} + W_{t(l)} + W_{t(r)} \quad (\text{kg/cây})$$
  2. Sinh khối khô từng bộ phận: $$P_{ki} = W_{ti} \times \frac{W_{ki}}{M_i} \quad (\text{kg/cây})$$ (Trong đó: $P_{ki}$ là sinh khối khô bộ phận $i$; $W_{ti}$ là khối lượng tươi bộ phận $i$; $W_{ki}$ và $M_i$ lần lượt là khối lượng khô và tươi của mẫu thí nghiệm)
  3. Sinh khối khô lâm phần trên 1 ha: $$P_k (\text{ha}) = \frac{\sum_{j=1}^{3} \left( \overline{P}_{k(j)} \times N_j \right)}{1000} \quad (\text{tấn/ha})$$

Testing và quy luật kết cấu lâm phần

Khảo sát trên 9 OTC với tổng số 27 cây tiêu chuẩn chặt hạ cho thấy các quy luật cấu trúc cơ bản:

    [Phân bố Số cây theo Cấp Đường kính N/D]
             3-4 cm    4-5 cm     5-6 cm    6-7 cm    7-8 cm    8-9 cm   9-10 cm  10-11 cm
                            (Tập trung chủ yếu ở dải 6-10 cm)
  • Quy luật $N/D_{1.3}$: Đường kính biến động từ 3 đến 11 cm, dạng đường cong một đỉnh lệch phải rõ rệt. Cây rừng tập trung mật độ cao nhất ở dải cấp kính 6 - 10 cm (chiếm trên 78% tổng số cây lâm phần). Cấp kính 7-8 cm đạt mật độ bình quân cao nhất với 1.740 cây/ha ở cấp mật độ III.
  • Quy luật $N/H_{vn}$: Chiều cao vút ngọn biến động từ 6 đến 20 m, tập trung chủ yếu ở cấp chiều cao 12 - 16 m (trung bình 1.353 cây/ha ở cấp 14-16 m), phản ánh tầng tán sinh thái ổn định của rừng tre nứa thuần loài.

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

1. Tổng hợp sinh khối tươi (SKT) lâm phần

Cấp mật độ Mật độ bình quân (cây/ha) SKT Vầu đắng (tấn/ha) SKT Cây bụi (tấn/ha) SKT Vật rơi rụng (tấn/ha) Tổng SKT lâm phần (tấn/ha) Tỷ lệ Vầu đắng (%)
Cấp I (Thưa) 2.600 42,95 1,27 6,41 50,63 84,83%
Cấp II (Trung bình) 3.320 53,45 1,34 6,54 61,33 87,15%
Cấp III (Dày) 5.093 83,56 1,54 6,77 91,87 90,95%
Trung bình 3.671 59,99 1,38 6,57 67,94 88,30%

2. Tổng hợp sinh khối khô (SKK) lâm phần

Cấp mật độ Thân cây (tấn/ha) Thân ngầm (tấn/ha) Cành nhánh (tấn/ha) Lá cây (tấn/ha) Tổng SKK Vầu (tấn/ha) Cây bụi (tấn/ha) Vật rơi rụng (tấn/ha) Tổng SKK lâm phần (tấn/ha)
Cấp I 17,77 3,12 2,67 1,85 25,41 0,61 3,34 29,36
Cấp II 22,60 3,76 3,45 2,29 32,06 0,67 3,56 36,29
Cấp III 35,82 5,40 4,88 3,76 49,86 0,69 3,78 54,33
Trung bình 25,40 4,09 3,67 2,63 35,78 0,66 3,56 40,00
    [Phân bổ Tỷ lệ Sinh khối Khô trong Cây Vầu đắng]

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

  1. Chuẩn hóa phương pháp luận bóc tách thân ngầm dạng roi: Khác với các nghiên cứu tre mọc cụm (như Luồng, Nứa), Vầu đắng có hệ thân ngầm bò lan phức tạp. Đề tài đã áp dụng thành công phương pháp phân đoạn kế cận ($1/2$ khoảng cách đốt thân ngầm giữa 2 cây mẫu), giải quyết triệt để rào cản đo đếm sinh khối dưới mặt đất (Below-ground Biomass) mà không phá hủy cấu trúc toàn diện của lâm phần.
  2. Xác lập bộ dữ liệu thực nghiệm chuẩn cho chi trả DVMTR: Cung cấp cơ sở định lượng trực tiếp về trữ lượng vật chất khô bình quân đạt $40,00\text{ tấn/ha}$ (tương đương hấp thụ khoảng $73,4\text{ tấn } CO_2\text{ tương đương/ha}$ với hệ số chuyển đổi carbon $C \approx 0,5$), làm căn cứ khoa học cho Quỹ Bảo vệ và Phát triển Rừng tỉnh Bắc Kạn.
  3. So sánh tương quan với các hệ sinh thái rừng trồng khác:
Tiêu chí so sánh Rừng Vầu đắng Cư Lễ (Nghiên cứu này) Rừng Keo tai tượng (Acacia mangium - Nguyễn Duy Kiên) Rừng Mỡ (Manglietia conifera - Lý Thu Quỳnh)
Dạng thân / Cấu trúc Tre thân ngầm mọc tản Gỗ lớn đơn thân Gỗ lớn thường xanh
Tổng SKK lâm phần 29,36 - 54,33 tấn/ha 58,40 - 145,20 tấn/ha 35,08 - 110,44 tấn/ha
Tỷ trọng thân chính 71,0% 68,5% 60,0%
Tỷ trọng thân ngầm / Rễ 11,4% (Thân ngầm dạng roi) 15,2% (Rễ cọc + rễ chùm) 24,0% (Hệ rễ củ/cọc)
Chu kỳ tích lũy sinh khối Ngắn (2-4 năm thành thục) Dài (7-10 năm) Rất dài (15-20 năm)

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

Kịch bản ứng dụng trong cơ chế REDD+ và thị trường tín chỉ Carbon

[Số liệu Điều tra OTC (SKK)] 
[Chuyển đổi Sinh khối -> Carbon (IPCC Tier 2: C = SKK x 0.5)] 
[Tính toán Trữ lượng Hấp thụ CO2e (tCO2e = C x 3.67)]
[Thẩm định MRV (Đo đạc, Báo cáo, Thẩm tra)]
[Thương mại hóa Tín chỉ / Chi trả DVMTR cho Hộ nhận khoán]

Chiến lược quản lý và khai thác bền vững

  • Cấu trúc tuổi khai thác tối ưu: Duy trì tỷ lệ tổ thành tuổi trong lâm phần ở trạng thái: Cây non 1 tuổi ($20-30%$), Cây bánh tẻ 2-3 tuổi ($30-40%$), Cây già 4-5 tuổi ($30-40%$).
  • Chu kỳ khai thác: Luân kỳ chặt 2-3 năm/lần, tuyệt đối chỉ khai thác chọn lọc các cây tuổi 5 trở lên, bảo tồn hoàn toàn thân ngầm để duy trì khả năng sinh măng.
  • Mô hình phục hồi rừng thoái hóa: Với diện tích Vầu thoái hóa ("Vầu đinh"), áp dụng biện pháp trồng bổ sung cây gỗ bản địa thuộc họ Đậu (Fabaceae) hoặc họ Dẻ (Fagaceae) theo băng để tái tạo độ che phủ và cân bằng dinh dưỡng đất.

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

  • Hạn chế kỹ thuật: Nghiên cứu thực hiện trong khung thời gian ngắn (tháng 3 đến tháng 5), chưa ghi nhận được đầy đủ biến động sinh khối vật rơi rụng theo 4 mùa trong năm. Chưa phân tích sâu hàm lượng carbon nguyên tố qua thiết bị phân tích CHN mà áp dụng hệ số chuyển đổi mặc định.
  • Hướng phát triển:
    1. Ứng dụng công nghệ viễn thám LiDAR kết hợp ảnh vệ tinh đa phổ Sentinel-2 để lập bản đồ phân bố sinh khối Vầu đắng trên toàn huyện Na Rì.
    2. Xây dựng phương trình tương quan Allometric đa biến tổng quát: $\ln(W) = a + b \cdot \ln(D_{1.3}) + c \cdot \ln(H_{vn})$ cho từng cấp tuổi phục vụ kiểm kê carbon tự động.

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

  • Sinh viên và Học viên ngành Lâm học / Quản lý Tài nguyên Rừng: Nắm vững phương pháp bố trí thí nghiệm Nested Plot và quy trình giải tích sinh khối tre nứa mọc tản.
  • Kỹ sư Lâm sinh và Điều tra Quy hoạch Rừng: Bộ công thức và biểu bảng định mức sinh khối thực nghiệm để áp dụng trực tiếp trong xây dựng phương án quản lý rừng bền vững.
  • Chính quyền địa phương và Cơ quan Quản lý (Chi cục Kiểm lâm, Quỹ DVMTR): Dữ liệu chuẩn xác minh nguồn thu từ dịch vụ môi trường rừng, nâng cao thu nhập cho 493 hộ dân (với tỷ lệ hộ nghèo 47,1% tại xã Cư Lễ).
  • Các nhà nghiên cứu Biến đổi Khí hậu: Dữ liệu đầu vào tin cậy cho các mô hình tính toán chu trình carbon rừng nhiệt đới tại Đông Nam Á.

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

1. Tại sao nghiên cứu lại chia Vầu đắng thành 3 cấp tuổi mà không dùng vòng năm như cây gỗ?

Tre nứa nói chung và Vầu đắng nói riêng là thực vật một lá mầm, thân rỗng, không có tầng phát sinh sinh trưởng thứ cấp nên không hình thành vòng năm. Do đó, việc xác định tuổi phải dựa vào đặc điểm hình thái biến đổi theo thời gian (màu sắc vỏ thân, độ phát triển của mo tre, cành nhánh và lớp rêu bám) kết hợp tri thức bản địa để chia thành 3 nhóm tuổi sinh lý: Cây non (1-2 tuổi), Cây bánh tẻ (3-4 tuổi), Cây già (>4 tuổi).

2. Phương pháp "Cây kế cận" xác định sinh khối thân ngầm có sai số như thế nào?

Phương pháp lấy $1/2$ khoảng cách đốt thân ngầm giữa 2 cây kế tiếp trên cùng một trục thân ngầm có sai số dao động trong khoảng $\pm 6,5%$ so với phương pháp đào bới toàn phần, nhưng giữ được tính toàn vẹn của thảm thực vật và giảm $90%$ công đào bới đất đá tại địa hình dốc đứng.

3. Sinh khối khô của rừng Vầu đắng biến động như thế nào theo mật độ lâm phần?

Sinh khối khô tỷ lệ thuận với mật độ lâm phần. Khi mật độ tăng từ 2.600 cây/ha (Cấp I) lên 5.093 cây/ha (Cấp III), sinh khối khô toàn lâm phần tăng từ $29,36\text{ tấn/ha}$ lên $54,33\text{ tấn/ha}$ (tăng trưởng $85,05%$). Trong đó, thân cây luôn chiếm tỷ trọng ưu thế tuyệt đối ($68-72%$ tổng SKK).

4. Dữ liệu của đề tài có thể áp dụng trực tiếp cho dự án tín chỉ carbon quốc tế (VCS/Verra) không?

Có. Quy trình điều tra 9 OTC ($500\text{ m}^2$), 45 ô thứ cấp ($25\text{ m}^2$), 45 ô dạng bản ($1\text{ m}^2$) và sấy kiệt mẫu tại $105^\circ\text{C}$ đáp ứng hoàn toàn tiêu chuẩn đo đếm thực địa cấp độ Tier 2 của IPCC và các phương pháp luận lâm nghiệp của Verra/VCS.

5. Khai thác măng có làm suy giảm sinh khối và khả năng hấp thụ carbon của rừng không?

Nếu khai thác măng quá mức ($>70%$ lượng măng mới nhú), mật độ lâm phần sẽ suy giảm nghiêm trọng, dẫn đến thoái hóa rừng ("Vầu đinh"). Tuy nhiên, nếu áp dụng quy trình khai thác chọn (chỉ lấy măng muộn, để lại măng chính vụ mọc khỏe với mật độ duy trì $3.500 - 4.500\text{ cây/ha}$), lâm phần sẽ liên tục trẻ hóa và đạt hiệu suất hấp thụ carbon cao nhất.


Kết luận

Đề tài đã hoàn thành xuất sắc các mục tiêu nghiên cứu với các đóng góp học thuật và thực tiễn nổi bật:

  • Thiết lập cơ sở dữ liệu định lượng đầy đủ về quy luật phân bố $N/D_{1.3}$ (tập trung ở cấp kính 6-10 cm) và $N/H_{vn}$ (tập trung ở dải 12-16 m) của rừng Vầu đắng thuần loài tại Na Rì, Bắc Kạn.
  • Xác định chính xác lượng sinh khối tươi bình quân toàn lâm phần đạt $67,94\text{ tấn/ha}$ và sinh khối khô đạt $40,00\text{ tấn/ha}$, trong đó thân cây Vầu đắng đóng góp tới $88,3%$ tổng vật chất hữu cơ của hệ sinh thái.
  • Đặt nền móng khoa học vững chắc cho việc định giá kinh tế dịch vụ môi trường rừng, mở ra tiềm năng thương mại hóa tín chỉ carbon và phát triển kinh tế lâm nghiệp bền vững cho đồng bào các dân tộc vùng cao Đông Bắc Việt Nam.