Giới thiệu dự án

Biến đổi khí hậu toàn cầu đang là thách thức nghiêm trọng nhất đối với nhân loại trong thế kỷ 21. Theo báo cáo đánh giá của Ủy ban Liên Chính phủ về Biến đổi Khí hậu (IPCC), có tới 95% khả năng các hoạt động phát thải khí nhà kính (KNK) của con người là nguyên nhân chính gây ra hiện tượng ấm lên toàn cầu. Để đối phó với tình trạng này, các công cụ tài chính và cơ chế quốc tế như Công ước Khung của Liên Hợp Quốc về Biến đổi Khí hậu (UNFCCC), Cơ chế Phát triển Sạch (CDM - Clean Development Mechanism), chương trình REDD+ và Quỹ Khí hậu Xanh (GCF với cam kết 100 tỷ USD/năm) đã được thiết lập nhằm lượng hóa và thương mại hóa khả năng hấp thụ $\text{CO}_2$ của hệ sinh thái rừng.

Tại Việt Nam, các mô hình canh tác nông lâm kết hợp (NLKH) giữ vai trò then chốt trong việc bảo vệ đất dốc, chống xói mòn và tạo sinh kế bền vững. Tuy nhiên, phần lớn các công trình nghiên cứu trữ lượng carbon trước đây chỉ tập trung vào rừng tự nhiên hoặc rừng trồng thuần loài (Keo, Thông, Bạch đàn), bỏ ngỏ mảng định lượng hấp thụ $\text{CO}_2$ của cây gỗ trong hệ thống NLKH đa tầng tán.

                  HỆ SINH THÁI NÔNG LÂM KẾT HỢP (CHÈ - RỪNG)

Đề tài "Nghiên cứu khả năng hấp thụ $\text{CO}_2$ của một số loài cây gỗ trồng xen trong mô hình nông lâm kết hợp (chè – rừng) tại xã Yên Ninh, huyện Phú Lương, tỉnh Thái Nguyên" được thực hiện nhằm giải quyết bài toán thiếu hụt dữ liệu khoa học về trữ lượng carbon ở quy mô địa phương.

Mục tiêu dự án

  1. Xác định các chỉ tiêu sinh trưởng và cấu trúc lâm phần: Đo đếm đường kính ngang ngực ($D_{1.3}$), chiều cao vút ngọn ($H_{vn}$), mật độ ($N$) và trữ lượng gỗ ($M$) của 4 loài cây gỗ trồng xen: Keo tai tượng (Acacia mangium), Muồng (Cassia siamea), Lát hoa (Chukrasia tabularis) và Xoan ta (Melia azedarach).
  2. Thiết lập quy luật tương quan sinh trưởng: Xây dựng phương trình tương quan một nhân tố dạng $H_{vn} = f(D_{1.3})$ để phục vụ dự đoán sinh trưởng phi phá hủy.
  3. Định lượng cấu trúc sinh khối và hàm lượng Carbon: Phân tích sinh khối tươi (SKT), sinh khối khô kiệt (SKK) theo từng bộ phận (thân, cành, lá, rễ) và tỷ lệ carbon nguyên tố ($C_{%}$).
  4. Lượng hóa khả năng hấp thụ $\text{CO}_2$ và giá trị kinh tế môi trường: Tính toán tổng lượng $\text{CO}_2$ cố định quy đổi ($\text{tấn/ha}$) và ước tính giá trị thương mại thông qua chứng chỉ giảm phát thải CERs với mức giá tham chiếu quốc tế 20 USD/tấn $\text{CO}_2$.

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

  • Địa bàn: Xã Yên Ninh, huyện Phú Lương, tỉnh Thái Nguyên (tổng diện tích mô hình chè - rừng khoảng 330 ha).
  • Đối tượng: 4 loài cây gỗ trồng xen với mật độ dao động từ 73 đến 233 cây/ha.
  • Giới hạn kỹ thuật: Sinh khối rễ cây gỗ được xác định gián tiếp qua hệ số tương quan sinh thái ($B_{\text{dưới mặt đất}} = \frac{1}{4} B_{\text{trên mặt đất}}$) nhằm bảo vệ nguyên vẹn hệ thống rễ của cây chè kinh doanh bên dưới.

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

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

Trước khi nghiên cứu được triển khai, việc đánh giá giá trị môi trường của các mô hình canh tác tại địa phương gặp nhiều bất cập do áp dụng sai lệch hệ số từ rừng thuần loài.

