Mô Hình Hóa Động Thái Sinh Khối và Hấp Thu Carbon Rừng Đước (Rhizophora apiculata) Tại Kiên Giang Bằng Mô Hình 3-PG


Tóm tắt nghiên cứu (200-250 từ)

  • Câu hỏi nghiên cứu chính: Động thái tích lũy sinh khối và khả năng hấp thu $CO_2$ của rừng ngập mặn Đước đôi (Rhizophora apiculata Blume) tại Kiên Giang biến đổi như thế nào theo các cấp tuổi (10 đến 19 tuổi), và khả năng thích ứng, ứng dụng của mô hình động thái sinh lý 3-PG trong việc mô phỏng sản lượng loại rừng này ra sao?
  • Phương pháp nghiên cứu: Nghiên cứu thiết lập hệ thống các ô tiêu chuẩn định lượng ($500\text{ m}^2$) trên các lâm phần rừng Đước 10, 14, 15, 16 và 19 tuổi tại huyện An Minh và An Biên, tỉnh Kiên Giang. Kết hợp phương pháp chặt hạ cây giải tích tiêu chuẩn để phân tích sinh khối tươi/khô từng bộ phận (thân, cành, lá, rễ chống, hoa quả), ước tính sinh khối rễ ngầm qua mô hình tương quan chuyên biệt, tính toán lượng hấp thu carbon theo chuẩn IPCC và hiệu chỉnh bộ tham số sinh lý khí hậu đầu vào cho mô hình 3-PG.
  • Phát hiện chính: Sinh khối tổng số của quần thể rừng Đước dao động từ $136,12\text{ tấn/ha}$ (10 tuổi) đến $229,45\text{ tấn/ha}$ (19 tuổi). Tương ứng, lượng $CO_2$ hấp thu tích lũy đạt từ $249,56\text{ tấn/ha}$ lên đến $420,66\text{ tấn/ha}$. Sinh khối thân chiếm ưu thế tuyệt đối ($62,14% - 72,68%$). Mô hình 3-PG thể hiện tiềm năng vượt trội trong việc số hóa và dự báo động thái sinh trưởng rừng ngập mặn.
  • Ý nghĩa thực tiễn: Cung cấp cơ sở khoa học định lượng chuẩn xác phục vụ đánh giá tiềm năng bể chứa "Carbon xanh" (Blue Carbon), tạo tiền đề tham gia thị trường tín chỉ carbon rừng và tối ưu hóa các biện pháp lâm sinh cho vùng ven biển Đồng bằng sông Cửu Long.

Bối cảnh và tầm quan trọng (300-350 từ)

Rừng ngập mặn nhiệt đới là một trong những hệ sinh thái có năng suất sinh học cao nhất hành tinh, đóng vai trò "lá chắn xanh" bảo vệ bờ biển trước sóng gió, hạn chế xói lở, ngăn mặn và điều hòa vi khí hậu. Trong số các loài cây ngập mặn, Đước đôi (Rhizophora apiculata Blume) giữ vị trí ưu thế tuyệt đối về trữ lượng và cấu trúc sinh thái tại vùng ven biển Nam Bộ. Đặc biệt, các nghiên cứu quốc tế đã khẳng định rừng ngập mặn có khả năng lưu giữ và cô lập carbon hữu cơ cao hơn nhiều lần so với phần lớn các kiểu rừng trên cạn.

Tuy nhiên, tại Việt Nam, phần lớn các công trình nghiên cứu trước đây về sinh khối và hấp thu carbon của rừng ngập mặn chỉ dừng lại ở các phương pháp điều tra thực nghiệm truyền thống (đo đếm tĩnh tại một thời điểm cụ thể). Cách tiếp cận này chỉ phản ánh trạng thái hiện tại, mang tính cục bộ cho từng địa điểm điều tra và không thể ngoại suy hay dự báo biến động sinh trưởng khi điều kiện khí hậu, môi trường sống hoặc chế độ lâm sinh thay đổi.

