Tổng quan về luận án

Sự biến đổi sử dụng đất và thảm phủ (Land Use/Land Cover - LULC) dưới áp lực đô thị hóa, chuyển dịch cơ cấu nông - lâm nghiệp kết hợp cùng biến đổi khí hậu (BĐKH) toàn cầu đang làm thay đổi căn bản cơ chế hình thành dòng chảy trên các lưu vực sông nhiệt đới gió mùa. Lưu vực sông Cả — một trong những hệ thống lưu vực xuyên biên giới lớn nhất miền Trung Việt Nam với tổng diện tích tự nhiên $27.200\text{ km}^2$ (trong đó $9.470\text{ km}^2$, tương đương $34,8%$ diện tích thượng lưu nằm trên lãnh thổ Lào và $17.730\text{ km}^2$ thuộc Việt Nam) — đang đối mặt với nguy cơ an ninh nguồn nước nghiêm trọng. Dự báo đến năm 2030, tổng nhu cầu nước cho phát triển kinh tế - xã hội và dân sinh tại Việt Nam sẽ chạm ngưỡng 122 tỷ $\text{m}^3/\text{năm}$ (tăng 1,5 lần so với hiện nay), trong khi tổng lượng nước mặt có xu thế suy giảm còn khoảng $96%$ so với thời kỳ nền.

                                KHUNG NGHIÊN CỨU TÍCH HỢP ĐA PHƯƠNG PHÁP
                                
   +---------------------------------------------------------------------------------------------------+
   |                                     NGUỒN DỮ LIỆU ĐẦU VÀO                                         |
   | - Viễn thám: Landsat 8 OLI, Sentinel-2 (10-30m)     - Khí tượng - Thủy văn (5 trạm thực đo)       |
   | - Địa hình: DEM 30m, Thổ nhưỡng, HTSDĐ 2005-2015    - Kịch bản BĐKH: RCP4.5 & RCP8.5 (2020-2039)  |
   +-------------------------------------------------+-------------------------------------------------+
                                                     |
                                                     v
   +-------------------------------------------------+-------------------------------------------------+
   |               MODULE 1: MÔ PHỎNG LULC           |            MODULE 2: MÔ HÌNH SWAT               |
   | - Ma trận xác suất chuyển đổi chuỗi Markov     | - Phân chia 29 tiểu lưu vực                     |
   | - Hệ tự hành tế bào Cellular Automata (CA)      | - Phân tích đơn vị phản ứng thủy văn (HRUs)     |
   | - Trọng số MCE/WLC chuẩn hóa qua Delphi-KAMET   | - Hiệu chỉnh & kiểm định thông số (NSE, R2)     |
   +-------------------------------------------------+-------------------------------------------------+
                                                     |
                                                     v
   +---------------------------------------------------------------------------------------------------+
   |                              DỰ TÍNH ĐỘNG THÁI DÒNG CHẢY 2030                                     |
   |     Lũ cực đoan gia tăng (+)  <--------------------->  Kiệt sâu, hạn hán mùa cạn trầm trọng (-)  |
   +---------------------------------------------------------------------------------------------------+

Khoảng trống nghiên cứu (research gap) cốt lõi xuất phát từ thực tiễn: các nghiên cứu thủy văn trước đây thường cô lập tác động của BĐKH hoặc đánh giá thay đổi thảm phủ tĩnh mang tính định tính, thiếu sự liên kết động lực học không gian - thời gian trên quy mô toàn lưu vực liên quốc gia, nơi dữ liệu thảm phủ phía thượng nguồn ngoài lãnh thổ thường bị gián đoạn. Luận án giải quyết hai câu hỏi nghiên cứu (Research Questions):

  • RQ1: Phương pháp luận nào cho phép mô phỏng không gian và lượng hóa độ bất định trong dự tính biến đổi thảm phủ tương lai cho một lưu vực sông xuyên biên giới thiếu hụt dữ liệu?
  • RQ2: Dòng chảy mặt trên lưu vực sông Cả biến động ra sao dưới tác động đồng thời của biến động LULC và kịch bản BĐKH trong giai đoạn đến năm 2030?

Tương ứng với hai giả thuyết nghiên cứu (Hypotheses):

  • H1: Động lực thảm phủ tương lai có thể lượng hóa chính xác bằng mô hình xác suất không gian tích hợp Chuỗi Markov, Hệ tự hành dạng tế bào (Cellular Automata - CA) và kỹ thuật đồng thuận chuyên gia Delphi theo quy tắc KAMET.
  • H2: Tác động đồng thời của chuyển đổi thảm phủ và BĐKH tạo ra hiệu ứng cộng hưởng tiêu cực, làm gia tăng biên độ dao động dòng chảy mùa lũ và làm suy kiệt dòng chảy mùa cạn nghiêm trọng hơn so với kịch bản chỉ xét riêng rẽ yếu tố khí hậu.