Tiêu chí so sánh Rừng trồng thuần loài Độc canh nông nghiệp (Chè thuần) Mô hình NLKH Chè - Rừng (Giải pháp đề xuất)
Hiệu quả kinh tế ngắn hạn Thấp (chờ chu kỳ khai thác 5–10 năm) Cao (thu hoạch búp chè định kỳ) Rất cao (kết hợp thu nhập ngắn hạn từ chè và dài hạn từ gỗ)
Khả năng cố định $\text{CO}_2$ Cao ($100 - 180\text{ tấn/ha}$) Rất thấp ($< 15\text{ tấn/ha}$) Khá cao ($40 - 70\text{ tấn/ha}$ tích lũy bổ sung)
Xói mòn và rửa trôi đất Trung bình (kém khi mới trồng/khai thác) Cao (đặc biệt vùng đồi dốc trung du) Tối thiểu (cấu trúc đa tầng giữ đất và điều hòa ẩm)
Cơ sở dữ liệu chi trả PFES Đã có quy chuẩn cơ bản Chưa đủ điều kiện áp dụng Cung cấp luận cứ kỹ thuật hoàn thiện biểu phí

Phân tích yêu cầu kỹ thuật theo mô hình MoSCoW:

  • Must have (Bắt buộc): Số liệu đo đếm thực địa $D_{1.3}$, $H_{vn}$, sinh khối các bộ phận thân - cành - lá - rễ, tỷ lệ $C_{%}$ phân tích phòng thí nghiệm, phương trình $H_{vn} - D_{1.3}$.
  • Should have (Nên có): Bảng tra phân bố số cây theo cấp kính ($N/D$), hệ số tương quan hồi quy $R \ge 0.85$.
  • Could have (Có thể mở rộng): Bản đồ số hóa phân vùng trữ lượng $\text{CO}_2$ trên toàn xã Yên Ninh.
  • Won't have (Tạm thời loại trừ): Đào bới cơ học toàn bộ bộ rễ cây gỗ (tránh phá hủy đất chè).

Thiết kế hệ thống thu thập và xử lý số liệu

Quy trình nghiên cứu kết hợp giữa điều tra lâm học ngoại nghiệp và phân tích hóa lý nội nghiệp:

Danh mục thiết bị và công cụ sử dụng:

  • Thước kẹp kính (Caliper) & Thước dây sợi thủy tinh (50m): Đo $D_{1.3}$ chính xác đến $\pm 0.1\text{ cm}$.
  • Thước đo cao Blume-Leiss: Đo chiều cao vút ngọn $H_{vn}$.
  • Cân điện tử hiện trường & Cân phân tích phòng thí nghiệm: Độ chính xác $0.01\text{ g}$.
  • Tủ sấy đối lưu cưỡng bức: Dải nhiệt độ $50 - 250^\circ\text{C}$.
  • Phần mềm phân tích dữ liệu: Microsoft Excel (Analysis ToolPak, hàm SUM, AVERAGE, REGRESSION).

Implementation và kết quả

Quy trình tính toán và mô hình thuật toán

Dưới đây là mã giả thuật toán chuẩn hóa quy trình tính toán sinh khối, carbon và giá trị kinh tế môi trường:

def calculate_carbon_and_economic_value(trees_data, c_concentrations, cer_price_usd=20, exchange_rate=21000):
    """
    Tính toán sinh khối, Carbon, CO2 và giá trị kinh tế theo chuẩn ICRAF (2010)
    """
    results = []
    for tree in trees_data:
        # 1. Sinh khối trên mặt đất (AGB) và dưới mặt đất (BGB)
        dry_aboveground = tree['dry_stem'] + tree['dry_branch'] + tree['dry_leaf']
        dry_belowground = 0.25 * dry_aboveground  # Hệ số rễ chuẩn
        total_dry_biomass_kg = dry_aboveground + dry_belowground
        
        # 2. Quy đổi sinh khối khô trên 1 ha (tấn/ha)
        total_dry_biomass_ton_ha = (total_dry_biomass_kg * tree['density_n_ha']) / 1000.0
        
        # 3. Lượng tích lũy Carbon (tấn C/ha)
        c_accum_ton_ha = (
            (tree['dry_stem'] * c_concentrations[tree['species']]['stem'] +
             tree['dry_branch'] * c_concentrations[tree['species']]['branch'] +
             tree['dry_leaf'] * c_concentrations[tree['species']]['leaf'] +
             dry_belowground * c_concentrations[tree['species']]['root']) * tree['density_n_ha']
        ) / 1000.0
        
        # 4. Lượng CO2 hấp thụ tương đương (tấn CO2/ha) theo tỷ lệ mol 44/12
        co2_ton_ha = c_accum_ton_ha * (44.0 / 12.0)
        
        # 5. Lượng giá kinh tế (VND/ha)
        economic_val_vnd = co2_ton_ha * cer_price_usd * exchange_rate
        
        results.append({
            'species': tree['species'],
            'co2_ton_ha': round(co2_ton_ha, 2),
            'economic_val_vnd': round(economic_val_vnd, 0)
        })
    return results