Trong bối cảnh biến đổi khí hậu diễn biến phức tạp và mục tiêu giảm phát thải ròng (Net Zero) đang trở thành cam kết chiến lược, việc áp dụng các mô hình động thái quá trình (Process-based Models) là một khoảng trống nghiên cứu cấp thiết. Phần mềm 3-PG (Physiological Principles in Predicting Growth) – một mô hình sinh thái sinh lý tích hợp các quá trình quang hợp, hô hấp, bức xạ mặt trời, độ hụt hơi nước và cân bằng dinh dưỡng đất – đã chứng minh hiệu quả cao trên thế giới đối với rừng trồng trên cạn, nhưng hầu như chưa từng được thử nghiệm trên hệ sinh thái rừng ngập mặn tại Việt Nam. Nghiên cứu này tại Kiên Giang mang tính tiên phong, giải quyết trực tiếp bài toán dự báo năng suất sinh khối và trữ lượng carbon lâu dài, đóng góp thiết thực cho công tác quy hoạch lâm nghiệp thích ứng biến đổi khí hậu.


Methodology và approach (350-400 từ)

Nghiên cứu áp dụng thiết kế thực nghiệm chuỗi thời gian (chronosequence) kết hợp giữa điều tra lâm học thực địa, phân tích giải tích phòng thí nghiệm và mô hình hóa sinh lý học 3-PG:

  • Thiết kế lấy mẫu thực địa: Nghiên cứu được triển khai tại xã Thuận Hòa (huyện An Minh) và xã Nam Thái (huyện An Biên) thuộc Ban quản lý rừng phòng hộ An Minh – An Biên, tỉnh Kiên Giang. Các lâm phần rừng Đước thuần loài được chọn ở 5 cấp tuổi đại diện: 10, 14, 15, 16 và 19 tuổi. Tại mỗi tuổi rừng, bố trí 3 ô tiêu chuẩn ngẫu nhiên diện tích $500\text{ m}^2$ ($20\text{ m} \times 25\text{ m}$).

  • Phương pháp thu thập và phân tích sinh khối:

    • Đo đếm lâm phần: Đo đường kính ngang ngực ($D_{1.3}$) bằng thước vanh và chiều cao vút ngọn ($H_{vn}$) từ gốc sát trên rễ chống đến đỉnh ngọn cao nhất.
    • Chặt hạ giải tích: Mỗi ô tiêu chuẩn chặt hạ 01 cây cá thể có kích thước trung bình về tiết diện ngang. Phân chia cây thành 5 hợp phần: thân, cành, lá, rễ chống trên mặt đất và hoa/quả.
    • Xác định khối lượng khô: Cân toàn bộ khối lượng tươi tại hiện trường; lấy mẫu lặp lại 3 lần cho từng bộ phận (50–100g/mẫu), sấy ở nhiệt độ $105^\circ\text{C}$ đối với thân, cành, rễ và $75^\circ\text{C}$ đối với lá đến khối lượng không đổi để tính tỷ lệ độ ẩm và sinh khối khô kiệt.
    • Sinh khối dưới mặt đất ($W_R$): Do tính chất phức tạp của hệ rễ ngập triều, sinh khối rễ ngầm được định lượng qua phương trình allometric tiêu chuẩn:

    $$W_R = 0{,}00698 \times DBH^{2{,}61}$$

  • Lượng hóa Carbon và $CO_2$: Khối lượng carbon tích lũy được tính bằng $50%$ sinh khối khô theo chuẩn quốc tế IPCC ($C = 0{,}5 \times W_{khô}$), từ đó quy đổi ra lượng $CO_2$ hấp thu:

    $$Q_{CO_2} = C \times \frac{44}{12}$$

  • Ứng dụng mô hình 3-PG: Thu thập chuỗi số liệu khí tượng nhiều năm (nhiệt độ trung bình/tối cao/tối thấp, lượng mưa, bức xạ mặt trời, độ bão hòa hơi nước, số ngày sương) từ trạm khí tượng địa phương; thiết lập các tham số sinh lý loài Đước đôi (chỉ số diện tích lá LAI, hiệu suất lượng tử tán rừng $\alpha_x$, tỷ lệ phân chia quang hợp vào rễ/thân/lá) và các đặc tính lập địa để mô phỏng động thái năng suất sơ cấp tổng số (GPP) và năng suất sơ cấp thực (NPP).