Luận án thiết lập khung lý thuyết liên ngành kết hợp giữa Lý thuyết Hệ phức tạp (Complex Systems Theory), Đánh giá đa tiêu chí (Multi-Criteria Evaluation - MCE) và Mô hình thủy văn bán phân bố SWAT (Soil and Water Assessment Tool), mang lại đột phá phương pháp luận và cung cấp cơ sở dữ liệu định lượng cho công tác quản trị lưu vực sông Cả trong thời kỳ biến đổi môi trường toàn cầu.

Literature Review và Positioning

Nghiên cứu về tương tác giữa thảm phủ, khí hậu và tài nguyên nước trên thế giới được cấu thành từ hai dòng học thuật chính:

  1. Dòng mô hình hóa động thái thảm phủ: Bắt nguồn từ nền tảng giải thuật chuỗi Markov (Markov, 1906) nhằm xây dựng ma trận xác suất chuyển đổi trạng thái đất đai, kết hợp với lý thuyết Hệ tự hành dạng tế bào (Cellular Automata - Wolfram, 1984) để phân bổ không gian. Các công trình của Sayemuzzaman & Jha (2014) tại North Carolina (Hoa Kỳ) và Vázquez-Quintero et al. (2016) tại Durango (Mexico) đã chứng minh tính ưu việt của tích hợp Markov-CA trong việc dự báo suy giảm diện tích rừng thông và mở rộng đất xây dựng với hệ số thỏa thuận Kappa đạt trên $90%$. Tương tự, Isma’il & Abubakar (2015) ứng dụng Markov-CA tại Wudil (Nigeria), chỉ ra xác suất đất trống chuyển đổi thành đất xây dựng lên tới $80,3%$, trong khi thảm thực vật chỉ có $41,8%$ khả năng bảo tồn trạng thái.
  2. Dòng đánh giá phản ứng thủy văn lưu vực: Tồn tại tranh luận học thuật sâu sắc giữa các nhà nghiên cứu về vai trò điều tiết của thảm phủ rừng. Một trường phái (Bosch & Hewlett, 1982; Farley et al., 2005; Ge Sun, 2006) khẳng định việc tăng diện tích rừng trồng làm suy giảm từ $20 - 40%$ tổng lượng dòng chảy năm và làm giảm dòng chảy kiệt do gia tăng thoát hơi nước tầng tán (canopy evapotranspiration) từ $50 - 300\text{ mm/năm}$. Ngược lại, trường phái thủy văn truyền thống và sinh thái học cảnh quan (Costa et al., 2003; Zhang et al., 2007) lập luận rằng rừng tự nhiên duy trì cấu trúc xốp của đất, cải thiện hệ số thấm và bổ sung tầng chứa nước ngầm, giúp duy trì dòng chảy kiệt và hạ thấp đỉnh lũ. Cụ thể, Costa et al. (2003) phân tích chuỗi số liệu 50 năm trên lưu vực sông Tocantins (Brazil - $175.360\text{ km}^2$) đã phát hiện khi diện tích nông nghiệp tăng từ $30%$ lên $49%$, lưu lượng trung bình năm tăng tới $24%$ mặc dù lượng mưa không biến đổi có ý nghĩa thống kê.
                  BẢN ĐỒ ĐỊNH VỊ HỌC THUẬT (ACADEMIC POSITIONING)
                  
   Dự tính LULC            |
   Động lực học            |   Sayemuzzaman & Jha (2014)     LUẬN ÁN NGUYỄN THANH BẰNG (2023)
   (CA-Markov + Viễn thám) |   Vázquez-Quintero (2016)       [Tích hợp Viễn thám đa nguồn +
                           |                                  Markov-CA-Delphi + Bán phân bố SWAT
                           |                                  trên lưu vực liên quốc gia]
                           |
   Dự tính LULC Tĩnh/      |   Trần Thục & cs (2003)         Hoàng Minh Tuyển (2006)
   Biến đổi Khí hậu rời rạc|   Ngô Đình Quế (2005)           Nguyễn Ý Như & cs (2009)
                           +------------------------------------------------------------------>
                               Thủy văn Quy mô Cục bộ        Thủy văn Lưu vực Xuyên biên giới

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, hay Nguyễn Ý Như & Nguyễn Thanh Sơn (2009) trên sông Bến Hải đã bước đầu ứng dụng SWAT để mô phỏng tác động của rừng. Tuy nhiên, phần lớn các công trình này chỉ thiết lập kịch bản LULC giả định tĩnh hoặc nghiên cứu quy mô nhỏ. Đóng góp vượt trội của luận án nằm ở vị trí tiên phong tích hợp liên hoàn dữ liệu viễn thám đa thời gian (Landsat, Sentinel-2), mô hình hóa chuyển đổi LULC động bằng Markov-CA hiệu chỉnh qua kỹ thuật Delphi-KAMET, và đồng bộ với các kịch bản phát thải đại diện (RCP4.5, RCP8.5) của IPCC trên toàn bộ lưu vực xuyên biên giới sông 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 phong phú hệ thống lý thuyết thủy văn môi trường và địa tin học thông qua các đóng góp cụ thể:

  1. Bổ sung Lý thuyết Hệ phức tạp và Tự tổ chức Không gian: Minh chứng rằng quá trình biến động thảm phủ vùng nhiệt đới không tuân theo các hàm tuyến tính đơn giản mà là kết quả của quá trình tự hành phi tuyến tính, được chi phối bởi ma trận xác suất chuyển đổi Markov và các điều kiện biên lân cận (neighborhood configuration) của mô hình Cellular Automata.
  2. Lý thuyết Cân bằng Thủy văn Động: Phát triển luận điểm bảo vệ cốt lõi: "Tác động đồng thời của sự thay đổi thảm phủ và biến đổi khí hậu sẽ làm thay đổi dòng chảy lưu vực sông Cả trong tương lai và có xu hướng khắc nghiệt hơn so với chỉ chịu tác động riêng rẽ của biến đổi khí hậu". Công trình cung cấp bằng chứng thực nghiệm khẳng định sự suy thoái thảm phủ rừng đầu nguồn làm triệt tiêu khả năng đệm (buffering capacity) của lưu vực, khuếch đại tính cực đoan của dòng chảy do BĐKH gây ra.
  3. Chuẩn hóa Lý thuyết Quyết định Đa tiêu chí trong Địa không gian: Tích hợp quy tắc tri thức chuyên gia KAMET (Knowledge Acquisition for Multiple Experts with Time scale) vào hàm mờ (Fuzzy membership standardization), giải quyết triệt để tính chủ quan trong gán trọng số không gian cho các nhân tố động lực biến động đất đai.

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

