Tổng quan về luận án

Bảo tồn và phục hồi các hệ sinh thái đất ngập nước nội địa (Inland Wetlands) trên nền đất than bùn nhiệt đới là một trong những thách thức sinh thái - kỹ thuật phức tạp nhất tại khu vực hạ lưu lưu vực sông Mê Kông. Luận án tiến sĩ kỹ thuật chuyên ngành Môi trường đất và nước của nghiên cứu sinh Phạm Văn Tùng (2017) với đề tài "Nghiên cứu, đề xuất chế độ nước hợp lý để phát triển rừng tràm tái sinh Vườn quốc gia U Minh Thượng" giải quyết trực tiếp cuộc khủng hoảng quản lý thủy văn sinh thái sau thảm họa cháy lịch sử tháng 3/2002 thiêu rụi 3.212 ha rừng nguyên sinh trên đất than bùn.

Khoảng trống nghiên cứu then chốt (Research Gap) xuất phát từ nghịch lý quản lý nước đơn tầng: nhằm ngăn ngừa triệt để nguy cơ cháy lớp than bùn trong mùa khô, Ban quản lý Vườn quốc gia (VQG) đã áp dụng giải pháp đóng cống, đắp đập giữ nước ở cao trình cao quanh năm ($H_{nam} = 117\text{ cm}$, vượt cao độ địa hình trung bình $Z_{tb} = 105\text{ cm}$ đến 12 cm trong giai đoạn 2003–2009). Biện pháp này làm gia tăng diện tích ngập úng vĩnh viễn lên 50% (3.979 ha), triệt tiêu hoàn toàn chu kỳ bán ngập tự nhiên, dẫn đến suy thoái nghiêm trọng sinh khối, ức chế tái sinh chồi/hạt của cây tràm (Melaleuca cajuputi) và làm biến đổi sinh cảnh quần xã động thực vật đặc trưng.

Nghiên cứu tập trung giải quyết 3 câu hỏi nghiên cứu và hệ giả thuyết khoa học:

  1. $RQ_1$ & $H_1$: Mức độ và thời gian ngập nước tương tác như thế nào đối với các chỉ tiêu sinh trưởng lâm sinh (đường kính $D_{1,3}$, chiều cao vút ngọn $H_{vn}$, sinh khối $TSK$) của rừng tràm tái sinh? Giả thuyết $H_1$: Tồn tại ngưỡng ngập tối ưu mà tại đó cây tràm đạt sinh khối cực đại và tỷ lệ cây tốt cao nhất, vượt qua ngưỡng này sự ngập úng liên tục sẽ gây thoái hóa sinh học.
  2. $RQ_2$ & $H_2$: Chế độ vận hành thủy lực theo không gian và thời gian nào có thể dung hòa đồng thời mục tiêu phòng cháy rừng than bùn và phục hồi sinh cảnh ngập nước theo mùa? Giả thuyết $H_2$: Việc phân chia lưu vực nội đồng thành các tiểu vùng thủy lực đa bậc kết hợp xác định thời điểm tích nước theo xác suất mưa sẽ tối ưu hóa cán cân ẩm.
  3. $RQ_3$ & $H_3$: Quy mô công trình thủy lợi (cống, đập tràn, trạm bơm) cần thiết lập ra sao để kiểm soát chế độ thủy văn nội đồng khép kín? Giả thuyết $H_3$: Hệ thống đập tràn tự do kết hợp cống tiêu hai chiều và trạm bơm cơ động sẽ đảm bảo thoát lũ khẩn cấp trong mùa mưa và duy trì cao trình mực nước khống chế trong mùa khô.

Khung lý thuyết của công trình tích hợp giữa Lý thuyết Sinh thái Thủy văn (Eco-hydrology Theory) và Lý thuyết Phân tích Hệ thống Kỹ thuật Tài nguyên Nước (Water Resources Systems Engineering). Phạm vi nghiên cứu bao quát toàn bộ vùng lõi VQG U Minh Thượng với quy mô 8.003 ha (trong tổng thể 21.107 ha bao gồm cả vùng đệm) trong chuỗi thời gian điều tra thực địa liên tục 8 năm (2009–2016) kết hợp chuỗi số liệu khí tượng thủy văn 31 năm (1985–2015). Nghiên cứu mang lại đóng góp đột phá khi xác định chính xác lịch biểu tích nước theo tần suất mưa thiết kế: ngày 11/9 cho năm ít nước ($P = 75%$), ngày 01/10 cho năm trung bình ($P = 50%$) và ngày 21/10 cho năm nhiều nước ($P = 25%$), giúp khôi phục 93% diện tích ngập theo mùa tại các khu vực thấp trũng.

Literature Review và Positioning

Nghiên cứu về thủy văn đất ngập nước và phục hồi hệ sinh thái rừng tràm trên nền than bùn phân hóa thành hai dòng học thuật chính với các tranh luận đối lập sâu sắc:

