Giới thiệu dự án

Nghiên cứu sinh khối và chu trình carbon trong hệ sinh thái rừng đang là trọng tâm then chốt của các chiến lược thích ứng biến đổi khí hậu toàn cầu theo Công ước Khung của Liên Hợp Quốc về Biến đổi Khí hậu (UNFCCC) và cơ chế giảm phát thải từ mất rừng và suy thoái rừng (REDD+). Theo các báo cáo từ Ủy ban Liên chính phủ về Biến đổi Khí hậu (IPCC), sinh khối dưới mặt đất (Belowground Biomass - BGB), đặc biệt là rễ nhỏ (Fine Roots - đường kính $\le 2\text{ mm}$), đóng vai trò sống còn trong chu trình tuần hoàn dinh dưỡng và tích lũy carbon hữu cơ trong đất. Mặc dù rễ nhỏ chỉ chiếm dưới 1,5% tổng sinh khối thực vật trong lâm phần, nhưng năng suất sinh khối sơ cấp rễ nhỏ (Fine Root Production - FRP) có thể đóng góp từ 33% đến 75% tổng năng suất sơ cấp hữu hiệu (Net Primary Productivity - NPP) của cả hệ sinh thái rừng.

Tại Việt Nam, cây Mỡ (Manglietia conifera) là loài cây gỗ bản địa đa mục đích, có tốc độ sinh trưởng nhanh và giá trị kinh tế - phòng hộ cao, được trồng phổ biến tại các tỉnh miền núi phía Bắc như Bắc Kạn, Thái Nguyên, Tuyên Quang và Yên Bái. Tuy nhiên, các nghiên cứu định lượng lâm nghiệp truyền thống phần lớn tập trung vào sinh khối trên mặt đất (thân, cành, lá) và đường kính ngang ngực ($D_{1.3}$), trong khi động thái sinh khối rễ nhỏ và lượng carbon ngầm tích lũy trong đất chưa được lượng hóa chi tiết theo các cấp tuổi và tầng đất phát sinh.

Vấn đề nghiên cứu và các điểm nghẽn thực tiễn

Việc thiếu hụt cơ sở dữ liệu thực nghiệm về sinh khối rễ nhỏ dẫn đến các sai số hệ thống trong tính toán trữ lượng carbon rừng theo tiêu chuẩn Tier 3 của IPCC. Các khó khăn chính bao gồm:

  • Thiếu hụt dữ liệu phân tầng ngầm: Chưa xác định rõ quy luật phân bố không gian theo chiều sâu thẳng đứng (0–50 cm) của hệ rễ nhỏ theo từng giai đoạn khép tán và thành thục sinh học của cây Mỡ.
  • Đánh giá thấp tiềm năng bồi hoàn carbon: Bỏ sót nguồn carbon hữu cơ tuần hoàn nhanh do rễ nhỏ liên tục tái sinh và phân hủy, dẫn đến đánh giá chưa đầy đủ giá trị dịch vụ môi trường rừng.
  • Hạn chế trong mô hình hóa lâm phần: Thiếu các tham số đầu vào định lượng phục vụ xây dựng mô hình dự báo sinh trưởng và cân bằng dinh dưỡng đất vùng nhiệt đới gió mùa.

Mục tiêu nghiên cứu cụ thể

  1. Định lượng cấu trúc lâm phần: Đánh giá các chỉ số sinh trưởng lâm học gồm đường kính ngang ngực ($D_{1.3}$), chiều cao vút ngọn ($H_{vn}$), chiều cao dưới cành ($H_{dc}$) và mật độ cây ($N$) tại các cấp tuổi 3, 5, 7 và 9 năm tuổi ở xã Chu Hương, huyện Ba Bể, tỉnh Bắc Kạn.
  2. Xác định sinh khối tươi và sinh khối khô rễ nhỏ: Định lượng khối lượng rễ nhỏ ($g/\text{m}^2$ và $\text{tấn/ha}$) theo 5 tầng đất phát sinh ($0-10\text{ cm}$, $10-20\text{ cm}$, $20-30\text{ cm}$, $30-40\text{ cm}$, $40-50\text{ cm}$).
  3. Lượng hóa trữ lượng carbon ngầm: Xác định lượng carbon hữu cơ tích lũy trong hệ rễ nhỏ của các lâm phần rừng trồng Mỡ qua các giai đoạn tuổi, làm cơ sở khoa học cho công tác chi trả dịch vụ môi trường rừng và thẩm định tín chỉ carbon.