Khung phân tích của luận án được xây dựng dựa trên sự liên kết giữa 3 khối lý thuyết: Lý thuyết Viễn thám Nhận dạng Quang phổ, Lý thuyết Động lực học Không gian Rời rạc và Lý thuyết Thủy văn Bán phân bố HRU.

                           KHUNG PHÂN TÍCH LIÊN KẾT 3 TẦNG
                           
   +---------------------------------------------------------------------------------------+
   | TẦNG 1: KHÔNG GIAN ĐỊA LÝ & VIỄN THÁM (Landsat 8, Sentinel-2, DEM, Thổ nhưỡng)        |
   | - Phân loại ảnh có kiểm định (MLC) -> 5 lớp thảm phủ chính                            |
   | - Kiểm định Kappa, Producer's Accuracy (PA), User's Accuracy (UA)                     |
   +-------------------------------------------+-------------------------------------------+
                                               |
                                               v
   +---------------------------------------------------------------------------------------+
   | TẦNG 2: MÔ HÌNH HÓA ĐỘNG LỰC THẢM PHỦ (CA - Markov - Delphi KAMET)                   |
   | - Ma trận xác suất chuyển đổi P(i,j) & Ma trận diện tích A(i,j)                      |
   | - Bản đồ tính phù hợp (Suitability Maps) qua đánh giá đa tiêu chí MCE/WLC              |
   | - Mô phỏng dự tính kịch bản LULC 2030                                                 |
   +-------------------------------------------+-------------------------------------------+
                                               |
                                               v
   +---------------------------------------------------------------------------------------+
   | TẦNG 3: MÔ HÌNH HÓA THỦY VĂN & ĐÁNH GIÁ TÁC ĐỘNG (SWAT Engine)                       |
   | - Chia 29 tiểu lưu vực, chồng xếp LULC 2030 + Thổ nhưỡng + Độ dốc -> HRUs             |
   | - Tích hợp Kịch bản Khí hậu RCP4.5 & RCP8.5 (Nhiệt độ, Lượng mưa 2020-2039)           |
   | - Định lượng biến động dòng chảy: Dòng chảy Năm, Dòng chảy Mùa lũ, Dòng chảy Kiệt     |
   +---------------------------------------------------------------------------------------+

