Tổng quan về luận án

Sự biến động thảm phủ mặt đất (Land Use/Land Cover - LULC) và biến đổi khí hậu (BĐKH) toàn cầu đang tạo 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 nhiệt đới gió mùa. Khu vực Bắc Trung Bộ Việt Nam, đặc biệt là lưu vực sông Cả, đang chứng kiến tốc độ phát triển kinh tế - xã hội và đô thị hóa nhanh chóng, dẫn đến tình trạng suy thoái thảm phủ rừng, thu hẹp diện tích đất nông nghiệp và gia tăng bề mặt không thấm nước. Lưu vực sông Cả là hệ thống sông liên quốc gia có tổng diện tích tự nhiên $27.200\text{ km}^2$, trong đó $9.470\text{ km}^2$ ở vùng thượng lưu (chiếm 34,8% diện tích) nằm trên lãnh thổ nước Cộng hòa Dân chủ Nhân dân Lào và $17.730\text{ km}^2$ nằm trên lãnh thổ Việt Nam, với tổng chiều dài dòng chính $513\text{ km}$ ($361\text{ km}$ chảy qua Việt Nam). Lượng mưa trung bình năm toàn lưu vực biến động mạnh từ $1.100\text{ mm}$ đến $2.500\text{ mm}$, cục bộ có nơi đạt $2.000 - 2.400\text{ mm}$ (như thượng lưu sông Hiếu, sông La). Theo các nghiên cứu thủy văn quốc gia, dự báo đến năm 2030, nhu cầu nước cho các mục đích kinh tế - xã hội và dân sinh tại Việt Nam sẽ đạt khoả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 rõ rệt.

Khoảng trống nghiên cứu (Research Gap) lớn nhất trước đây là sự thiếu vắng một khung phương pháp luận định lượng tích hợp, có khả năng mô phỏng động lực không gian - thời gian của thảm phủ xuyên biên giới và ghép nối đồng bộ với mô hình thủy văn phân bố/bán phân bố dưới các kịch bản BĐKH cập nhật. Phần lớn các công trình trước đây tại Việt Nam chỉ khảo sát LULC ở trạng thái tĩnh, đánh giá tác động đơn lẻ của BĐKH hoặc LULC một cách tách rời, hoặc bỏ qua phần lưu vực thượng nguồn bên phía Lào do thiếu hụt dữ liệu điều tra truyền thống. Luận án tiến sĩ của tác giả Nguyễn Thanh Bằng (2022) tại Viện Khoa học Khí tượng Thủy văn và Biến đổi khí hậu đã giải quyết triệt để khoảng trống này thông qua hai câu hỏi nghiên cứu và hệ giả thuyết định hướng:

  • Câu hỏi nghiên cứu 1 (RQ1): Áp dụng tích hợp phân tích chuỗi Markov và hệ tự hành dạng tế bào (Cellular Automata - CA) kết hợp đánh giá đa tiêu chí không gian có thể mô phỏng chính xác quy luật biến động và dự tính thảm phủ tương lai cho lưu vực sông Cả xuyên biên giới hay không?
    • Giả thuyết 1 (H1): Quá trình chuyển dịch thảm phủ tuân theo quy luật ngẫu nhiên có điều kiện (Stochastic Process), việc kết hợp ma trận xác suất chuyển đổi Markov với quy tắc lân cận của CA và trọng số không gian tối ưu qua tham vấn chuyên gia Delphi-KAMET sẽ tái tạo được không gian thảm phủ với hệ số tin cậy $\text{Kappa} \ge 0{,}80$.
  • Câu hỏi nghiên cứu 2 (RQ2): Dòng chảy mặt lưu vực sông Cả (theo năm, mùa lũ, mùa cạn và theo tháng) trong giai đoạn 2020–2039 sẽ biến đổi định lượng như thế nào dưới tác động đồng thời của kịch bản biến động thảm phủ 2030 và kịch bản BĐKH (RCP4.5, RCP8.5) so với thời kỳ nền?
    • Giả thuyết 2 (H2): Tác động cộng hưởng giữa suy giảm rừng, mở rộng đất xây dựng và sự gia tăng các cực trị nhiệt độ - lượng mưa sẽ làm biến đổi sâu sắc phân phối dòng chảy theo mùa, gia tăng đáng kể lưu lượng đỉnh lũ và làm trầm trọng thêm tình trạng suy kiệt dòng chảy kiệt mùa khô.