[Dòng 1: Quản lý Bảo tồn Sinh thái]          [Dòng 2: Kỹ thuật Lâm nghiệp & PCCC]
(Hoàng Văn Thắng, Lê Diên Dực, 2004)        (Phạm Trọng Thịnh, 2005; Safford, 2002)
   Ưu tiên duy trì nhịp điệu ngập mùa           Ưu tiên giữ nước mức tối đa quanh năm
   Nhằm bảo tồn đa dạng sinh học                Để chống cháy tầng than bùn triệt để
                     \                               /
                      \                             /
                   [Bế tắc: Ngập úng làm chết rừng tràm]
                                    |
            [POSITIONING CỦA LUẬN ÁN PHẠM VĂN TÙNG, 2017]
  Tích hợp Phân tích Hệ thống Thủy lực Đa bậc + Sinh thái Thủy văn
   -> Thiết lập Chế độ Nước Tối ưu Động theo Tần suất Mưa Thiết kế

Dòng thứ nhất, đại diện bởi Hoàng Văn Thắng và Lê Diên Dực (2004), Lê Diên Dực (1989), Phan Nguyên Hồng và cộng sự (1996), Nguyễn Chu Hồi (1999), tiếp cận theo hướng phân loại đất ngập nước của IUCN (với các mã sinh cảnh $X_p, T_s, 1, 3$ tương ứng phân loại Việt Nam $A\text{ }II\text{ }11, A\text{ }II\text{ }16, B\text{ }III\text{ }1, B\text{ }IV\text{ }2$). Quan điểm này nhấn mạnh tầm quan trọng cốt lõi của nhịp điệu ngập nước theo mùa (Seasonal Pulsing) đối với chu trình sinh địa hóa và sự đa dạng của các quần xã chim nước, thực vật thủy sinh. Dòng thứ hai, đại diện bởi các nhà kỹ thuật lâm nghiệp và phòng cháy chữa cháy rừng (Phạm Trọng Thịnh, 2005; Thái Thành Lượm, 2002), khẳng định than bùn U Minh là vật liệu cháy đặc biệt nguy hiểm khi độ ẩm ranh giới tụt dưới ngưỡng bén lửa; do đó bắt buộc phải dâng cao mực nước kênh để mao dẫn làm ẩm tầng đất mặt.

Sự bế tắc của y văn nằm ở chỗ chưa có công trình nào định lượng hóa được mối quan hệ hàm số giữa các mức ngập thủy lực thực tế và động thái tăng trưởng sinh khối của rừng tràm tái sinh sau thảm họa cháy, cũng như thiếu một mô hình toán thủy văn điều tiết đa vùng chi tiết.

Trên trường quốc tế, công trình được định vị tương đương với các nghiên cứu phục hồi đất than bùn quy mô lớn:

  • So sánh với nghiên cứu phục hồi than bùn nhiệt đới vùng Central Kalimantan - Indonesia của Ritzema et al. (2014) và Page et al. (2009): Các tác giả Indonesia tập trung vào kỹ thuật chặn kênh thoát nước tự nhiên (Canal Blocking) để nâng mực nước ngầm diện rộng chống cháy. Tuy nhiên, giải pháp của Page et al. chưa giải quyết được bài toán phân bậc địa hình phức tạp, dẫn đến hiện tượng úng cục bộ ở hạ lưu kênh chặn. Luận án của Phạm Văn Tùng đã giải quyết triệt để hạn chế này bằng cách kết hợp phân vùng thủy lực độc lập và kiểm soát xả tràn linh hoạt.
  • So sánh với Chương trình Phục hồi Thủy văn Vườn quốc gia Everglades, Hoa Kỳ (CERP - Comprehensive Everglades Restoration Plan) của Sklar et al. (2005) và Light & Dineen (1994): CERP áp dụng mô hình điều tiết lưu vực mở với quy mô hàng trăm nghìn hecta chịu ảnh hưởng của thủy triều và dòng chảy tràn lục địa. Ngược lại, U Minh Thượng là hệ sinh thái "thủy văn nội đồng khép kín" độc nhất vô nhị, hoàn toàn cách ly với triều biển Tây và chỉ lấy nguồn cấp từ mưa tại chỗ. Công trình của tác giả đã xác lập chuẩn mực mới cho kỹ thuật điều tiết tài nguyên nước cục bộ trong điều kiện cô lập thủy văn tuyệt đối.

Đó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 Lý thuyết Sinh thái Thủy văn của Zalewski et al. (1997) và Rodríguez-Iturbe (2000) bằng cách chứng minh rằng mối quan hệ giữa chế độ ngập ($H$) và sinh khối tích lũy ($TSK$) của loài cây ngập nước nhiệt đới Melaleuca cajuputi không phải là một hàm tuyến tính đơn điệu mà tuân theo quy luật đường cong parabol cực trị. Cụ thể, việc duy trì mực nước ngập sâu liên tục ($>0,8\text{ m}$ so với mặt đất) làm suy giảm nghiêm trọng khả năng hô hấp rễ, ức chế tạo tầng rễ khí sinh thứ cấp, làm giảm $34,8%$ lượng tăng trưởng trữ lượng gỗ hàng năm so với sinh cảnh bán ngập.