Phương pháp tiếp cận và phạm vi giới hạn

  • Phương pháp tiếp cận: Sử dụng phương pháp ống dung trọng khoan lõi đất phân tầng (Stratified Soil Coring Method) kết hợp phân tích trọng lượng nhiệt khô (Thermo-gravimetric method) trong phòng thí nghiệm.
  • Phạm vi nghiên cứu: 12 ô tiêu chuẩn (OTC) diện tích $500\text{ m}^2$ ($20\text{ m} \times 25\text{ m}$) đại diện cho 4 cấp tuổi (3, 5, 7, 9 tuổi) tại 3 vị trí địa hình (chân, sườn, đỉnh đồi) với độ dốc $6^\circ - 12^\circ$ tại xã Chu Hương, huyện Ba Bể, tỉnh Bắc Kạn; thu thập 108 điểm mẫu lõi đất (9 ô dạng bản $1\text{ m}^2$ mỗi OTC).
  • Giới hạn: Đề tài chỉ tập trung vào rễ nhỏ ($d \le 2\text{ mm}$) của cây Mỡ; không phân định động thái giữa rễ sống và necromass (rễ chết) và không đo lường hệ rễ thô ($d > 2\text{ mm}$).

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

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

Phương pháp đánh giá Ưu điểm kỹ thuật Hạn chế kỹ thuật Độ chính xác Mức độ phù hợp với rừng trồng Mỡ
Khai quật toàn phần (Whole-tree Excavation) Thu thập toàn diện rễ thô và rễ cọc. Phá hủy hoàn toàn kết cấu đất; làm đứt gãy và mất rễ nhỏ ($<2\text{ mm}$); chi phí nhân công rất cao. Thấp đối với rễ nhỏ Không khả thi cho điều tra mở rộng
Hộp quan sát rễ ngầm (Rhizotron / Minirhizotron) Quan sát trực tiếp động thái sinh trưởng theo thời gian thực (in-situ). Chi phí trang thiết bị cao; khó triển khai trên địa hình đồi núi dốc $6^\circ-12^\circ$; diện tích quan sát hạn hẹp. Cao về động thái sinh trưởng Hạn chế trong điều kiện thực địa Bắc Kạn
Ống sinh trưởng (Ingrowth Core) Đo lường chính xác tốc độ phát triển rễ mới hàng năm (FRP). Mất thời gian ủ mẫu kéo dài (6–12 tháng); làm xáo trộn vi môi trường đất ban đầu. Trung bình - Cao Phù hợp cho nghiên cứu chu kỳ dài
Khoan lõi đất phân tầng (Soil Coring Method - Lựa chọn) Định lượng chính xác mật độ sinh khối theo từng lớp đất thẳng đứng; giảm thiểu biến dạng cấu trúc mẫu; tính lặp lại cao. Đòi hỏi quy trình đãi mẫu và tách rễ tỉ mỉ; phụ thuộc vào độ chính xác của ống lấy mẫu dung trọng. Rất cao cho mật độ rễ tại thời điểm Tối ưu nhất cho điều kiện thực địa

So sánh với các nghiên cứu tiền nghiệm trên thế giới và Việt Nam

  • Thế giới: Ammer & Wagner (2005) chứng minh mô hình hóa sinh khối rễ nhỏ đạt 62–72% độ chính xác so với phương pháp ống dung trọng; Katrin Heinsoo & cs. (2009) chỉ ra rễ nhỏ chiếm 39–54% tổng sinh khối rễ ở tầng $0-10\text{ cm}$; Knute & Raich (1992) ghi nhận FRP đạt từ $52$ đến $1708\text{ g/m}^2/\text{năm}$ trên các hệ sinh thái rừng toàn cầu.
  • Việt Nam: Các nghiên cứu của Đỗ Hoàng Chung (2012), Võ Đại Hải (2008) và Vũ Tấn Phương (2006) tại rừng tự nhiên và rừng phục hồi chủ yếu đo lường thảm mục hoặc carbon trên mặt đất; số liệu tại Tây Nguyên (Kon Hà Nừng) cho thấy $>50%$ khối lượng rễ tập trung ở tầng $0-30\text{ cm}$. Nghiên cứu này bổ sung lỗ hổng thực nghiệm trên loài Mỡ tại Đông Bắc Bộ.