Định nghĩa phân loại 5 lớp thảm phủ chủ đạo có tác động quyết định đến dòng chảy bề mặt bao gồm: Đất rừng (rừng tự nhiên, rừng phòng hộ, rừng sản xuất), Đất nông nghiệp (lúa nước, hoa màu, cây lâu năm), Đất xây dựng (đô thị, điểm dân cư, công trình), Vùng nước (sông suối, hồ chứa), và Đất trống (đất đồi núi trọc, bãi bồi). Điều kiện biên phân tích được xác lập chặt chẽ: lưu vực sông Cả với ranh giới tự nhiên bao trùm $9.470\text{ km}^2$ thuộc lãnh thổ Lào và chuỗi thời gian nền tảng giai đoạn 2005 - 2015 phục vụ huấn luyện, hướng tới mốc dự báo 2030 (chu kỳ 2020 - 2039).

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ủ thế giới quan thực chứng luận kết hợp chủ nghĩa hiện thực phản biện (Positivism & Critical Realism), vận dụng phương pháp tiếp cận đa ngành, phân tích nhân - quả và mô hình hóa bán phân bố đa cấp độ. Thiết kế nghiên cứu xử lý bài toán mô phỏng không gian đa thời gian trên diện tích $27.200\text{ km}^2$ với độ phân giải không gian $30\text{ m}$.

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu được thực hiện qua các giai đoạn chuẩn hóa:

  1. Xử lý dữ liệu viễn thám và phân loại LULC: Thu thập tư liệu vệ tinh Landsat TM/ETM+/OLI và Sentinel-2 các thời kỳ 2005, 2010, 2015. Áp dụng thuật toán phân loại giám sát Xác suất Cực đại (Maximum Likelihood Classification - MLC) kết hợp giải đoán chỉ số thực vật sai khác (NDVI) và mô hình số độ cao (DEM Aster $30\text{ m}$). Khảo sát thực địa thu thập tọa độ mẫu phục vụ huấn luyện và giải đoán.
  2. Mô hình hóa CA-Markov kết hợp Delphi: Thiết lập ma trận xác suất chuyển đổi $P$ và ma trận diện tích chuyển đổi $S$ theo công thức: $$P = \begin{bmatrix} P_{11} & P_{12} & \dots & P_{1n} \ P_{21} & P_{22} & \dots & P_{2n} \ \vdots & \vdots & \ddots & \vdots \ P_{n1} & P_{n2} & \dots & P_{nn} \end{bmatrix} \quad \left(0 \le P_{ij} \le 1, , \sum_{j=1}^n P_{ij} = 1\right)$$ Áp dụng kỹ thuật Delphi qua 3 vòng tham vấn chuyên gia liên ngành theo quy tắc KAMET để sàng lọc các nhân tố động lực (khoảng cách tới đường giao thông, sông suối, trung tâm đô thị, độ dốc, cao độ). Chuẩn hóa các nhân tố bằng hàm mờ (Fuzzy logic) 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) nhằm tạo bản đồ mức độ phù hợp (Suitability Maps). Chạy lặp mô hình CA không gian 2 chiều để mô phỏng bản đồ thảm phủ năm 2015, kiểm định độ chính xác trước khi dự tính thảm phủ năm 2030.
  3. Mô phỏng Thủy văn SWAT: Phân chia lưu vực sông Cả thành 29 tiểu lưu vực (Sub-basins). Tiếp tục phân tách thành hàng trăm Đơn vị Phản ứng Thủy văn (Hydrologic Response Units - HRUs) thông qua việc chồng xếp các lớp bản đồ LULC, bản đồ thổ nhưỡng (đất feralit, phù sa...) và phân cấp độ dốc.
             MA TRẬN HIỆU CHỈNH & KIỂM ĐỊNH MÔ HÌNH (CALIBRATION & VALIDATION)
             
   +----------------------------------+----------------------------------+
   |   PHÂN LOẠI & DỰ TÍNH LULC       |     MÔ PHỎNG THỦY VĂN SWAT       |
   +----------------------------------+----------------------------------+
   | - Chỉ số Kappa tổng thể: > 0,85  | - Trạm Nghĩa Khánh, Hòa Duyệt,   |
   | - Producer's Accuracy: > 85%     |   Dừa, Quỳ Châu, Sơn Diệm        |
   | - User's Accuracy: > 85%         | - Hiệu chỉnh (2007-2010):        |
   | - Kiểm định chéo với hiện trạng   |   NSE: 0,76 - 0,88 | R2: 0,78 - 0,89 |
   |   thực địa và QHSDĐ 2020         | - Kiểm định (2011-2014):         |
   |                                  |   NSE: 0,72 - 0,85 | R2: 0,75 - 0,87 |
   +----------------------------------+----------------------------------+
   |   ĐỘ TIN CẬY: ĐẠT MỨC TỐT ĐẾN RẤT TỐT (EXCELLENT HYDROLOGICAL FIT)  |
   +---------------------------------------------------------------------+

Độ tin cậy của mô phỏng thủy văn được lượng hóa qua Hệ số Hiệu quả Nash-Sutcliffe ($NSE$) và Hệ số Xác định ($R^2$): $$NSE = 1 - \frac{\sum_{t=1}^T (Q_{o}^t - Q_{m}^t)^2}{\sum_{t=1}^T (Q_{o}^t - \overline{Q_o})^2}; \quad R^2 = \frac{\left[\sum (Q_o - \overline{Q_o})(Q_m - \overline{Q_m})\right]^2}{\sum (Q_o - \overline{Q_o})^2 \sum (Q_m - \overline{Q_m})^2}$$ Kết quả hiệu chỉnh (2007 - 2010) và kiểm định (2011 - 2014) tại 5 trạm thủy văn chủ lực (Hòa Duyệt, Nghĩa Khánh, Dừa, Quỳ Châu, Sơn Diệm) đều đạt $NSE \ge 0,75$ và $R^2 \ge 0,80$, chứng minh mô hình đủ năng lực mô phỏng quá trình thủy văn thực tế.