Công trình xác lập mô hình lý thuyết điều tiết nước sinh thái thông qua 3 mệnh đề khoa học (Theoretical Propositions):

  • Mệnh đề $P_1$: Động lực tăng trưởng sinh khối rừng tràm tái sinh đạt trạng thái tối ưu sinh học khi thời gian rút nước khô mặt đất ($T_{kh.ngap}$) đạt từ 3,5 đến 6 tháng/năm, trùng khớp với pha phát triển chồi ngọn và phát tán hạt giống tự nhiên.
  • Mệnh đề $P_2$: Rủi ro bùng phát cháy ngầm tầng than bùn tỷ lệ nghịch phi tuyến với độ ẩm mặt đất, trong đó độ ẩm mao dẫn được duy trì an toàn khi mực nước ngầm ($H_{n.ngam}$) cách mặt đất không quá 30–40 cm vào đỉnh điểm mùa khô (tháng 4).
  • Mệnh đề $P_3$: Cân bằng sinh thái trong đầm lầy than bùn khép kín chỉ đạt được khi tỷ lệ sinh cảnh được cấu trúc theo tỷ lệ vàng: $0%$ ngập quanh năm tại các tiểu vùng phục hồi, $\ge 80%$ diện tích ngập theo mùa ($S_{ng.mua}$), và diện tích có nguy cơ cháy cao trong tháng kiệt nhất được kiểm soát dưới $20%$.

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

Khung phân tích của luận án là sự hợp nhất liên ngành giữa 3 nền tảng: (1) Thổ nhưỡng - Địa mạo học than bùn, (2) Động lực học Lâm sinh thực nghiệm, và (3) Thủy lực công trình tài nguyên nước.

+-------------------------------------------------------------------------------+
|                    KHUNG PHÂN TÍCH THỦY VĂN - SINH THÁI ĐA TẦNG               |
+-------------------------------------------------------------------------------+
                                        |
        +-------------------------------+-------------------------------+
        |                               |                               |
        v                               v                               v
[1. Địa mạo & Đất]              [2. Lâm sinh Học]             [3. Cân bằng Nước]
- 2.106 ha than bùn             - Đo 16 cây tiêu chuẩn        - Chuỗi mưa 31 năm (Pearson III)
- Phẫu diện than/sét            - Tương quan D1,3 & SK        - Bốc thoát hơi: 2,77 mm/ngày
- Bản đồ cao độ Ztb             - Kiểm định ANOVA mức ngập    - Rò rỉ qua đê bao: 11%
        |                               |                               |
        +-------------------------------+-------------------------------+
                                        |
                                        v
                    +---------------------------------------+
                    |    MÔ HÌNH HÓA THỦY LỰC ĐA BẬC        |
                    |    Phân chia Khu A, B, C độc lập      |
                    +---------------------------------------+
                                        |
                                        v
                    +---------------------------------------+
                    |    CHẾ ĐỘ NƯỚC ĐIỀU TIẾT TỐI ƯU       |
                    |    Lịch tích nước P = 25%, 50%, 75%   |
                    +---------------------------------------+

Cách tiếp cận phân tích mới nằm ở việc chuyển hóa bài toán quản lý nước từ tĩnh (một mực nước cho toàn VQG) sang quản lý động đa bậc không gian. Luận án định nghĩa tường minh khái niệm "Chế độ nước hợp lý cho rừng tràm than bùn" là tập hợp các đường quá trình mực nước $H(t)$ theo thời gian trong năm, thỏa mãn đồng thời 3 điều kiện biên (Boundary Conditions):

  • Điều kiện sinh học: Cung cấp đủ tầng thoáng khí cho rễ tràm trong mùa khô ít nhất 90–120 ngày.
  • Điều kiện an toàn nhiệt: Duy trì độ ẩm vật liệu cháy dưới tán rừng ngăn ngừa bén lửa trong 5 tháng mùa khô.
  • Điều kiện cân bằng trữ lượng ẩm: Tổng dung tích trữ cuối mùa mưa $W_{tru}$ bù đắp đủ lượng hao hụt bốc thoát hơi mặt thoáng, thoát hơi tầng tán và rò rỉ đê bao ($2,77\text{ mm/ngày}$).

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 dụng thế giới quan Thực chứng (Positivism) kết hợp Phân tích Hệ thống Kỹ thuật (Systems Engineering), bảo đảm tính khách quan tuyệt đối thông qua dữ liệu định lượng đo đạc trực tiếp. Thiết kế nghiên cứu kết hợp giữa khảo sát thực địa lâm sinh học cắt ngang nhiều thời kỳ (2009, 2011, 2014, 2016) và quan trắc chuỗi thời gian liên tục (Time-series) về khí tượng, mực nước mặt và mực nước ngầm.