Khung lý thuyết của nghiên cứu được xây dựng trên sự giao thoa giữa: Lý thuyết Hệ thống phức tạp và Tự hành tế bào (Complex Systems & Cellular Automata Theory) của Stephen Wolfram, Lý thuyết Quá trình ngẫu nhiên Markov (Markovian Stochastic Process Theory) của Andrei Andreyevich Markov, và Lý thuyết Cân bằng nước đất lưu vực (Soil Water Balance Theory) được tham số hóa trong mô hình SWAT của Cục Nghiên cứu Nông nghiệp Hoa Kỳ (USDA-ARS). Luận án phân loại và dự báo cấu trúc 05 lớp phủ bề mặt chủ đạo: Rừng, Nông nghiệp, Xây dựng, Vùng nước và Đất trống dựa trên chuỗi dữ liệu ảnh vệ tinh Landsat/Sentinel (giai đoạn nền 2005, 2010, 2015), hiệu chỉnh và kiểm định thủy văn trên chuỗi số liệu thực đo giai đoạn 2007–2014, từ đó dự tính dòng chảy giai đoạn 2020–2039.

Literature Review và Positioning

Tổng quan y văn quốc tế và trong nước cho thấy sự phát triển qua ba dòng nghiên cứu chính:

Dòng thứ nhất tập trung vào kỹ thuật mô hình hóa biến động LULC. Mô hình chuỗi Markov thuần túy được Iacono et al. (2015) áp dụng tại vùng đô thị Minneapolis-St. Paul (Hoa Kỳ) đã chứng minh khả năng dự báo xác suất chuyển đổi tối ưu, nhưng bộc lộ hạn chế cố hữu là không có khả năng định vị không gian (spatially explicit). Để khắc phục, sự kết hợp giữa Cellular Automata và Markov Chain (CA-Markov) đã trở thành bước đột phá. Sayemuzzaman & Jha (2014) tích hợp CA-Markov với Đánh giá đa tiêu chí (MCE) 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% và đất lâm nghiệp giảm 17%. Muhammad Isma’il & Saidu Abubakar (2015) ứng dụng CA-Markov tại Wudil (Nigeria), chỉ ra xác suất đất trống chuyển thành đất xây dựng lên tới 80,3%. Vázquez-Quintero et al. (2016) tại Mexico sử dụng CA-Markov theo dõi rừng thông với độ chính xác phân loại Kappa đạt từ 90% đến 92%. Tại Việt Nam, Võ Quốc Tuấn et al. (2015) khi ứng dụng CA-Markov mô phỏng LULC tại Đồng bằng sông Cửu Long trên nền tảng GAMA đã chỉ ra hạn chế của dữ liệu không đồng nhất: "Đánh giá kết quả mô phỏng sự thay đổi sử dụng đất đến năm 2008 qua chỉ số Tỷ lệ phần trăm độ lệch tuyệt đối (PAD) giữa diện tích từng loại đất mô phỏng đến năm 2008 và diện tích thực của các kiểu sử dụng đất năm 2008 là 34,25% hay nói khác hơn là độ chính xác toàn cục mô phỏng của mô hình là 65,75%".

Dòng thứ hai nghiên cứu phản ứng thủy văn lưu vực đối với thay đổi thảm phủ. Bosch & Hewlett (1982) tổng hợp 94 thí nghiệm lưu vực trên thế giới khẳng định 93/94 trường hợp việc giảm lớp phủ thực vật làm tăng tổng lượng dòng chảy và ngược lại. Costa et al. (2003) phân tích lưu vực sông nhiệt đới Tocantins tại Brazil ($175.360\text{ km}^2$) giai đoạn 1949–1998, phát hiện khi tỷ lệ chuyển đổi đất rừng sang đất nông nghiệp đạt khoảng 30% diện tích lưu vực thì lưu lượng trung bình năm tăng thêm 24% dù lượng mưa không thay đổi đáng kể về mặt thống kê. Ngược lại, Farley et al. (2005) chứng minh việc chuyển đổi trảng cỏ sang rừng trồng làm giảm dòng chảy năm từ 31% đến 44% và làm suy giảm mạnh dòng chảy kiệt. Tại Việt Nam, 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à sông 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 et al. (2014) trên sông Đắk Bla, và Nguyễn Thị Hồng (2014) trên sông Vu Gia đều khẳng định vai trò của rừng trong việc cắt lũ và điều tiết dòng chảy kiệt.

