Giới thiệu dự án

Sự gia tăng nồng độ khí nhà kính trong khí quyển, đặc biệt là khí carbon dioxide ($CO_2$), đang là nguyên nhân hàng đầu dẫn đến hiện tượng ấm lên toàn cầu và biến đổi khí hậu cực đoan. Theo các báo cáo của Ban Liên chính phủ về Biến đổi Khí hậu (IPCC) và Tổ chức Nông Lương Liên Hợp Quốc (FAO), các hệ sinh thái thực vật trên cạn đóng vai trò là "bể hấp thụ carbon" (Carbon Sink) tự nhiên với chi phí kinh tế thấp nhất và tính bền vững lâu dài. Tại khu vực trung du miền núi Bắc Bộ, canh tác chè (Camellia sinensis) thuần loài truyền thống đang bộc lộ nhiều điểm hạn chế: nguy cơ xói mòn đất mặt vào mùa mưa bão, khả năng giữ ẩm kém vào mùa khô và tính tổn thương sinh thái cao.

Mô hình Nông Lâm Kết Hợp (NLKH) "Chè - Rừng" tại xã Vô Tranh, huyện Phú Lương, tỉnh Thái Nguyên với diện tích canh tác chè kinh doanh lên tới 565 ha (trong tổng số 613 ha đất mô hình NLKH toàn xã) đã giải quyết được bài toán canh tác bền vững. Việc trồng xen các loài cây gỗ bản địa và cây sinh trưởng nhanh giúp tạo tầng tán che bóng, cản gió bão, điều hòa vi khí hậu và bảo vệ đất. Tuy nhiên, việc lượng hóa khả năng hấp thụ khí $CO_2$ và trữ lượng carbon tích lũy của các loài cây gỗ trồng xen trong mô hình NLKH này vẫn còn thiếu cơ sở dữ liệu thực nghiệm chuẩn hóa, gây khó khăn cho việc định giá rừng và xây dựng cơ chế Chi trả Dịch vụ Môi trường Rừng (PFES/DVMTR) hay tham gia Cơ chế Phát triển Sạch (CDM) và REDD+.

+------------------------------------------------------------------------------------+
|                         KHUNG CANH TÁC NÔNG LÂM KẾT HỢP                            |
|                                                                                    |
|  [ Tầng Cây Gỗ Vượt Tán: Keo lai, Keo tai tượng, Xoan ta, Mỡ ]                    |
|    - Cản gió, che bóng, tích lũy sinh khối lớn                                     |
|    - Hấp thụ CO2 khí quyển (Tích lũy Carbon hữu cơ vào Thân, Cành, Lá, Rễ)         |
|                                                                                    |
|  [ Tầng Canh Tác Nông Nghiệp: Nương Chè Kinh Doanh ]                               |
|    - Giữ ẩm vi khí hậu dưới tán, chống xói mòn bề mặt đất                         |
|    - Tạo nguồn thu nhập tiền mặt định kỳ cho nông hộ                               |
|                                                                                    |
|  [ Bể Chứa Dưới Mặt Đất: Hệ rễ đa tầng & Hữu cơ đất ]                             |
|    - Tương tác rễ cây gỗ + rễ chè, cố định đạm (đối với họ Đậu)                    |
|    - Lưu giữ 20% tổng sinh khối cây gỗ dưới mặt đất                                |
+------------------------------------------------------------------------------------+

Mục tiêu nghiên cứu của dự án

  1. Mục tiêu lý luận: Xây dựng cơ sở dữ liệu thực nghiệm và luận cứ khoa học để định lượng giá trị hấp thụ carbon của các hệ thống NLKH đặc trưng vùng trung du miền núi phía Bắc.
  2. Mục tiêu thực tiễn 1: Đánh giá các chỉ tiêu sinh trưởng ($D_{1.3}$, $H_{vn}$, $N$, $M$) và xác định cấu trúc sinh khối tươi (SKT), sinh khối khô (SKK) theo từng bộ phận (thân, cành, lá, rễ) của 4 loài cây gỗ trồng xen: Xoan ta (Melia azedarach), Keo tai tượng (Acacia mangium), Keo lai (Acacia hybrid), và Mỡ (Manglietia conifera).
  3. Mục tiêu thực tiễn 2: Xác định chính xác trữ lượng carbon tích lũy ($tC/ha$) và khả năng hấp thụ $CO_2$ ($tCO_2/ha$) thông qua phân tích mẫu khô kiệt và hệ số chuyển đổi hóa học lượng giác.
  4. Mục tiêu ứng dụng: Đề xuất phương pháp định giá kinh tế dịch vụ môi trường hấp thụ $CO_2$ dựa trên mức giá thị trường carbon tự nguyện quốc tế.