MoSCoW Requirements Framework:

Thiết kế hệ thống thu mẫu và phân tích

Ngăn xếp công cụ và thiết bị đo đạc (Technology & Tool Stack)

  • Thiết bị ngoại nghiệp: Thước dây chia vạch, thước sào đo cao chuyên dụng, địa bàn địa hình Silva v4.0, ống khoan đất dung trọng hình trụ bằng thép không gỉ ($d = 5{,}7\text{ cm}$, tiết diện bề mặt $s = 25{,}50\text{ cm}^2$, chiều dài $10\text{ cm}$).
  • Thiết bị nội nghiệp: Tủ sấy đối lưu Memmert UN55 (độ ổn định $\pm 0{,}5^\circ\text{C}$), Cân phân tích điện tử chuẩn xác Ohaus Pioneer PA214 ($4$ số lẻ, độ chia $0{,}0001\text{ g}$).
  • Phần mềm xử lý dữ liệu: Microsoft Excel 2013, Minitab v17.0, Python v3.10 (Pandas v2.0, NumPy v1.24, SciPy v1.10) cho tính toán thống kê và xử lý phân phối dữ liệu.

Mô hình thuật toán tính toán sinh khối và Carbon

Công thức xác định đường kính thân cây tại vị trí $1{,}3\text{ m}$: $$D_{1.3} = \frac{C}{\pi} \quad (\text{cm})$$ Trong đó: $C$ là chu vi thân cây ($\text{cm}$), $\pi \approx 3{,}14159$.

Công thức tính toán mật độ sinh khối rễ nhỏ ($SK_{\text{Fineroot}}$): $$SK_{\text{Fineroot}} = \frac{\sum_{i=1}^{n} P_i}{s \times n} \times a \quad (\text{g/m}^2)$$ Chuyển đổi sang đơn vị $\text{tấn/ha}$: $$SK_{\text{Fineroot (tấn/ha)}} = SK_{\text{Fineroot (g/m}^2)} \times 10^{-2}$$ Trong đó:

  • $P_i$: Khối lượng rễ nhỏ thu được trong lõi đất thứ $i$ ($\text{g}$).
  • $n$: Số lõi đất lấy mẫu trên mỗi OTC ($n = 9$).
  • $s$: Tiết diện mặt cắt ngang của ống khoan dung trọng ($s = 25{,}50\text{ cm}^2$).
  • $a$: Hệ số chuẩn hóa diện tích ($a = 10000\text{ cm}^2/\text{m}^2$).

Công thức tính trữ lượng Carbon tích lũy ngầm ($C_{\text{stock}}$): $$C_{\text{stock}} = SK_{\text{Fineroot (tấn/ha)}} \times CF$$ Trong đó: $CF = 0{,}50$ là hệ số chuyển đổi hàm lượng Carbon theo tiêu chuẩn IPCC.

import numpy as np
import pandas as pd

def calculate_fine_root_metrics(sample_weights_g: list[float], core_diameter_cm: float = 5.7, cf: float = 0.5) -> dict:
    """
    Tính toán sinh khối tươi/khô rễ nhỏ và trữ lượng Carbon tích lũy.
    """
    n = len(sample_weights_g)
    radius_cm = core_diameter_cm / 2.0
    surface_area_cm2 = np.pi * (radius_cm ** 2) # 25.5176 cm2 ~ 25.50 cm2
    
    total_weight_g = sum(sample_weights_g)
    # Tính toán g/m2
    biomass_g_per_m2 = (total_weight_g / (surface_area_cm2 * n)) * 10000.0
    # Chuyển đổi sang tấn/ha
    biomass_ton_per_ha = biomass_g_per_m2 * 0.01
    # Tính lượng Carbon (tấn C/ha)
    carbon_ton_per_ha = biomass_ton_per_ha * cf
    
    return {
        "biomass_g_m2": round(biomass_g_per_m2, 2),
        "biomass_ton_ha": round(biomass_ton_ha, 2),
        "carbon_ton_ha": round(carbon_ton_ha, 2)
    }

Phương pháp luận (Methodology) và Kế hoạch thực hiện