Kết quả kiểm định tương quan sinh trưởng

Phân tích hồi quy đơn biến giữa $H_{vn}$ và $D_{1.3}$ cho thấy mối quan hệ tương quan thuận rất chặt chẽ ở tất cả các loài nghiên cứu:

Loài cây Phương trình tương quan thực nghiệm Hệ số tương quan ($R$) Sai số phương sai ($S%$) Mức ý nghĩa ($F$)
Muồng $H_{vn} = -1.459 + 0.911 \times D_{1.3}$ 0.98 0.21 0.00 ($p < 0.01$)
Keo tai tượng $H_{vn} = -4.252 + 1.069 \times D_{1.3}$ 0.85 – 0.95 0.20 – 0.48 0.00 ($p < 0.01$)
Lát hoa $H_{vn} = -1.489 + 0.991 \times D_{1.3}$ 0.85 0.20 0.00 ($p < 0.01$)
Xoan ta $H_{vn} = -1.776 + 0.753 \times D_{1.3}$ 0.98 0.14 0.00 ($p < 0.01$)

[!NOTE] Giá trị $R \ge 0.85$ và sai số phương sai cực thấp ($S% \le 0.48%$) chứng minh các phương trình trên có độ tin cậy rất cao, cho phép suy toán trực tiếp chiều cao lâm phần từ đường kính mà không cần chặt hạ cây.

Tổng hợp số liệu sinh khối, Carbon và hấp thụ $\text{CO}_2$

Kết quả điều tra thực nghiệm tại 3 Ô tiêu chuẩn (OTC) được tổng hợp chi tiết trong bảng dưới đây:

OTC / Loài cây Mật độ $N$ (cây/ha) $\bar{D}_{1.3}$ (cm) $\bar{H}_{vn}$ (m) Sinh khối tươi (tấn/ha) Sinh khối khô (tấn/ha) Tích lũy Carbon (tấn C/ha) Hấp thụ $\text{CO}_2$ (tấn/ha) Giá trị kinh tế (triệu VNĐ/ha)
OTC 1 - Muồng 233 6.95 5.00 7.77 3.09 1.40 5.13 2.154
OTC 1 - Keo tai tượng 73 12.07 9.00 6.79 3.15 1.41 5.18 2.175
OTC 2 - Lát hoa 200 4.46 6.00 4.64 2.21 1.06 3.90 1.638
OTC 2 - Keo tai tượng 173 12.36 8.73 15.27 6.13 2.96 10.86 4.561
OTC 3 - Keo tai tượng 173 12.32 8.92 15.68 6.49 3.03 11.09 4.657
OTC 3 - Xoan ta 207 10.24 6.00 16.53 7.35 3.54 12.97 5.447
                       TRỮ LƯỢNG CO2 HẤP THỤ THEO LOÀI CÂY (TẤN/HA)
                        0    2    4    6    8    10   12   14 (tấn CO2/ha)