Phương pháp tiếp cận và phạm vi nghiên cứu

  • Phương pháp tiếp cận: Nghiên cứu thực địa kết hợp rút mẫu thực nghiệm chặt hạ cây tiêu chuẩn bình quân giải tích, sấy khô kiệt trong phòng thí nghiệm đến khối lượng không đổi và phân tích toán học hồi quy tương quan.
  • Phạm vi không gian: 03 Ô tiêu chuẩn (OTC) đại diện, mỗi ô có diện tích $1.500\text{ m}^2$ ($50\text{ m} \times 30\text{ m}$) tại xã Vô Tranh, huyện Phú Lương, tỉnh Thái Nguyên.
  • Phạm vi thời gian: Thu thập và xử lý mẫu thực địa từ tháng 01 đến tháng 05 năm 2014.

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

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

Trước khi mô hình NLKH "Chè - Cây gỗ" được áp dụng rộng rãi, thực trạng sản xuất nông lâm nghiệp tại địa phương đối mặt với các thách thức lớn được tổng hợp trong ma trận phân tích:

Phương thức canh tác Ưu điểm Nhược điểm / Rủi ro Khả năng tích lũy Carbon Giá trị kinh tế dài hạn
Chè độc canh Thu hoạch búp định kỳ, quay vòng vốn nhanh. Rửa trôi xói mòn đất, dễ bùng phát sâu bệnh, chịu hạn kém khi biến đổi khí hậu. Thấp ($< 15\text{ tC/ha}$), chủ yếu tích lũy ở tầng cây bụi thấp. Trung bình, phụ thuộc hoàn toàn vào giá nông sản.
Rừng trồng thuần loài (Keo/Bạch đàn) Tích lũy sinh khối tập trung, kỹ thuật lâm sinh đơn giản. Chu kỳ kinh doanh dài (5–7 năm) không có thu nhập ngắn hạn; sau khai thác trắng đất bị thoái hóa nặng. Trung bình – Cao ($50 - 120\text{ tC/ha}$), nhưng bị mất sạch khi khai thác trắng. Chỉ thu hoạch một lần cuối chu kỳ, rủi ro cháy rừng cao.
Mô hình NLKH Chè - Rừng (Giải pháp đề xuất) Kết hợp nguồn thu ngắn ngày (búp chè) và dài ngày (gỗ); tán cây che bóng giảm thoát hơi nước, tăng độ phì đất. Kỹ thuật điều tiết mật độ phức tạp, cần kiểm soát cạnh tranh ánh sáng và dinh dưỡng. Ổn định và đa tầng ($70 - 150\text{ tC/ha}$ tính cả chè, cây gỗ và thảm mục). Tối ưu hóa trên cùng đơn vị diện tích; hưởng lợi từ thương mại carbon.

Yêu cầu nghiên cứu được phân loại theo mô hình chuẩn hóa MoSCoW:

  • Must have (Bắt buộc): Đo đếm chính xác đường kính ngang ngực ($D_{1.3}$), chiều cao vút ngọn ($H_{vn}$), mật độ cây ($N$), trữ lượng thân cây ($M$); xác định tỷ lệ sấy khô kiệt từng bộ phận (thân, cành, lá) và sinh khối rễ.
  • Should have (Nên có): Xây dựng phương trình tương quan toán học $H_{vn} = f(D_{1.3})$ cho từng loài với hệ số tương quan $R > 0.90$.
  • Could have (Có thể có): Lượng hóa giá trị kinh tế theo chứng chỉ giảm phát thải CERs với đơn giá thị trường cập nhật.
  • Won't have (Tạm thời chưa làm): Khai đào toàn bộ hệ rễ chè do tránh phá hủy vùng sinh trưởng cây chè đang thu hoạch của nông dân; giải pháp thay thế là sử dụng phương trình allometric sinh khối dưới mặt đất chuẩn quốc tế.

Thiết kế hệ thống đo đếm và phân tích toán học