Phạm vi không gian bao quát toàn bộ 8.003 ha vùng lõi, phân tầng theo độ sâu ngập nước và độ dày tầng than bùn:

  • Lớp than bùn mỏng ($0 - <30\text{ cm}$): 646 ha.
  • Lớp than bùn trung bình ($30 - <60\text{ cm}$): 546 ha.
  • Lớp than bùn dày ($60 - <90\text{ cm}$): 657 ha.
  • Lớp than bùn rất dày ($\ge 90\text{ cm}$): 257 ha.
  • Đất không than bùn/cháy triệt để dạng da quy: 1.106 ha.

Quy trình nghiên cứu rigorous

Quy trình thu thập và xử lý dữ liệu được thiết lập theo các giao thức tiêu chuẩn quốc tế:

  • Khảo sát lâm sinh: Thiết lập hệ thống các Ô tiêu chuẩn (OTC) diện tích $500\text{ m}^2$ ($20\text{ m} \times 25\text{ m}$) và $1000\text{ m}^2$ định vị tọa độ GPS chính xác (Phụ lục 17). Tiến hành đo đếm toàn diện đường kính thân tại vị trí $1,3\text{ m}$ ($D_{1,3}$), chiều cao vút ngọn ($H_{vn}$), chiều cao dưới cành ($H_{dc}$), đường kính tán ($D_{tan}$) và mật độ cây ($N\text{ cây/ha}$).
  • Giải tích sinh khối phá hủy (Destructive Sampling): Chặt hạ và giải tích chi tiết 16 cây tiêu chuẩn đại diện cho các cấp đường kính để xác định sinh khối tươi các bộ phận: thân gỗ ($SKTt$), vỏ ($SKVt$), cành ($SKCt$), lá ($SKLt$). Lấy mẫu sấy khô kiệt ở nhiệt độ $105^\circ\text{C}$ đến khối lượng không đổi để xác định hệ số chuyển đổi khô/tươi ($k = p/P$), từ đó tính toán sinh khối khô thực tế ($SKGkvk, SKVk, SKCk, SKLk, TSKk$).
  • Mạng lưới quan trắc thủy văn: Thu thập số liệu mực nước mặt hàng ngày trên hệ thống kênh chính và 7 tuyến đo thủy văn chuyên biệt cắt ngang vùng lõi (gồm 82 ống đo mực nước ngầm $H_{n.ngam}$ sâu từ 1,5 đến 3,0 m xuyên qua tầng than bùn vào tầng khoáng).
  • Tam giác đạc dữ liệu (Data Triangulation): Đối chiếu chéo giữa số liệu ảnh vệ tinh GIS giải đoán thảm phủ (năm 2002, 2006, 2011), bản đồ cao độ địa hình kỹ thuật số DEM với mắt lưới chi tiết, và số đọc thực tế từ các trạm thủy chí nội đồng.

Data và phân tích

Số liệu khí tượng thủy văn được trích xuất từ chuỗi số liệu 31 năm (1985–2015) của Trạm Khí tượng Rạch Giá. Tác giả áp dụng giải thuật phân phối xác suất thống kê Pearson loại III (Pearson Type III Distribution) để tính toán lượng mưa thiết kế cực trị $X_p$ ứng với các tần suất xuất hiện $P = 25%, 50%, 75%$ và lượng mưa lũ lớn nhất các thời đoạn 1, 3, 5, 7 ngày:

$$X_p = X_{tb} \cdot K_p(C_v, C_s)$$

Trong đó:

  • $X_{tb}$ là lượng mưa trung bình nhiều năm ($2.068\text{ mm/năm}$).
  • $C_v$ là hệ số phân tán; $C_s$ là hệ số thiên lệch ($C_s = 2C_v$).
  • $K_p$ là hệ số mô-đun phụ thuộc tần suất.

Tính toán tổn thất ẩm được mô hình hóa theo phương trình cân bằng ẩm thực nghiệm:

$$q_{hao} = E_{tan} + E_{bui} + E_{dat} + E_{mat} + Q_{ro}$$

Kết quả xác định tổng lượng nước thất thoát bình quân mùa mưa là $2,77\text{ mm/ngày}$, trong đó lượng bốc thoát hơi mặt đệm chiếm $89%$ và lượng rò rỉ qua thân đê bao chiếm $11%$ ($1.085.666\text{ m}^3/49\text{ ngày}$).

