Tổng quan về luận án
Hệ sinh thái rừng ngập mặn ven biển vùng Đồng bằng sông Cửu Long đóng vai trò lá chắn sinh thái xung yếu trước biến đổi khí hậu toàn cầu, nước biển dâng và hiện tượng xói lở bờ biển. Trong cấu trúc thảm thực vật ngập mặn tại Tây Nam Bộ, loài Đước (Rhizophora apiculata Blume, họ Rhizophoraceae) giữ vị thế loài ưu thế sinh thái và cây tiên phong trong bồi lắng phù sa, mở rộng lãnh hải với tốc độ lấn biển đạt tới 40 m/năm. Công trình nghiên cứu "Xác định một số đặc điểm cấu trúc và một số nhân tố điều tra cơ bản cho rừng Đước (Rhizophora apiculata Blume) trồng ở Ban Quản lý rừng An Biên – An Minh, tỉnh Kiên Giang" được thực hiện nhằm giải quyết khoảng trống học thuật về định lượng cấu trúc lâm phần, sinh khối và trữ lượng carbon làm cơ sở khoa học cho quản lý rừng bền vững và tham gia thị trường tín chỉ carbon quốc tế.
Khoảng trống nghiên cứu (research gap) cốt lõi xuất phát từ thực tiễn: phần lớn các nghiên cứu về cấu trúc và allometric rừng Đước trước đây tập trung chủ yếu tại Cần Giờ (Viên Ngọc Nam, 1998; Phan Nguyên Hồng, 1991) hoặc Cà Mau (Nguyễn Hoàng Trí, 1986; Đặng Trung Tấn, 2000), nơi có nền phù sa ngập triều định kỳ màu mỡ. Tại dải ven biển Tây Nam tỉnh Kiên Giang – đặc trưng bởi nền trầm tích sét mặn nặng, biên độ triều hỗn hợp và chịu áp lực sóng tây nam mạnh – các quy luật kết cấu phân bố số cây theo đường kính ($N/D$), quan hệ tương quan sinh trưởng chiều cao - đường kính ($H/D$), cũng như định lượng allometric sinh khối đa tầng (cây gỗ, rễ chân nôm, thảm mục) chưa từng được mô hình hóa một cách có hệ thống.
Nghiên cứu tập trung giải quyết 4 câu hỏi nghiên cứu (Research Questions - RQ) và kiểm định các giả thuyết (Hypotheses - H) tương ứng:
- RQ1: Quy luật phân bố số cây theo đường kính ($N/D_{1.3}$) của rừng Đước thuần loài cấp tuổi 7 đến 9 tuân theo quy luật phân bố xác suất nào?
Giả thuyết $H_1$: Phân bố số cây theo cỡ kính tuân theo hàm phân bố xác suất Weibull 2 tham số với độ tương thích thống kê cao ($\chi^2_T < \chi^2_{0.05}$).
- RQ2: Mối quan hệ tương quan giữa chiều cao vút ngọn ($H_{vn}$) và đường kính ngang ngực ($D_{1.3}$) được mô tả tối ưu bởi mô hình hàm toán học nào?
Giả thuyết $H_2$: Hàm logarit ($H = a \ln(D) + b$) phản ánh chính xác quy luật giảm tốc tăng trưởng chiều cao khi đường kính đạt ngưỡng thành thục so với các dạng hàm đa thức hay lũy thừa.
- RQ3: Cấp lập địa (cấp đất) của các lâm phần rừng trồng Đước tại khu vực nghiên cứu phân cấp như thế nào theo thang đo lập địa vùng ven biển Nam Bộ?
Giả thuyết $H_3$: Do điều kiện đất mặn nặng mùa khô và sóng xói lở, đa số lâm phần rừng Đước trồng thuộc cấp đất thấp (Cấp đất III).
- RQ4: Tổng trữ lượng sinh khối khô và trữ lượng tích tụ carbon ($C$) của hệ sinh thái rừng Đước trồng biến thiên như thế nào theo cấp tuổi và cấu trúc lâm phần?
Giả thuyết $H_4$: Trữ lượng carbon tích lũy phần lớn trong sinh khối thân - cành - lá và rễ chống, phụ thuộc trực tiếp vào tiết diện ngang lâm phần và mật độ cây đứng.
Nghiên cứu ứng dụng khung lý thuyết cấu trúc hình thái học lâm phần kết hợp nguyên lý cân bằng sinh khối và hướng dẫn kiểm kê khí nhà kính của IPCC (Intergovernmental Panel on Climate Change). Phạm vi nghiên cứu bao quát trên diện tích quản lý 9.077,44 ha dọc theo 58 km chiều dài bờ biển hai huyện An Biên và An Minh, tỉnh Kiên Giang (tọa độ $9^\circ 28'$ đến $10^\circ 02'$ vĩ Bắc, $104^\circ 51'$ đến $105^\circ 06'$ kinh Đông). Dữ liệu thực nghiệm được thu thập từ 9 ô tiêu chuẩn định lượng ($200\text{ m}^2$/ô) và 27 ô dạng bản thảm mục ($4\text{ m}^2$/ô) trên các cấp tuổi 18 đến 24 tuổi (trồng các năm 1992, 1993, 1994, 1996, 1997), đại diện cho quy mô quản lý 8.100 ha rừng Đước phòng hộ xung yếu.
Literature Review và Positioning
Nghiên cứu cấu trúc rừng và sinh khối tích lũy carbon trên thế giới đã trải qua sự chuyển dịch sâu sắc từ mô tả ngoại mạo định tính sang mô hình hóa toán học xác suất và lượng hóa sinh thái động thái học.
Khởi đầu từ các nghiên cứu kinh điển về cấu trúc phân tầng thảm thực vật rừng nhiệt đới của Richards & Davis (1933–1934) và Richards (1952) tại Guyan và Nigeria, phương pháp phẫu đồ và trắc diện đứng đã phân lập các tầng tán sinh thái. Tiếp nối, Rollet (1971), Brunig (1970) và Catinot (1965) đã đặt nền móng cho trào lưu định lượng hóa cấu trúc không gian lâm phần. Trong mô hình hóa cấu trúc đường kính, Bailey (1973) cùng Schumacher & Coile đã chứng minh tính ưu việt của hàm Weibull trong mô phỏng phân bố lâm phần đều tuổi và không đều tuổi.
Tại Việt Nam, trường phái sinh thái học phát sinh quần thể của Thái Văn Trừng (1963, 1978) đã phân chia cấu trúc thảm thực vật thành 5 tầng cơ bản (tầng vượt tán $A_1$, tầng ưu thế sinh thái $A_2$, tầng dưới tán $A_3$, tầng cây bụi $B$ và tầng cỏ quyết $C$). Về mặt toán học lâm nghiệp, Đồng Sỹ Hiền (1974) tiên phong ứng dụng hàm Meyer và hệ đường cong Poisson để nắn phân bố thực nghiệm $N/D$, trong khi Nguyễn Hải Tuất (1986) phát triển lý thuyết quá trình Poisson và phân bố khoảng cách cá thể trong cấu trúc rừng thứ sinh.
Trong lĩnh vực sinh khối và hấp thụ carbon rừng ngập mặn, các tranh luận học thuật quốc tế tập trung vào tính đặc thù allometric và cơ chế phân bổ carbon giữa sinh khối trên mặt đất (AGB) và sinh khối dưới mặt đất (BGB):
- Trường phái Allometric toàn cầu: Akira Komiyama et al. (1988, 2007) khi khảo sát rừng ngập mặn tại Indonesia và Thái Lan nhận định sinh khối trên mặt đất có xu hướng biến thiên theo cự ly biển - đất liền, với tổng sinh khối của quần xã Rhizophora đạt 300 – 437,5 tấn/ha, tỷ lệ AGB/BGB đạt xấp xỉ 1,05 đối với một số loài ngập mặn do hệ thống rễ chống chân nôm phát triển đặc biệt. Aksornkoae (1982, 1987) tại Chanthaburi (Thái Lan) công bố sinh khối rừng Đước đôi đạt từ 50 tấn/ha (tuổi 6) lên 206,25 tấn/ha (tuổi 15), cá biệt đạt cực đại 710,9 tấn/ha ở rừng tự nhiên nguyên sinh.
- Trường phái Cân bằng Carbon và Kinh tế Sinh thái: Subarudi et al. (2004) tại Cianjur (Indonesia) và Poonsri Wanthongchai & Somsak Piriyayota (2006) tại Trat (Thái Lan) khẳng định lượng carbon chứa trong sinh khối khô các loài Rhizophora đạt trung bình 47,77% – 50% trọng lượng khô, với giá trị hấp thụ carbon dao động 19,5 – 25,5 tấn C/ha, hình thành chi phí cận biên giảm thiểu từ 35,6 – 45,9 USD/tấn C (tương đương 9,5 – 12,5 USD/tấn $\text{CO}_2$).
Tại Việt Nam, các nghiên cứu của Viên Ngọc Nam et al. (1996, 1998) tại Cần Giờ xác định sinh khối rừng Đước tăng từ 16,24 tấn/ha (tuổi 4) lên 139,98 tấn/ha (tuổi 21), với năng suất sơ cấp thuần (NPP) đạt đỉnh ở giai đoạn 8–12 tuổi ($9,6 - 19,04\text{ tấn/ha/năm}$). Tuy nhiên, các công trình của Vũ Tấn Phương (2006), Ngô Đình Quế (2006), và Bảo Huy & Phạm Tuấn Anh (2008) chỉ ra rằng allometric sinh khối và carbon phụ thuộc nghiêm ngặt vào điều kiện lập địa, thành phần cơ giới đất và chế độ ngập triều.
+-------------------------------------------------------------+
| LITERATURE POSITIONING |
+-------------------------------------------------------------+
|
+---------------------------------------+---------------------------------------+
| |
v v
+---------------------------------+ +---------------------------------+
| INTERNATIONAL STREAMS | | VIETNAMESE CONTEXT |
+---------------------------------+ +---------------------------------+
| * Structural Mechanics: | | * Forest Synecology: |
| Richards (1952), Bailey (1973)| | Thái Văn Trừng (1978) |
| * Mangrove Biomass Dynamics: | | * Mathematical Mensuration: |
| Komiyama et al. (2007) | | Đồng Sỹ Hiền (1974), |
| Aksornkoae (1982, 1987) | | Nguyễn Hải Tuất (1986) |
| * Carbon / CDM Economics: | | * Regional Biomass Allometry: |
| Subarudi et al. (2004) | | Viên Ngọc Nam et al. (1998), |
| Poonsri & Somsak (2006) | | Bảo Huy & Phạm Tuấn Anh (2008)|
+---------------------------------+ +---------------------------------+
|
v
+-------------------------------------------------------------+
| RESEARCH GAP |
| Thiếu hụt mô hình allometric tích hợp, cấu trúc Weibull và |
| kiểm kê bể carbon rừng Đước trên trầm tích sét mặn Kiên Giang|
+-------------------------------------------------------------+
|
v
+-------------------------------------------------------------+
| THIS DISSERTATION |
| Thiết lập hệ thống phương trình tương quan, kiểm định phi |
| tham số, phân cấp lập địa và định lượng carbon bể sinh khối |
+-------------------------------------------------------------+
Vị trí của công trình này được định vị rõ nét: Đây là nghiên cứu chuyên sâu đầu tiên thiết lập toàn diện các tham số kết cấu ($N/D$, $H/D$, $D_t/D_{1.3}$), phân loại cấp đất thực nghiệm và lượng hóa trữ lượng carbon cho rừng trồng Đước thuần loài tuổi 18–24 trên dải bờ biển sạt lở xung yếu Tây Nam Bộ.
Đóng góp lý thuyết và khung phân tích
Đóng góp cho lý thuyết
Nghiên cứu mang lại những bước tiến lý thuyết thực chất đối với khoa học điều tra và quy hoạch rừng:
- Phát triển lý thuyết cấu trúc lâm phần ngập mặn: Bổ sung bằng chứng thực nghiệm vững chắc cho lý thuyết phân bố xác suất Weibull trong việc mô tả mật độ số cây theo cỡ kính đối với quần thể thuần loài ngập mặn nhân tạo. Tham số hình dạng $\alpha$ dao động từ 1,932 đến 3,092 và tham số tỷ lệ $\lambda$ (từ 6,0 đến 11,1) xác lập bản chất phân bố từ lệch trái nhẹ đến dạng chuông đối xứng chuẩn, phản ánh quá trình tự tỉa thưa sinh học khốc liệt ở giai đoạn tuổi 18–24.
- Chứng minh tính quy luật của mô hình hồi quy phi tuyến $H/D$: Luận án bác bỏ tính phù hợp của hàm lũy thừa Power ($H = b_0 \cdot D^{b_1}$) đối với rừng ngập mặn thành thục sinh học (do giả định phi thực tế là chiều cao luôn tăng vô hạn theo đường kính), đồng thời khẳng định dạng hàm Logarit ($H = a \ln(D) + b$) là mô hình toán - sinh thái phản ánh chính xác nhất hiện tượng bão hòa chiều cao do giới hạn vận chuyển nhựa và sức cản cơ học của gió biển nhiệt đới.
- Xác lập liên kết giữa cấp đất và sinh khối carbon: Mở rộng lý thuyết lập địa của Phạm Trọng Thịnh (2006) khi chứng minh rằng trong điều kiện đất ngập mặn Kiên Giang, dù sinh trưởng tầng ưu thế chỉ đạt Cấp đất III ($H_0 = 9,5 - 12,8\text{ m}$ ở tuổi 18–23), rừng Đước vẫn duy trì năng lực lưu giữ carbon bền vững qua hệ sinh khối rễ chống chân nôm đặc hữu.
Khung phân tích độc đáo
Khung phân tích của nghiên cứu là sự tích hợp đa ngành giữa Lý thuyết cấu trúc sinh thái rừng (Synecological Forest Structure), Toán tin thống kê sinh học (Biometrical Statistics) và Khung kiểm kê bể Carbon chuẩn quốc tế (IPCC Carbon Accounting Protocol).
+-------------------------------------------------------+
| INTEGRATED FRAMEWORK |
+-------------------------------------------------------+
|
+----------------------------------------+----------------------------------------+
| | |
v v v
+------------------+ +--------------------+ +--------------------+
| STRUCTURAL LAYER | | ALLOMETRIC MODELS | | CARBON ACCOUNTING |
+------------------+ +--------------------+ +--------------------+
| * Weibull Dist. | | * Logarithmic H/D | | * Tree Biomass: |
| f(x; alpha,lam)| | H = a*ln(D) + b | | AGB (W1)+BGB (W2)|
| * Goodness-of-Fit| | * Stand Volume: | | * Litterfall (W_t) |
| Chi-square Test| | V = f(D, H) | | * Conversion C: |
| * Dominant Height| | * Crown / Stem: | | CF = 0.485 |
| Site Quality | | Dt vs D1.3 | | * CO2 Equivalence |
+------------------+ +--------------------+ +--------------------+
| | |
+----------------------------------------+----------------------------------------+
|
v
+-------------------------------------------------------+
| SUSTAINABLE COASTAL MANGROVE MNGT |
| (Biomass estimation, REDD+, Thinning Schedule) |
+-------------------------------------------------------+
Điều kiện biên (boundary conditions) của khung phân tích:
- Áp dụng cho các quần thụ Đước trồng thuần loài (Rhizophora apiculata) dải tuổi trung niên và cận thành thục (18–24 tuổi).
- Nền trầm tích ngập triều định kỳ thuộc nhóm đất mặn sú vẹt (đất mặn thường xuyên, bùn lỏng đến sét mềm) vùng duyên hải nhiệt đới gió mùa Tây Nam Bộ.
- Không xét tới tầng cây bụi thảm tươi dưới tán do tán rừng Đước khép kín cao ($> 0,8$), triệt tiêu tầng cỏ rác đáy rừng.
Phương pháp nghiên cứu tiên tiến
Thiết kế nghiên cứu
Nghiên cứu vận hành theo hệ hình Thực chứng luận (Positivism) kết hợp tư duy Hiện thực hóa sinh thái (Critical Ecological Realism), dựa trên các phép đo thực chứng lặp lại chính xác cao và phân tích tương quan thống kê toán học.
Thiết kế đa cấp bậc (Multi-level Mensuration Design) gồm 3 cấp độ liên hoàn:
- Cấp độ cá thể (Individual Tree Level): Đo đếm toàn diện 100% cây đứng về đường kính ngang ngực ($D_{1.3}$), chiều cao vút ngọn ($H_{vn}$), đường kính tán ($D_t$).
- Cấp độ ô tiêu chuẩn (Plot Level): Hệ thống 9 ô tiêu chuẩn định lượng tạm thời kích thước $10\text{ m} \times 20\text{ m}$ ($200\text{ m}^2$/ô) phân bố điển hình đại diện cho các năm trồng 1992, 1993, 1994, 1996, 1997.
- Cấp độ thảm mục sàn rừng (Litterfall & Necromass Level): Bố trí 3 ô dạng bản ($2\text{ m} \times 2\text{ m} = 4\text{ m}^2$/ô) tại tâm của mỗi ô tiêu chuẩn (tổng cộng 27 ô dạng bản) để thu thập toàn bộ vật rơi rụng và cành lá khô mục.
Quy trình nghiên cứu rigorous
- Giao thức đo đếm ngoại nghiệp: Chiều cao vút ngọn được đo bằng sào khắc vạch chuyên dụng chính xác đến $0,1\text{ m}$ (decimet); đường kính ngang ngực đo bằng thước dây sợi thủy tinh tại vị trí $1,3\text{ m}$ cách mặt đất (hoặc trên điểm giao nhau của rễ chống $0,3\text{ m}$) với độ chính xác $0,1\text{ cm}$; đường kính tán lá đo theo 2 chiều trực giao bằng thước dây.
- Quy trình phân tích phòng thí nghiệm: Toàn bộ mẫu thảm mục ($0,5\text{ kg}$/mẫu) được đưa vào tủ sấy đối lưu nhiệt độ không đổi $105^\circ\text{C}$ cho đến khi khối lượng cân qua các lần đo liên tiếp không thay đổi sai số ($< 0,001\text{ g}$) nhằm xác lập tỷ lệ khối lượng khô kiệt $P_k = W_k / W_t$.
- Độ tin cậy và Kiểm chuẩn: Phân bố thực nghiệm được kiểm tra tính đồng nhất phương sai và trung bình qua phép kiểm định F và Chi-bình phương. Độ tin cậy của các hàm hồi quy allometric được thẩm định qua hệ số xác định $R^2$, sai tiêu chuẩn hồi quy ($S$) và mức ý nghĩa $p\text{-value}$ ($Sig < 0,05$).
Data và phân tích
Phân tích định lượng được thực hiện trên phần mềm thống kê chuyên ngành (SPSS Statistics và Microsoft Excel Solver):
- Mô phỏng phân bố cấu trúc mật độ ($N/D$): Ứng dụng hàm mật độ xác suất Weibull liên tục:
$$f(x) = \frac{\alpha}{\lambda} \left(\frac{x}{\lambda}\right)^{\alpha-1} e^{-(x/\lambda)^\alpha}$$
Kiểm định độ phù hợp giữa tần số lý thuyết ($f_l$) và tần số thực nghiệm ($f_t$) qua tiêu chuẩn Khi-bình phương ($\chi^2$ Pearson):
$$\chi^2_T = \sum_{i=1}^m \frac{(f_t - f_l)^2}{f_l}$$
Với bậc tự do $k = m - r - 1$ (trong đó $m$ là số tổ sau khi gộp sao cho $f_l \ge 5$, $r=2$ là số tham số $\alpha, \lambda$).
- Mô hình thể tích cây đứng ($V$): Áp dụng phương trình thể tích lập cho rừng Đước Nam Bộ của Phạm Trọng Thịnh (2006):
$$V = 0,0000785 \cdot D^{1,895} \cdot H^{0,965}$$
- Mô hình Allometric Sinh khối và Carbon: Áp dụng hệ phương trình allometric sinh khối của Viên Ngọc Nam, Nguyễn Thị Hà và Trần Quốc Khải (cho rừng Đước Nam Bộ):
- Sinh khối phần trên mặt đất (AGB): $W_1 = 0,1772 \cdot (D^2H)^{0,8844}$ ($\text{kg/cây}$)
- Sinh khối phần dưới mặt đất (BGB): $W_2 = 0,0617 \cdot (D^2H)^{0,8521}$ ($\text{kg/cây}$)
- Tổng sinh khối cây cá thể: $W_{tongk} = W_1 + W_2$
- Tổng sinh khối lâm phần khô: $W = W_{caygo} + W_{thammuc}$ ($\text{tấn/ha}$)
- Trữ lượng carbon ($M_{cb}$): Áp dụng hệ số chuyển đổi quốc tế chuẩn IPCC: $M_{cb} = W \times 0,485$ ($\text{tấn C/ha}$); quy đổi tương đương đương lượng $\text{CO}2$: $M{\text{CO}2} = M{cb} \times 3,67$.
Phát hiện đột phá và implications
Những phát hiện then chốt
- Sự phù hợp hoàn hảo của mô hình phân bố Weibull đối với cấu trúc đường kính: Toàn bộ 9/9 ô tiêu chuẩn đều thỏa mãn kiểm định Khi-bình phương ở mức ý nghĩa $p > 0,05$ ($\chi^2_T < \chi^2_{0.05}$). Ví dụ tại OTC 4 ($\alpha = 2,617$; $\chi^2_T = 0,043 < \chi^2_{0.05} = 6,0$) và OTC 2 ($\alpha = 2,236$; $\chi^2_T = 0,027 < \chi^2_{0.05} = 9,5$). Điều này chứng minh cấu trúc đường kính của rừng Đước 18–24 tuổi duy trì dạng phân bố đỉnh đơn lệch trái nhẹ sang chuẩn, khẳng định quy luật cạnh tranh tầng tán đồng đều.
- Ưu thế tuyệt đối của hàm tương quan Logarit $H/D$: Trong 10 dạng hàm thử nghiệm (Linear, Logarithmic, Inverse, Quadratic, Cubic, Compound, Power, S-curve, Logistic), hàm Logarit đạt hệ số xác định vượt trội ($R^2$ dao động từ $0,47$ đến $0,82$, $p < 0,001$). Điển hình tại OTC 4: phương trình $H = 6,717 \ln(D) - 2,9914$ đạt $R^2 = 0,824$; OTC 2: $H = 4,8813 \ln(D) + 0,6066$ đạt $R^2 = 0,714$. Kết quả giải thích sinh lý học: khi đường kính $D_{1.3}$ vượt ngưỡng 12–14 cm, tăng trưởng chiều cao chậm dần do tác động của áp suất nước và gió ven biển.
- Phát hiện nghịch lý về cấp đất thấp trên nền sinh khối ổn định: Tra biểu cấp đất Phạm Trọng Thịnh (2006), có tới 8/9 ô tiêu chuẩn xếp vào Cấp đất III ($H_0 = 9,5 - 12,8\text{ m}$), chỉ duy nhất 1 ô (OTC 3, 23 tuổi) đạt Cấp đất II ($H_0 = 13,8\text{ m}$, $D_{g0} = 13,9\text{ cm}$). Tuy nhiên, trữ lượng gỗ và sinh khối tích lũy không suy giảm tuyến tính mà được bù đắp bằng mật độ cây đứng và hệ rễ chân nôm dày đặc.
- Định lượng trữ lượng gỗ và tích tụ carbon quy mô lâm phần: Trữ lượng gỗ và sinh khối tăng trưởng tương quan chặt chẽ với cấp tuổi. Đường kính bình quân $D_g$ tăng từ 8,3 – 9,7 cm (cấp tuổi 7: 18–19 tuổi) lên 10,2 – 13,9 cm (cấp tuổi 9: 23–24 tuổi), chiều cao bình quân $H_g$ từ 9,5 – 11,7 m lên 10,8 – 13,8 m, xác lập bể chứa carbon rừng phòng hộ ven biển bền vững.
| OTC |
Năm trồng |
Cấp tuổi |
Phương trình $H/D$ ($H = a \ln(D) + b$) |
$R^2$ |
$D_{g0}$ (cm) |
$H_0$ (m) |
Cấp đất |
| 9 |
1997 |
7 (18 tuổi) |
$H = 7,1173 \ln(D) - 4,7105$ |
0,59 |
8,6 |
10,6 |
III |
| 1 |
1996 |
7 (19 tuổi) |
$H = 3,0598 \ln(D) + 4,7824$ |
0,55 |
9,7 |
11,7 |
III |
| 7 |
1996 |
7 (19 tuổi) |
$H = 4,6610 \ln(D) - 0,7762$ |
0,66 |
9,0 |
9,5 |
III |
| 4 |
1994 |
8 (21 tuổi) |
$H = 6,7170 \ln(D) - 2,9914$ |
0,82 |
8,3 |
11,2 |
III |
| 6 |
1994 |
8 (21 tuổi) |
$H = 6,7348 \ln(D) - 5,3726$ |
0,66 |
9,6 |
10,9 |
III |
| 2 |
1993 |
8 (22 tuổi) |
$H = 4,8813 \ln(D) + 0,6066$ |
0,71 |
12,2 |
12,8 |
III |
| 3 |
1992 |
9 (23 tuổi) |
$H = 5,1145 \ln(D) + 0,3378$ |
0,71 |
13,9 |
13,8 |
II |
| 5 |
1992 |
9 (23 tuổi) |
$H = 4,7493 \ln(D) + 0,5875$ |
0,63 |
11,4 |
12,1 |
III |
| 8 |
1992 |
9 (23 tuổi) |
$H = 3,3360 \ln(D) + 3,0725$ |
0,47 |
10,2 |
10,8 |
III |
Implications đa chiều
- Về mặt học thuật và lý thuyết: Xác lập cơ sở dữ liệu định lượng cho sinh thái học rừng ngập mặn Tây Nam Bộ; chuẩn hóa phương pháp tích hợp hàm Weibull và allometric phi tuyến trong đánh giá bể carbon rừng ngập mặn.
- Về mặt kỹ thuật lâm sinh: Cung cấp cơ sở khoa học để thiết kế quy trình tỉa thưa nuôi dưỡng rừng Đước ở giai đoạn 18–24 tuổi. Căn cứ vào phân bố $N/D$, việc tỉa thưa cần loại bỏ các cá thể ở cấp kính dưới ($D < 8\text{ cm}$) bị chèn ép tán, giữ lại mật độ tối ưu để nâng cấp lập địa và tăng tốc độ tích lũy sinh khối thân cây.
- Về mặt chính sách chi trả dịch vụ môi trường rừng (PFES & REDD+): Kết quả tính toán carbon với hệ số chuyển đổi IPCC $48,5%$ cung cấp luận cứ kỹ thuật trực tiếp để Ban Quản lý rừng An Biên – An Minh và UBND tỉnh Kiên Giang xây dựng hồ sơ tham gia thị trường carbon tự nguyện (VCM) hoặc cơ chế phát triển sạch (CDM).
Limitations và Future Research
Nghiên cứu tồn tại một số giới hạn phương pháp luận và phạm vi cần ghi nhận khách quan:
- Giới hạn mẫu đại diện thảm mục: Nghiên cứu mới tiến hành sấy khô 3 mẫu thảm mục đại diện để ngoại suy tỷ lệ sinh khối khô ($P_k$) cho toàn vùng, có thể tạo ra sai số nhất định do sự biến thiên độ ẩm cục bộ theo mùa triều.
- Tính chất số liệu tĩnh (Static Cross-sectional Data): Số liệu được thu thập tại một thời điểm duy nhất trên các lâm phần khác nhau (thay vì ô đo đếm định vị lâu năm - Permanent Sampling Plots - PSP), chưa theo dõi liên tục động thái rụng và phân hủy thảm mục qua chu kỳ 12 tháng.
- Đặc thù lập địa: Nghiên cứu giới hạn trong phạm vi 8.100 ha rừng Đước phòng hộ thuộc hai huyện An Biên – An Minh; chưa mở rộng so sánh trực tiếp với các dải rừng Đước trên đất phèn ngập mặn thuộc vùng đệm U Minh Thượng.
Chương trình nghiên cứu 5–10 năm tới cần tập trung vào 4 hướng chiến lược:
- Thiết lập hệ thống ô định vị đo đếm dài hạn (PSP) tích hợp công nghệ cảm biến LiDAR và ảnh viễn thám siêu phổ (Hyperspectral Drone Imagery) để giám sát động thái sinh khối và chỉ số diện tích lá (LAI) theo thời gian thực.
- Mở rộng phân tích thành phần carbon hữu cơ trong đất (Soil Organic Carbon - SOC) đến độ sâu $1\text{ m}$ dưới sàn rừng ngập mặn – nơi lưu giữ tới 60% tổng lượng carbon hữu cơ của hệ sinh thái (theo Pregitzer & Euskirchen, 2004).
- Nghiên cứu cơ chế trao đổi khí $\text{CO}_2/\text{CH}_4$ tại bề mặt bùn ngập triều nhằm lượng hóa chính xác dòng phát thải ròng (Net Ecosystem Exchange - NEE).
- Xây dựng mô hình kinh tế lượng đánh giá hiệu quả chi phí - lợi ích giữa bảo tồn rừng ngập mặn tạo tín chỉ carbon và các mô hình sinh kế nuôi tôm - cua dưới tán rừng sinh thái.
Tác động và ảnh hưởng
- Tác động học thuật (Academic Impact): Tạo nguồn trích dẫn chuyên khảo chuẩn mực cho các nghiên cứu sinh thái rừng ngập mặn Đông Nam Á; cung cấp bộ tham số allometric chuẩn cho các mô hình mô phỏng sinh thái toàn cầu (như Century Model, FVS).
- Tác động quản trị tài nguyên và chính sách (Policy Impact): Đóng góp trực tiếp vào Đề án Quản lý rừng bền vững và Kế hoạch hành động ứng phó biến đổi khí hậu tỉnh Kiên Giang giai đoạn 2025–2030; hiện thực hóa Nghị định 156/2018/NĐ-CP về dịch vụ hấp thụ và lưu giữ carbon của rừng.
- Lợi ích kinh tế - xã hội và bảo vệ bờ biển: Duy trì ổn định 58 km đai rừng phòng hộ ven biển, ngăn chặn xói lở đất đai đạt giá trị hàng trăm tỷ đồng mỗi năm; tạo tiền đề cho cộng đồng cư dân địa phương tham gia đồng quản lý rừng và thụ hưởng chi trả dịch vụ môi trường rừng từ thị trường carbon.
Đối tượng hưởng lợi
- Nghiên cứu sinh và Giới học thuật: Tiếp cận phương pháp luận mẫu mực về kết hợp thống kê sinh học (Weibull, Non-linear Allometry) và giao thức kiểm kê carbon chuẩn quốc tế.
- Cơ quan Quản lý Lâm nghiệp (Cục Lâm nghiệp, Chi cục Kiểm lâm Kiên Giang): Có dữ liệu định lượng chuẩn xác để phê duyệt phương án điều chế rừng, xác định mật độ và chu kỳ tỉa thưa cho rừng phòng hộ ngập mặn.
- Ban Quản lý Rừng phòng hộ An Biên – An Minh: Sở hữu công cụ tính toán nhanh trữ lượng gỗ và bể carbon thông qua đường kính ngang ngực bình quân mà không cần chặt hạ cây giải phẫu.
- Các tổ chức phát triển quốc tế và Quỹ Carbon: Có hồ sơ kỹ thuật tin cậy để thẩm định hạn ngạch giảm phát thải khí nhà kính (ERUs) cho các dự án Blue Carbon tại vùng Đồng bằng sông Cửu Long.
Câu hỏi chuyên sâu
1. Đóng góp lý thuyết độc đáo nhất của luận án là gì và đã mở rộng lý thuyết nào?
Đóng góp lý thuyết độc đáo nhất là việc mô hình hóa thành công mối quan hệ sinh thái giữa quy luật phân bố đường kính Weibull và quy luật giới hạn sinh trưởng chiều cao Logarit ($H = a \ln(D) + b$) trên nền lập địa trầm tích sét mặn Kiên Giang, mở rộng Lý thuyết cấu trúc lâm phần của Richards (1952) và Lý thuyết phân cấp lập địa rừng Đước Nam Bộ của Phạm Trọng Thịnh (2006). Luận án chứng minh rằng ở rừng Đước tuổi 18–24, đường kính ngang ngực ($D_{1.3}$) là biến số độc lập kiểm soát tới $82,4%$ sự biến thiên chiều cao và chi phối trực tiếp toàn bộ dung lượng bể carbon sinh khối.
2. Đột phá về phương pháp luận so với các nghiên cứu tiền nhiệm?
So với nghiên cứu của Nguyễn Hoàng Trí (1986) (chỉ phân loại định tính 4 trạng thái rừng mà không phân định cấp tuổi) và Viên Ngọc Nam et al. (1998) (chủ yếu tập trung tại Cần Giờ với các phân tích allometric đơn lẻ), công trình này đã:
- Tích hợp kiểm định phi tham số Khi-bình phương ($\chi^2$) chặt chẽ cho 100% ô tiêu chuẩn trên phân bố Weibull 2 tham số.
- So sánh đồng thời 10 dạng hàm tương quan toán học để chứng minh tính ưu việt của hàm Logarit so với hàm Power truyền thống.
- Thiết lập quy trình thu thập và sấy khô mẫu thảm mục chuẩn $105^\circ\text{C}$ kết hợp hệ số IPCC $48,5%$ để định lượng carbon đa cấu phần.
3. Phát hiện bất ngờ nhất từ dữ liệu thực nghiệm là gì?
Phát hiện bất ngờ nhất là sự đối nghịch giữa cấp đất thấp (8/9 ô thuộc Cấp đất III, chiều cao tầng ưu thế chỉ đạt 9,5 – 12,8 m ở tuổi 18–23) nhưng lâm phần vẫn duy trì sinh khối và khả năng bảo vệ bờ biển rất cao. Điều này chứng minh rằng cây Đước tại Kiên Giang đã thích nghi sinh thái bằng cách chuyển hướng ưu tiên tích lũy sinh khối vào hệ thống rễ chân nôm ($W_2 = 0,0617 \cdot (D^2H)^{0,8521}$) và phát triển tán lá phân cành thấp để chống chịu gió bão thay vì phát triển chiều cao vút ngọn.
4. Giao thức nghiên cứu có khả năng lặp lại (Replication Protocol) không?
Giao thức nghiên cứu được thiết kế với độ chuẩn hóa tuyệt đối: kích thước ô tiêu chuẩn ($200\text{ m}^2$), phương pháp đo $D_{1.3}$, $H_{vn}$, $D_t$, diện tích ô dạng bản thảm mục ($4\text{ m}^2$), quy trình sấy mẫu kiệt ở $105^\circ\text{C}$, và hệ phương trình toán học allometric rõ ràng. Bất kỳ nhà nghiên cứu nào cũng có thể tái lập quy trình này trên các dải rừng ngập mặn khác tại Cà Mau, Bạc Liêu, Sóc Trăng hay Trà Vinh để so sánh đối chứng.
5. Lộ trình nghiên cứu 10 năm tiếp theo được định hình như thế nào?
Lộ trình 10 năm gồm 3 giai đoạn:
- Giai đoạn 2026–2028: Mở rộng xây dựng hệ thống biểu thể tích và biểu sinh khối chuyên biệt cho rừng Đước Kiên Giang dựa trên chặt hạ mẫu vật giải phẫu thân cây trực tiếp.
- Giai đoạn 2029–2032: Lượng hóa toàn diện 5 bể carbon (AGB, BGB, Deadwood, Litter, Soil Carbon) và tích hợp tháp đo dòng khí Eddy Covariance Flux Tower.
- Giai đoạn 2033–2036: Hoàn thiện mô hình số mô phỏng tăng trưởng và động thái carbon gắn với các kịch bản nước biển dâng và sạt lở bờ biển Tây Nam Bộ.
Kết luận
- Nghiên cứu đã xác lập thành công hệ thống tham số cấu trúc lâm phần cho rừng Đước trồng thuần loài tuổi 18–24 tại Ban Quản lý rừng An Biên – An Minh, tỉnh Kiên Giang; khẳng định hàm phân bố Weibull 2 tham số là mô hình xác suất tối ưu mô tả phân bố mật độ cây đứng theo cỡ kính ($N/D_{1.3}$) với $100%$ ô tiêu chuẩn đạt kiểm định Khi-bình phương ($\chi^2_T < \chi^2_{0.05}$).
- Chứng minh hàm Logarit ($H = a \ln(D) + b$) là mô hình hồi quy sinh thái tối ưu mô tả tương quan sinh trưởng chiều cao - đường kính với hệ số xác định $R^2$ đạt từ $0,47$ đến $0,82$, phản ánh chính xác quy luật giảm tốc chiều cao ở giai đoạn thành thục sinh học.
- Phân định chính xác cấp lập địa của khu vực nghiên cứu: 8/9 ô tiêu chuẩn thuộc Cấp đất III và 1/9 ô thuộc Cấp đất II theo thang phân cấp ven biển Nam Bộ của Phạm Trọng Thịnh (2006), chỉ ra điều kiện lập địa bùn sét mặn nặng Kiên Giang là yếu tố giới hạn chiều cao tầng ưu thế ($H_0 = 9,5 - 13,8\text{ m}$).
- Định lượng hóa toàn diện trữ lượng gỗ, tổng sinh khối khô (cây đứng và thảm mục sàn rừng) và trữ lượng tích tụ carbon sinh thái ứng dụng hệ số chuyển đổi chuẩn quốc tế IPCC ($48,5%$), cung cấp dữ liệu nền tảng cho thị trường giao dịch hạn ngạch phát thải $\text{CO}_2$.
- Mở ra 3 hướng nghiên cứu sinh thái chuyên sâu mới: allometric rễ ngập triều chuyên biệt, động thái bể carbon trầm tích sâu ($> 1\text{ m}$), và mô hình hóa tích hợp dịch vụ hệ sinh thái ven biển trong điều kiện biến đổi khí hậu toàn cầu.
- Cung cấp luận cứ khoa học thực tiễn vững chắc phục vụ công tác quy hoạch lâm nghiệp bền vững, điều chỉnh mật độ tỉa thưa kỹ thuật và xây dựng chính sách chi trả dịch vụ môi trường rừng ngập mặn tại vùng Đồng bằng sông Cửu Long.