Phát hiện chính (400-450 từ)

Kết quả phân tích từ hệ thống ô tiêu chuẩn thực nghiệm và dữ liệu giải tích đã chỉ ra các quy luật quan trọng về động thái sinh khối và hấp thu carbon của rừng Đước Kiên Giang:

+-----------------------------------------------------------------------------------+
| Lâm phần      | Mật độ     | Sinh khối cá thể | Sinh khối quần thể | CO2 hấp thu   |
| (Tuổi rừng)   | (thân/ha)  | (kg khô/cây)     | (tấn khô/ha)       | (tấn/ha)      |
+---------------+------------+------------------+--------------------+---------------+
| 10 tuổi       | 5.047      | 25,94            | 136,12             | 249,56        |
| 14 tuổi       | 3.847      | 39,20            | 158,07             | 289,79        |
| 15 tuổi       | 4.327      | 51,10            | 228,98             | 419,79        |
| 16 tuổi       | 3.820      | 40,56            | 162,69             | 298,26        |
| 19 tuổi       | 3.127      | 70,17            | 229,45             | 420,66        |
+-----------------------------------------------------------------------------------+

1. Xu thế tăng trưởng sinh khối cá thể và quần thể

  • Sinh khối khô trung bình của một cây cá thể tăng trưởng tuyến tính theo tuổi, từ $25,94\text{ kg/cây}$ ở tuổi 10 lên đến $70,17\text{ kg/cây}$ ở tuổi 19 (tăng gấp gần 2,7 lần).
  • Tổng sinh khối quần thể (bao gồm cả trên và dưới mặt đất) tăng từ $136,12\text{ tấn/ha}$ (tuổi 10) lên $229,45\text{ tấn/ha}$ (tuổi 19).

2. Cấu trúc phân bố sinh khối giữa các bộ phận

  • Thân cây chiếm tỷ trọng lớn nhất trong toàn bộ sinh khối cây ở mọi lứa tuổi, dao động từ $62,14%$ đến $72,68%$, và tỷ lệ này có xu hướng tăng dần khi rừng già đi do quá trình khép tán và tỉa cành tự nhiên.
  • Tổ hợp thân và rễ chống trên mặt đất chiếm tỷ lệ áp đảo từ $75,11%$ đến $86,28%$, thể hiện đặc tính thích nghi cơ học vượt trội giúp cây đứng vững trên nền đất bùn lầy ven biển.
  • Cành chiếm $8,48% - 16,27%$, rễ chống chiếm $11,14% - 15,84%$, lá chiếm $5,08% - 10,84%$, trong khi hoa/quả chiếm tỷ lệ không đáng kể ($0% - 0,37%$).

3. Khả năng hấp thu $CO_2$ vượt bậc so với rừng trên cạn

  • Lượng $CO_2$ tích lũy của quần thể rừng Đước đạt mức cao: từ $249,56\text{ tấn } CO_2\text{/ha}$ (10 tuổi) lên $420,66\text{ tấn } CO_2\text{/ha}$ (19 tuổi).
  • Giai đoạn tăng tốc hấp thu mạnh nhất diễn ra từ tuổi 10 đến tuổi 15 (tăng thêm hơn $170\text{ tấn } CO_2\text{/ha}$), sau đó tốc độ tích lũy chậm lại ở các tuổi cao hơn. Khi so sánh với các loài cây trồng lâm nghiệp nội địa cùng nhóm tuổi như Keo lai ($60 - 407\text{ tấn/ha}$), Keo tai tượng ($57,63 - 281,40\text{ tấn/ha}$), rừng Đước ngập mặn thể hiện năng lực cố định carbon vượt trội.