+-----------------------------------------------------------------------------------------+
|                    QUY TRÌNH ĐO ĐẾM VÀ TÍNH TOÁN SINH KHỐI - CARBON                     |
|                                                                                         |
| [1. Khảo sát thực địa]  -->  Lập 3 OTC (50m x 30m = 1.500 m2)                           |
|                                    |                                                    |
| [2. Đo đếm lâm sinh]    -->  Đo D1.3 (thước kẹp), Hvn (máy đo cao), tính D_tb, H_tb     |
|                                    |                                                    |
| [3. Chặt cây tiêu chuẩn]-->  Phân đoạn thân (L=1m), cành (<5cm), lá                     |
|                              Cân sinh khối tươi tại chỗ (P_tươi)                        |
|                                    |                                                    |
| [4. Xử lý mẫu khô kiệt] -->  Lấy mẫu 0.5kg -> Rút gọn 30g                               |
|                              Sấy tủ sấy (Thân/Cành: 80-105°C; Lá: 70-85°C)              |
|                              Kiểm tra 4 lần cân liên tiếp (delta_m = 0)                 |
|                                    |                                                    |
| [5. Mô hình hóa toán học]--> Tính SKK -> Tính Trữ lượng Carbon (C) -> Chuyển đổi CO2    |
|                              Xây dựng hồi quy H/D trên Excel ToolPak                    |
+-----------------------------------------------------------------------------------------+

Công cụ và Thiết bị Kỹ thuật Áp dụng

  • Phần mềm xử lý dữ liệu: Microsoft Excel 2013 với tiện ích mở rộng Analysis ToolPak (thực thi hồi quy OLS - Ordinary Least Squares Regression).
  • Phần mềm đối chuẩn: Mô hình CO2Fix V3.1 (Wageningen University) dùng làm cơ sở tham chiếu dòng carbon trong lâm phần.
  • Thiết bị đo đếm thực địa: Thước dây chuyên dụng đo chu vi/đường kính Haglof (độ chính xác 0.1 cm), thước đo cao Blume-Leiss (độ chính xác 0.5 m), cân đồng hồ Nhơn Hòa 5kg và 100kg.
  • Thiết bị phòng lab: Tủ sấy đối lưu Memmert UN110 (dải nhiệt độ hoạt động $30 - 250^\circ\text{C}$), cân phân tích điện tử Sartorius (độ chính xác $\pm 0.01\text{ g}$).

Methodology

Quy trình nghiên cứu áp dụng các phương pháp luận lâm học và hóa sinh thực vật chuẩn của ICRAF (2010) và IPCC (2006):

  1. Phương pháp xác định Cây tiêu chuẩn: Cây tiêu chuẩn đại diện cho ô tiêu chuẩn được lựa chọn thỏa mãn điều kiện: $$D_{1.3\text{ (cây tc)}} \approx \bar{D}{1.3} \quad \text{và} \quad H{vn\text{ (cây tc)}} \approx \bar{H}_{vn}$$

  2. Phương pháp phân đoạn sinh khối thân và các bộ phận:

    • Thân cây: Chặt sát gốc, chia thành từng lóng $L = 1.0\text{ m}$, các phần ngọn có đường kính $< 5\text{ cm}$ được quy về sinh khối cành.
    • Cành: Tách toàn bộ lá, chặt nhỏ cành và cân toàn bộ khối lượng tươi tại chỗ.
    • Lá: Gom toàn bộ lá tươi, trộn đều và cân khối lượng tổng.
    • Rễ (Sinh khối dưới mặt đất - $SK_{DMĐ}$): Áp dụng phương trình tương quan chuẩn allometric quốc tế của ICRAF đối với hệ thống canh tác không đào rễ: $$SK_{DMĐ} = TSK_{TMĐ} \times \frac{1}{4} = 0.25 \times TSK_{TMĐ} \quad (\text{tấn/ha})$$
  3. Phương pháp sấy khô kiệt xác định tỷ lệ khô/tươi ($K_i$):

    • Rút mẫu đại diện: Thân lấy tại vị trí $D_{1.3}$ ($0.5\text{ kg}$), cành băm nhỏ trộn đều ($0.5\text{ kg}$), lá trộn đều ($0.5\text{ kg}$).
    • Xử lý mẫu rút gọn trong phòng thí nghiệm: Lấy chính xác $m_0 = 30\text{ g}$ mẫu tươi.
    • Chế độ sấy: Thân và cành sấy ở $80 - 105^\circ\text{C}$ trong $8 - 12\text{ giờ}$ (kiểm tra lặp lại sau mỗi 2h). Lá sấy ở $70 - 85^\circ\text{C}$ trong $4 - 8\text{ giờ}$. Mẫu đạt độ khô kiệt tuyệt đối khi khối lượng không đổi qua 4 lần cân liên tiếp.
    • Tỷ lệ chuyển đổi khô/tươi bộ phận $i$: $$K_i = \frac{m_{\text{khô kiệt}}}{m_{\text{tươi}}} = \frac{m_{\text{khô kiệt}}}{30\text{ g}}$$
  4. Thuật toán quy đổi trữ lượng Carbon và $CO_2$ hấp thụ:

    • Sinh khối khô bộ phận $i$ trên 1 ha: $$Pk_i = \frac{Pt_i \times K_i \times N}{1000} \quad (\text{tấn/ha})$$
    • Trữ lượng Carbon tích lũy ($C$): $$C = \sum C_i = \sum \left( Pk_i \times c_i \right) \quad (\text{tấn C/ha})$$ (Trong đó $c_i$ là hàm lượng carbon phân tích hóa học trong vật chất khô từng bộ phận).
    • Lượng $CO_2$ hấp thụ tương đương (theo tỷ lệ khối lượng phân tử $M_{CO_2} / M_C = 44 / 12$): $$W_{CO_2} = C \times \frac{44}{12} = C \times 3.667 \quad (\text{tấn } CO_2/\text{ha})$$