Timeline Đề tài (Tháng 08/2014 - Tháng 12/2014):
  • Quản trị rủi ro & Kiểm soát chất lượng (QA/QC):
    • Rủi ro lẫn tạp chất đất và rễ thô: Sử dụng rây mắt lưới $0{,}5\text{ mm}$ ngâm nước để loại bỏ hoàn toàn hạt sét, mùn hữu cơ bám dính; dùng thước kẹp panme kiểm tra loại bỏ rễ $>2\text{ mm}$.
    • Rủi ro sai số khối lượng do ẩm kế: Toàn bộ mẫu sấy đẳng nhiệt $105^\circ\text{C}$ đủ $8\text{ giờ}$ đến khi khối lượng không đổi, chuyển ngay vào bình hút ẩm silicagel trước khi cân trên cân phân tích Ohaus PA214.

Implementation và kết quả

Đặc điểm cấu trúc lâm phần rừng Mỡ qua các cấp tuổi

Dữ liệu điều tra từ 12 OTC ($500\text{ m}^2/\text{OTC}$) phân bố đồng đều tại các vị trí Chân, Sườn và Đỉnh đồi phản ánh rõ quá trình sinh trưởng và tỉa thưa tự nhiên:

OTC Tuổi (Năm) Độ dốc ($^\circ$) Vị trí địa hình $\bar{D}_{1.3}$ ($\text{cm}$) $\bar{H}_{vn}$ ($\text{m}$) $\bar{H}_{dc}$ ($\text{m}$) Mật độ ($N/\text{OTC}$) Mật độ quy đổi ($N/\text{ha}$)
01 3 7 Chân 6,66 5,01 2,39 39 780
02 3 8 Sườn 6,21 4,84 2,15 41 820
03 3 9 Đỉnh 6,29 4,88 2,10 41 820
04 5 7 Chân 11,23 8,10 4,74 34 680
05 5 8 Sườn 11,04 7,93 3,72 37 740
06 5 10 Đỉnh 10,77 7,83 3,76 32 640
07 7 6 Chân 16,23 14,34 8,27 11 220
08 7 8 Sườn 18,31 16,69 6,71 19 380
09 7 9 Đỉnh 17,17 15,58 6,78 16 320
10 9 9 Chân 19,47 17,43 10,56 9 180
11 9 10 Sườn 24,72 18,57 11,23 11 220
12 9 12 Đỉnh 18,86 15,84 9,54 14 280

Động thái sinh khối tươi rễ nhỏ ($g/\text{m}^2$)

Sinh khối tươi của hệ rễ nhỏ biến động mạnh theo cấp tuổi lâm phần và suy giảm nghiêm ngặt theo độ sâu tầng đất:

Nhóm tuổi Tầng 0–10 cm Tầng 10–20 cm Tầng 20–30 cm Tầng 30–40 cm Tầng 40–50 cm Tổng trung bình ($g/\text{m}^2$)
Tuổi 3 $1827{,}16$ ($57{,}4%$) $519{,}97$ ($16{,}3%$) $386{,}35$ ($12{,}1%$) $281{,}77$ ($8{,}9%$) $168{,}48$ ($5{,}3%$) $3183{,}73$
Tuổi 5 $2800{,}29$ ($57{,}8%$) $944{,}02$ ($19{,}5%$) $554{,}83$ ($11{,}5%$) $424{,}11$ ($8{,}7%$) $122{,}00$ ($2{,}5%$) $4845{,}32$
Tuổi 7 $2649{,}27$ ($72{,}4%$) $360{,}20$ ($9{,}8%$) $540{,}30$ ($14{,}8%$) $110{,}39$ ($3{,}0%$) $0{,}00$ ($0%$) $3660{,}13$
Tuổi 9 $1693{,}54$ ($70{,}6%$) $458{,}97$ ($19{,}2%$) $244{,}01$ ($10{,}2%$) $0{,}00$ ($0%$) $0{,}00$ ($0%$) $2396{,}50$

Sinh khối khô và Trữ lượng Carbon tích lũy ngầm

Khối lượng khô của rễ nhỏ và lượng Carbon tích lũy đạt đỉnh ở giai đoạn rừng Mỡ khép tán mạnh nhất (5 năm tuổi), sau đó giảm dần khi lâm phần bước vào giai đoạn thành thục công nghệ và tỉa thưa tự nhiên mạnh:

