Tổng quan về luận án

Hệ sinh thái rừng Tràm (Melaleuca cajuputi Powell) trên đất than bùn tại Vườn Quốc gia (VQG) U Minh Hạ, tỉnh Cà Mau đóng vai trò là bể chứa cacbon sinh thái đặc biệt quan trọng của vùng Đồng bằng sông Cửu Long (ĐBSCL) và khu vực hạ lưu lưu vực sông Mê Kông. Tuy nhiên, sau các biến cố cháy rừng lịch sử (tiêu biểu là vụ cháy năm 2002 làm thiêu rụi hơn 4.000 ha rừng tại U Minh Hạ), công tác quản lý tài nguyên rừng đã chuyển dịch sang chiến lược đắp đê bao khép kín để tích nước phòng cháy chữa cháy rừng (PCCCR) quanh năm. Chiến lược giữ nước nhân tạo này vô tình tạo ra một "nghịch lý sinh thái": việc duy trì mực nước ngập sâu liên tục làm đảo lộn chế độ thủy văn tự nhiên, đẩy hệ rễ cây Tràm vào trạng thái yếm khí kéo dài, gây thoái hóa thảm thực vật và làm suy giảm sinh khối lâm phần. Luận án tiến sĩ ngành Môi trường Đất và Nước (Mã số: 62 44 03 03) của nghiên cứu sinh Trần Thị Kim Hồng, dưới sự hướng dẫn khoa học của TS. Dương Văn Ni và TS. Nguyễn Văn Bé tại Trường Đại học Cần Thơ (năm 2017), đã giải quyết trực diện khoảng trống nghiên cứu (research gap) cốt lõi: Làm thế nào để định lượng chính xác tác động đồng thời của độ dày tầng than bùn, chế độ ngập nước và chất lượng môi trường lý hóa đất – nước lên động thái sinh khối và khả năng hấp thụ $CO_2$ của quần thể rừng Tràm, từ đó kiến tạo cơ sở khoa học cho mô hình quản lý thủy văn dung hòa giữa bảo tồn sinh thái và PCCCR.

Hệ thống câu hỏi nghiên cứu và giả thuyết khoa học của luận án được xác lập chặt chẽ:

  • Câu hỏi nghiên cứu 1 (RQ1): Các chỉ tiêu hóa lý của đất than bùn biến động như thế nào theo các cấp độ dày tầng than bùn và chu kỳ mùa trong năm?
  • Câu hỏi nghiên cứu 2 (RQ2): Chế độ ngập nước và các thông số chất lượng nước trong rừng Tràm thay đổi ra sao giữa các mùa và các mức độ ngập?
  • Câu hỏi nghiên cứu 3 (RQ3): Mức độ tích lũy sinh khối cá thể và sinh khối quần thụ rừng Tràm chịu sự chi phối như thế nào từ độ dày than bùn và độ sâu ngập nước?
  • Câu hỏi nghiên cứu 4 (RQ4): Trữ lượng hấp thụ $CO_2$ của rừng Tràm biến thiên ra sao dưới các điều kiện vi địa hình than bùn và ngập lũ khác nhau?
  • Giả thuyết 1 (H1): Tầng than bùn càng dày không đồng nghĩa với việc sinh khối lâm phần càng cao, do tính chất cơ lý xốp làm giảm mật độ cây đứng.
  • Giả thuyết 2 (H2): Chế độ ngập nước sâu kéo dài quá 6 tháng kìm hãm sự tăng trưởng sinh khối và giảm trữ lượng tích lũy $CO_2$ của quần thụ Tràm so với chế độ ngập nông.
  • Giả thuyết 3 (H3): Tồn tại mối quan hệ tương tác đa biến có ý nghĩa thống kê giữa các thông số hóa lý môi trường ($BOD_5$, $N\text{-}NO_3^-$, dung trọng đất, chất hữu cơ) với tổng sinh khối rừng Tràm.

Khung lý thuyết của nghiên cứu được xây dựng trên nền tảng Thuyết Thích ứng Sinh thái Vùng ngập nước (Ecological Wetland Adaptation Theory), Thuyết Chịu đựng Yếm khí Thực vật (Plant Anaerobic Stress Theory của Crawford, 1992; Armstrong et al., 1994) và Thuyết Sinh trắc học Allometric Rừng ngập phèn (Allometric Scaling Theory). Nghiên cứu được triển khai trên quy mô thực địa gồm 18 ô tiêu chuẩn cố định ($100\text{ m}^2$/ô) phân bố tại các phân khu bảo vệ nghiêm ngặt (Tiểu khu II, Tiểu khu IV và Tiểu khu 75) thuộc VQG U Minh Hạ trong khung thời gian 24 tháng liên tục (từ tháng 01/2014 đến tháng 12/2015), trải qua 2 mùa mưa và 2 mùa nắng hoàn chỉnh. Đóng góp đột phá của luận án là đã chứng minh bằng số liệu thực nghiệm: sinh khối rừng Tràm đạt giá trị cực đại ($95,9\text{ tấn/ha}$, hấp thụ $147\text{ tấn } CO_2\text{/ha}$) ở độ dày than bùn $20\text{--}40\text{ cm}$ và mức ngập nông $< 30\text{ cm}$ ($91\text{ tấn/ha}$, hấp thụ $143\text{ tấn } CO_2\text{/ha}$), đồng thời chỉ ra việc ngập sâu $> 60\text{ cm}$ làm suy giảm sinh khối nghiêm trọng xuống mức $75\text{ tấn/ha}$.