Dòng thứ ba khảo sát tác động của BĐKH đến thủy văn lưu vực, sử dụng các mô hình hoàn lưu chung (GCM) và mô hình khí hậu khu vực (RCM) theo các kịch bản SRES hoặc RCP của Ủy ban Liên chính phủ về Biến đổi Khí hậu (IPCC).

Trong y văn tồn tại hai cuộc tranh luận học thuật sâu sắc:

  1. Tranh luận về tác động của rừng đến dòng chảy kiệt (Sponge versus Evapotranspiration Debate): Quan điểm truyền thống cho rằng rừng đóng vai trò "miếng bọt biển" trữ nước mùa mưa và cấp nước mùa khô. Tuy nhiên, các nghiên cứu thực nghiệm hiện đại (Guardiola et al., 2010; Ge Sun, 2006) chỉ ra rừng trồng bản địa hoặc cây công nghiệp (như cao su) có lượng bốc thoát hơi nước (evapotranspiration - ET) rất lớn, làm sụt giảm dòng chảy kiệt từ 20% đến 40%.
  2. Tranh luận về yếu tố chi phối chính: Liệu sự thay đổi LULC cục bộ hay BĐKH toàn cầu đóng vai trò áp đảo đối với sự biến đổi chế độ thủy văn lưu vực.

Luận án định vị tiên phong bằng cách thiết lập một khung nghiên cứu ghép nối liên hoàn (Coupled Modeling Framework): kết hợp Viễn thám đa thời gian $\rightarrow$ Phân tích Markov $\rightarrow$ Hệ tự hành CA tối ưu hóa bằng Delphi-KAMET và Đại số mờ (Fuzzy Logic) $\rightarrow$ Mô hình thủy văn bán phân bố SWAT. Khung phân tích này cho phép bóc tách và định lượng riêng rẽ cũng như đồng thời tác động của LULC và BĐKH trên một lưu vực sông liên quốc gia phức tạp.

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

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

Luận án đóng góp vào hệ thống lý thuyết thủy văn - sinh thái và khoa học thông tin địa lý thông qua các luận điểm khoa học vững chắc:

  1. Mở rộng Lý thuyết Tự hành Tế bào (Cellular Automata) và Quá trình Ngẫu nhiên: Luận án chứng minh rằng ma trận xác suất chuyển đổi bậc nhất Markov kết hợp với luật lân cận không gian CA (Moore/Von Neumann neighborhood) có thể biểu diễn chính xác động lực thay đổi cảnh quan phức tạp khi được ràng buộc bởi các yếu tố sinh - lý - hóa và nhân sinh.
  2. Tham số hóa cơ chế phản ứng thủy văn đa cấp độ (HRUs): Bằng việc phân chia lưu vực thành các Đơn vị Phản ứng Thủy văn (Hydrologic Response Units - HRUs) đồng nhất về thổ nhưỡng, địa hình và thảm phủ, luận án làm sáng tỏ cơ chế truyền dẫn nước mặt và nước ngầm theo phương trình cân bằng nước đất của Arnold et al. (1998): $$SW_t = SW_0 + \sum_{i=1}^t (R_{day} - Q_{surf} - E_a - w_{seep} - Q_{gw})$$ (trong đó $SW_t$ là độ ẩm đất cuối cùng, $SW_0$ là độ ẩm đất ban đầu, $R_{day}$ là lượng mưa, $Q_{surf}$ là dòng chảy bề mặt, $E_a$ là lượng bốc thoát hơi nước, $w_{seep}$ là lượng nước thấm vào tầng ngấm, $Q_{gw}$ là dòng chảy ngầm ngày $i$).
  3. Khung đề xuất về sự tương tác phi tuyến giữa LULC và BĐKH: Luận án đưa ra mô hình khái niệm chứng minh rằng BĐKH đóng vai trò là "động lực ngoại sinh" chi phối tổng lượng đầu vào (nhiệt - mưa), trong khi LULC đóng vai trò là "bộ điều biến nội sinh" quyết định hệ số dòng chảy, tốc độ tập trung lũ và khả năng duy trì dòng chảy cơ bản (baseflow).

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

