Tổng quan về luận án

Sự biến động phức tạp của các chu trình thủy văn dưới tác động kép của nhân sinh học và biến đổi khí hậu (BĐKH) toàn cầu đang đặt ra những thách thức nghiêm trọng đối với an ninh nguồn nước tại các lưu vực sông xuyên biên giới. Trong bối cảnh đó, công trình nghiên cứu của Nghiên cứu sinh Nguyễn Thanh Bằng (2022) với đề tài "Nghiên cứu đánh giá tác động của thay đổi thảm phủ và biến đổi khí hậu đến dòng chảy trên lưu vực sông Cả" (Chuyên ngành: Biến đổi khí hậu, Mã số: 9440221, Viện Khoa học Khí tượng Thủy văn và Biến đổi khí hậu) mang tính tiên phong trong việc thiết lập khung phân tích tích hợp mô hình không gian - thời gian phục vụ quản trị tài nguyên nước bền vững.

Khoảng trống nghiên cứu (Research Gap) cốt lõi mà luận án tập trung giải quyết xuất phát từ thực tiễn: hầu hết các công trình nghiên cứu thủy văn trước đây tại Việt Nam thường đánh giá tác động đơn lẻ của BĐKH hoặc giả định thảm phủ tĩnh/thay đổi mang tính định tính mà chưa định lượng hóa được sự biến đổi thảm phủ động trong tương lai xuyên suốt toàn lưu vực, đặc biệt là phần diện tích ngoài biên giới quốc gia. Trích dẫn nguyên văn từ nguồn tư liệu luận án: "Sông Cả dài 513 km với diện tích lưu vực 27.200 km² – một trong những lưu vực sông dài và lớn ở Việt Nam... trong đó 9470 km² ở thượng lưu (chiếm 34,8%) nằm trong lãnh thổ Lào". Sự thiếu hụt số liệu đo đạc xuyên biên giới tại vùng thượng nguồn tỉnh Xiêng Khoảng (Lào) đã tạo ra rào cản lớn cho việc mô phỏng chính xác chế độ dòng chảy hạ du thuộc hai tỉnh Nghệ An và Hà Tĩnh.

Luận án xác lập 02 câu hỏi nghiên cứu (Research Questions) và 02 giả thuyết khoa học (Hypotheses) tương ứng:

  • Câu hỏi nghiên cứu 1 (RQ1): Phương pháp phân tích chuỗi Markov kết hợp Hệ tự hành dạng tế bào (Cellular Automata - CA) có khả năng mô phỏng và dự tính chính xác không gian - thời gian sự biến động thảm phủ trên một lưu vực xuyên biên giới phức tạp như lưu vực sông Cả hay không?
  • Câu hỏi nghiên cứu 2 (RQ2): Chế độ dòng chảy mặt (dòng chảy năm, mùa lũ, mùa cạn và dòng chảy tháng) trên lưu vực sông Cả sẽ biến đổi như thế nào dưới tác động đồng thời của biến động thảm phủ năm 2030 và các kịch bản BĐKH (RCP4.5 và RCP8.5) giai đoạn 2020–2039?
  • Giả thuyết khoa học 1 (H1): Việc tích hợp chuỗi Markov, Cellular Automata cùng các ràng buộc không gian xác định qua phương pháp Delphi/AHP cho phép tái tạo ma trận chuyển đổi và mô phỏng chính xác phân bố 05 lớp thảm phủ chính với độ tin cậy thống kê cao (Hệ số Kappa > 0,85).
  • Giả thuyết khoa học 2 (H2): Tác động hiệp đồng giữa suy giảm diện tích rừng/gia tăng đất xây dựng và biến động cực đoan của nhiệt độ - lượng mưa theo các kịch bản RCP sẽ làm gia tăng tính cực đoan của dòng chảy: khuếch đại đỉnh lũ mùa mưa và làm trầm trọng thêm tình trạng suy kiệt dòng chảy mùa khô.

Khung lý thuyết của luận án được xây dựng dựa trên sự giao thoa giữa Lý thuyết Hệ thống Phức tạp (Complex Systems Theory), Quá trình Ngẫu nhiên Markovian (Markovian Stochastic Process) và Lý thuyết Cân bằng Thủy văn Vật lý Bán phân bố (Semi-distributed Physical Hydrological Theory). Phạm vi không gian nghiên cứu bao trùm 27.200 km² lưu vực, với chuỗi dữ liệu viễn thám nền tảng 2005–2015, dữ liệu khí tượng thủy văn đo đạc và kịch bản dự tính đến năm 2030 (thời kỳ 2020–2039). Công trình giải quyết bài toán nhu cầu nước dự báo tăng gấp 1,5 lần (đạt khoảng 122 tỷ m³/năm vào năm 2030 trên toàn quốc), đóng góp cơ sở dữ liệu quan trọng cho các quy hoạch phát triển kinh tế - xã hội vùng Bắc Trung Bộ.

Literature Review và Positioning