Implementation và kết quả

Development process

Quá trình điều tra thực địa và xử lý dữ liệu được tiến hành qua 5 giai đoạn:

+---------------------------------------------------------------------------------------+
|                                TIẾN ĐỘ THỰC HIỆN DỰ ÁN                                |
|                                                                                       |
|  Tháng 01/2014: Khảo sát sơ bộ địa bàn xã Vô Tranh, phân loại mô hình NLKH.           |
|  Tháng 02/2014: Thiết lập 03 OTC (50m x 30m), đo đếm lâm sinh tổng thể D1.3 và Hvn.   |
|  Tháng 03/2014: Chọn và chặt hạ cây tiêu chuẩn; cân tươi tại chỗ, phân lập mẫu 0.5kg. |
|  Tháng 04/2014: Sấy khô kiệt mẫu tại phòng thí nghiệm, kiểm tra cân lặp 4 lần.        |
|  Tháng 05/2014: Xây dựng phương trình tương quan trên Excel, tính Carbon và viết báo cáo|
+---------------------------------------------------------------------------------------+

Đoạn mã giả dưới đây mô tả quá trình xử lý thống kê lâm sinh và tính toán carbon tự động:

import numpy as np

def calculate_carbon_and_co2(N_density, fresh_weights, dry_ratio, carbon_content):
    """
    Tính toán sinh khối khô, trữ lượng Carbon và CO2 hấp thụ cho 1 ha
    fresh_weights: dict {'than': kg, 'canh': kg, 'la': kg}
    dry_ratio: dict {'than': k_th, 'canh': k_c, 'la': k_la}
    carbon_content: dict {'than': c_th, 'canh': c_c, 'la': c_la}
    """
    # 1. Tính sinh khối khô trên mặt đất của cây tiêu chuẩn (kg/cây)
    pk_aboveground = {}
    for organ in ['than', 'canh', 'la']:
        pk_aboveground[organ] = fresh_weights[organ] * dry_ratio[organ]
    
    total_pk_aboveground_tree = sum(pk_aboveground.values())
    
    # 2. Sinh khối rễ dưới mặt đất (Allometric proxy 25%)
    pk_root_tree = total_pk_aboveground_tree * 0.25
    
    # 3. Quy đổi ra 1 ha (tấn/ha)
    pk_ha = {organ: (pk_aboveground[organ] * N_density) / 1000.0 for organ in pk_aboveground}
    pk_ha['re'] = (pk_root_tree * N_density) / 1000.0
    total_skk_ha = sum(pk_ha.values())
    
    # 4. Trữ lượng Carbon (tấn C/ha)
    c_ha = {organ: pk_ha[organ] * carbon_content.get(organ, 0.50) for organ in pk_ha}
    total_carbon_ha = sum(c_ha.values())
    
    # 5. Hấp thụ CO2 tương đương (tấn CO2/ha)
    total_co2_ha = total_carbon_ha * (44.0 / 12.0)
    
    return {
        'total_skk_ha': total_skk_ha,
        'total_carbon_ha': total_carbon_ha,
        'total_co2_ha': total_co2_ha,
        'pk_ha_breakdown': pk_ha
    }