Phân bổ không gian tích lũy:

  • Trên mặt đất (Thân, cành, lá): Chiếm ưu thế tuyệt đối với tổng cộng 39.30 tấn $\text{CO}_2$/ha (trong đó thân cây chiếm $60 - 70%$ tổng lượng hấp thụ).
  • Dưới mặt đất (Hệ rễ): Đóng góp 9.83 tấn $\text{CO}_2$/ha (chiếm $20 - 25%$ tổng sinh khối).

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

  1. Định lượng hóa sinh khối vi mô trong mô hình chè - rừng: Đây là công trình thực nghiệm đầu tiên tại huyện Phú Lương xác định tỷ lệ carbon giải tích thực tế của từng bộ phận cho 4 loài cây gỗ trồng xen trên nương chè, thay vì giả định hệ số 0.5 cố định của IPCC.
  2. Phương pháp bảo tồn rễ chè: Ứng dụng giải pháp tính toán sinh khối rễ gián tiếp thông qua hệ số tương quan sinh thái ($B_{\text{rễ}} = 0.25 B_{\text{thân cành lá}}$), giải quyết triệt để rủi ro làm hư hại hệ thống canh tác chè đặc sản của nông dân.
  3. So sánh với các nghiên cứu tiền nhiệm:
    • So với nghiên cứu của Bảo Huy và cs (2009) trên mô hình Bời lời đỏ - Sắn ($3 - 84\text{ tấn }\text{CO}_2\text{/ha}$ tùy tuổi), mô hình Chè - Keo/Xoan tại Yên Ninh đạt mức hấp thụ ổn định $10.86 - 12.97\text{ tấn }\text{CO}_2\text{/ha}$ ngay ở giai đoạn cây con và trung niên.
    • So với nghiên cứu của Nguyễn Thanh Tiến (2012) trên rừng phục hồi IIB tại Thái Nguyên ($106.91\text{ tấn }\text{CO}_2\text{/ha}$ tầng cây gỗ), tầng cây gỗ trồng xen trong mô hình NLKH đạt mật độ thấp hơn ($73 - 233\text{ cây/ha}$) nhưng mang lại tốc độ tăng trưởng đường kính nhanh hơn nhờ tận dụng phân bón từ chè.

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

Tiềm năng mở rộng và kịch bản kinh tế

Với tổng diện tích 330 ha chè - rừng tại xã Yên Ninh, tiềm năng kinh tế môi trường mang lại rất đáng kể:

graph TD
    A[Mô hình Chè - Rừng 330 ha] --> B[Sản phẩm Nông nghiệp: Búp chè thương phẩm]
    A --> C[Sản phẩm Lâm nghiệp: Gỗ lớn gia dụng / công nghiệp]
    A --> D[Dịch vụ Môi trường Rừng: Tích lũy Carbon]
    
    B --> E[Doanh thu thường xuyên hàng tháng]
    C --> F[Tích lũy tài sản dài hạn 5-10 năm]
    D --> G[Nguồn thu từ CERs / Chi trả DVMTR PFES]
    
    E & F & G --> H[Nâng cao sinh kế & Xóa đói giảm nghèo bền vững]

Lộ trình triển khai khuyến nghị

  1. Giai đoạn 1 (0 – 6 tháng): Tuyên truyền nâng cao nhận thức cho các hộ dân tại 16 xóm thuộc xã Yên Ninh về lợi ích kép của cây gỗ che bóng (giảm thoát hơi nước cho chè và tạo tín chỉ carbon).
  2. Giai đoạn 2 (6 – 18 tháng): Chuẩn hóa mật độ trồng xen lý tưởng ($150 - 200\text{ cây/ha}$ đối với Keo tai tượng hoặc Xoan ta) để tối ưu hóa khả năng hấp thụ $\text{CO}_2$ mà không làm giảm năng suất búp chè.
  3. Giai đoạn 3 (18 – 36 tháng): Liên kết với Quỹ Bảo vệ và Phát triển Rừng tỉnh Thái Nguyên nhằm tích hợp mô hình vào đề án chi trả dịch vụ môi trường rừng (PFES).

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

Hạn chế kỹ thuật

  • Dung lượng mẫu khảo sát còn giới hạn ở 03 OTC đại diện ($4.500\text{ m}^2$) do kinh phí và thời gian có hạn.
  • Chưa tính toán trực tiếp lượng carbon tích lũy trong sinh khối của chính tán chè bên dưới và bể carbon trong đất ($SOC - \text{Soil Organic Carbon}$).
  • Biến động giá tín chỉ carbon trên thị trường quốc tế có thể ảnh hưởng đến doanh thu thực tế.