Tổng quan y văn quốc tế và trong nước về động lực thảm phủ và phản ứng thủy văn cho thấy sự phát triển của ba trường phái học thuật chủ đạo:

Trường phái thứ nhất tập trung vào kỹ thuật mô hình hóa biến động sử dụng đất/thảm phủ (LULC). Xuất phát từ nền tảng toán học của Andrei Markov (1906) và lý thuyết hệ phức tạp của Stephen Wolfram (1980s), sự kết hợp giữa chuỗi Markov và Cellular Automata (CA-Markov) đã trở thành chuẩn mực quốc tế. Griselda Vázquez-Quintero và cộng sự (2016) đã ứng dụng thành công CA-Markov tại vùng rừng ôn đới Pueblo Nuevo (Durango, Mexico) đạt chỉ số Kappa từ 90% đến 92%, chứng minh xu hướng suy giảm liên tục của diện tích rừng thông giai đoạn 2014–2028. Tương tự, Sayemuzzaman và Manoj K. Jha (2014) tích hợp Đánh giá Đa tiêu chí (MCE) và CA-Markov tại bang North Carolina (Hoa Kỳ), dự báo đến năm 2030 diện tích đất xây dựng tăng 20% trong khi đất lâm nghiệp giảm 17% so với năm 1992. Tại châu Á và châu Phi, các công trình của Arak (Iran, 2015) và Wudil, Kano (Nigeria - Muhammad Isma’il & Saidu Abubakar, 2015) cũng khẳng định tính ưu việt của CA-Markov khi kết hợp dữ liệu Landsat/GIS để mô phỏng mở rộng đô thị với độ chính xác trên 80–95%.

Trường phái thứ hai khảo sát mối quan hệ định lượng giữa biến động thảm phủ và chế độ dòng chảy. Luận án trích xuất tổng hợp kinh điển từ y văn: "Đánh giá 94 thí nghiệm lưu vực cho thấy sự thay đổi hướng dòng chảy sau các hoạt động lâm nghiệp có thể được dự đoán với độ chính xác hợp lý vì 93/94 kết quả thí nghiệm đều chỉ ra rằng việc giảm sản lượng nước liên quan tới việc giảm độ che phủ hoặc tăng sản lượng nước liên quan tới sự gia tăng của độ che phủ" (Bosch & Hewlett, 1982; Bosch J., 1982). Tuy nhiên, tồn tại một cuộc tranh luận học thuật sâu sắc (academic debate) giữa hai quan điểm đối nghịch:

  • Quan điểm "Hiệu ứng miếng bọt biển" (Sponge effect hypothesis): Cho rằng tán rừng và lớp thảm mục làm tăng khả năng thấm, bổ cập nước ngầm và duy trì dòng chảy mùa kiệt (Nisbet, 2001; Vũ Văn Tuấn & Phạm Thị Hương Lan, 2003 tại lưu vực Mù Căng Chải).
  • Quan điểm "Bốc thoát hơi nước vượt trội" (Evapotranspiration trade-off hypothesis): Các nghiên cứu thực nghiệm của Farley và cộng sự (2005) chỉ ra khi chuyển cỏ/cây bụi sang rừng trồng, dòng chảy năm giảm từ 31% đến 44% do lượng bốc thoát hơi nước tăng mạnh. Ge Sun và cộng sự (2006) chứng minh trồng rừng tại Trung Quốc làm giảm dòng chảy từ 50 mm/năm (vùng khô) đến 300 mm/năm (vùng ẩm), tương đương mức sụt giảm 20–40% tổng lượng nước năm.

So sánh với các nghiên cứu quốc tế điển hình, nghiên cứu của Costa và cộng sự (2003) trên lưu vực sông nhiệt đới Tocantins (Porto Nacional, Brazil - 175.360 km²) cho thấy khi diện tích rừng bị chuyển đổi 30% sang đất nông nghiệp thì lưu lượng trung bình năm tăng 24% dù lượng mưa không đổi về mặt thống kê. Nghiên cứu mô hình Markov-Kalman tại lưu vực XiuHe (Trung Quốc, 2010) cũng khẳng định vai trò chi phối của chuyển dịch nông nghiệp tới thủy văn bề mặt.

Tại Việt Nam, các nghiên cứu của Trần Thục & Huỳnh Thị Lan Hương (2003 trên sông Trà Khúc), Ngô Đình Quế (2005 trên sông Rào Nậy và Pơ Cô), Nguyễn Ý Như & Nguyễn Thanh Sơn (2009 trên sông Bến Hải), Nguyễn Thị Tịnh Ấu và cs (2014 trên sông Đắk Bla), Nguyễn Thị Hồng (2014 trên sông Vu Gia) đã ứng dụng mô hình SWAT. Tuy nhiên, điểm hạn chế của các công trình trong nước là chỉ phân tích kịch bản thảm phủ quá khứ hoặc kịch bản giả định tĩnh. Luận án của Nguyễn Thanh Bằng định vị chính xác vị trí học thuật tiên phong khi thiết lập cầu nối: Tích hợp dữ liệu viễn thám đa thời gian (Landsat, Sentinel-2), tham vấn chuyên gia Delphi/KAMET, mô phỏng không gian CA-Markov cho tương lai 2030 xuyên biên giới, và nạp vào mô hình thủy văn SWAT dưới các kịch bản BĐKH của IPCC AR5 (RCP4.5, RCP8.5).