Khung phân tích của luận án tích hợp liên hoàn 3 lý thuyết nền tảng:

  • Lý thuyết Chuỗi Markov (Markov Chain): Xác định xác suất chuyển đổi trạng thái giữa 5 lớp phủ (Rừng, Nông nghiệp, Xây dựng, Vùng nước, Đất trống).
  • Lý thuyết Đánh giá Đa tiêu chí Không gian (Spatial Multi-Criteria Evaluation - SMCE): Ứng dụng hàm thuộc tính mờ (Fuzzy Membership Functions) và Phân tích Thứ bậc (AHP của Thomas Saaty) kết hợp phương pháp Kết hợp Tuyến tính có Trọng số (Weighted Linear Combination - WLC) để tạo lập bản đồ tính phù hợp (Suitability Maps).
  • Lý thuyết Mô hình hóa Thủy văn Bán phân bố (SWAT): Tích hợp dữ liệu số độ cao (DEM), bản đồ đất và khí tượng phân bố để tính toán dòng chảy.

Điều kiện biên (Boundary conditions) của khung phân tích được xác lập rõ: Áp dụng cho các lưu vực sông nhiệt đới gió mùa có độ chia cắt địa hình lớn, tồn tại tính chất lưu vực liên quốc gia với sự khác biệt về chế độ quản lý thảm phủ giữa các vùng lãnh thổ.

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

Thiết kế nghiên cứu

Nghiên cứu được thiết kế theo trường phái Thực chứng luận (Positivism) kết hợp Hiện thực phê phán (Critical Realism), sử dụng phương pháp tiếp cận hỗn hợp đa cấp độ và đa thời gian (Multi-level & Multi-temporal):

  • Cấp độ không gian: Toàn bộ lưu vực sông Cả ($27.200\text{ km}^2$), phân chia thành các tiểu lưu vực (sub-basins) và hệ thống hàng trăm HRUs.
  • Cấp độ thời gian: Chuỗi quan sát viễn thám và thống kê quá khứ (2005, 2010, 2015); giai đoạn hiệu chỉnh - kiểm định thủy văn (2007–2014); giai đoạn dự tính tương lai 2030 (chuỗi khí hậu 2020–2039).

Quy trình nghiên cứu rigorous

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

flowchart TD
    A[Dữ liệu Viễn thám Landsat/Sentinel & Bản đồ HTSDĐ] --> B[Phân loại MLC & Đánh giá độ chính xác Kappa]
    B --> C[Tính toán Ma trận Xác suất Chuyển đổi Markov]
    D[Tham vấn Chuyên gia Delphi 3 vòng theo Quy tắc KAMET] --> E[Chuẩn hóa Hàm mờ & Trọng số AHP/WLC]
    C & E --> F[Mô hình CA-Markov: Mô phỏng LULC 2015 & Dự tính LULC 2030]
    G[Dữ liệu DEM 30m, Thổ nhưỡng FAO, Khí tượng - Thủy văn ngày] --> H[Thiết lập Mô hình SWAT & Phân chia Tiểu lưu vực / HRUs]
    F & H --> I[Hiệu chỉnh & Kiểm định SWAT 2007-2014 qua NSE, R2, PBIAS]
    J[Kịch bản BĐKH RCP4.5 & RCP8.5 giai đoạn 2020-2039] --> K[Mô phỏng Dòng chảy Tương lai dưới Tác động Kép LULC & BĐKH]
    I & K --> L[Định lượng Biến đổi Dòng chảy Năm, Mùa Lũ, Mùa Cạn và Theo Tháng]
  1. Phân loại thảm phủ: Ứng dụng kỹ thuật Phân loại Hợp lý Cực đại (Maximum Likelihood Classification - MLC) trên ảnh vệ tinh Landsat 8, Sentinel 2 và Landsat TM (độ phân giải $30\text{m} \times 30\text{m}$). Khái niệm thảm phủ tuân thủ định nghĩa của FAO: "thảm phủ là lớp vỏ vật lý sinh học quan sát được trên bề mặt Trái đất" [68]. Kết quả phân loại được kiểm chứng chéo với dữ liệu điều tra thực địa và bản đồ hiện trạng sử dụng đất có tính pháp lý.
  2. Quy trình Delphi 3 vòng theo Quy tắc KAMET (Knowledge Acquisition for Multiple Experts with Time scale): Tham vấn nhóm chuyên gia liên ngành (Khí tượng, Thủy văn, Lâm nghiệp, Viễn thám, Quy hoạch) qua 3 vòng lấy ý kiến nghiêm ngặt nhằm xác định và đạt sự đồng thuận cao về các nhân tố động lực và ràng buộc không gian (khoảng cách tới đường giao thông, sông suối, trung tâm đô thị, độ dốc DEM).
  3. Mô hình hóa CA-Markov: Tích hợp ma trận xác suất chuyển đổi từ phân tích Markov với bản đồ tính phù hợp MCE, chạy lặp qua các chu kỳ tế bào CA để mô phỏng không gian LULC năm 2015 (dùng kiểm định mô hình) và dự tính không gian LULC năm 2030.
  4. Mô phỏng thủy văn SWAT: Thiết lập mô hình bán phân bố, nạp các trường dữ liệu địa hình (DEM $30\text{m}$), phân loại đất thổ nhưỡng, thuộc tính vật lý đất và chuỗi số liệu khí tượng ngày tại các trạm trong lưu vực.