Literature Review và Positioning

Tổng quan y văn quốc tế và trong nước về sinh thái rừng đầm lầy than bùn (peat swamp forest) phản ánh ba dòng nghiên cứu chính:

Dòng nghiên cứu thứ nhất tập trung vào phương pháp ước tính sinh khối và trữ lượng carbon thông qua quan hệ tương quan sinh trắc học allometric. Các tác giả tiên phong như Wang et al. (1981), Reddy & Sugur (1992), Finlayson (1993), Rayachhetry et al. (2001) và Van et al. (2000, 2002) đã chứng minh đường kính ngang ngực (DBH hay $D_{1.3}$) và chiều cao vút ngọn ($H_{vn}$) là hai biến số chỉ thị tin cậy nhất để dự đoán sinh khối thân, cành, lá của chi Melaleuca, trong đó thành phần gỗ chiếm $83%\text{--}96%$ tổng sinh khối. Tại Việt Nam, các mô hình hồi quy của Phạm Xuân Quý (2008), Bảo Huy (2008), Viên Ngọc Nam & Lư Ngọc Trâm Anh (2011), Viên Ngọc Nam & Lê Hoàng Long (2012) tại Tràm Chim và Gáo Giồng đã chuẩn hóa công thức sinh khối cho rừng Tràm ĐBSCL. Tuy nhiên, các mô hình này chủ yếu áp dụng cho đất phèn khoáng mà chưa tính toán chuyên biệt trên nền đất than bùn có tầng mùn dày.

Dòng nghiên cứu thứ hai phân tích đặc tính lý hóa đất than bùn nhiệt đới và động thái lưu trữ carbon. Andriesse (1988), Siregar & Sambas (1999), Đỗ Đình Sâm et al. (2011) chỉ ra rằng đất than bùn nhiệt đới tuy giàu hữu cơ nhưng rất nghèo dinh dưỡng khả tiêu và nhạy cảm với sụt lún, suy thoái do cháy. Alpian et al. (2013) khi nghiên cứu rừng Tràm tại Central Kalimantan (Indonesia) ghi nhận sinh khối ở độ dày than bùn $51\text{--}100\text{ cm}$ đạt $144,1\text{ tấn/ha}$ và $101\text{--}200\text{ cm}$ đạt $127,8\text{ tấn/ha}$, nhưng chưa phân tích sâu ở các tầng than bùn nông ($< 50\text{ cm}$).

Dòng nghiên cứu thứ ba xoay quanh cuộc tranh luận học thuật về quản lý thủy văn rừng Tràm. Một trường phái (đại diện bởi Vương Văn Quỳnh et al., 2005; Thái Thành Lượm, 2005; Nguyễn Văn Hiệp, 2005) ủng hộ giải pháp đắp đê, trữ nước sâu quanh năm ($0,5\text{--}1,0\text{ m}$) để ngăn chặn triệt để nguy cơ cháy than bùn. Ngược lại, trường phái sinh thái học (đại diện bởi Dương Văn Ni, 2005; Nguyễn Ngọc Anh, 2005; Trần Văn Thắng & Trần Quang Bảo, 2011 tại VQG U Minh Thượng) lập luận rằng việc ngập nước liên tục làm ức chế hệ rễ khí sinh, giảm khả năng hô hấp, làm chết cây và tích tụ vật liệu cháy dạng bổi mục chưa phân hủy, gây nguy cơ đại hỏa hoạn trong tương lai.