Đó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 ứng dụng của Lý thuyết Quá trình Ngẫu nhiên Markovian và Lý thuyết Hệ tự hành dạng tế bào trong khoa học thủy văn lưu vực xuyên biên giới. Bằng việc lượng hóa ma trận xác suất chuyển đổi trạng thái và tích hợp hàm lân cận không gian 2 chiều, nghiên cứu đã vượt qua giới hạn phi không gian của chuỗi Markov truyền thống. Về mặt thủy văn học, công trình củng cố Lý thuyết Đơn vị Phản ứng Thủy văn (Hydrologic Response Unit - HRU), chứng minh rằng sự phân mảnh vi mô của thảm phủ bề mặt tạo ra tác động phi tuyến tính lên quá trình hình thành dòng chảy mặt và dòng chảy ngầm sát mặt.

+-----------------------------------------------------------------------------+
|                      KHUNG PHÂN TÍCH TÍCH HỢP ĐA MÔ HÌNH                     |
+-----------------------------------------------------------------------------+
|  [Dữ liệu Viễn thám: Landsat 8 / Sentinel-2 (2005-2010-2015)]               |
|  [Phân loại Thảm phủ MLC -> 05 Lớp: Rừng, Nông nghiệp, Đô thị, Nước, Trống] |
|  [Chuỗi Markov: Ma trận Xác suất & Diện tích Chuyển đổi P(ij)]              |
|  [Phương pháp Delphi / KAMET]                           [Phân tích Đa Tiêu   |
|  (Tham vấn Chuyên gia 3 Vòng)                           Chí MCE / Hàm Mờ]   |
|  [Hệ tự hành dạng Tế bào (Cellular Automata) -> Bản đồ LULC 2030]            |
|  [Kịch bản Biến đổi Khí hậu: RCP4.5 & RCP8.5 (Giai đoạn 2020-2039)]         |
|  [Dữ liệu Địa hình (DEM 30m)]                           [Dữ liệu Thổ nhưỡng |
|                                                         (Bản đồ Đất FAO)]   |
|  [Mô hình Thủy văn Bán phân bố SWAT: 29 Tiểu lưu vực -> Phân tích HRU]     |
|  [Định lượng Tác động Hiệp đồng: Dòng chảy Năm, Mùa Lũ, Mùa Cạn, Tháng]    |
+-----------------------------------------------------------------------------+

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

Khung phân tích của luận án tích hợp chặt chẽ 3 trụ cột lý thuyết và mô hình:

  1. Lý thuyết Viễn thám & Nhận dạng Không gian: Sử dụng thuật toán Phân loại Hợp lý Cực đại (Maximum Likelihood Classification - MLC) để bóc tách 05 lớp phủ chính: Rừng, Đất nông nghiệp, Đất xây dựng, Vùng nước, và Đất trống.
  2. Lý thuyết Quyết định Đa tiêu chí (MCE) & Quy tắc KAMET (Knowledge Acquisition for Multiple Experts with Time scale): Chuẩn hóa các biến nhân tố (khoảng cách tới đường giao thông, sông suối, trung tâm đô thị, độ dốc, độ cao DEM) thông qua hàm mờ (Fuzzy membership) và xác định trọng số bằng phương pháp Kết hợp Tuyến tính có Trọng số (Weighted Linear Combination - WLC) sau 3 vòng tham vấn chuyên gia Delphi.
  3. Mô hình Thủy văn Vật lý Bán phân bố SWAT: Thiết lập liên kết hữu cơ giữa các yếu tố khí tượng (nhiệt độ, lượng mưa), thổ nhưỡng, địa hình và LULC thông qua các phương trình bảo toàn khối lượng và cân bằng nước lưu vực:

$$SW_t = SW_0 + \sum_{i=1}^{t} (R_{day} - Q_{surf} - E_a - w_{seep} - Q_{gw})$$

Điều kiện biên (Boundary Conditions): Khung phân tích được xác lập tối ưu cho vùng khí hậu nhiệt đới gió mùa Bắc Trung Bộ chịu ảnh hưởng của gió mùa Tây Nam khô nóng (gió Phơn) và bão/áp thấp nhiệt đới tập trung từ tháng 8 đến tháng 11, với cấu trúc lưu vực kéo dài từ độ cao trên 2.000m tại Lào xuống đồng bằng ven biển Cửa Hội.

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

Thiết kế nghiên cứu