Data và phân tích

Độ tin cậy của mô hình được đánh giá qua các chỉ số thống kê tiêu chuẩn quốc tế:

  • Kiểm định phân loại và mô phỏng thảm phủ: Sử dụng Hệ số $\text{Kappa}$ ($\text{Kappa standard}$, $\text{Kappa location}$, $\text{Kappa quantity}$), Độ chính xác Toàn cục (Overall Accuracy), Độ chính xác Nhà sản xuất (Producer Accuracy - PA) và Độ chính xác Người dùng (User Accuracy - UA). Kết quả phân loại thảm phủ và mô phỏng CA-Markov năm 2015 đều đạt hệ số $\text{Kappa} \ge 0{,}85$, chứng minh tính tương thích rất cao giữa mô phỏng và thực tế.
  • Hiệu chỉnh và kiểm định mô hình SWAT: Thực hiện hiệu chỉnh bộ thông số (Calibration: 2007–2010) và kiểm định độc lập (Validation: 2011–2014) tại 5 trạm thủy văn chính trên lưu vực sông Cả: Yên Thượng, Hòa Duyệt, Dừa, Quỳ Châu và Sơn Diệm.
    • Các chỉ số đánh giá bao gồm: Hệ số hiệu quả mô hình Nash–Sutcliffe Efficiency ($\text{NSE}$), Hệ số xác định ($R^2$) và Phần trăm sai số khối lượng ($\text{PBIAS}$).
    • Kết quả: Trong cả hai giai đoạn hiệu chỉnh và kiểm định, tại hầu hết các trạm, chỉ số $\text{NSE}$ đều đạt từ $0{,}75$ đến $0{,}88$, và $R^2$ đạt từ $0{,}78$ đến $0{,}90$. Theo tiêu chuẩn phân loại mức độ chính xác của Moriasi et al., kết quả mô phỏng đạt cấp "Tốt" đến "Rất tốt" (Good to Very Good).

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

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

  1. Quy luật biến động thảm phủ giai đoạn 2005–2030: Kết quả mô phỏng CA-Markov chỉ ra sự phân hóa mạnh mẽ giữa các lớp phủ. Đến năm 2030, diện tích đất xây dựng gia tăng nhanh chóng (tập trung mạnh tại vùng hạ lưu, TP Vinh và các trục giao thông huyết mạch). Diện tích đất nông nghiệp có xu hướng thu hẹp và chuyển đổi mục đích. Thảm phủ rừng thượng nguồn thuộc lãnh thổ Lào duy trì độ ổn định tương đối, tuy nhiên tại khu vực miền núi Tây Nghệ An và Hà Tĩnh, chất lượng rừng và độ che phủ có sự suy giảm cục bộ do áp lực kinh tế và trồng rừng sản xuất luân canh ngắn ngày.
  2. Sự phân cực khốc liệt của chế độ dòng chảy (Flow Regime Polarization): Tác động đồng thời của LULC 2030 và BĐKH (theo cả hai kịch bản RCP4.5 và RCP8.5 giai đoạn 2020–2039) làm gia tăng rõ rệt tính cực đoan của dòng chảy mặt lưu vực sông Cả:
    • Dòng chảy mùa lũ: Lưu lượng dòng chảy mùa lũ tại các trạm thủy văn tăng từ $8{,}5%$ đến $18{,}2%$ so với thời kỳ nền. Đỉnh lũ có xu hướng xuất hiện sớm hơn, biên độ lũ dâng cao nhanh hơn do diện tích đất xây dựng tăng và rừng phòng hộ bị phân mảnh, làm giảm hệ số thấm và rút ngắn thời gian tập trung dòng chảy.
    • Dòng chảy mùa kiệt: Lưu lượng dòng chảy mùa cạn suy giảm nghiêm trọng từ $11{,}8%$ đến $24{,}6%$ so với thời kỳ nền. Tình trạng hạn hán, thiếu nước cục bộ vào các tháng 3 đến tháng 6 trở nên trầm trọng hơn tại các vùng đồng bằng hạ du Nghệ An và Hà Tĩnh.
  3. Phát hiện phản trực giác về tương tác Thảm phủ - Dòng chảy kiệt: Nghiên cứu phát hiện tại một số tiểu lưu vực có sự gia tăng diện tích rừng trồng đơn loài thuần loại, dòng chảy kiệt không tăng như kỳ vọng lý thuyết truyền thống mà lại có xu thế giảm nhẹ vào các tháng cao điểm mùa khô. Nguyên nhân là do mức độ bốc thoát hơi nước (ET) của cây trồng ngoại lai sinh trưởng nhanh lớn hơn lượng nước ngầm tích lũy bổ sung, củng cố thêm bằng chứng thực nghiệm của Guardiola et al. (2010) tại Đông Nam Á.
  4. Bóc tách định lượng đóng góp riêng rẽ giữa LULC và BĐKH: BĐKH (lượng mưa và nhiệt độ) đóng vai trò chi phối chính đến sự biến thiên tổng lượng dòng chảy năm ($65% - 75%$ mức độ tác động), trong khi sự thay đổi LULC đóng vai trò then chốt trong việc làm thay đổi phân phối nội năm, làm gia tăng hệ số dòng chảy mặt mùa mưa và triệt tiêu dòng chảy ngầm mùa khô ($25% - 35%$ mức độ tác động).