Quan hệ sinh trưởng sinh khối được phân tích hồi quy phi tuyến đa biến trên phần mềm thống kê chuyên dụng. Sự sai khác sinh khối giữa các nghiệm thức ngập nước được kiểm định thông qua phân tích phương sai một nhân tố (One-way ANOVA) tại mức ý nghĩa thống kê $\alpha = 0,05$ ($p < 0,001$). Tính toán thủy lực kích thước công trình áp dụng phương trình dòng chảy qua cống ngập và đập tràn đỉnh rộng:

$$Q_C = \mu \cdot \omega_d \cdot \sqrt{2g(H_0 - h_h)}$$

$$Q_{tran} = m \cdot B \cdot \sqrt{2g} \cdot H_0^{3/2}$$

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

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

Luận án công bố 5 phát hiện cốt lõi mang tính bước ngoặt khoa học:

Chỉ tiêu lâm sinh / Thủy văn Chế độ ngập sâu liên tục (2003–2009) Chế độ ngập theo mùa tối ưu (Đề xuất) Mức độ cải thiện / Giá trị kiểm định
Đường kính bình quân ($D_{1,3}$) $6,82 \pm 0,45\text{ cm}$ $9,45 \pm 0,38\text{ cm}$ Tăng $+38,5%$ ($p < 0,01$)
Chiều cao vút ngọn ($H_{vn}$) $7,12 \pm 0,52\text{ m}$ $10,85 \pm 0,61\text{ m}$ Tăng $+52,4%$ ($p < 0,01$)
Trữ lượng gỗ bình quân ($M$) $32,15\text{ m}^3/\text{ha}$ $58,40\text{ m}^3/\text{ha}$ Tăng $+81,6%$ ($p < 0,001$)
Tổng sinh khối khô ($TSK_k$) $24,80\text{ tấn/ha}$ $45,62\text{ tấn/ha}$ Tăng $+83,9%$ ($F = 42,6; p < 0,001$)
Tỷ lệ diện tích ngập theo mùa $29%\text{ (2.376 ha)}$ $93%\text{ tại Khu A, B (2.535 ha)}$ Mở rộng $+64%$ diện tích sinh thái tối ưu

Thứ nhất, luận án phát hiện hiện tượng sinh lý nghịch thường: cây tràm sống trong điều kiện úng nước kinh niên hình thành lớp giác gỗ xốp, mật độ rễ cọc tiêu giảm nghiêm trọng và chỉ phát triển hệ rễ thở bề mặt yếu ớt. Khi gặp gió xoáy mùa mưa, tỷ lệ đổ ngã gãy đổ lên tới $18,4%$. Ngược lại, tại chế độ bán ngập (rút nước 3,5–4 tháng), tổng sinh khối khô đạt đỉnh $45,62\text{ tấn/ha}$, hệ rễ ăn sâu bám chắc vào tầng khoáng dưới than bùn.

Thứ hai, chứng minh sự thất bại của quy hoạch quản lý nước một bậc giai đoạn 2002–2009. Việc đóng kín cống giữ mực nước $140\text{ cm}$ đã biến $50%$ diện tích VQG thành đầm lầy chết, tiêu diệt sự tái sinh của thảm tươi bản địa (Phragmites australis, Eleocharis dulcis) và ngăn cản quá trình khoáng hóa than bùn hiếu khí.

Thứ ba, xác lập bản đồ phân vùng thủy lực tối ưu 3 khu vực chức năng dựa trên cao trình địa hình tự nhiên:

  • Khu A & B (2.723 ha, $Z_{tb} = +0,83\text{ m}$): Vùng địa hình trũng, khống chế mực nước mùa khô ở mức $+0,63\text{ m}$ đến $+0,70\text{ m}$, giải phóng $93%$ diện tích khỏi ngập úng quanh năm.
  • Khu C (5.280 ha, $Z_{tb} = +1,22\text{ m}$): Vùng địa hình gò cao chứa phần lớn diện tích than bùn dày, khống chế mực nước kiệt tháng 4 ở mức $+1,17\text{ m}$, bảo đảm diện tích có nguy cơ cháy cao tuyệt đối không vượt quá $20%$.

Thứ tư, xác định chính xác thời điểm đóng cống tích nước hàng năm căn cứ vào lượng mưa tích lũy cuối mùa mưa:

  • Năm ít nước ($P = 75%$): Bắt đầu tích nước từ ngày 11 tháng 9.
  • Năm nước trung bình ($P = 50%$): Bắt đầu tích nước từ ngày 01 tháng 10.
  • Năm nhiều nước ($P = 25%$): Bắt đầu tích nước từ ngày 21 tháng 10.

Thứ năm, phát hiện về cơ chế rò rỉ thủy lực: Thân đê bao vùng lõi dài 38 km bị tổn thất thẩm lậu qua thân đê và nền đất sét pha cát trung bình $1,085\text{ triệu m}^3$ trong mùa mưa, đòi hỏi mặt bằng tích nước thiết kế phải tính dôi dư $11%$ dung tích hữu ích.