Luận án vận dụng thế giới quan thực chứng (Positivism) kết hợp phương pháp mô hình hóa tất định - ngẫu nhiên. Thiết kế nghiên cứu đa cấp độ (Multi-level hierarchical design):

  • Cấp toàn lưu vực (Macro-level): Lưu vực sông Cả (27.200 km²) được số hóa và thiết lập mạng lưới quan trắc khí tượng thủy văn chuẩn mực.
  • Cấp tiểu lưu vực (Meso-level): Toàn bộ lưu vực được phân chia thành 29 tiểu lưu vực (sub-basins) dựa trên dữ liệu mô hình số độ cao (DEM 30m) từ USGS/NASA.
  • Cấp đơn vị phản ứng thủy văn (Micro-level): Sự giao thoa đồng nhất giữa 05 lớp sử dụng đất, các loại thổ nhưỡng chính (đất feralit, phù sa, xám xói mòn...) và 05 cấp độ dốc hình thành nên hàng trăm HRU đặc thù.
+-------------------------------------------------------------------------+
|                    QUY TRÌNH THỰC THI NGHIÊN CỨU                        |
+-------------------------------------------------------------------------+
|  GIAI ĐOẠN 1: TIỀN XỬ LÝ & GIẢI ĐOÁN VIỄN THÁM                          |
|  - Thu thập ảnh Landsat 8 (15-30m) & Sentinel-2 (10-20m) 2005, 2010, 2015|
|  - Hiệu chỉnh khí quyển, tính toán chỉ số thực vật NDVI, NDWI, NDBI     |
|  - Phân loại MLC -> Đánh giá độ chính xác qua Ma trận Nhầm lẫn          |
|  GIAI ĐOẠN 2: MÔ PHỎNG THẢM PHỦ CA-MARKOV 2030                          |
|  - Xây dựng ma trận xác suất chuyển đổi P(ij) & ma trận diện tích       |
|  - Tham vấn Delphi 3 vòng theo quy tắc KAMET xác định yếu tố động lực  |
|  - Chuẩn hóa Fuzzy & lập bản đồ phù hợp bằng WLC                       |
|  - Chạy mô phỏng CA-Markov -> Kiểm định độ chính xác bằng hệ số Kappa   |
|  GIAI ĐOẠN 3: HIỆU CHỈNH & KIỂM ĐỊNH MÔ HÌNH THỦY VĂN SWAT              |
|  - Nhập DEM, Bản đồ đất, Bản đồ LULC mô phỏng, Chuỗi số liệu KTTV       |
|  - Phân tích độ nhạy thông số (CN2, ALPHA_BF, GW_DELAY, SOL_AWC...)     |
|  - Hiệu chỉnh (2007-2010) & Kiểm định (2011-2014) tại 06 trạm thủy văn  |
|  - Đánh giá chỉ số: Hệ số Nash-Sutcliffe (NSE), R², PBIAS               |
|  GIAI ĐOẠN 4: ĐÁNH GIÁ TÁC ĐỘNG TƯƠNG LAI (2020-2039)                   |
|  - Tích hợp Kịch bản LULC 2030 x Kịch bản BĐKH (RCP4.5 & RCP8.5)        |
|  - Mô phỏng biến đổi dòng chảy: Dòng chảy năm, mùa lũ, mùa cạn, tháng   |
+-------------------------------------------------------------------------+

Quy trình nghiên cứu rigorous

Quy trình tam giác đạc phương pháp và dữ liệu (Data & Method Triangulation) được kiểm soát nghiêm ngặt:

  • Dữ liệu thực địa: Thu thập hàng trăm điểm GPS kiểm chứng thực địa trên các sinh cảnh điển hình tại Nghệ An và Hà Tĩnh nhằm đối soát kết quả giải đoán ảnh.
  • Kỹ thuật Delphi/KAMET: Thực hiện qua 3 vòng độc lập với sự tham gia của các chuyên gia đầu ngành trong lĩnh vực viễn thám, lâm nghiệp, khí tượng và quy hoạch đất đai để triệt tiêu thiên kiến chủ quan trong việc chọn trọng số nhân tố.
  • Kiểm tra độ tin cậy và giá trị (Validity & Reliability): Độ chính xác phân loại viễn thám được kiểm định qua Ma trận Nhầm lẫn (Error Matrix) với các chỉ số Độ chính xác Người sản xuất (Producer Accuracy - PA), Độ chính xác Người dùng (User Accuracy - UA) và Độ chính xác Tổng thể (Overall Classification Accuracy - CA) đều vượt ngưỡng 85%.

Data và phân tích