Cấp tuổi Chỉ tiêu định lượng Tầng 0–10 cm Tầng 10–20 cm Tầng 20–30 cm Tầng 30–40 cm Tầng 40–50 cm Tổng cộng
Tuổi 3 Sinh khối khô ($\text{tấn/ha}$) 4,71 1,51 1,23 0,90 0,52 8,87
Lượng Carbon ($\text{tấn C/ha}$) 2,35 ($53{,}1%$) 0,76 ($17{,}0%$) 0,61 ($13{,}9%$) 0,45 ($10{,}1%$) 0,26 ($5{,}9%$) 4,43
Tuổi 5 Sinh khối khô ($\text{tấn/ha}$) 7,16 2,56 1,77 1,25 0,52 13,41
Lượng Carbon ($\text{tấn C/ha}$) 3,58 ($54{,}0%$) 1,28 ($19{,}3%$) 0,89 ($13{,}4%$) 0,62 ($9{,}4%$) 0,26 ($3{,}9%$) 6,63
Tuổi 7 Sinh khối khô ($\text{tấn/ha}$) 7,20 0,93 1,42 0,29 0,00 9,85
Lượng Carbon ($\text{tấn C/ha}$) 3,60 ($73{,}2%$) 0,47 ($9{,}4%$) 0,71 ($14{,}4%$) 0,15 ($3{,}0%$) 0,00 ($0%$) 4,92
Tuổi 9 Sinh khối khô ($\text{tấn/ha}$) 4,24 1,19 0,64 0,00 0,00 6,07
Lượng Carbon ($\text{tấn C/ha}$) 2,12 ($69{,}9%$) 0,60 ($19{,}6%$) 0,32 ($10{,}5%$) 0,00 ($0%$) 0,00 ($0%$) 3,04
Trữ lượng Carbon tích lũy rễ nhỏ qua các cấp tuổi (tấn C/ha):
Tuổi 3: [████████████░░░░░░░░] 4.43 t C/ha
Tuổi 5: [███████████████████] 6.63 t C/ha  <-- Cực đại (Khép tán hoàn toàn)
Tuổi 7: [██████████████░░░░░] 4.92 t C/ha
Tuổi 9: [█████████░░░░░░░░░░] 3.04 t C/ha

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

  1. Xác lập bộ dữ liệu thực nghiệm chuẩn xác về rễ nhỏ loài Mỡ (Manglietia conifera): Cung cấp dữ liệu định lượng đầu tiên về sinh khối rễ nhỏ phân tầng cho loài cây bản địa này tại vùng núi đá vôi Đông Bắc Bộ, thay thế các hệ số giả định chung chung trước đây.
  2. Làm rõ quy luật tập trung không gian tầng mặt: Chứng minh rằng từ $53{,}1%$ đến $73{,}2%$ tổng lượng carbon ngầm của rễ nhỏ tập trung tại tầng đất mặt $0-10\text{ cm}$. Khi tuổi rừng tăng lên (từ 7–9 tuổi), hệ rễ nhỏ hầu như tiêu biến ở độ sâu $>30\text{ cm}$ do rễ cái hóa gỗ và nhu cầu hấp thu tập trung ở tầng thảm mục mặt đất.
  3. Phát hiện đỉnh sinh khối tại tuổi 5 (Giai đoạn khép tán cực đại): Lượng sinh khối khô rễ nhỏ đạt cực đại $13{,}41\text{ tấn/ha}$ (tương đương $6{,}63\text{ tấn C/ha}$) ở tuổi 5, cao hơn $51{,}2%$ so với tuổi 3 ($8{,}87\text{ tấn/ha}$) và cao hơn $120{,}9%$ so với tuổi 9 ($6{,}07\text{ tấn/ha}$). Đây là bằng chứng quan trọng chứng minh cường độ trao đổi chất ngầm mạnh nhất diễn ra khi tán rừng bắt đầu giao nhau hoàn toàn.
  4. Cải tiến độ chính xác thẩm định Carbon cho cơ chế REDD+: Giúp các đơn vị kiểm toán lâm nghiệp nâng cấp mức độ tính toán bể carbon dưới mặt đất từ Tier 1 (mặc định IPCC) lên Tier 3 (dữ liệu thực nghiệm đo đạc địa phương), tăng độ chính xác ước tính bể carbon ngầm thêm $25-40%$.

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