Luận án của Trần Thị Kim Hồng định vị ở giao điểm của ba dòng nghiên cứu trên, tạo bước đột phá khi giải quyết tranh luận bằng thực chứng: so sánh đối sánh trực tiếp với nghiên cứu của Tanit Nuyim (2000, 2005) tại đầm lầy than bùn To Daeng (Thái Lan) – nơi mực nước ngập tối ưu được duy trì $10\text{--}40\text{ cm}$ và mật độ sống giảm còn $41%$ ở tuổi 14; và công trình của Takashi Yamanoshita et al. (2001) về cơ chế chịu ngập 6 tháng của M. cajuputi. Luận án chứng minh rằng việc ngập sâu trên $60\text{ cm}$ kéo dài trong hệ thống đê bao khép kín ở U Minh Hạ là nguyên nhân trực tiếp làm suy giảm mật độ lâm phần xuống còn $1.150\text{--}1.200\text{ cây/ha}$ và sinh khối chỉ còn $75\text{ tấn/ha}$, mở ra hướng điều chỉnh chính sách quản lý thủy văn động mang tầm khu vực.

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án mở rộng và làm sâu sắc thêm Thuyết Ứng suất Môi trường Ngập úng (Plant Flooding Stress Theory) của Crawford (1992) và Armstrong et al. (1994) khi áp dụng trên loài cây tiên phong ngập mặn/ngập phèn Melaleuca cajuputi. Mặc dù Tràm sở hữu hệ mô xốp dẫn khí (aerenchyma) và rễ khí sinh biến dưỡng thích nghi cao, nghiên cứu đã chứng minh ngưỡng chống chịu sinh lý của loài này có giới hạn xác định. Khi bị ngập sâu liên tục ($> 60\text{ cm}$) trong môi trường nước phèn yếm khí (oxy hòa tan DO cực thấp, $1,33\text{--}3,00\text{ mg/L}$), chức năng trao đổi chất của hệ rễ bị tê liệt, dẫn đến hiện tượng rễ biến dạng, mất khả năng neo bám cơ học trong nền đất mùn xốp, thúc đẩy tỷ lệ đổ ngã sinh học.

Đồng thời, luận án đóng góp vào Lý thuyết Động thái Sinh khối Lâm phần (Stand Biomass Dynamics Theory) thông qua việc phát hiện "nghịch lý sinh khối giữa cá thể và quần thụ":

  • Ở cấp độ cá thể: Cây Tràm sinh trưởng trên tầng than bùn dày ($60\text{--}80\text{ cm}$) đạt kích thước đường kính và sinh khối cá thể lớn nhất ($77,4\text{--}78,9\text{ kg/cây}$), cao hơn rõ rệt so với cây ở tầng than bùn mỏng $20\text{--}40\text{ cm}$ ($60,7\text{ kg/cây}$).
  • Ở cấp độ quần thụ: Tổng sinh khối trên một đơn vị diện tích lại diễn biến theo chiều ngược lại: đạt cực đại tại tầng than bùn $20\text{--}40\text{ cm}$ ($95,9\text{ tấn/ha}$) và giảm thấp nhất ở tầng than bùn $60\text{--}80\text{ cm}$ ($72,3\text{ tấn/ha}$).

Nguyên nhân lý thuyết được luận giải là do nền than bùn dày có dung trọng cực thấp ($0,19\text{--}0,37\text{ g/cm}^3$) kết hợp ngập nước làm giảm sức tải cơ học của đất, khiến mật độ cây đứng bị suy giảm mạnh ($1.150\text{ cây/ha}$ so với $1.564\text{ cây/ha}$ ở tầng than bùn mỏng), dẫn đến tổng sinh khối tích lũy của cả lâm phần bị sụt giảm.

Khung phân tích độc đáo

Khung phân tích của luận án tích hợp liên ngành ba trụ cột: Khoa học Đất (Pedology), Thủy văn Môi trường nước (Hydrology & Water Quality) và Lâm học Đo học Rừng (Forest Biometrics). Cách tiếp cận phân tích đa biến cho phép bóc tách ảnh hưởng độc lập cũng như tương tác chéo giữa các thông số hóa lý môi trường lên sinh khối thực vật thông qua phương trình hàm tổng quát:

$$TSK = f(\text{Độ dày than bùn}, \text{Độ sâu ngập}, \text{Dung trọng}, \text{CHC}, TN, TP, pH, DO, BOD_5, N\text{-}NO_3^-, N\text{-}NH_4^+, \text{Mùa})$$

Điều kiện biên (boundary conditions) của khung phân tích được xác lập rõ ràng: Áp dụng cho kiểu phụ thổ nhưỡng than bùn nhiệt đới hình thành trên nền trầm tích sét đầm lầy nội địa, thuộc hệ sinh thái rừng úng phèn có đê bao nhân tạo kiểm soát dòng chảy, chịu chế độ khí hậu nhiệt đới gió mùa bán đảo với mùa mưa - khô tương phản sâu sắc.

Phương pháp nghiên cứu tiên tiến

Thiết kế nghiên cứu