Hệ thống thông số hiệu chỉnh mô hình thủy văn SWAT được tối ưu hóa qua thuật toán phân tích độ nhạy. Các chỉ số thống kê đánh giá hiệu quả mô hình tại 06 trạm thủy văn quan trắc đại diện:

  • Trạm Yên Thượng (Dòng chính sông Cả): Giai đoạn hiệu chỉnh (2007–2010) đạt $NSE \ge 0,82$, $R^2 \ge 0,85$; Giai đoạn kiểm định (2011–2014) đạt $NSE \ge 0,80$, $R^2 \ge 0,83$.
  • Trạm Hòa Duyệt (Sông La): $NSE$ dao động từ $0,78$ đến $0,84$, chỉ số tương quan $R^2 > 0,81$.
  • Các trạm phụ lưu (Dừa, Quỳ Châu, Sơn Diệm, Nghĩa Khánh): Tất cả đều đáp ứng tiêu chuẩn khắt khe của Moriasi et al., với $NSE > 0,75$, $R^2 > 0,78$, chỉ số sai số tương đối $PBIAS < \pm 15%$, khẳng định mô hình SWAT đã phản ánh xuất sắc cơ chế thủy văn thực tế của lưu vực.

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

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

1. Xu thế chuyển dịch thảm phủ sâu sắc giai đoạn 2005–2030: Phân tích ma trận chuyển đổi diện tích xác nhận sự thu hẹp liên tục của đất rừng tự nhiên và sự bùng nổ của đất xây dựng, nông nghiệp. Trích dẫn nguyên văn dữ liệu lịch sử từ luận án: "Năm 1943 có khoảng 1,2 triệu ha rừng (~67% diện tích lưu vực trong lãnh thổ VN), đến nay diện tích rừng chiếm khoảng 40% diện tích lưu vực, so với diện tích đất của các huyện miền núi và Hương Khê, Hương Sơn thì diện tích đất có rừng chiếm đến 43%". Đến năm 2030, mô hình CA-Markov dự báo diện tích đất xây dựng tiếp tục tăng trưởng mạnh tại khu vực đô thị hóa dọc hạ lưu (Vinh, Cửa Lò, Nghi Lộc) và các trục hành lang kinh tế, trong khi diện tích đất rừng phòng hộ đầu nguồn phía Việt Nam đối mặt áp lực suy thoái cục bộ nếu không có giải pháp bảo vệ nghiêm ngặt.

2. Sự biến đổi cực đoan của các yếu tố khí hậu tương lai (2020–2039): Dưới kịch bản RCP4.5 và RCP8.5, nhiệt độ trung bình năm trên toàn lưu vực tăng từ 1,1°C đến 1,8°C so với thời kỳ nền. Đáng chú ý, lượng mưa có sự phân hóa mùa sâu sắc: lượng mưa mùa mưa (tháng 8–11) tăng từ 8,5% đến 18,2%, trong khi lượng mưa mùa khô (tháng 12 đến tháng 5 năm sau) giảm từ 6,4% đến 15,3%, làm gia tăng đáng kể biên độ dao động ẩm của lưu vực.

+-------------------------------------------------------------------------------------+
|         BẢNG TỔNG HỢP BIẾN ĐỔI DÒNG CHẢY GIAI ĐOẠN 2020-2039 SO VỚI THỜI KỲ NỀN      |
+----------------------+--------------------+--------------------+--------------------+
| Kịch bản Mô phỏng    | Dòng chảy Năm (%)  | Dòng chảy Mùa Lũ   | Dòng chảy Mùa Cạn  |
|                      |                    | (Tháng 8-11) (%)   | (Tháng 12-5) (%)   |
+----------------------+--------------------+--------------------+--------------------+
| Kịch bản BĐKH Đơn lẻ | +3,5% đến +6,8%    | +8,2% đến +14,5%   | -7,5% đến -13,2%   |
| (RCP4.5 / RCP8.5)    |                    |                    |                    |
+----------------------+--------------------+--------------------+--------------------+
| Kịch bản Tích hợp:   | +5,8% đến +11,4%   | +14,8% đến +24,6%  | -12,4% đến -21,8%  |
| LULC 2030 + RCP8.5   |                    |                    |                    |
+----------------------+--------------------+--------------------+--------------------+

3. Tác động khuếch đại dòng chảy đỉnh lũ (Flood Amplification): Khi tích hợp đồng thời biến động thảm phủ 2030 và kịch bản BĐKH cực đoan (RCP8.5), lưu lượng đỉnh lũ tại các trạm Yên Thượng và Hòa Duyệt tăng đột biến từ 18,5% đến 27,3%. Sự gia tăng lớp phủ không thấm (bề mặt bê tông hóa, đô thị hóa) kết hợp với giảm mật độ che phủ rừng ở các tiểu lưu vực miền núi làm giảm thời gian tập trung dòng chảy ($T_c$), gia tăng hệ số dòng chảy mặt ($C_{runoff}$), khiến lũ lên nhanh hơn và đạt đỉnh cao hơn.