Testing và validation

Phân tích hồi quy tương quan chiều cao vút ngọn ($H_{vn}$) theo đường kính ngang ngực ($D_{1.3}$) trên phần mềm Microsoft Excel (Analysis ToolPak) đạt độ tin cậy thống kê rất cao:

$$H_{VN} = a + b \times D_{1.3}$$

  • Xoan ta (Melia azedarach): $H_{VN} = 0.378 + 0.590 \times D_{1.3}$ ($R = 0.92, S% = 0.295, F = 0.00$)
  • Keo tai tượng (Acacia mangium): $H_{VN} = -3.063 + 1.594 \times D_{1.3}$ ($R = 0.93, S% = 0.431, F = 0.00$)
  • Keo lai (Acacia hybrid): $H_{VN} = 3.308 + 0.964 \times D_{1.3}$ ($R = 0.92, S% = 0.298, F = 0.00$)
  • Mỡ (Manglietia conifera): Phương trình tương quan tuyến tính đạt $R = 0.92$.

Tất cả các phương trình đều có hệ số tương quan $R > 0.92$, mức ý nghĩa kiểm định Fisher $F < 0.01$, sai số tiêu chuẩn $S% < 0.5%$, chứng minh đường kính $D_{1.3}$ là biến số tin cậy tuyệt đối để dự đoán chiều cao và sinh khối cây cá thể.

                  BIỂU ĐỒ TƯƠNG QUAN Hvn THEO D1.3
  Hvn (m)
   14 |                                          * (Keo lai: R=0.92)
   12 |                                  *   * (Keo tai tượng: R=0.93)
   10 |                              *
    8 |                      *
    6 |              *   * (Mỡ & Xoan ta: R=0.92)
    4 |          *
    2 |
    0 +-------------------------------------------------
      0     2     4     6     8     10    12    14  D1.3 (cm)

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

Dưới đây là bảng tổng hợp đầy đủ các chỉ tiêu sinh trưởng, sinh khối tươi, sinh khối khô, trữ lượng carbon và lượng $CO_2$ hấp thụ của các loài cây gỗ trong 3 ô tiêu chuẩn thực nghiệm:

OTC Loài cây Mật độ $N$ (cây/ha) $D_{1.3}$ (cm) $H_{vn}$ (m) Trữ lượng $M$ ($\text{m}^3/\text{ha}$) Sinh khối tươi SKT (tấn/ha) Sinh khối khô SKK (tấn/ha) Trữ lượng Carbon (tC/ha) Lượng $CO_2$ hấp thụ ($\text{t}CO_2/\text{ha}$)
OTC 1 Xoan ta 273 9.52 5.96 5.19 19.45 8.69 4.55 16.69
Keo tai tượng 113 9.03 11.13 3.50 7.77 3.03 1.54 5.61
OTC 2 Mỡ 213 10.36 6.30 5.11 9.32 3.87 2.23 8.19
Keo tai tượng 133 9.25 11.80 4.52 9.98 3.76 1.86 6.83
OTC 3 Keo lai 213 9.80 12.70 8.73 12.51 5.24 2.63 9.63
Xoan ta 127 8.80 5.66 1.91 8.10 3.54 1.83 6.70

Phân bố hấp thụ $CO_2$ theo cấu trúc cơ quan cây gỗ:

  • Thân cây: Đóng vai trò là bể chứa carbon lớn nhất, chiếm $53.2% - 63.9%$ tổng lượng $CO_2$ hấp thụ toàn cây.
  • Cành cây: Chiếm $13.7% - 21.3%$.
  • Hệ rễ dưới mặt đất: Chiếm cố định $20.0%$ tổng lượng carbon tích lũy.
  • Lá cây: Chiếm tỷ trọng nhỏ nhất, dao động từ $4.1% - 6.4%$.

Tỷ trọng đóng góp hấp thụ $CO_2$ bình quân giữa các loài cây trồng xen:

  1. Xoan ta (Melia azedarach): Đạt tỷ lệ hấp thụ bình quân cao nhất với $35.95%$ (nhờ ưu thế phân cành rộng và mật độ cá thể cao trong mô hình nương chè).
  2. Keo lai (Acacia hybrid): Chiếm $25.61%$ (tốc độ tăng trưởng thể tích thân nhanh vượt trội đạt $M = 8.73\text{ m}^3/\text{ha}$).
  3. Mỡ (Manglietia conifera): Chiếm $21.78%$.
  4. Keo tai tượng (Acacia mangium): Chiếm $16.67%$.

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