4. Phát hiện bất thường về mối tương quan giữa tuổi 15 và 16

Nghiên cứu phát hiện sinh khối quần thể ở tuổi 16 ($162,69\text{ tấn/ha}$) thấp hơn đáng kể so với tuổi 15 ($228,98\text{ tấn/ha}$). Hiện tượng này không phản ánh sự suy thoái sinh trưởng sinh lý mà bắt nguồn từ sự khác biệt về mật độ lâm phần ($3.820\text{ thân/ha}$ ở tuổi 16 so với $4.327\text{ thân/ha}$ ở tuổi 15) và tác động của các biện pháp tỉa thưa, chặt cành cục bộ của các hộ dân nhận khoán bảo vệ rừng theo mô hình kinh tế Lâm - Ngư 7/3.


Đóng góp khoa học (250-300 từ)

  • Về mặt lý thuyết sinh thái: Công trình bổ sung hệ thống số liệu định lượng chi tiết về quy luật phân bổ sinh khối các bộ phận của loài Đước đôi (R. apiculata) theo cấp tuổi thành thục tại vùng sinh thái Tây Nam Bộ. Nghiên cứu làm sáng tỏ mối quan hệ giữa sinh trưởng sinh lý cây cá thể và mật độ quần thể trong việc kiến tạo trữ lượng bể chứa carbon ven biển.
  • Về mặt đổi mới phương pháp luận: Đây là một trong những nghiên cứu tiên phong tại Việt Nam thực hiện chuyển giao và hiệu chỉnh mô hình động thái sinh lý 3-PG cho hệ sinh thái rừng ngập mặn. Nghiên cứu đã xác lập quy trình xử lý các giả định sinh lý đặc thù (như ảnh hưởng không đáng kể của độ hụt nước khả dụng trong đất ngập triều - ASW) để mô hình 3-PG hoạt động chính xác trong môi trường thủy văn mặn.
  • Về mặt ứng dụng và chính sách: Cung cấp bộ công cụ dự báo định lượng sản lượng gỗ và sinh khối theo kịch bản thời gian, hỗ trợ đắc lực cho các cơ quan quản lý lâm nghiệp tỉnh Kiên Giang. Kết quả trữ lượng $CO_2$ là căn cứ khoa học pháp lý phục vụ việc xây dựng "Đường phát thải cơ sở" (Baseline) cho các dự án tín chỉ Carbon rừng ngập mặn, đóng góp vào cam kết Đóng góp do quốc gia tự quyết định (NDC) của Việt Nam.

Đối tượng quan tâm (200-250 từ)

  • Các nhà nghiên cứu học thuật & Nhà sinh thái học: Cung cấp dữ liệu thực nghiệm chuẩn hóa về cấu trúc sinh thái rừng ngập mặn, bộ tham số chuyển đổi sinh học cho các mô hình tính toán chu trình carbon toàn cầu và mô hình hóa sinh quyển.
  • Cán bộ quản lý lâm nghiệp & Ban quản lý rừng phòng hộ: Giúp Ban quản lý rừng phòng hộ An Minh – An Biên và Sở Nông nghiệp & PTNT Kiên Giang tối ưu hóa mật độ trồng, xác định luân kỳ tỉa thưa hợp lý và nâng cao hiệu quả kỹ thuật lâm sinh.
  • Các tổ chức phát triển dự án Carbon & Nhà hoạch định chính sách: Cung cấp số liệu kiểm chứng độc lập cho các dự án chi trả dịch vụ môi trường rừng (PES), thị trường tín chỉ carbon tự nguyện (VCM) và các cơ chế tài chính khí hậu quốc tế.
  • Cộng đồng cư dân nhận khoán đất rừng ven biển: Cung cấp định hướng cân đối bền vững giữa nuôi trồng thủy sản và bảo tồn độ che phủ rừng trong mô hình sinh kế kết hợp Lâm – Ngư 70/30.