Kịch bản ứng dụng lâm nghiệp thực tiễn

  • Thương mại hóa tín chỉ Carbon (Forest Carbon Offsets): Tích hợp số liệu tích lũy Carbon rễ nhỏ ($3{,}04 - 6{,}63\text{ tấn C/ha}$) vào hồ sơ đăng ký dự án tín chỉ carbon theo tiêu chuẩn Quốc tế Verra (VCS - Verified Carbon Standard) hoặc Plan Vivo cho các diện tích rừng trồng Mỡ tại Đông Bắc Bộ.
  • Tối ưu hóa kỹ thuật lâm sinh và bón phân: Dữ liệu cho thấy rễ nhỏ hấp thụ dinh dưỡng tập trung $>70%$ ở tầng $0-10\text{ cm}$, do đó các biện pháp chăm sóc, bón thúc hoặc xới xáo nông chỉ nên tác động ở tầng mặt $0-10\text{ cm}$ để tránh gây tổn thương hệ rễ tơ của cây Mỡ.
  • Xác định thời điểm tỉa thưa hợp lý: Sự sụt giảm mật độ sinh khối rễ nhỏ sau 5 năm tuổi ($13{,}41\text{ tấn/ha}$ xuống $9{,}85\text{ tấn/ha}$ ở tuổi 7) chỉ ra sự cạnh tranh không gian ngầm khốc liệt; thời điểm $5-6\text{ tuổi}$ là cột mốc tối ưu để can thiệp tỉa thưa nuôi dưỡng rừng.
Lộ trình triển khai ứng dụng quản lý rừng bền vững:

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

Các hạn chế kỹ thuật thừa nhận

  • Chưa phân lập rễ sống và rễ chết (Live/Dead Root Partitioning): Do thu mẫu đồng loạt qua ống khoan dung trọng, nghiên cứu chưa tách biệt được tỷ lệ giữa rễ sống có chức năng sinh lý và rễ chết đang phân hủy (Necromass).
  • Tính chất thời điểm (Snapshot sampling): Thu thập số liệu tập trung vào mùa khô (tháng 08 đến tháng 12/2014), chưa phản ánh được biến động động thái rễ nhỏ qua các mùa trong năm (mùa sinh trưởng mưa nhiều vs mùa nghỉ khô lạnh).
  • Mối tương quan allometric: Chưa xây dựng phương trình hồi quy tương quan đa biến giữa sinh khối rễ nhỏ với diện tích tán lá ($S_{\text{tán}}$) và hàm lượng dinh dưỡng hữu cơ (N, P, K) trong đất.

Hướng nghiên cứu phát triển

  • Áp dụng kỹ thuật nhuộm rễ (Tetrazolium staining) hoặc phương pháp đồng vị phóng xạ carbon ($^{13}\text{C}/^{14}\text{C}$) để theo dõi tốc độ chu chuyển rễ (Fine root turnover rate).
  • Thiết lập hệ thống ống quan sát rễ ngầm kỹ thuật số (Minirhizotron system) để ghi nhận động thái sinh trưởng rễ 24/7 theo mùa vụ.
  • Mở rộng nghiên cứu trên các mô hình rừng trồng hỗn giao giữa Mỡ với các loài cây bản địa khác (Re gừng, Lát hoa, Trám).

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

Ma trận phân bổ giá trị cho các bên liên quan:

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

1. Yêu cầu kỹ thuật trang thiết bị tối thiểu để triển khai phương pháp khoan lõi đất là gì?

Cần chuẩn bị ống khoan sắt chuyên dụng hình trụ ($d = 5{,}7\text{ cm}$, $h = 10\text{ cm}$, tiết diện $25{,}50\text{ cm}^2$), búa cao su đóng ống lấy mẫu, rây đãi mẫu mắt lưới $0{,}5\text{ mm}$, khay nhôm, tủ sấy đẳng nhiệt điều chỉnh được nhiệt độ $105^\circ\text{C}$ ($\pm 1^\circ\text{C}$) và cân phân tích điện tử độ chính xác tối thiểu $0{,}001\text{ g}$.

2. Tại sao sinh khối rễ nhỏ lại đạt đỉnh ở tuổi 5 và giảm dần ở tuổi 7 và tuổi 9?