Các điểm mới mang tính đột phá của đề tài

  • Lượng hóa thực nghiệm đầu tiên trên mô hình NLKH Chè - Rừng: Khác với phần lớn các công trình trước đây chỉ tập trung vào rừng trồng thuần loài (Keo, Thông, Bạch đàn), nghiên cứu này cung cấp bộ số liệu thực nghiệm gốc đầu tiên về khả năng tích lũy sinh khối của các loài cây gỗ khi trồng xen trong hệ thống nông lâm kết hợp chè tại miền Bắc Việt Nam.
  • Giải pháp thu mẫu không xâm lấn rễ chè: Ứng dụng thành công hệ số tương quan allometric $SK_{DMĐ} = 0.25 \times TSK_{TMĐ}$ giải quyết triệt để xung đột giữa yêu cầu nghiên cứu khoa học và bảo tồn nguyên trạng nương chè kinh doanh của nông hộ.
  • Xây dựng hệ phương trình sinh trưởng chuẩn hóa: Cung cấp 3 phương trình hồi quy $H/D$ với độ tin cậy $R > 0.92$, cho phép các nhà quản lý rừng tính nhanh sinh khối cây đứng mà không cần chặt hạ.
+------------------------------------------------------------------------------------+
|                SO SÁNH KHẢ NĂNG HẤP THỤ CO2 VỚI CÁC NGHIÊN CỨU KHÁC                |
|                                                                                    |
| [1] Mô hình NLKH Chè - Rừng (Đề tài này):                                          |
|     - Hấp thụ: 12.33 - 22.30 tCO2/ha (Tầng cây gỗ trồng xen)                       |
|     - Điểm nổi bật: Kết hợp tạo thu nhập kép (búp chè + gỗ + carbon)               |
|                                                                                    |
| [2] Rừng trồng Keo lai thuần loài (Nguyễn Viết Khoa, 2007):                        |
|     - Trữ lượng Carbon: 49.6 - 113.8 tC/ha (Toàn lâm phần thuần loài)              |
|     - Điểm yếu: Độc canh, nguy cơ suy thoái đất cao sau chu kỳ khai thác           |
|                                                                                    |
| [3] Rừng tự nhiên phục hồi IIB tại Thái Nguyên (Nguyễn Thanh Tiến, 2012):          |
|     - Lượng CO2 tầng cây gỗ: 106.91 tCO2/ha; Tổng lâm phần: 460.69 tCO2/ha        |
|     - Đặc thù: Cần bảo vệ nghiêm ngặt, không thể khai thác thương mại nông sản     |
+------------------------------------------------------------------------------------+

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

Kịch bản ứng dụng trong đời sống và quản lý

  1. Thương mại Carbon và Chi trả Dịch vụ Môi trường Rừng (PFES):
    • Cơ sở dữ liệu của đề tài cho phép các cơ quan chức năng (Chi cục Kiểm lâm, Quỹ Bảo vệ và Phát triển Rừng) định mức chi trả tiền dịch vụ môi trường rừng cho các hộ nông dân trồng chè có xen cây lâm nghiệp.
  2. Chứng chỉ Nông nghiệp Bền vững (Rainforest Alliance / VietGAP):
    • Số liệu hấp thụ $CO_2$ là bằng chứng khoa học phục vụ việc dán nhãn sinh thái "Chè phát thải thấp" (Low Carbon Tea), nâng cao giá trị xuất khẩu chè Vô Tranh - Phú Lương sang thị trường EU và Hoa Kỳ.

Hạch toán giá trị kinh tế môi trường từ hấp thụ $CO_2$

Áp dụng đơn giá giao dịch trên thị trường carbon quốc tế thời điểm nghiên cứu là $20\text{ USD/tấn } CO_2$ (tương đương $420.000\text{ VNĐ/tấn } CO_2$ với tỷ giá $1\text{ USD} = 21.000\text{ VNĐ}$):

Mô hình thực nghiệm Thành phần loài cây gỗ Lượng $CO_2$ tích lũy ($\text{t}CO_2/\text{ha}$) Giá trị môi trường quy đổi (USD/ha) Giá trị môi trường quy đổi (VNĐ/ha)
OTC 1 Xoan ta + Keo tai tượng 22.30 $446.00 9.366.000
OTC 2 Mỡ + Keo tai tượng 15.02 $300.40 6.308.400
OTC 3 Keo lai + Xoan ta 16.33 $326.60 6.858.600