Implications đa chiều

  • Ý nghĩa Lý thuyết: Cung cấp bằng chứng thực nghiệm chuẩn xác khẳng định quy luật nhịp điệu ngập nước trong phục hồi rừng ngập nước nội địa nhiệt đới; làm cầu nối hoàn hảo giữa lý thuyết lâm sinh và thủy lực công trình.
  • Đổi mới Phương pháp luận: Thiết lập quy trình 4 bước tích hợp: Phân tích DEM địa hình $\rightarrow$ Mô hình hóa cân bằng ẩm Pearson III $\rightarrow$ Giải tích sinh khối thực nghiệm $\rightarrow$ Thiết kế thủy lực đa bậc. Phương pháp này có thể chuyển giao nguyên vẹn cho VQG Tràm Chim, VQG U Minh Hạ và Khu bảo tồn Láng Sen.
  • Ứng dụng Thực tiễn: Cung cấp cuốn cẩm nang vận hành quy chuẩn cho Ban quản lý VQG U Minh Thượng, chấm dứt hoàn toàn tình trạng đóng mở cống cảm tính, loại trừ nguy cơ cháy rừng trong khi vẫn tăng gấp đôi tốc độ tăng trưởng rừng tràm.
  • Khuyến nghị Chính sách: Kiến nghị Bộ Nông nghiệp và PTNT cùng UBND tỉnh Kiên Giang phê duyệt đề án hiện đại hóa hệ thống 3 cống hộp tự động hóa, 2 đập tràn điều tiết phân khu và nạo vét $38\text{ km}$ kênh bao nhằm chủ động ứng phó biến đổi khí hậu.

Limitations và Future Research

Mặc dù đạt được những kết quả đột phá, luận án thẳng thắn ghi nhận 4 giới hạn nghiên cứu mang tính khách quan:

  1. Giới hạn số liệu bốc thoát hơi: Do điều kiện thực địa thiếu tháp đo vi khí tượng Eddy Covariance chuyên dụng, số liệu bốc thoát hơi mặt đệm phải tính chuyển gián tiếp thông qua bốc hơi ống Piche tại trạm Rạch Giá (cách vùng nghiên cứu 30 km) kết hợp hệ số thực nghiệm, chưa phản ánh đầy đủ độ biến thiên vi khí hậu dưới tán rừng dày.
  2. Giới hạn về mô hình thấm: Nghiên cứu giả định hệ số thấm qua thân đê bao và nền than bùn là đồng nhất theo không gian; chưa mô phỏng được dòng thấm cục bộ không bão hòa qua các khe nứt than bùn khi bị co ngót mạnh trong mùa khô hạn cực đoan.
  3. Giới hạn chuỗi điều tra sinh học: Số liệu sinh khối mới chỉ tập trung vào cấu trúc tầng cây gỗ tràm (Melaleuca cajuputi), chưa định lượng chi tiết sinh khối ngầm của hệ rễ mịn dưới độ sâu $>1,0\text{ m}$ và chưa lượng hóa đầy đủ chỉ số đa dạng sinh học của quần xã thủy sinh vật đáy.
  4. Điều kiện biên về biến đổi khí hậu: Các kịch bản tính toán thủy văn mới chỉ dựa trên chuỗi số liệu quá khứ 31 năm, chưa tích hợp kịch bản nước biển dâng và hiện tượng El Niño/La Niña cực đoan gia tăng tần suất khô hạn tại Bán đảo Cà Mau.

Chương trình nghiên cứu 10 năm tiếp theo (Future Research Agenda) cần tập trung vào 4 hướng ưu tiên:

  • Ứng dụng công nghệ ảnh viễn thám không gian LiDAR và Drone siêu phổ để quan trắc độ ẩm tầng than bùn thời gian thực.
  • Thiết lập mô hình toán thủy văn - sinh thái 3 chiều (3D Eco-hydrological Model) mô phỏng dòng thấm đa pha trong môi trường than bùn nứt nẻ.
  • Nghiên cứu cơ chế hấp thụ và phát thải khí nhà kính ($CO_2, CH_4$) của đất than bùn dưới các chế độ ngập nước điều tiết khác nhau nhằm phục vụ thị trường tín chỉ Carbon rừng ngập nước.
  • Đánh giá toàn diện tác động của xâm nhập mặn từ biển Tây qua hệ thống kênh rạch vùng đệm đến sự ổn định hóa học của tầng sinh phèn ($C_p$) nằm sâu dưới lớp than bùn.

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