4. Suy kiệt nghiêm trọng dòng chảy mùa cạn (Dry-Season Baseflow Depletion): Phát hiện mang tính cảnh báo cao nhất của luận án là sự sụt giảm dòng chảy kiệt trong các tháng 3, 4 và 5. Lưu lượng trung bình tháng mùa kiệt tại các trạm chính giảm từ 15,2% đến 26,4% dưới kịch bản kết hợp. Hiện tượng này bắt nguồn từ việc bốc thoát hơi nước tiềm năng ($PET$) gia tăng do nhiệt độ tăng cao, kết hợp với suy giảm khả năng trữ ẩm và bổ cập tầng chứa nước ngầm do diện tích rừng bị chia cắt.

Implications đa chiều

+------------------------------------------------------------------------------------+
|                   HỆ THỐNG HÀM Ý VÀ ỨNG DỤNG ĐA CHIỀU                              |
+------------------------------------------------------------------------------------+
| 1. HÀM Ý LÝ THUYẾT & HỌC THUẬT                                                     |
| - Hoàn thiện khung phân tích thủy văn không gian cho lưu vực sông xuyên biên giới. |
| - Làm sáng tỏ cơ chế phản ứng phi tuyến tính giữa thảm phủ và khí hậu cực đoan.   |
|                                                                                    |
| 2. HÀM Ý QUẢN TRỊ TÀI NGUYÊN NƯỚC (IWRM)                                           |
| - Điều chỉnh quy trình vận hành liên hồ chứa (Bản Vẽ, Hủa Na, Khe Bố...) trên      |
|   lưu vực sông Cả nhằm ứng phó với lũ sớm và hạn hán kéo dài.                      |
| - Quy hoạch hệ thống công trình trữ nước phân tán tại các tiểu lưu vực nhạy cảm.   |
|                                                                                    |
| 3. CHÍNH SÁCH QUY HOẠCH SỬ DỤNG ĐẤT & LÂM NGHIỆP                                   |
| - Thiết lập hành lang bảo vệ rừng phòng hộ đầu nguồn nghiêm ngặt tại miền Tây      |
|   Nghệ An (Kỳ Sơn, Tương Dương, Con Cuông) và Hà Tĩnh (Hương Sơn, Hương Khê).      |
| - Tăng cường tỷ lệ cây bản địa có độ tán che dày, hạn chế trồng rừng sản xuất độc   |
|   canh keo/bạch đàn chu kỳ ngắn vốn làm giảm dòng chảy kiệt.                       |
|                                                                                    |
| 4. NGOẠI GIAO NGUỒN NƯỚC & HỢP TÁC QUỐC TẾ                                        |
| - Cung cấp cơ sở khoa học để Việt Nam và Lào xây dựng hiệp định chia sẻ dữ liệu     |
|   thủy văn và phối hợp bảo vệ rừng đầu nguồn tại tỉnh Xiêng Khoảng.                |
+------------------------------------------------------------------------------------+

Limitations và Future Research

Mặc dù đạt được những kết quả xuất sắc, luận án thẳng thắn thừa nhận 04 giới hạn nghiên cứu (Limitations):

  1. Độ phân giải và tính liên tục của dữ liệu xuyên biên giới: Số liệu quan trắc khí tượng thủy văn bề mặt tại vùng thượng nguồn phía Lào (chiếm 34,8% diện tích lưu vực) còn rất thưa thớt, phụ thuộc chủ yếu vào dữ liệu viễn thám và tái phân tích toàn cầu, dẫn đến sai số tiềm ẩn trong mô hình hóa quá trình mưa - dòng chảy cục bộ phía thượng lưu.
  2. Bất định trong hạ quy mô kịch bản BĐKH (Downscaling Uncertainty): Việc sử dụng các mô hình khí hậu khu vực (RCM) để hạ quy mô từ các GCM toàn cầu vẫn mang độ bất định nhất định đối với các hiện tượng thời tiết cực đoan quy mô vi khí hậu (như mưa cục bộ do địa hình chắn gió dãy Trường Sơn).
  3. Phân loại thảm phủ dừng lại ở 05 lớp khái quát: Mô hình chưa bóc tách chi tiết các cấu trúc rừng nội vi (rừng giàu, rừng nghèo, rừng hỗn giao, rừng keo lai chu kỳ ngắn) vốn có các hệ số thoát hơi nước và giữ nước tầng đất rất khác nhau.
  4. Giả định tĩnh về hành vi con người: Mô hình CA-Markov chưa tích hợp đầy đủ các kịch bản can thiệp chính sách đột xuất, biến động kinh tế vĩ mô hoặc các dự án công trình thủy điện - thủy lợi mới có thể phát sinh trong tương lai.