Trung bình mỗi hecta mô hình NLKH Chè - Rừng mang lại nguồn thu tiềm năng từ 6.3 đến 9.4 triệu VNĐ/ha từ dịch vụ lưu giữ carbon, bên cạnh nguồn thu chính từ sản lượng chè búp khô.


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

Hạn chế kỹ thuật

  • Ước tính sinh khối rễ gián tiếp: Việc sử dụng hệ số tương quan $0.25$ cho sinh khối dưới mặt đất tuy phù hợp thực tế bảo vệ chè nhưng chưa phản ánh chính xác cấu trúc giải phẫu rễ của từng loài cụ thể.
  • Tính thời điểm: Dữ liệu điều tra thu thập trong 5 tháng đầu năm 2014, chưa theo dõi được biến động sinh khối và dòng carbon theo chu kỳ 4 mùa trong nhiều năm liên tục.

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

  • Ứng dụng Viễn thám và LiDAR: Kết hợp dữ liệu ảnh vệ tinh Sentinel-2 và thiết bị bay không người lái (UAV LiDAR) để xây dựng bản đồ số hóa trữ lượng carbon cho toàn bộ diện tích 613 ha mô hình NLKH xã Vô Tranh.
  • Phân tích hàm lượng Carbon chi tiết bằng máy phân tích nguyên tố (CHNS Analyzer): Thay thế việc quy đổi hệ số trung bình bằng phân tích trực tiếp hàm lượng C trong từng mô thực vật.

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

+---------------------------------------------------------------------------------------+
|                             MA TRẬN ĐỐI TƯỢNG HƯỞNG LỢI                               |
|                                                                                       |
|  [ Sinh Viên & Học Viên Lâm Nghiệp ]                                                  |
|    - Cung cấp quy trình thực nghiệm chuẩn từ lập OTC, giải tích thân cây đến sấy mẫu. |
|    - Tài liệu tham khảo phương pháp luận phân tích hồi quy tương quan sinh học.        |
|                                                                                       |
|  [ Nhà Khoa Học & Nghiên Cứu Viên ]                                                   |
|    - Bộ số liệu gốc về sinh khối và carbon của Xoan ta, Keo lai, Keo tai tượng, Mỡ.   |
|    - Dữ liệu đầu vào tin cậy để hiệu chỉnh mô hình phần mềm CO2Fix và C-BM CFS3.      |
|                                                                                       |
|  [ Nông Dân & Hợp Tác Xã Trồng Chè ]                                                  |
|    - Định hướng kỹ thuật chọn loài cây trồng xen tối ưu (Xoan ta: hấp thụ CO2 tốt nhất|
|      và tán thưa ít cạnh tranh ánh sáng với chè).                                     |
|    - Cơ hội gia tăng thu nhập từ 6 - 9 triệu VNĐ/ha qua chi trả môi trường rừng.      |
|                                                                                       |
|  [ Cơ Quan Quản Lý & Hoạch Định Chính Sách ]                                          |
|    - Cơ sở khoa học để UBND tỉnh Thái Nguyên xây dựng biểu giá dịch vụ môi trường rừng.|
|    - Luận cứ phục vụ đàm phán các dự án tín chỉ carbon CDM/REDD+ quốc tế.             |
+---------------------------------------------------------------------------------------+

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

1. Tại sao loài Xoan ta lại có tỷ lệ hấp thụ $CO_2$ cao nhất ($35.95%$) trong mô hình?

Xoan ta (Melia azedarach) là cây gỗ bản địa có khả năng thích ứng cao, hệ cành nhánh phát triển phân tầng rộng tạo diện tích quang hợp lớn. Trong mô hình nương chè tại OTC 1, mật độ Xoan ta đạt $273\text{ cây/ha}$ kết hợp với đường kính bình quân lớn ($9.52\text{ cm}$) giúp tổng sinh khối tươi đạt tới $19.45\text{ tấn/ha}$, vượt trội so với các loài khác.

2. Việc tính sinh khối rễ bằng phương trình $SK_{DMĐ} = 0.25 \times TSK_{TMĐ}$ có đảm bảo tính khoa học?