Data và phân tích

Nguồn cơ sở dữ liệu đầu vào bao gồm:

  • Mô hình số độ cao (DEM) độ phân giải $30\text{ m}$ từ USGS/NASA.
  • Bản đồ hiện trạng và quy hoạch sử dụng đất tỉnh Nghệ An, Hà Tĩnh năm 2010, 2020.
  • Số liệu quan trắc khí tượng ngày (lượng mưa, nhiệt độ tối cao, nhiệt độ tối thấp, độ ẩm, bức xạ, gió) tại 9 trạm đại diện lưu vực giai đoạn 1980 - 2015.
  • Dữ liệu lưu lượng dòng chảy thực đo ngày tại 5 trạm khống chế thủy văn.
  • Kịch bản BĐKH phiên bản cập nhật của Bộ Tài nguyên và Môi trường (kịch bản RCP4.5 và RCP8.5 giai đoạn 2020 - 2039).

Phương pháp nghiên cứu khắc phục triệt để sai số thường gặp trong các nghiên cứu trước (vốn có sai lệch $PAD = 34,25%$ do không đồng bộ bản đồ), đảm bảo tính nhất quán topology và tính chuẩn xác của các ma trận chuyển đổi.

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

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

  1. Động lực chuyển dịch LULC giai đoạn 2005 - 2030: Kết quả dự tính chỉ ra sự phân hóa không gian mạnh mẽ. Đất xây dựng tại lưu vực sông Cả có tốc độ tăng trưởng nhanh nhất, mở rộng tập trung quanh vùng đồng bằng hạ lưu (TP. Vinh, Nghi Lộc, Đức Thọ). Đất nông nghiệp có xu hướng thu hẹp do đô thị hóa và chuyển đổi mục đích sử dụng. Đất rừng tự nhiên tại khu vực miền núi phía Tây Nghệ An có xu hướng bảo tồn hoặc tăng nhẹ nhờ các chính sách lâm nghiệp, nhưng chất lượng thảm phủ và độ giàu sinh khối có dấu hiệu suy thoái phân tán.
  2. Hiệu ứng cộng hưởng làm gia tăng đỉnh lũ mùa mưa: Theo kịch bản kết hợp biến động LULC năm 2030 và kịch bản BĐKH RCP8.5, lưu lượng dòng chảy đỉnh lũ trung bình tăng từ $12,5%$ đến $28,4%$ tại các tiểu lưu vực thượng nguồn và trung lưu (các trạm Nghĩa Khánh, Dừa, Sơn Diệm). Sự gia tăng diện tích mặt đệm không thấm (đất xây dựng) cùng với suy giảm thảm phủ tự nhiên làm rút ngắn thời gian tập trung dòng chảy ($T_c$), đẩy nhanh quá trình tập trung lũ, dẫn tới nguy cơ ngập lụt gia tăng đột biến ở hạ du.
  3. Khủng hoảng suy kiệt dòng chảy mùa cạn: Dưới tác động đồng thời của gia tăng nhiệt độ (làm tăng bốc thoát hơi nước tiềm năng $ET_0$) và sự suy giảm lượng mưa mùa khô theo kịch bản RCP4.5/RCP8.5 kết hợp với biến đổi LULC, lưu lượng dòng chảy mùa cạn (tháng 1 đến tháng 4) tại các trạm giảm sâu từ $8,6%$ đến $22,3%$. Hiện tượng "lũ tăng, kiệt giảm" trở nên gay gắt hơn so với kịch bản chỉ xét riêng rẽ yếu tố BĐKH từ $3,5 - 6,2%$.
  4. Tính phi đồng nhất không gian giữa các lưu vực nhánh: Lưu vực sông Ngàn Sâu, Ngàn Phố (trạm Hòa Duyệt, Sơn Diệm) có mức độ nhạy cảm thủy văn cao nhất trước biến động thảm phủ, với biên độ dao động dòng chảy năm vượt trội so với lưu vực nhánh sông Hiếu (trạm Quỳ Châu).
                      BIẾN ĐỘNG DÒNG CHẢY THEO KỊCH BẢN TƯƠNG LAI
                      
   Lưu lượng (Q)
        ^
        |                                  Kịch bản Tương lai (LULC 2030 + RCP8.5)
        |                 /\               ---------------------------------------
        |                /  \  <--- Lũ tăng mạnh (+12.5% đến +28.4%)
        |    /\         /    \             Thời kỳ nền (Baseline)
        |   /  \       /      \            .......................................
        |  /    \_____/        \
        | /                     \_________
        |/                                \___  <--- Kiệt sâu (-8.6% đến -22.3%)
        +----------------------------------------> Thời gian (Tháng trong năm)
          Tháng 1       Tháng 8-10       Tháng 12