FAQ - Các câu hỏi thường gặp (250-300 từ)

1. Phát hiện quan trọng nhất của nghiên cứu này là gì?

Khẳng định rừng Đước đôi tại Kiên Giang có tốc độ hấp thu carbon và tích lũy sinh khối cực kỳ cao (đạt $420,66\text{ tấn } CO_2\text{/ha}$ ở tuổi 19), vượt trội so với hầu hết các loại rừng trồng trên cạn, với thân cây là kho lưu giữ carbon chính yếu (chiếm trên $62-72%$).

2. Điểm đặc biệt trong phương pháp tiếp cận của nghiên cứu là gì?

Kết hợp đồng thời giữa phương pháp cân đo giải tích thực nghiệm phá hủy truyền thống với mô hình động thái quá trình 3-PG, cho phép liên kết dữ liệu sinh lý thực vật với yếu tố khí hậu để dự báo sản lượng dài hạn thay vì chỉ đo đếm tĩnh.

3. Tại sao sinh khối rừng 16 tuổi lại thấp hơn rừng 15 tuổi?

Do mật độ cây thực tế tại lâm phần 16 tuổi thấp hơn ($3.820\text{ thân/ha}$ so với $4.327\text{ thân/ha}$) cùng với sự khác biệt về lịch sử tác động tỉa thưa, chăm sóc của các hộ dân nhận khoán, chứng minh mật độ là yếu tố chi phối mạnh tới sinh khối quần thể.

4. Kết quả nghiên cứu có thể khái quát hóa (generalize) cho các vùng khác không?

Các quy luật về cấu trúc tỷ lệ bộ phận và bộ tham số 3-PG đã hiệu chỉnh có thể mở rộng áp dụng cho các vùng rừng ngập mặn Đước đôi tại Cà Mau, Bạc Liêu, Sóc Trăng hoặc TP.HCM (Cần Giờ) khi được hiệu chỉnh theo các biến khí tượng địa phương tương ứng.

5. Ứng dụng thực tiễn lớn nhất của đề tài là gì?

Làm cơ sở kỹ thuật xây dựng hồ sơ tín chỉ "Carbon xanh" thương mại, tính toán giá trị dịch vụ hấp thu $CO_2$ của hệ sinh thái rừng ngập mặn và hỗ trợ quy hoạch rừng thích ứng biến đổi khí hậu.


Kết luận (150 từ)

Nghiên cứu đã định lượng hóa một cách toàn diện động thái sinh khối và tiềm năng hấp thu carbon của rừng Đước đôi (Rhizophora apiculata) qua các giai đoạn tuổi từ 10 đến 19 tại Kiên Giang, đồng thời khẳng định tính khả thi của mô hình sinh lý 3-PG trong mô phỏng sản lượng rừng ngập mặn. Kết quả cho thấy rừng Đước ven biển là bể chứa carbon sinh học khổng lồ, đóng vai trò then chốt trong giảm thiểu tác động khí nhà kính.

Hướng nghiên cứu tiếp theo cần tập trung tích hợp mô hình không gian 3-PGS với dữ liệu viễn thám vi mô, đồng thời mở rộng mô phỏng dưới các kịch bản nước biển dâng và xâm nhập mặn cực đoan. Các cơ quan quản lý và doanh nghiệp cần nhanh chóng phối hợp ứng dụng bộ chỉ số này để thúc đẩy các dự án thị trường carbon rừng ngập mặn thực chất tại Đồng bằng sông Cửu Long.