Phương pháp này tuân thủ đúng tài liệu hướng dẫn kỹ thuật của Trung tâm Nghiên cứu Nông Lâm Kết Hợp Thế Giới (ICRAF, 2010) và Hướng dẫn kiểm kê khí nhà kính quốc gia của IPCC (2006). Trong điều kiện nông lâm kết hợp, việc đào bới toàn bộ hệ rễ sẽ phá hủy bầu rễ chè kinh doanh, do đó áp dụng hệ số allometric quốc tế $25%$ là giải pháp chuẩn mực được chấp nhận toàn cầu.

3. Căn cứ nào để quy đổi từ khối lượng Carbon ($C$) sang khối lượng $CO_2$?

Dựa trên định luật bảo toàn khối lượng và phương trình hóa học quang hợp: $$6CO_2 + 6H_2O \xrightarrow{\text{Ánh sáng, Diệp lục}} C_6H_{12}O_6 + 6O_2$$ Khối lượng phân tử của $CO_2$ là $44\text{ g/mol}$ (gồm 1 nguyên tử $C = 12$ và 2 nguyên tử $O = 16$). Do đó, $1\text{ đơn vị Carbon}$ trong sinh khối tương đương với $\frac{44}{12} \approx 3.667\text{ đơn vị } CO_2$ được hấp thụ từ khí quyển.

4. Người dân trồng chè có bị giảm năng suất búp khi trồng xen cây gỗ không?

Không, nếu kiểm soát mật độ hợp lý ($100 - 250\text{ cây/ha}$). Cây gỗ tầng cao che bớt bức xạ gay gắt vào mùa hè, giảm thoát hơi nước mặt đất, giữ ẩm cho nương chè, đồng thời lá rụng từ các loài cây họ Đậu (như Keo tai tượng, Keo lai) cung cấp lượng đạm hữu cơ tự nhiên cải tạo đất chè.

5. Dự án tín chỉ carbon từ mô hình NLKH này cần những thủ tục pháp lý gì để giải ngân?

Cần thiết lập: (1) Đường phát thải cơ sở (Baseline), (2) Đề án thiết kế dự án (PDD) theo tiêu chuẩn VCS (Verified Carbon Standard) hoặc Gold Standard, (3) Hệ thống giám sát - đo đạc - báo cáo - thẩm tra (MRV) định kỳ bằng ô tiêu chuẩn cố định, và (4) Giấy chứng nhận quyền sử dụng đất hợp pháp của các hộ tham gia.


Kết luận

Công trình khóa luận tốt nghiệp của tác giả Phạm Đức Cảnh đã giải quyết trọn vẹn mục tiêu định lượng khả năng hấp thụ $CO_2$ của các loài cây gỗ trồng xen trong mô hình NLKH Chè - Rừng tại xã Vô Tranh, huyện Phú Lương, tỉnh Thái Nguyên.

Các kết quả nghiên cứu nổi bật gồm:

  1. Xác lập hệ thống dữ liệu lâm sinh chi tiết cho 4 loài cây gỗ chủ lực: đường kính $D_{1.3}$ từ $8.80 - 10.36\text{ cm}$, chiều cao $H_{vn}$ từ $5.66 - 12.70\text{ m}$, và trữ lượng thân cây đạt từ $1.91 - 8.73\text{ m}^3/\text{ha}$.
  2. Xây dựng thành công 3 phương trình hồi quy allometric $H/D$ với hệ số tương quan cao ($R = 0.92 - 0.93$).
  3. Định lượng chính xác khả năng hấp thụ $CO_2$ dao động từ $5.61\text{ đến } 16.69\text{ tấn } CO_2/\text{ha}$, trong đó Xoan ta đóng góp tỷ trọng lớn nhất ($35.95%$), theo sau là Keo lai ($25.61%$), Mỡ ($21.78%$) và Keo tai tượng ($16.67%$).
  4. Xác định giá trị kinh tế môi trường quy đổi đạt từ $6.3\text{ đến } 9.4\text{ triệu VNĐ/ha}$, khẳng định tính ưu việt toàn diện cả về mặt kinh tế, xã hội và bảo vệ môi trường sinh thái của mô hình Nông Lâm Kết Hợp.

Kết quả này là tiền đề khoa học vững chắc để các địa phương vùng trung du nhân rộng mô hình canh tác chè sinh thái, đồng thời tham gia hiệu quả vào thị trường tín chỉ carbon trong nước và quốc tế.