Luận án tuân thủ triết lý thực chứng (positivism) với thiết kế nghiên cứu định lượng thực nghiệm sinh thái (ecological empirical design). Thiết kế thí nghiệm phân khối theo mô hình đa nhân tố:

  • Nhân tố A: Cấp độ dày tầng than bùn (3 nghiệm thức: $20\text{--}40\text{ cm}$, $40\text{--}60\text{ cm}$ và $60\text{--}80\text{ cm}$).
  • Nhân tố B: Cấp độ sâu ngập nước (3 nghiệm thức: ngập thấp $< 30\text{ cm}$, ngập trung bình $30\text{--}60\text{ cm}$ và ngập sâu $> 60\text{ cm}$).
  • Nhân tố C: Biến thiên thời gian theo chu kỳ mùa (Mùa mưa: tháng 5 đến tháng 10; Mùa nắng: tháng 11 đến tháng 4 năm sau).

Cỡ mẫu nghiên cứu gồm 18 ô tiêu chuẩn (OTC) diện tích $100\text{ m}^2$ ($10\text{ m} \times 10\text{ m}$) được bố trí ngẫu nhiên đại diện trong các phân khu rừng tự nhiên nguyên sinh thuần loài trên 15 năm tuổi tại Khoảnh 13, 14 (Tiểu khu II), Khoảnh 23, 24, 26, 27 (Tiểu khu IV) và Khoảnh 3 (Tiểu khu 75) thuộc VQG U Minh Hạ.

                             THIẾT KẾ NGHIÊN CỨU THỰC ĐỊA TẠI VQG U MINH HẠ
                               (18 Ô Tiêu Chuẩn OTC - 100 m² / Ô - 2014-2015)

Quy trình nghiên cứu rigorous

Quy trình thu thập và phân tích mẫu đất, nước và lâm sinh được chuẩn hóa nghiêm ngặt:

  • Đo đạc chỉ tiêu sinh thái rừng: Toàn bộ cây Tràm trong 18 OTC được kiểm kê, đánh số cố định. Đo đường kính ngang ngực ($D_{1.3}$) bằng thước kẹp kính chuyên dụng và thước dây chu vi sai số $\pm 0,1\text{ cm}$; đo chiều cao vút ngọn ($H_{vn}$) và chiều cao dưới cành ($H_{dc}$) bằng thước đo cao lâm nghiệp Blume-Leiss.
  • Phân tích lý hóa đất than bùn: Lấy mẫu đất tại các tầng chỉ định bằng ống đai thể tích hình trụ để xác định dung trọng ($g/\text{cm}^3$). Đo $pH_{\text{H}2\text{O}}$ và $pH{\text{KCl}}$ (tỷ lệ 1:2,5) bằng máy đo pH điện tử để bàn. Hàm lượng chất hữu cơ (CHC %) được phân tích theo phương pháp Walkley-Black oxy hóa bằng $K_2Cr_2O_7 + H_2SO_4$ đậm đặc; Tổng Nitơ (TN %) xác định theo phương pháp Micro-Kjeldahl; Tổng Photpho (TP % $P_2O_5$) xác định bằng phương pháp so màu Molybdate sau khi công phá bằng hỗn hợp $H_2SO_4 + HClO_4$; Đạm khoáng ($N\text{-}NH_4^+$, $N\text{-}NO_3^-$) trích ly bằng dung dịch $KCl\text{ 2M}$ và định lượng theo phương pháp chưng cất Kjeldahl.
  • Phân tích chất lượng môi trường nước: Mẫu nước được lấy định kỳ tại các OTC theo tiêu chuẩn TCVN. Thông số nhiệt độ, pH, Oxy hòa tan (DO) đo trực tiếp tại hiện trường bằng máy đo đa chỉ tiêu cầm tay Hanna HI98194. Thông số nhu cầu oxy sinh hóa ($BOD_5$) phân tích theo phương pháp ủ mẫu 5 ngày ở $20^\circ\text{C}$ và đo suy giảm DO (TCVN 6001-1:1995). Nồng độ $N\text{-}NH_4^+$ và $N\text{-}NO_3^-$ xác định theo Tiêu chuẩn Phương pháp Chuẩn Phân tích Nước và Nước thải (SMEWW).

Tam giác hóa phương pháp (Methodological Triangulation) được thực hiện thông qua việc kết hợp điều tra đo đếm thực địa, phân tích hóa nghiệm phòng thí nghiệm (ISO/IEC 17025) và đối chiếu chéo với mô hình hồi quy sinh khối allometric chuẩn khu vực ĐBSCL.

Data và phân tích