Ở tuổi 5, rừng Mỡ đạt độ khép tán cao nhất, tán lá tiếp xúc tối đa với ánh sáng và cây cần lượng lớn dinh dưỡng/nước để hoàn thiện sinh khối thân lá, thúc đẩy hệ rễ tơ phát triển mạnh mẽ. Sau tuổi 5, quá trình tỉa thưa tự nhiên làm giảm mật độ cây/ha (từ $680-740\text{ cây/ha}$ ở tuổi 5 xuống còn $180-280\text{ cây/ha}$ ở tuổi 9), đồng thời các rễ nhỏ chuyển hóa thành rễ thô hóa gỗ ($d > 2\text{ mm}$), khiến mật độ rễ nhỏ trên một đơn vị diện tích giảm xuống.

3. Tỷ lệ carbon $0{,}50$ áp dụng cho rễ nhỏ có độ tin cậy khoa học như thế nào?

Hệ số chuyển đổi $CF = 0{,}50$ (tương đương Carbon chiếm 50% khối lượng sinh khối khô) là tiêu chuẩn quốc tế được công nhận bởi IPCC Good Practice Guidance for LULUCF cho các bộ phận gỗ và rễ cây rừng thân gỗ nhiệt đới.

4. Phương pháp thu mẫu này có tích hợp được vào phần mềm GIS/MRV để giám sát rừng không?

Hoàn toàn có thể. Tọa độ 12 OTC và dữ liệu sinh khối rễ nhỏ theo tầng đất có thể được nhập vào cơ sở dữ liệu không gian (Geodatabase) trong ArcGIS hoặc QGIS, từ đó liên kết với ảnh viễn thám đa phổ (Sentinel-2, Landsat 8/9) để xây dựng bản đồ số trữ lượng Carbon toàn diện cho toàn khu vực huyện Ba Bể.

5. Chi phí triển khai đo đạc sinh khối rễ nhỏ trên quy mô diện rộng là bao nhiêu?

Chi phí lấy mẫu và phân tích phòng thí nghiệm dao động khoảng $1.500.000 - 2.500.000\text{ VNĐ}$ cho mỗi OTC (gồm 9 điểm mẫu phân tầng 5 lớp đất), tiết kiệm hơn $80%$ so với phương pháp khai quật toàn phần và rẻ hơn $95%$ so với lắp đặt hệ thống Minirhizotron.


Kết luận

Công trình nghiên cứu của tác giả Hoàng Văn Tuấn đã hoàn thành xuất sắc việc định lượng sinh khối rễ nhỏ ($d \le 2\text{ mm}$) và lượng Carbon ngầm tích lũy tại rừng trồng Mỡ (Manglietia conifera) qua 4 cấp tuổi (3, 5, 7, 9 tuổi) tại xã Chu Hương, huyện Ba Bể, tỉnh Bắc Kạn. Nghiên cứu mang lại những kết luận then chốt:

  • Tính quy luật phân tầng: Sinh khối rễ nhỏ tập trung cao nhất ở tầng đất mặt $0-10\text{ cm}$ (chiếm $53{,}1% - 73{,}2%$) và giảm mạnh theo độ sâu phẫu diện đất.
  • Động thái theo tuổi: Sinh khối khô và lượng Carbon tích lũy đạt cực đại ở giai đoạn 5 năm tuổi ($13{,}41\text{ tấn/ha}$ sinh khối khô; $6{,}63\text{ tấn C/ha}$), trước khi giảm dần ở tuổi 7 ($4{,}92\text{ tấn C/ha}$) và tuổi 9 ($3{,}04\text{ tấn C/ha}$).
  • Giá trị thực tiễn: Cung cấp cơ sở khoa học trực tiếp cho các biện pháp kỹ thuật lâm sinh tầng mặt, xác định thời điểm tỉa thưa và bổ sung dữ liệu định lượng giá trị dịch vụ môi trường rừng phục vụ thị trường tín chỉ carbon.

Quý độc giả, các nhà nghiên cứu lâm nghiệp và các nhà phát triển dự án carbon quan tâm đến việc ứng dụng bộ dữ liệu hoặc mở rộng mô hình điều tra sinh khối ngầm cho các loài cây bản địa khác có thể tham khảo phương pháp luận chuẩn hóa này để triển khai tại đơn vị.