Implications đa chiều

  • Về mặt học thuật và lý thuyết: Công trình xác lập bằng chứng thực nghiệm vững chắc về quy luật tương tác ghép bộ (coupled interaction) giữa biến động mặt đệm cục bộ và biến đổi khí hậu quy mô khu vực, mở rộng hiểu biết về thủy văn học lưu vực xuyên biên giới tại Đông Nam Á.
  • Về mặt phương pháp luận: Cung cấp quy trình tích hợp mô hình chuẩn hóa (Remote Sensing -> Markov-CA-Delphi -> SWAT) có thể nhân rộng trực tiếp cho các lưu vực sông tương tự tại khu vực miền Trung (sông Mã, sông Thu Bồn, sông Ba) và các quốc gia đang phát triển.
  • Về ứng dụng thực tiễn và chính sách:
    • Cảnh báo các nhà quy hoạch tỉnh Nghệ An và Hà Tĩnh về nguy cơ vỡ cấu trúc cân bằng nước nếu tiếp tục chuyển đổi đất rừng đầu nguồn sang đất phát triển cơ sở hạ tầng thiếu kiểm soát.
    • Cung cấp luận cứ khoa học để điều chỉnh quy trình vận hành liên hồ chứa trên lưu vực sông Cả (Bản Vẽ, Khe Bố, Chi Khê), yêu cầu nâng dung tích phòng lũ mùa mưa và gia tăng lưu lượng xả duy trì dòng chảy môi trường trong mùa kiệt.
    • Đặt nền móng cho cơ chế hợp tác song phương Việt - Lào trong việc chia sẻ dữ liệu thủy văn và đồng quản lý bảo vệ thảm phủ rừng xuyên biên giới tại tỉnh Xiêng Khoảng.

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 3 giới hạn nghiên cứu chính:

  1. Độ phân giải và mật độ quan trắc thượng nguồn: Phần diện tích $34,8%$ thuộc lãnh thổ Lào thiếu mạng lưới trạm quan trắc khí tượng thủy văn mặt đất dày đặc; số liệu khí tượng chủ yếu dựa trên tái phân tích và nội suy không gian, tạo ra độ bất định nhất định đối với các tiểu lưu vực vùng biên.
  2. Giả định tĩnh về công trình nhân tạo: Mô hình chưa mô phỏng chi tiết quy trình vận hành đóng - mở thời gian thực của toàn bộ hệ thống hồ đập thủy điện, công trình điều tiết thủy lợi vừa và nhỏ xây dựng mới trong tương lai.
  3. Chưa tích hợp mô hình kinh tế - tác tử (Agent-Based Models): Mô hình CA-Markov phản ánh quy luật chuyển đổi không gian theo chuỗi lịch sử nhưng chưa phản ánh đầy đủ các cú sốc chính sách bất thường hoặc biến động thị trường đất đai vi mô.
                            CHƯƠNG TRÌNH NGHIÊN CỨU 10 NĂM TỚI
                            
   [2024 - 2026]              [2027 - 2029]              [2030 - 2033+]
   Tích hợp CMIP6/SSPs        Mô hình Đa tác tử (ABM)    Bản sao số Lưu vực (Digital Twin)
   + Downscaling độ phân      + Tương tác Hành vi nông   + Cập nhật Viễn thám Real-time
     giải cao 1-5 km            hộ & Thị trường đất      + Tối ưu hóa Hồ chứa bằng AI/ML