Công trình của Tiến sĩ Phạm Văn Tùng tạo ra những xung lực biến đổi mạnh mẽ trên 5 phương diện:

  • Tác động Học thuật: Thiết lập một điển cứu mẫu mực (Benchmark Case Study) cho ngành kỹ thuật môi trường đất và nước tại Việt Nam; mở ra hướng nghiên cứu liên ngành Kỹ thuật Thủy lợi - Bảo tồn Đa dạng sinh học. Dự kiến mang lại chỉ số trích dẫn cao trong các nghiên cứu về phục hồi đất than bùn Đông Nam Á.
  • Chuyển đổi Ngành Kỹ thuật Thủy lợi: Thay đổi căn bản tư duy thiết kế công trình thủy lợi truyền thống (vốn chỉ thuần túy phục vụ tưới tiêu nông nghiệp) sang thiết kế công trình thủy lợi sinh thái đa mục tiêu phục vụ bảo tồn tài nguyên rừng.
  • Hiệu quả Xã hội & Vùng đệm: Việc điều tiết xả nước hợp lý từ vùng lõi trong mùa khô cung cấp nguồn nước ngọt quý giá cứu trợ cho 13.069 ha đất sản xuất nông nghiệp và nuôi trồng thủy sản của người dân hai xã An Minh Bắc và Minh Thuận, giảm thiểu mâu thuẫn sinh kế giữa cộng đồng địa phương và khu bảo tồn.
  • Bảo tồn Đa dạng Sinh học Toàn cầu: Bảo vệ an toàn nguồn gen quý hiếm của VQG U Minh Thượng - Vườn Di sản ASEAN (ASEAN Heritage Park) và Khu Ramsar thứ 2.228 của thế giới, bảo tồn sinh cảnh sống của 188 loài chim, 32 loài thú (trong đó có Rái cá lông mũi Lutra sumatrana đặc hữu).
  • Lợi ích Kinh tế Lượng hóa: Ngăn ngừa tổn thất do cháy rừng ước tính hàng trăm tỷ đồng mỗi năm, đồng thời gia tăng tốc độ tích lũy sinh khối gỗ của 8.003 ha rừng tràm thêm $26,25\text{ m}^3/\text{ha}$, tương đương gia tăng hàng trăm nghìn mét khối gỗ dự trữ sinh thái.

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

  • Nghiên cứu sinh và Giới học giả Trẻ: Tiếp cận được một phương pháp luận mẫu mực về giải tích sinh khối thực địa kết hợp mô hình hóa thủy văn Pearson III; kế thừa hệ thống cơ sở dữ liệu đo đạc chi tiết về thổ nhưỡng và thủy văn VQG U Minh Thượng.
  • Các Nhà Khoa học Lâm nghiệp & Môi trường: Sở hữu bộ công cụ phương trình tương quan sinh khối allometric chính xác cho cây tràm tái sinh nhiệt đới, làm nền tảng tính toán trữ lượng carbon rừng ngập nước.
  • Cơ quan Quản lý Vườn Quốc gia & Khu Bảo tồn: Ban Quản lý VQG U Minh Thượng, U Minh Hạ, Tràm Chim, Láng Sen có được cẩm nang quy trình vận hành đóng/mở công trình thủy lợi chuẩn hóa theo từng tuần và từng tháng trong năm.
  • Các Cơ quan Hoạch định Chính sách (Bộ NN&PTNT, Tổng cục Thủy lợi, Chi cục Kiểm lâm): Cơ sở dữ liệu khoa học chuẩn xác để xây dựng các định mức kỹ thuật, quy chuẩn quốc gia (QCVN) về phòng cháy chữa cháy rừng kết hợp bảo tồn đất ngập nước than bùn.
  • Cộng đồng Cư dân Vùng đệm: Hưởng lợi trực tiếp từ nguồn nước ngọt bổ sung trong mùa khô, giảm thiểu thiệt hại do hạn mặn và nâng cao năng suất canh tác lúa - tôm bền vững.

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 độc đáo nhất là việc luận án đã bác bỏ định kiến quản lý giữ nước tĩnh sâu quanh năm để chống cháy, đồng thời mở rộng Lý thuyết Sinh thái Thủy văn (Eco-hydrology Theory) thông qua việc chứng minh quy luật "Nhịp điệu ngập nước mùa" là động lực bắt buộc cho sự phát triển cực đỉnh của rừng tràm than bùn. Tác giả đã mô hình hóa thành công mối quan hệ phi tuyến giữa cao trình mực nước khống chế và sinh khối tích lũy, xác lập ngưỡng rút nước sinh thái $3,5 - 6\text{ tháng/năm}$ để tái cấu trúc hệ rễ cọc và tối đa hóa trữ lượng gỗ.

2. Đổi mới phương pháp luận của công trình so với các nghiên cứu trước đây?
So với nghiên cứu của Phạm Trọng Thịnh (2005) (chỉ khảo sát tĩnh về thổ nhưỡng than bùn) và Hoàng Văn Thắng (2004) (chỉ phân loại định tính sinh cảnh đất ngập nước), luận án của Phạm Văn Tùng tạo bước đột phá phương pháp luận khi tích hợp tam giác nghiên cứu: Số hóa địa hình DEM chi tiết $\rightarrow$ Đo đạc giải tích 16 cây tiêu chuẩn xác định sinh khối thực tế $\rightarrow$ Mô hình hóa cân bằng ẩm đa biến theo phân phối xác suất Pearson III. Đặc biệt, phương pháp phân vùng thủy lực đa bậc độc lập đã thay thế hoàn toàn cách tiếp cận quản lý cào bằng đơn vùng trước đó.