Implications đa chiều

  • Về mặt Lý thuyết: Cung cấp cơ sở khoa học định lượng khẳng định mối quan hệ tương tác phi tuyến giữa LULC và BĐKH trên lưu vực nhiệt đới liên quốc gia, hoàn thiện khung nhận thức về tính dễ bị tổn thương tài nguyên nước.
  • Về mặt Phương pháp luận: Thiết lập quy trình chuẩn hóa kết hợp viễn thám, Delphi-KAMET, CA-Markov và SWAT, có khả năng chuyển giao và áp dụng hiệu quả cho các lưu vực sông khác tại Việt Nam và khu vực Đông Nam Á có điều kiện tương đồng hoặc thiếu hụt số liệu quan trắc.
  • Về mặt Thực tiễn và Quản lý:
    • Cung cấp luận cứ cốt lõi để các nhà hoạch định chính sách tại Nghệ An và Hà Tĩnh điều chỉnh Quy hoạch Sử dụng Đất giai đoạn 2021–2030, tầm nhìn 2050.
    • Là căn cứ kỹ thuật để tối ưu hóa quy trình vận hành liên hồ chứa trên lưu vực sông Cả (như hồ Bản Vẽ, Khe Bố, Hủa Na), cân bằng giữa mục tiêu cắt lũ mùa mưa và tích trữ dung tích hữu ích chống hạn mùa kiệt.
    • Đặt nền móng khoa học cho việc đàm phán và thiết lập cơ chế hợp tác quản lý tài nguyên nước xuyên biên giới giữa Việt Nam và Lào tại thượng nguồn sông Cả.

Limitations và Future Research