Các hướng nghiên cứu tiếp nối bao gồm:

  • Nâng cấp kịch bản BĐKH sang thế hệ mô hình CMIP6 với các kịch bản kinh tế - xã hội Shared Socioeconomic Pathways (SSPs).
  • Tích hợp mô hình kinh tế tác tử (Agent-Based Modeling) để lượng hóa hành vi ra quyết định chuyển đổi đất đai của cộng đồng dân cư địa phương.
  • Phát triển hệ thống bản sao số (Digital Twin) lưu vực sông Cả, kết nối dữ liệu viễn thám vi sóng (SAR Sentinel-1) thời gian thực với mô hình thủy lực 2 chiều mô phỏng ngập lụt chi tiết.

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

  • Tác động học thuật (Academic Impact): Tạo lập một khung phương pháp luận tiêu chuẩn, có khả năng trích dẫn cao trên các tạp chí quốc tế chuyên ngành (như Journal of Hydrology, Science of The Total Environment, Land Use Policy), đóng góp bộ dữ liệu tham chiếu chuẩn cho các nghiên cứu thủy văn xuyên biên giới Đông Nam Á.
  • Tác động hoạch định chính sách (Policy Influence): Trực tiếp hỗ trợ 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, Hà Tĩnh trong việc lập Quy hoạch tài nguyên nước thời kỳ 2021 - 2030, tầm nhìn 2050 và Kế hoạch hành động ứng phó BĐKH cấp tỉnh.
  • Lợi ích kinh tế - xã hội và môi trường: Bảo vệ sinh kế bền vững cho hơn 3,5 triệu dân cư trên lưu vực; tối ưu hóa hiệu quả phát điện và cấp nước tưới tiêu; giảm thiểu thiệt hại kinh tế hàng trăm tỷ đồng mỗi năm do thiên tai lũ lụt và hạn hán cực đoan gây ra.
  • Ý nghĩa quốc tế: Tăng cường quan hệ hợp tác bảo vệ môi trường xuyên biên giới Việt Nam - Lào, đóng góp vào các mục tiêu phát triển bền vững của Liên Hợp Quốc (SDG 6: Nước sạch và Vệ sinh, SDG 13: Hành động vì Khí hậu, SDG 15: Tài nguyên trên Cạn).

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

                                  ĐỐI TƯỢNG HƯỞNG LỢI CHÍNH
                                  
   +------------------------------+-------------------------------+------------------------------+
   |   NGHIÊN CỨU SINH & HỌC THUẬT|   CƠ QUAN QUẢN LÝ NHÀ NƯỚC   |   KHU VỰC KINH TẾ & DOANH NGHIỆP|
   | - Kế thừa khung Markov-CA-   | - Bộ TN&MT, Cục Quản lý TNN   | - Các nhà máy thủy điện      |
   |   Delphi-SWAT chuẩn tắc      | - UBND tỉnh Nghệ An & Hà Tĩnh | - Công ty Khai thác Thủy lợi |
   | - Bộ dữ liệu hiệu chỉnh chuẩn| - Quy hoạch không gian & rừng | - Nông nghiệp công nghệ cao  |
   +------------------------------+-------------------------------+------------------------------+
  • Nghiên cứu sinh và Giới học thuật: Kế thừa bộ tham số mô hình SWAT đã được hiệu chuẩn ($NSE > 0,75$), tiếp cận khung quy trình mô phỏng thảm phủ kết hợp Delphi-KAMET để phát triển các đề tài chuyên sâu.
  • Nhà hoạch định chính sách và Quản lý Lưu vực: Sở TN&MT, Ủy ban Lưu vực sông Cả sở hữu công cụ định lượng khoa học để cấp phép khai thác tài nguyên nước, khoanh vùng bảo vệ hành lang nguồn nước và phân bổ hạn ngạch sử dụng nước mùa khô.
  • Đơn vị Quản lý Rừng và Doanh nghiệp Thủy điện/Thủy lợi: Ban Quản lý rừng phòng hộ, Công ty Thủy điện Bản Vẽ có cơ sở dữ liệu chính xác để xác định lưu vực trọng điểm chi trả dịch vụ môi trường rừng (PES) và lập lịch tích - xả nước tối ưu.
  • Cộng đồng dân cư địa phương: Hưởng lợi từ hệ thống cảnh báo sớm rủi ro thiên tai và nguồn nước sinh hoạt, sản xuất ổn định, giảm thiểu tổn thương sinh kế.

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 mở rộng Lý thuyết Động lực học Không gian Phức tạp (Complex Spatial Dynamics) và Lý thuyết Cân bằng Nước Lưu vực thông qua việc chứng minh cơ chế tương tác phi tuyến tính giữa quá trình chuyển đổi mặt đệm cục bộ và biến đổi khí hậu toàn cầu. Luận án đã phá vỡ giả định truyền thống xem thảm phủ là yếu tố tĩnh trong các kịch bản khí hậu dài hạn, cung cấp công thức tích hợp đa nhân tố lượng hóa được độ trễ và tính khuếch đại của dòng chảy mặt.

2. Đột phá phương pháp luận so với các nghiên cứu tiền nhiệm quốc tế và trong nước là gì?

So với các nghiên cứu sử dụng Markov-CA đơn thuần (như Sayemuzzaman & Jha, 2014) hoặc các mô hình nội suy tĩnh tại Việt Nam (Nguyễn Ý Như, 2009; Trần Thục, 2003), luận án tạo đột phá bằng việc thiết lập quy trình khép kín: Viễn thám đa nguồn (Landsat 8 + Sentinel-2) $\rightarrow$ Mô hình Markov $\rightarrow$ Hệ tự hành CA $\rightarrow$ Kỹ thuật Delphi theo quy tắc KAMET $\rightarrow$ Mô hình bán phân bố SWAT. Giải pháp này loại bỏ sai số không gian (giảm tỷ lệ sai lệch $PAD$ từ mức $34,25%$ của các nghiên cứu trước xuống dưới $10%$), xử lý thành công vùng dữ liệu trống xuyên biên giới.

   TIÊU CHÍ PHƯƠNG PHÁP       CÁC NGHIÊN CỨU TRƯỚC ĐÂY            LUẬN ÁN NGUYỄN THANH BẰNG (2023)
   -----------------------    ---------------------------------   ----------------------------------
   Xử lý biến động LULC       Giả định kịch bản tĩnh / nội suy    Mô hình động CA-Markov + Viễn thám
   Xác định trọng số không gian Gán chủ quan hoặc AHP đơn giản    Delphi 3 vòng chuẩn hóa KAMET-Fuzzy
   Quy mô không gian          Lưu vực nội địa / Tiểu vùng nhỏ     Toàn lưu vực xuyên biên giới (27.200 km2)
   Tích hợp kịch bản          Tách rời LULC hoặc chỉ xét BĐKH     Ghép bộ đồng thời LULC 2030 + RCP4.5/8.5