Hướng nghiên cứu tiếp theo

  • Mở rộng phạm vi lập ô tiêu chuẩn trên toàn bộ các huyện vùng chè trọng điểm của Thái Nguyên (Đại Từ, Định Hóa, Đồng Hỷ).
  • Ứng dụng công nghệ viễn thám GIS và ảnh chụp UAV đa phổ để giám sát tự động sinh khối tán cây theo thời gian thực.
  • Nghiên cứu chu trình dinh dưỡng khoáng và tương tác vi khí hậu giữa tầng cây che bóng và cây chè.

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

  • Sinh viên & Học viên ngành Lâm nghiệp/Quản lý tài nguyên: Tiếp cận tài liệu tham khảo chuẩn mực về phương pháp giải tích thân cây và quy trình phân tích carbon phòng thí nghiệm.
  • Kỹ sư lâm nghiệp & Cán bộ khuyến nông: Sở hữu các phương trình tương quan $H_{vn} = f(D_{1.3})$ thực nghiệm để kiểm kê nhanh trữ lượng rừng mà không cần chặt hạ.
  • Hộ nông dân trồng chè: Nhận thức rõ giá trị kinh tế bổ sung từ cây gỗ che bóng ($1.64 - 5.45\text{ triệu VNĐ/ha}$ từ dịch vụ môi trường bên cạnh nguồn thu từ gỗ và chè).
  • Các nhà hoạch định chính sách & Cơ quan quản lý môi trường: Có căn cứ khoa học vững chắc để xây dựng cơ chế chi trả dịch vụ hấp thụ carbon rừng cấp tỉnh và quốc gia.

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

1. Yêu cầu kỹ thuật để triển khai trồng xen cây gỗ vào nương chè là gì?

Cần chọn các loài cây đa mục đích có bộ rễ cọc ăn sâu (như Keo tai tượng, Xoan ta, Muồng) để không cạnh tranh dinh dưỡng tầng mặt với rễ chè. Mật độ trồng xen khuyến nghị từ $150 - 200\text{ cây/ha}$, cự ly hàng cách hàng $8 - 10\text{ m}$.

2. Giới hạn hấp thụ $\text{CO}_2$ của mô hình là bao nhiêu và làm sao để tối ưu?

Khả năng hấp thụ $\text{CO}_2$ tăng dần theo cấp tuổi và cấp kính của cây gỗ, đạt đỉnh sinh trưởng sinh khối ở giai đoạn 6–10 tuổi. Để tối ưu hóa, cần tỉa cành định kỳ nhằm giữ độ tàn che tán cây gỗ ở mức $0.2 - 0.3$, vừa đủ che bóng cho chè vừa tối đa hóa quang hợp.

3. Phương pháp tích hợp số liệu này vào hệ thống kiểm kê KNK quốc gia?

Số liệu sinh khối và tỷ lệ $C_{%}$ thực nghiệm của đề tài có thể cung cấp trực tiếp các hệ số phát thải/hấp thụ đặc thù quốc gia (Tier 2/Tier 3 theo hướng dẫn của IPCC LULUCF), thay thế cho hệ số mặc định toàn cầu (Tier 1).

4. Chi phí duy trì và quản lý mô hình NLKH chè - rừng có cao không?

Chi phí tăng thêm gần như bằng không do cây gỗ tận dụng trực tiếp công chăm sóc, làm cỏ và lượng phân bón dư thừa trong quá trình thâm canh chè hàng năm của người dân.

5. Khả năng hoàn vốn (ROI) từ dịch vụ carbon được tính toán như thế nào?

Với mức giá 20 USD/tấn $\text{CO}_2$, mỗi hecta chè - rừng mang lại thêm từ $1.64 - 5.45\text{ triệu VNĐ}$ từ tín chỉ carbon. Khi cây gỗ đạt tuổi khai thác (sau 7–10 năm), nguồn thu từ gỗ thương phẩm sẽ mang lại lợi nhuận ròng vượt trội so với độc canh cây chè.


Kết luận

Đề tài đã hoàn thành xuất sắc các mục tiêu nghiên cứu, chứng minh rằng mô hình nông lâm kết hợp chè – rừng tại xã Yên Ninh, huyện Phú Lương không chỉ đem lại hiệu quả kinh tế nông nghiệp thuần túy mà còn là bể hấp thụ carbon vô cùng hiệu quả. Với trữ lượng hấp thụ từ $3.90\text{ đến }12.97\text{ tấn }\text{CO}_2\text{/ha}$ (tương đương giá trị môi trường $1.638 - 5.447\text{ triệu VNĐ/ha}$), việc nhân rộng mô hình này là giải pháp chiến lược giúp bảo vệ môi trường sinh thái, ứng phó biến đổi khí hậu và hiện thực hóa các cam kết giảm phát thải ròng của Việt Nam trên trường quốc tế.