3. Phát hiện bất ngờ nhất được hỗ trợ bởi dữ liệu thực nghiệm là gì?
Phát hiện bất ngờ nhất là việc ngập nước sâu liên tục không hề giúp bảo vệ rừng tràm mà ngược lại làm suy giảm $83,9%$ tốc độ tăng trưởng sinh khối khô ($TSK_k$ chỉ đạt $24,80\text{ tấn/ha}$ so với $45,62\text{ tấn/ha}$ ở vùng bán ngập tối ưu, $F = 42,6; p < 0,001$) và làm gia tăng tỷ lệ đổ ngã cơ học lên $18,4%$ do thối rễ cọc. Thêm vào đó, việc rút nước có kiểm soát trong mùa khô không làm gia tăng nguy cơ cháy nếu duy trì mực nước ngầm mao dẫn trong tầng than bùn cách mặt đất không quá 30–40 cm.

4. Luận án có cung cấp quy trình nhân bản (Replication Protocol) hoàn chỉnh không?
Có. Luận án cung cấp đầy đủ và chi tiết toàn bộ quy trình nhân bản gồm: (1) Bản đồ tọa độ hệ thống ô tiêu chuẩn GPS (Phụ lục 17); (2) Bộ phương trình hồi quy tương quan sinh khối allometric giữa đường kính $D_{1,3}$ và các thành phần sinh khối tươi/khô; (3) Bảng thông số kích thước và bản vẽ kỹ thuật 3D hệ thống cống tiêu, cống lấy nước và đập tràn điều tiết; (4) Lịch biểu tích nước chi tiết theo 3 kịch bản tần suất mưa $P = 25%, 50%, 75%$.

5. Chương trình nghiên cứu 10 năm được phác thảo ra sao?
Chương trình nghiên cứu dài hạn được định hình tập trung vào 3 trục chính: (1) Xây dựng hệ thống giám sát tự động IoT quan trắc mực nước ngầm và độ ẩm than bùn truyền dữ liệu thời gian thực; (2) Thiết lập mô hình toán thủy văn - chất lượng nước 3D dự báo biến đổi tầng phèn hoạt động ($B_j$) và tầng sinh phèn ($C_p$) dưới tác động của biến đổi khí hậu; (3) Lượng giá kinh tế dịch vụ hệ sinh thái và nghiên cứu cơ chế bồi hoàn carbon của rừng tràm U Minh Thượng phục vụ thị trường carbon quốc tế.

Kết luận

  1. Luận án đã giải quyết trọn vẹn và triệt để bài toán hóc búa về mâu thuẫn giữa phòng chống cháy rừng than bùn và bảo tồn sinh thái tại VQG U Minh Thượng sau thảm họa cháy lịch sử năm 2002.
  2. Xác định chính xác chế độ nước tối ưu sinh thái cho rừng tràm tái sinh với thời gian rút nước mặt từ 3,5 đến 6 tháng/năm, giúp nâng tổng sinh khối rừng từ $24,80\text{ tấn/ha}$ lên $45,62\text{ tấn/ha}$ và trữ lượng gỗ tăng từ $32,15\text{ m}^3/\text{ha}$ lên $58,40\text{ m}^3/\text{ha}$ ($p < 0,001$).
  3. Phát triển thành công giải pháp phân vùng thủy lực đa bậc không gian, chia tách vùng lõi 8.003 ha thành 3 phân khu A, B, C độc lập để điều tiết nước theo địa hình, khôi phục $93%$ diện tích ngập theo mùa tại Khu A và B.
  4. Đưa ra lịch biểu đóng cống tích nước chuẩn xác theo xác suất mưa thiết kế Pearson III: ngày 11/9 ($P = 75%$), ngày 01/10 ($P = 50%$), và ngày 21/10 ($P = 25%$), bảo đảm cân đối hoàn hảo lượng hao hụt bốc thoát hơi $2,77\text{ mm/ngày}$.
  5. Thiết kế hoàn chỉnh hệ thống công trình thủy lợi sinh thái gồm cống hộp tiêu thoát lũ, trạm bơm cơ động và đập tràn tự do với đầy đủ thông số thủy lực và bản vẽ kỹ thuật 3D trực quan.
  6. Mở ra 3 hướng nghiên cứu học thuật tiên phong về thủy văn than bùn nứt nẻ, cân bằng phát thải khí nhà kính đầm lầy nhiệt đới, và chuyển đổi số trong quản lý sinh thái đất ngập nước phục vụ bảo tồn bền vững tài nguyên thiên nhiên quốc gia.