Mặc dù đạt được những kết quả đột phá, luận án thẳng thắn thừa nhận 4 hạn chế khoa học:

  1. Độ phân giải và sai số viễn thám: Do đặc thù mây mù dày đặc tại khu vực miền núi Bắc Trung Bộ, một số cảnh ảnh Landsat/Sentinel trong mùa mưa có độ phủ mây nhất định, đòi hỏi kỹ thuật lọc và nội suy phổ có thể tạo ra sai số biên nhỏ trong phân loại thảm phủ.
  2. Độ bất định của kịch bản khí hậu: Các kịch bản BĐKH chi tiết hóa động lực từ RCM vẫn chứa đựng độ bất định nội tại về lượng mưa cực trị ngắn hạn (mưa trận cục bộ).
  3. Tác động phi tuyến của công trình nhân tạo: Mô hình chưa xét hết toàn bộ tác động điều tiết phức tạp của hệ thống đập thủy điện vừa và nhỏ được xây dựng mới liên tục ở vùng thượng lưu Lào và Việt Nam trong những năm gần đây.
  4. Giả định ma trận chuyển đổi: Phân tích chuỗi Markov giả định xác suất chuyển đổi giữ nguyên tính dừng (stationary) tương đối, chưa dự báo hết những biến động đột biến về chính sách vĩ mô.

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

  • Tích hợp mạng nơ-ron tích chập (CNN) và học sâu (Deep Learning) với dữ liệu viễn thám siêu phổ (Hyperspectral) để nâng cao độ chính xác phân loại LULC.
  • Mở rộng mô hình SWAT kết hợp mô hình thủy lực 2D (như HEC-RAS 2D/MIKE 21) để mô phỏng ngập lụt chi tiết vùng hạ lưu ven biển đô thị.
  • Đánh giá bồi lắng bùn cát, xói mòn đất và lan truyền ô nhiễm nguồn nước nông nghiệp dưới tác động kép của LULC và BĐKH.
  • Lượng giá kinh tế dịch vụ hệ sinh thái nước (Ecosystem Services Valuation) phục vụ cơ chế chi trả dịch vụ môi trường rừng xuyên biên giới.

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

Công trình tạo ra những tác động lan tỏa mạnh mẽ trên nhiều phương diện:

  • Ảnh hưởng Học thuật: Đóng góp nguồn tư liệu khoa học chuẩn mực với chuỗi số liệu mô phỏng tin cậy, thúc đẩy các công bố quốc tế chuyên sâu về thủy văn lưu vực xuyên biên giới tại Đông Nam Á.
  • Chuyển đổi Ngành và Ứng dụng Quản lý: Giúp các cơ quan khí tượng thủy văn, các ban quản lý lưu vực sông và các công ty thủy điện nâng cao năng lực dự báo nguồn nước, lập kế hoạch phân bổ hạn ngạch khai thác nước mặt cho các ngành nông - lâm - thủy sản và công nghiệp.
  • Tác động Chính sách: Cung cấp cơ sở khoa học vững chắc phục vụ việc rà soát các quy chuẩn kỹ thuật quốc gia, kế hoạch hành động ứng phó BĐKH cấp tỉnh và Đề án An ninh nguồn nước quốc gia.
  • Lợi ích Xã hội và Môi trường: Giúp cộng đồng dân cư dọc lưu vực sông Cả giảm thiểu rủi ro sinh mạng và tài sản trước thiên tai lũ lụt cực đoan, đồng thời bảo vệ nguồn sinh kế nông nghiệp bền vững trong các đợt hạn hán khốc liệt.

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

  • Nghiên cứu sinh và Giới nghiên cứu Thủy văn - Môi trường: Tiếp cận khung phương pháp luận kết hợp CA-Markov - Delphi - SWAT hoàn chỉnh để giải quyết các bài toán thiếu số liệu và lưu vực xuyên biên giới.
  • Các Nhà khoa học cao cấp và Giảng viên Đại học: Kế thừa hệ thống luận điểm lý thuyết về tương tác thủy văn - sinh thái cảnh quan để phát triển các chương trình nghiên cứu chuyên sâu.
  • Các Ban Quản lý Lưu vực sông và Đơn vị Vận hành Công trình Thủy lợi - Thủy điện: Sử dụng bộ dữ liệu dự báo dòng chảy đến năm 2039 để lập phương án vận hành hồ chứa an toàn và hiệu quả kinh tế.
  • Các Cơ quan Quản lý Nhà nước và Hoạch định Chính sách (Bộ TN&MT, UBND tỉnh Nghệ An, Hà Tĩnh): Có được căn cứ khoa học trực quan, định lượng phục vụ công tác quy hoạch không gian sử dụng đất và chiến lược an ninh tài nguyên 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ì? Luận án đã mở rộng Lý thuyết Quá trình ngẫu nhiên Markov và Hệ tự hành tế bào trong việc giải bài toán động lực cảnh quan xuyên biên giới, chứng minh rằng sự kết hợp giữa quy tắc chuyển đổi Markovian và tối ưu hóa không gian MCE-AHP cho phép lượng hóa chính xác cơ chế phản hồi thủy văn (Hydrologic Feedback) theo phương trình cân bằng nước đất SWAT.
  2. Điểm đổi mới phương pháp luận nổi bật khi so sánh với các nghiên cứu trước đây? So với nghiên cứu của Sayemuzzaman & Jha (2014) tại Bắc Carolina hay Võ Quốc Tuấn et al. (2015) tại ĐBSCL, luận án đã loại bỏ tính chủ quan trong việc thiết lập ma trận phù hợp thông qua quy trình Delphi 3 vòng theo quy tắc KAMET và hàm thuộc tính mờ. Đồng thời, luận án giải quyết trọn vẹn "khoảng trống thượng nguồn bên ngoài biên giới" bằng ảnh vệ tinh độ phân giải cao kết hợp số liệu khí tượng tái phân tích, tạo ra chuỗi dữ liệu thủy văn liên tục và đồng nhất đạt $\text{NSE} > 0{,}80$.
  3. Phát hiện nào mang tính bất ngờ nhất từ dữ liệu thực nghiệm? Phát hiện về sự suy giảm dòng chảy kiệt tại các tiểu lưu vực có diện tích rừng trồng đơn loại tăng lên. Dữ liệu thực nghiệm bác bỏ quan điểm trực giác thông thường rằng cứ tăng diện tích rừng là tăng dòng chảy kiệt, chỉ ra rằng mức độ bốc thoát hơi nước cao của các loài cây ngoại lai sinh trưởng nhanh đã lấy đi phần lớn lượng ẩm trong đất vào mùa khô.
  4. Luận án có cung cấp quy trình tái lập (Replication Protocol) rõ ràng không? Toàn bộ quy trình từ tiền xử lý ảnh viễn thám, thuật toán MLC, lập ma trận chuyển đổi Markov, bảng câu hỏi Delphi-KAMET, tiêu chuẩn hóa hàm mờ, thiết lập các HRUs, xác định độ nhạy và bộ thông số hiệu chỉnh SWAT đều được chuẩn hóa chi tiết qua các sơ đồ, bảng biểu và công thức toán học, đảm bảo khả năng tái lập hoàn toàn trên các lưu vực tương tự.
  5. Chương trình nghiên cứu 10 năm tới được phác thảo như thế nào? Chương trình tập trung vào 3 trụ cột: (1) Ứng dụng trí tuệ nhân tạo (AI/Deep Learning) trong viễn thám thời gian thực để dự báo biến động thảm phủ; (2) Mô hình hóa động lực học bùn cát và chất lượng nước tích hợp thủy lực 2D dòng tràn; (3) Nghiên cứu cơ chế chính sách chi trả dịch vụ hệ sinh thái nước xuyên biên giới Việt - Lào.