Chương trình nghiên cứu tương lai (Future Research Agenda):

  • Hướng 1: Tích hợp dữ liệu viễn thám độ phân giải siêu cao (Radar Sentinel-1 SAR, ảnh vệ tinh VNREDSat-1) để theo dõi động học độ ẩm đất và cấu trúc tầng tán rừng theo thời gian thực.
  • Hướng 2: Ứng dụng các thuật toán máy học/học sâu (Machine Learning/Deep Learning như LSTM, Random Forest) kết hợp mô hình vật lý thủy văn (Physically-based AI) để nâng cao độ chính xác dự báo dòng chảy thời gian thực.
  • Hướng 3: Mở rộng đánh giá các tác động thứ cấp: xói mòn đất, vận chuyển bùn cát, bồi lắng hồ chứa thủy điện và xâm nhập mặn vùng cửa sông Lam dưới tác động của nước biển dâng.
  • Hướng 4: Phát triển khung hỗ trợ ra quyết định liên kết kinh tế - sinh thái - thủy văn (Hydro-Economic Modeling) phục vụ chi trả dịch vụ môi trường rừng (PES) xuyên biên giới.

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

  • Tác động học thuật (Academic Impact): Thiết lập một quy trình chuẩn mực (Benchmark protocol) về phương pháp luận tích hợp CA-Markov và SWAT cho các lưu vực sông có tính chất xuyên biên giới tại Đông Nam Á. Công trình mở ra tiềm năng trích dẫn cao trong các nghiên cứu chuyên ngành biến đổi khí hậu, thủy văn môi trường và địa tin học.
  • Chuyển đổi thực tiễn và công nghiệp (Industry & Practical Transformation): Cung cấp công cụ tính toán chính xác cho các tập đoàn năng lượng, công ty thủy điện (Thủy điện Bản Vẽ, Hủa Na, Nậm Nơn...) trong việc tối ưu hóa quy trình điều tiết hồ chứa, vừa gia tăng hiệu suất phát điện vừa đảm bảo xả dòng chảy môi trường phục vụ hạ du.
  • Tác động chính sách công (Policy Influence): Là tài liệu khoa học nền tảng phục vụ Sở Tài nguyên và Môi trường, Sở Nông nghiệp và Phát triển Nông thôn hai tỉnh Nghệ An và Hà Tĩnh trong việc lập "Quy hoạch tài nguyên nước lưu vực sông Cả thời kỳ 2021–2030, tầm nhìn 2050" và xây dựng Kế hoạch hành động ứng phó BĐKH cấp tỉnh.
  • Lợi ích xã hội (Societal Benefits): Giúp giảm thiểu thiệt hại kinh tế và bảo vệ tính mạng cho hơn 3,5 triệu người dân sinh sống tại lưu vực sông Cả trước các thảm họa thiên tai lũ quét, sạt lở đất trong mùa mưa và hạn hán, xâm nhập mặn khốc liệt trong mùa khô.

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

+-------------------------------------------------------------------------------------+
| Đối tượng Hưởng lợi        | Giá trị Đóng góp và Ứng dụng Cụ thể                    |
+----------------------------+--------------------------------------------------------+
| Nghiên cứu sinh &          | Tiếp cận bộ quy tắc hiệu chỉnh thông số SWAT, phương    |
| Giới Viện - Trường         | pháp luận Delphi/KAMET và ma trận chuyển đổi CA-Markov.|
+----------------------------+--------------------------------------------------------+
| Cơ quan Quản lý Nhà nước   | Căn cứ khoa học phục vụ phê duyệt quy hoạch sử dụng    |
| (Bộ TN&MT, UBND Tỉnh)      | đất và cấp phép phân bổ hạn ngạch khai thác tài nguyên |
|                            | nước mặt bền vững.                                     |
+----------------------------+--------------------------------------------------------+
| Ban Quản lý Lưu vực Sông & | Bộ dữ liệu kịch bản dòng chảy chi tiết đến cấp tiểu    |
| Đơn vị Vận hành Hồ chứa    | lưu vực phục vụ xây dựng biểu đồ điều tiết phòng lũ và |
|                            | cấp nước mùa cạn tối ưu.                               |
+----------------------------+--------------------------------------------------------+
| Cộng đồng Địa phương &     | Nâng cao khả năng thích ứng, chủ động chuyển đổi cơ cấu|
| Doanh nghiệp Nông nghiệp   | cây trồng, mùa vụ thích ứng với tình trạng thiếu nước. |
+----------------------------+--------------------------------------------------------+

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 thiết lập thành công mô hình tích hợp không gian - ngẫu nhiên liên kết giữa Lý thuyết Quá trình Ngẫu nhiên Markovian và Mô hình Cân bằng Nước Vật lý Bán phân bố (SWAT) trong điều kiện lưu vực xuyên biên giới thiếu hụt dữ liệu đo đạc. Nghiên cứu đã mở rộng Lý thuyết Đơn vị Phản ứng Thủy văn (HRU) bằng cách chứng minh rằng: sự bất đồng nhất về mặt không gian của quá trình chuyển đổi thảm phủ (LULC) đóng vai trò điều biến phi tuyến tính đối với độ nhạy của lưu vực trước các cú sốc khí hậu (climate shocks).

2. Điểm cải tiến phương pháp luận cốt lõi so với các nghiên cứu trước đây là gì?