Dữ liệu thực nghiệm được xử lý thống kê bằng phần mềm chuyên dụng SPSS Statistics (phiên bản 20.0). Kiểm định phân phối chuẩn và kiểm định tính đồng nhất của phương sai (Levene's test) được thực hiện trước khi phân tích phương sai hai nhân tố (Two-way ANOVA) để đánh giá tương tác giữa độ dày than bùn $\times$ mùa và độ sâu ngập $\times$ mùa. Phép thử hậu kiểm Duncan ($p < 0,05$) được dùng để so sánh sự khác biệt giữa các giá trị trung bình nghiệm thức.

Mô hình hồi quy đa biến từng bước (Stepwise Multiple Linear Regression) được áp dụng để xác định mức độ đóng góp của từng thông số môi trường lên sinh khối lâm phần, kiểm soát hiện tượng đa cộng tuyến qua hệ số phóng đại phương sai (VIF < 2,5).

Phát hiện đột phá và implications

Những phát hiện then chốt

Thứ nhất, xác lập bộ cơ sở dữ liệu chi tiết về đặc tính hóa lý độc đắc của đất than bùn VQG U Minh Hạ. Kết quả phân tích khẳng định đất than bùn nơi đây có độ xốp đặc biệt cao, dung trọng cực thấp dao động trong khoảng $0,19\text{--}0,37\text{ g/cm}^3$, tính chất phèn hóa mạnh với pH dao động từ $3,41\text{--}4,84$. Hàm lượng chất hữu cơ (CHC) áp đảo, chiếm tỷ lệ từ $83,71%\text{--}94,00%$. Về mặt dinh dưỡng khoáng, đất than bùn rất giàu đạm tổng số ($0,58\text{--}1,23%\text{ N}$) nhưng nghèo đến trung bình về lân tổng số ($0,03\text{--}0,12%\text{ P}_2\text{O}_5$). Chất lượng đất có sự biến động rõ rệt giữa mùa mưa và mùa khô ($p < 0,05$), chịu sự chi phối mạnh của quá trình hòa tan và rửa trôi axít hữu cơ.

Thứ hai, phát hiện mối tương quan nghịch giữa độ dày tầng than bùn và tổng sinh khối quần thụ. Ở độ dày than bùn $20\text{--}40\text{ cm}$, sinh khối khô cây cá thể đạt mức thấp nhất ($60,7\text{ kg/cây}$), khác biệt có ý nghĩa thống kê ($p < 0,05$) so với tầng $40\text{--}60\text{ cm}$ ($78,9\text{ kg/cây}$) và tầng $60\text{--}80\text{ cm}$ ($77,4\text{ kg/cây}$). Tuy nhiên, tổng sinh khối rừng Tràm lại đạt giá trị cao nhất ở tầng than bùn $20\text{--}40\text{ cm}$ ($95,9\text{ tấn/ha}$), kế tiếp là tầng $40\text{--}60\text{ cm}$ ($81,1\text{ tấn/ha}$) và thấp nhất ở tầng than bùn dày $60\text{--}80\text{ cm}$ ($72,3\text{ tấn/ha}$). Khả năng hấp thụ $CO_2$ tương ứng đạt đỉnh ở tầng $20\text{--}40\text{ cm}$ với $147\text{ tấn/ha}$ (sai số $\pm 66\text{ tấn/ha}$), vượt trội so với tầng $60\text{--}80\text{ cm}$ chỉ đạt $110\text{ tấn/ha}$ ($\pm 51\text{ tấn/ha}$).

Thứ ba, giải mã tác động tiêu cực của chế độ ngập úng sâu đối với sinh khối và mật độ rừng. Trong môi trường nước ngập, pH nước biến thiên từ $4,21\text{--}4,83$ (mùa mưa chua hơn mùa khô do nước mưa kéo axít từ thảm mục); lượng oxy hòa tan trong nước cực kỳ thiếu hụt (DO chỉ đạt $1,33\text{--}3,00\text{ mg/L}$); trong khi chỉ số $BOD_5$ dao động rất lớn từ $8,0\text{--}53,2\text{ mg/L}$, phản ánh mức độ phân hủy yếm khí mạnh. Mật độ cây Tràm ở nghiệm thức ngập thấp ($< 30\text{ cm}$) duy trì cao nhất ($1.564 \pm 306\text{ cây/ha}$), cho sinh khối rừng đạt $91\text{ tấn/ha}$ và hấp thụ $143\text{ tấn } CO_2\text{/ha}$. Khi độ sâu ngập tăng lên $> 60\text{ cm}$, sinh khối lâm phần tụt giảm xuống $75\text{ tấn/ha}$ do cây bị suy thoái và đổ ngã rễ.

Thứ tư, phương trình hồi quy đa biến nhận diện các nhân tố kiểm soát chính. Mô hình hồi quy tuyến tính chứng minh: các chỉ tiêu mức độ ngập ($p < 0,1$), nồng độ $BOD_5$ và $N\text{-}NO_3^-$ trong nước ($p < 0,1$), cùng với dung trọng đất, chất hữu cơ đất và chu kỳ mùa vụ ($p < 0,05$) là những biến số có tác động thúc đẩy hoặc kìm hãm trực tiếp đến sinh khối tích lũy của rừng Tràm. Ngược lại, các yếu tố như pH, DO, $N\text{-}NH_4^+$ trong nước và TN, TP trong đất chưa thể hiện mối tương quan độc lập rõ rệt trong mô hình đa biến.

Nghiệm thức nghiên cứu Độ dày than bùn / Độ sâu ngập Mật độ cây đứng (cây/ha) Sinh khối cá thể (kg/cây) Sinh khối lâm phần (tấn/ha) Khả năng tích lũy $CO_2$ (tấn/ha)
Cấp than bùn 1 $20\text{--}40\text{ cm}$ $1.564 \pm 306$ $60,7 \pm 12,4$ 95,9 147 ± 66
Cấp than bùn 2 $40\text{--}60\text{ cm}$ $1.320 \pm 215$ $78,9 \pm 15,1$ 81,1 $124 \pm 56$
Cấp than bùn 3 $60\text{--}80\text{ cm}$ $1.150 \pm 180$ $77,4 \pm 14,8$ 72,3 $110 \pm 51$
Cấp ngập nước 1 $< 30\text{ cm}$ $1.564 \pm 306$ $64,2 \pm 11,8$ 91,0 143
Cấp ngập nước 2 $30\text{--}60\text{ cm}$ $1.380 \pm 240$ $68,5 \pm 13,2$ 85,0 123
Cấp ngập nước 3 $> 60\text{ cm}$ $1.200 \pm 195$ $70,1 \pm 14,0$ 75,0 136

Implications đa chiều

  • Về mặt lý thuyết: Cung cấp cơ sở khoa học định lượng bổ sung vào sinh thái học đất ngập nước nhiệt đới, chứng minh loài Melaleuca cajuputi là loài cây "chịu đựng được ngập úng nhưng không tối ưu trong môi trường ngập úng sâu quanh năm".
  • Về mặt phương pháp luận: Chuẩn hóa quy trình tích hợp giữa quan trắc hóa lý đất - nước và đo đạc trắc lượng rừng trên nền đất than bùn nhiệt đới, làm tài liệu tham khảo chuẩn mực cho các nghiên cứu sinh thái đầm lầy tại khu vực Đông Nam Á.
  • Về mặt thực tiễn và chính sách: Đề xuất một lộ trình thay đổi mang tính cách mạng trong quản lý thủy văn tại VQG U Minh Hạ: Chấm dứt việc đóng cống giữ nước ngập sâu liên tục quanh năm. Thay vào đó, vận hành hệ thống cống đập linh hoạt để duy trì mực nước ngập nông tối ưu dưới $30\text{ cm}$ trong mùa khô nhằm vừa bảo vệ lớp than bùn không bị khô kiệt gây cháy, vừa cho phép hệ rễ Tràm tiếp xúc với oxy không khí, đảm bảo cây sinh trưởng mạnh, đạt năng suất sinh khối và hiệu quả hấp thụ $CO_2$ cao nhất.

Limitations và Future Research

Luận án ghi nhận một số giới hạn nghiên cứu khách quan:

  1. Số lượng ô tiêu chuẩn (18 OTC) tuy đáp ứng chuẩn mực thống kê nhưng tập trung chủ yếu tại Phân khu bảo vệ nghiêm ngặt, chưa bao phủ toàn bộ diện tích rừng trồng tái sinh ở các vùng đệm và phân khu phục hồi sinh thái.
  2. Nghiên cứu chỉ tập trung định lượng sinh khối trên mặt đất (Above-Ground Biomass - AGB) của cây gỗ Tràm mà chưa định lượng trực tiếp sinh khối rễ ngầm dưới mặt đất (Below-Ground Biomass - BGB) cũng như sinh khối thảm tươi, dây leo (Stenochlaena palustris, Diplazium esculentum) do tính chất phức tạp và hạn chế về phá hủy sinh cảnh rừng đặc dụng.
  3. Chuỗi số liệu thủy văn và sinh khối được đo đạc trong 2 năm (2014 - 2015), chưa theo dõi liên tục trong một chu kỳ biến đổi khí hậu dài hạn (10 - 20 năm) gắn với các hiện tượng El Niño/La Niña cực đoan.

Chương trình nghiên cứu tương lai cần tập trung vào các định hướng cụ thể:

  • Ứng dụng công nghệ viễn thám siêu phổ (Hyperspectral Remote Sensing) và LiDAR gắn trên thiết bị bay không người lái (UAV) để giải đoán sinh khối và lập bản đồ trữ lượng carbon toàn diện cho toàn bộ $8.527\text{ ha}$ VQG U Minh Hạ.
  • Quan trắc động thái phát thải khí nhà kính thực địa ($CO_2$, $CH_4$, $N_2O$) từ đất than bùn theo các chế độ ngập nước khác nhau bằng phương pháp buồng đo kín (Closed-chamber method) kết hợp Eddy Covariance.
  • Nghiên cứu quần xã vi sinh vật đất và nấm rễ nội cộng sinh (Mycorrhizae) trong việc hỗ trợ cây Tràm hấp thụ lân và cố định đạm trong điều kiện ngập phèn than bùn.

Tác động và ảnh hưởng

Luận án tạo ra tác động học thuật và thực tiễn sâu rộng:

  • Tác động học thuật: Kết quả nghiên cứu là nguồn trích dẫn nền tảng cho các công trình nghiên cứu về bể chứa carbon đất ngập nước (Blue/Teal Carbon) tại ĐBSCL, đóng góp trực tiếp vào hệ thống dữ liệu kiểm kê khí nhà kính quốc gia thuộc khuôn khổ Thỏa thuận Paris và Công ước RAMSAR về đất ngập nước.
  • Tác động chuyển đổi quản lý ngành: Cung cấp luận cứ khoa học trực tiếp để Ban Quản lý VQG U Minh Hạ, Chi cục Kiểm lâm tỉnh Cà Mau và Sở Nông nghiệp & PTNT điều chỉnh "Đề án Quản lý Thủy văn và Phòng chống Cháy rừng", chuyển từ phương thức trữ nước tĩnh cực đoan sang cơ chế điều tiết mực nước động theo mùa.
  • Tác động kinh tế - xã hội: Góp phần nâng cao năng suất sinh khối gỗ rừng Tràm vùng đệm, bảo tồn đa dạng sinh học các loài thủy sản nước ngọt bản địa, bảo vệ sinh kế nuôi ong lấy mật truyền thống dưới tán rừng Tràm của cộng đồng cư dân địa phương.

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

  • Nghiên cứu sinh và Giới học thuật: Tiếp cận khung phương pháp luận đo đạc sinh khối chuẩn xác, bộ dữ liệu thực nghiệm chi tiết về hóa lý than bùn nhiệt đới và các hàm hồi quy tương quan allometric phục vụ nghiên cứu biến đổi khí hậu.
  • Ban Quản lý các Vườn Quốc gia và Khu Bảo tồn: (VQG U Minh Hạ, VQG U Minh Thượng, VQG Tràm Chim) Sở hữu quy trình kỹ thuật vận hành mực nước tối ưu ($< 30\text{ cm}$) để vừa bảo vệ than bùn vừa chống suy thoái lâm phần.
  • Các Cơ quan Quản lý Nhà nước và Hoạch định Chính sách: Bộ Tài nguyên & Môi trường, Bộ Nông nghiệp & Phát triển Nông thôn có căn cứ vững chắc để xây dựng quy chuẩn kỹ thuật quốc gia về quản lý bền vững đất than bùn và định giá dịch vụ môi trường rừng gắn với thị trường tín chỉ carbon ($Carbon\text{ Credits}$).

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 phát hiện và chứng minh "Nghịch lý phân bổ sinh khối cá thể - quần thụ" trên nền đất than bùn nhiệt đới, mở rộng trực tiếp Thuyết Ứng suất Môi trường Ngập úng (Plant Flooding Stress Theory) của Crawford (1992). Luận án xác lập rằng đất than bùn dày tuy cung cấp lượng dinh dưỡng và mùn lớn giúp cây cá thể phát triển chu vi thân to hơn ($78,9\text{ kg/cây}$), nhưng tính chất cơ lý xốp của nền than bùn (dung trọng chỉ $0,19\text{--}0,37\text{ g/cm}^3$) kết hợp ngập úng kéo dài đã triệt tiêu lực bám của rễ, gây đổ ngã hàng loạt và làm giảm mật độ cây đứng, dẫn đến tổng sinh khối lâm phần sụt giảm $24,6%$ so với tầng than bùn mỏng ($72,3\text{ tấn/ha}$ so với $95,9\text{ tấn/ha}$).

2. Phương pháp nghiên cứu của luận án có tính đổi mới như thế nào khi so sánh với các nghiên cứu quốc tế trước đây? So với nghiên cứu của Tanit Nuyim (2000, 2005) tại Thái Lan vốn chỉ quan sát đơn biến mực nước và sinh trưởng, hay công trình của Alpian et al. (2013) tại Indonesia chỉ khảo sát sinh khối tĩnh trên 2 cấp than bùn ($50\text{--}100\text{ cm}$ và $101\text{--}200\text{ cm}$), luận án đã đổi mới bằng thiết kế ma trận tương tác 3 tầng: Tích hợp đồng thời 3 cấp độ dày than bùn ($20\text{--}40$, $40\text{--}60$, $60\text{--}80\text{ cm}$) $\times$ 3 mức ngập nước ($<30$, $30\text{--}60$, $>60\text{ cm}$) $\times$ 4 chu kỳ mùa (mưa - nắng qua 2 năm), kết hợp phân tích hồi quy đa biến giữa các thông số hóa lý môi trường ($pH, DO, BOD_5, N\text{-}NH_4^+, N\text{-}NO_3^-, \text{dung trọng}, CHC, TN, TP$) với sinh khối lâm 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à nồng độ đạm tổng số trong đất than bùn U Minh Hạ đạt mức rất cao ($0,58\text{--}1,23%\text{ N}$), hàm lượng chất hữu cơ lên tới $94%$, nhưng đây không phải là môi trường tối ưu cho rừng Tràm phát triển mật độ cao khi bị ngập nước liên tục. Nước giữ lại trong rừng chủ yếu là nước mưa có tính axít cao ($pH = 4,21\text{--}4,83$), oxy hòa tan cực thấp ($DO = 1,33\text{--}3,00\text{ mg/L}$), khiến lượng đạm hữu cơ này không thể khoáng hóa hiệu quả để cây hấp thu, biến môi trường ngập sâu thành "vùng nghèo dinh dưỡng sinh lý".

4. Luận án có cung cấp quy trình nhân rộng (Replication Protocol) cụ thể không? Có. Luận án cung cấp bộ quy chuẩn đo đếm thực địa chi tiết gồm: thiết lập 18 ô mẫu chuẩn $100\text{ m}^2$, quy trình lấy mẫu dung trọng bằng ống đai không phá hủy cấu trúc đất, phương pháp phân tích 7 chỉ tiêu đất và 5 chỉ tiêu nước theo TCVN/SMEWW, cùng hệ thống phương trình allometric hồi quy tương quan sinh khối khô từ đường kính $D_{1.3}$ chuẩn hóa cho hệ sinh thái rừng Tràm ĐBSCL.

5. Chương trình nghị sự nghiên cứu 10 năm (10-year Research Agenda) được phác thảo ra sao? Chương trình nghiên cứu 10 năm gồm 4 trụ cột chiến lược:

  • Xây dựng tháp quan trắc Eddy Covariance tự động đo đạc phát thải ròng carbon ($NEE$) và thông lượng khí nhà kính ($CH_4, CO_2$) trên đất than bùn U Minh Hạ.
  • Ứng dụng công nghệ LiDAR viễn thám kết hợp trí tuệ nhân tạo (AI) giám sát biến động sinh khối và mực nước thời gian thực.
  • Thiết lập các ô thí nghiệm thao túng thủy văn (hydrological manipulation experiments) mô phỏng các kịch bản El Niño hạn hán khốc liệt.
  • Xây dựng cơ chế kinh tế tài chính carbon sinh thái, tích hợp rừng Tràm than bùn vào thị trường tín chỉ carbon quốc tế.

Kết luận

Luận án tiến sĩ của Trần Thị Kim Hồng đã đạt được 5 đóng góp học thuật và thực tiễn mang tính then chốt:

  1. Định lượng toàn diện đặc tính thổ nhưỡng than bùn: Xác lập bản chất đất than bùn VQG U Minh Hạ có dung trọng cực thấp ($0,19\text{--}0,37\text{ g/cm}^3$), độ xốp cao ($81,4%\text{--}87,2%$), giàu chất hữu cơ ($83,71%\text{--}94,00%$) và tổng đạm ($0,58\text{--}1,23%$), nhưng có độ phèn cao ($pH = 3,41\text{--}4,84$).
  2. Khám phá ngưỡng sinh thái than bùn tối ưu: Chứng minh tầng than bùn $20\text{--}40\text{ cm}$ cho sinh khối lâm phần cao nhất ($95,9\text{ tấn/ha}$) và khả năng tích lũy $CO_2$ cao nhất ($147\text{ tấn/ha}$); tầng than bùn càng dày sinh khối lâm phần càng suy giảm.
  3. Xác lập giới hạn ngập nước sinh lý: Chứng minh chế độ ngập nông $< 30\text{ cm}$ là điều kiện lý tưởng nhất cho rừng Tràm phát triển sinh khối ($91\text{ tấn/ha}$) và hấp thụ $CO_2$ ($143\text{ tấn/ha}$), trong khi chế độ ngập sâu $> 60\text{ cm}$ kìm hãm sinh thái và làm giảm sinh khối xuống còn $75\text{ tấn/ha}$.
  4. Mô hình hóa tác động đa biến: Thiết lập thành công mô hình hồi quy chứng minh tác động đồng thời có ý nghĩa thống kê của mức ngập, $BOD_5$, $N\text{-}NO_3^-$ trong nước và dung trọng, chất hữu cơ trong đất lên sinh khối rừng.
  5. Đổi mới tư duy quản trị tài nguyên rừng: Cung cấp cơ sở khoa học xác đáng để chuyển đổi mô hình quản lý thủy văn từ "trữ nước ngập sâu chống cháy" sang "điều tiết thủy văn động thích ứng sinh thái", kiến tạo di sản khoa học phục vụ bảo tồn đa dạng sinh học và phát triển bền vững tài nguyên rừng đất ngập nước Việt Nam.