Kết luận

  1. Luận án đã giải quyết trọn vẹn bài toán khoa học về mô phỏng biến động thảm phủ và đánh giá tác động đồng thời của LULC và BĐKH đến dòng chảy trên một lưu vực sông liên quốc gia quy mô lớn ($27.200\text{ km}^2$).
  2. Thiết lập thành công khung phương pháp luận tích hợp tiên tiến: Viễn thám đa thời gian $\rightarrow$ Phân tích chuỗi Markov $\rightarrow$ Hệ tự hành CA tối ưu hóa bằng Delphi-KAMET và Đại số mờ $\rightarrow$ Mô hình thủy văn bán phân bố SWAT.
  3. Dự tính định lượng bức tranh thảm phủ lưu vực sông Cả năm 2030 với 05 lớp phủ chủ đạo, đạt độ chính xác kiểm định xuất sắc ($\text{Kappa} \ge 0{,}85$).
  4. Chỉ ra quy luật phân cực sâu sắc của chế độ dòng chảy giai đoạn 2020–2039: dòng chảy lũ tăng từ $8{,}5%$ đến $18{,}2%$, dòng chảy cạn suy giảm từ $11{,}8%$ đến $24{,}6%$, làm sáng tỏ cơ chế tác động phản trực giác của rừng trồng đối với dòng chảy kiệt.
  5. Cung cấp hệ thống luận cứ khoa học thực tiễn vững chắc phục vụ quy hoạch sử dụng đất, vận hành liên hồ chứa và hoạch định chiến lược an ninh nguồn nước bền vững cho khu vực Bắc Trung Bộ trong bối cảnh BĐKH toàn cầu.