So với nghiên cứu của Sayemuzzaman & Jha (2014) tại Mỹ hay Nguyễn Thị Hồng (2014) trên sông Vu Gia, luận án không chỉ dùng ma trận xác suất chuyển đổi quá khứ đơn thuần mà đã áp dụng Phương pháp Delphi 3 vòng theo quy tắc KAMET để tích hợp tri thức chuyên gia vào việc xác định các ràng buộc và nhân tố động lực không gian (độ dốc, giao thông, khoảng cách đô thị). Đồng thời, so với các nghiên cứu thủy văn thông thường vốn cố định thảm phủ khi chạy kịch bản BĐKH, luận án đã cho phép cả hai động lực (LULC tương lai 2030 và Khí hậu RCPs) cùng biến động đồng thời.

3. Phát hiện bất ngờ hoặc phản trực giác nhất từ dữ liệu thực nghiệm là gì?

Phát hiện phản trực giác (counter-intuitive finding) nằm ở phản ứng của dòng chảy kiệt: mặc dù tổng lượng mưa năm theo các kịch bản RCP có xu hướng tăng nhẹ (3,5% đến 6,8%), nhưng lưu lượng dòng chảy mùa cạn lại giảm rất sâu (lên tới hơn 21%). Điều này chứng minh rằng lượng mưa tăng chỉ tập trung vào một số đợt mưa bão cường suất lớn trong mùa mưa, trong khi sự gia tăng bốc thoát hơi nước ($PET$) do nhiệt độ tăng kết hợp với việc suy giảm rừng tự nhiên đã triệt tiêu hoàn toàn khả năng tích trữ nước của lưu vực trong mùa khô.

4. Luận án có cung cấp quy trình nhân bản (Replication Protocol) rõ ràng không?

Có. Luận án mô tả chi tiết từng bước: từ kỹ thuật tiền xử lý ảnh viễn thám (NDVI, MLC), bảng câu hỏi tham vấn Delphi, hàm mờ tiêu chuẩn hóa không gian, cấu trúc ma trận chuyển đổi, đến bộ 16 thông số hiệu chuẩn độ nhạy của SWAT (CN2, ALPHA_BF, GWQMN, ESCO...) kèm theo các ngưỡng đánh giá thống kê ($NSE$, $R^2$, $PBIAS$). Quy trình này hoàn toàn có thể nhân bản và áp dụng cho tất cả các lưu vực sông tương tự tại khu vực Miền Trung và Tây Nguyên.

5. Chương trình nghiên cứu 10 năm tiếp theo được định hình ra sao?

Chương trình nghiên cứu 10 năm tới hướng tới việc xây dựng "Bản sao Số Thủy văn Lưu vực" (Digital Twin for River Basin Hydrology), tích hợp mô hình SWAT với các hệ thống cảm biến IoT quan trắc tự động và ảnh vệ tinh thời gian thực, đồng thời mở rộng nghiên cứu về an ninh nguồn nước xuyên biên giới Việt - Lào dưới tác động tổng hợp của BĐKH và các dự án thủy điện thượng nguồn.

Kết luận

  1. Luận án đã xây dựng thành công phương pháp luận tích hợp tiên tiến giữa viễn thám, Hệ thông tin địa lý (GIS), chuỗi Markov, Hệ tự hành dạng tế bào (CA) và mô hình thủy văn SWAT bán phân bố để giải quyết bài toán đánh giá tác động kép của thảm phủ và BĐKH trên lưu vực sông Cả.
  2. Đã lượng hóa và dự tính chính xác bản đồ thảm phủ lưu vực sông Cả năm 2030 với 05 lớp phủ chính, đạt hệ số tin cậy thống kê Kappa cao, phản ánh rõ nét xu hướng suy giảm diện tích rừng và gia tăng diện tích đất xây dựng, nông nghiệp.
  3. Hiệu chỉnh và kiểm định thành công bộ thông số mô hình SWAT tại 06 trạm thủy văn chủ lực (Yên Thượng, Hòa Duyệt, Dừa, Quỳ Châu, Sơn Diệm, Nghĩa Khánh) với các chỉ số tối ưu ($NSE > 0,78$, $R^2 > 0,81$), bảo đảm độ tin cậy khoa học tuyệt đối cho các mô phỏng dòng chảy.
  4. Phát hiện quy luật tác động hiệp đồng mang tính cực đoan: sự biến đổi thảm phủ kết hợp kịch bản BĐKH (RCP4.5 và RCP8.5) làm gia tăng mạnh mẽ lưu lượng lũ mùa mưa (tăng tới 24,6%) và làm suy kiệt sâu sắc dòng chảy kiệt mùa cạn (giảm tới 21,8%), đe dọa trực tiếp đến an ninh nguồn nước hạ du.
  5. Đề xuất hệ thống giải pháp công trình và phi công trình mang tính chiến lược, cung cấp luận cứ khoa học nền tảng cho công tác quy hoạch tài nguyên nước, quản lý lưu vực bền vững và thúc đẩy hợp tác quản trị nguồn nước xuyên biên giới Việt Nam - Lào.