3. Phát hiện bất ngờ và 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 là: Mức độ suy kiệt dòng chảy mùa cạn trong kịch bản kết hợp LULC và BĐKH lớn hơn đáng kể so với tổng đại số của hai tác động riêng rẽ. Hiện tượng này xảy ra do sự suy giảm thảm thực vật tầng cao làm suy thoái cấu trúc tầng ngậm nước mặt, kết hợp với nhiệt độ tăng cao làm gia tăng đột biến bốc thoát hơi nước thực tế ($ET$), triệt tiêu khả năng nạp nước ngầm duy trì dòng chảy kiệt, dẫn đến suy giảm dòng chảy mùa kiệt lên tới $22,3%$.

4. Luận án có cung cấp đầy đủ giao thức để tái lặp nghiên cứu (Replication Protocol) không?

Có. Luận án cung cấp tường minh toàn bộ quy trình: từ bảng danh mục thông số SWAT được hiệu chỉnh ($CN2$, $SOL_AWC$, $ALPHA_BF$, $ESCO...$), ma trận xác suất chuyển đổi Markov, bảng câu hỏi tham vấn Delphi 3 vòng kèm ngưỡng đồng thuận KAMET, đến danh mục phân loại 5 lớp LULC và tiêu chuẩn kiểm định thống kê ($NSE$, $R^2$, $Kappa$), cho phép cộng đồng khoa học tái lập hoàn toàn trên các lưu vực khác.

5. Chương trình nghiên cứu 10 năm tới được phác thảo ra sao?

Chương trình nghiên cứu mở rộng bao gồm:

  • Tích hợp mô hình khí hậu CMIP6 độ phân giải siêu cao ($1 - 5\text{ km}$).
  • Ứng dụng Trí tuệ nhân tạo (AI/Deep Learning - mạng LSTM/Transformer) trong dự báo dòng chảy tức thời kết hợp mô hình vật lý SWAT.
  • Xây dựng mô hình tương tác Thủy văn - Xã hội học (Socio-hydrology) để lượng hóa phản hồi qua lại giữa chính sách sử dụng đất và nguồn tài nguyên nước.

Kết luận

Luận án tiến sĩ của NCS. Nguyễn Thanh Bằng đã hoàn thành xuất sắc các mục tiêu khoa học đề ra, mang lại 5 đóng góp cốt lõi:

  1. Xây dựng và kiểm định thành công khung phương pháp luận tích hợp Chuỗi Markov, Hệ tự hành dạng tế bào (CA) và kỹ thuật Delphi theo quy tắc KAMET để dự tính không gian thảm phủ tương lai với độ chính xác tin cậy ($Kappa > 0,85$).
  2. Thiết lập bộ thông số mô hình thủy văn bán phân bố SWAT tối ưu cho toàn bộ $27.200\text{ km}^2$ lưu vực sông Cả, đạt chỉ số kiểm định thủy văn xuất sắc ($NSE$ từ $0,72 - 0,88$; $R^2$ từ $0,75 - 0,89$).
  3. Định lượng hóa thành công tác động đồng thời của biến động thảm phủ năm 2030 và kịch bản BĐKH (RCP4.5, RCP8.5) đến chế độ dòng chảy: làm gia tăng lưu lượng lũ mùa mưa từ $12,5 - 28,4%$ và làm suy giảm dòng chảy mùa cạn từ $8,6 - 22,3%$.
  4. Khẳng định quy luật khuếch đại rủi ro thủy văn khi kết hợp suy thoái mặt đệm cục bộ với BĐKH toàn cầu, cung cấp cơ sở khoa học để định hình lại chiến lược an ninh nguồn nước vùng Bắc Trung Bộ.
  5. Đề xuất hệ thống giải pháp ứng phó đồng bộ từ quy hoạch không gian đất đai, bảo vệ rừng đầu nguồn xuyên biên giới Việt - Lào, đến vận hành liên hồ chứa thông minh.

Công trình đánh dấu bước tiến quan trọng trong chuyên ngành Biến đổi khí hậu và Thủy văn môi trường, để lại di sản khoa học thực chứng có giá trị ứng dụng lâu dài cho công tác quản trị tài nguyên và phát triển bền vững lưu vực sông tại Việt Nam.