Tổng quan về luận án

Tài nguyên nước (TNN) lưu vực thượng nguồn sông Đồng Nai giữ vai trò sống còn đối với an ninh nguồn nước, điều tiết thủy văn và phát triển kinh tế - xã hội cho toàn bộ vùng hạ lưu, đặc biệt là siêu đô thị TP. Hồ Chí Minh và vùng kinh tế trọng điểm phía Nam. Tuy nhiên, áp lực kép từ biến đổi khí hậu (BĐKH) và chuyển dịch sử dụng đất do nhân sinh đã làm suy giảm nghiêm trọng lớp phủ thực vật, gia tăng xói mòn, biến động dòng chảy mùa kiệt và suy thoái chất lượng nước. Luận án tiến sĩ chuyên ngành Quản lý Tài nguyên và Môi trường (Mã số: 62850101) của nghiên cứu sinh Phạm Hùng, dưới sự hướng dẫn của PGS. Lê Văn Trung và PGS. Võ Lê Phú tại Trường Đại học Bách khoa – ĐHQG TP. Hồ Chí Minh, là công trình tiên phong thiết lập khung tích hợp Hệ thống Thông tin Địa lý (GIS), Viễn thám (RS) và Mô hình Đánh giá Đất và Nước (SWAT) nhằm xây dựng Hệ thống Hỗ trợ Ra quyết định (DSS) trong Quản lý Tổng hợp Tài nguyên Nước (IWRM).

Khoảng trống nghiên cứu (research gap) cốt lõi mà luận án giải quyết xuất phát từ thực trạng quản lý TNN tại Việt Nam: dù thể chế đã định hướng IWRM (Luật Tài nguyên nước 2012, Quyết định 81/2006/QĐ-TTg), việc triển khai thực tế vẫn mang tính phân tán, thiếu công cụ kỹ thuật định lượng đa chiều kết nối động lực học bề mặt với phản hồi thủy văn. Đánh giá của Chương trình Môi trường Liên hợp quốc (UNEP, 2018) chỉ ra chỉ số thực thi IWRM của Việt Nam chỉ đạt mức trung bình thấp 38/100 (trong đó hợp phần công cụ quản lý chỉ đạt 36/100). Luận án giải quyết câu hỏi nghiên cứu trung tâm: "Làm thế nào để ứng dụng tích hợp GIS, RS và SWAT nhằm tạo lập công cụ hỗ trợ IWRM ba hợp phần (Đánh giá TNN, Thông tin TNN, Phân bổ TNN) thích ứng với BĐKH trên lưu vực sông xuyên vùng?"

Hai giả thuyết nghiên cứu chính được kiểm định:

  • H1: Sự suy thoái lớp phủ rừng trong giai đoạn 1994–2020 làm gia tăng phi tuyến tính lưu lượng dòng chảy mùa lũ, suy giảm dòng chảy kiệt và khuếch đại rủi ro xói mòn đất vượt ngưỡng tự điều tiết của lưu vực.
  • H2: Tích hợp mô hình SWAT với dữ liệu viễn thám đa thời gian và kịch bản BĐKH RCP4.5 cho phép định lượng chính xác sự mất cân bằng cung - cầu nước và chỉ số tổn thương nguồn nước (VI) tại từng tiểu lưu vực (WMU) đến năm 2030.

Nghiên cứu được triển khai trên phạm vi không gian toàn bộ thượng lưu sông Đồng Nai với chuỗi số liệu khí tượng - thủy văn 30 năm (1984–2014) từ 5 trạm khí tượng và 3 trạm thủy văn (Thanh Bình, Đại Nga, Tà Lài), kết hợp chuỗi dữ liệu giám sát chất lượng nước 8 năm (2012–2020) tại 49 vị trí quan trắc định kỳ.

Literature Review và Positioning

Khung lý thuyết IWRM được Tổ chức Hợp tác Nước Toàn cầu (GWP, 2000) định nghĩa: "Quản lý tổng hợp tài nguyên nước là một quá trình thúc đẩy sự phối hợp phát triển và quản lý nguồn nước, đất đai và tài nguyên liên quan, nhằm tối đa hoá lợi ích kinh tế và phúc lợi xã hội một cách công bằng mà không phương hại đến tính bền vững của các hệ sinh thái thiết yếu". Khung này dựa trên tam giác quản trị bền vững gồm: Môi trường thuận lợi (Enabling environment), Khung thể chế (Institutional framework) và Công cụ quản lý (Management instruments) [GWP, 2000; UNEP, 2018].

Trong dòng nghiên cứu quốc tế, việc số hóa công cụ IWRM đã chuyển dịch từ các mô hình thủy văn đơn lẻ sang các khung không gian tích hợp. Dòng học thuật thứ nhất tập trung vào việc mô phỏng phản hồi dòng chảy từ biến động lớp phủ (LULC) bằng SWAT (như Mango và cs., 2011 tại lưu vực sông Mara, Kenya; Wagner và cs., 2013 tại lưu vực Mula & Mutha, Ấn Độ). Dòng học thuật thứ hai tập trung vào tích hợp mô hình bề mặt và nước ngầm đa tỷ lệ (như Abbas và cs., 2018 tích hợp SWAT và MODFLOW-NWT cho lưu vực sông Dee, Vương quốc Anh). Tuy nhiên, các tranh biện học thuật lớn vẫn tồn tại giữa hai quan điểm đối nghịch: một bên xem biến động khí hậu là động lực chính chi phối suy giảm an ninh nước (IPCC, 2007; Can và cs., 2015), trong khi nhóm đối lập chứng minh sự chuyển đổi mục đích sử dụng đất do nhân sinh (LULC) mới là nhân tố gây biến dạng sâu sắc nhất đến chu trình thủy văn và gia tăng tải lượng trầm tích (Wischmeier & Smith, 1978; Renard và cs., 1997; Loveland và cs., 2000).

Luận án định vị nghiên cứu tại điểm giao thoa của hai trường phái này, kế thừa và phát triển vượt bậc so với các nghiên cứu tiền nhiệm tại Việt Nam (như Trần Thị Vân, 2011; Lê Quang Toàn, 2014; Vũ Minh Tuấn, 2014). Công trình không dừng lại ở việc quan trắc biến động rừng đơn thuần hay mô phỏng dòng chảy lý thuyết, mà trực tiếp kết nối chuỗi ảnh Landsat đa thời gian (1994–2020), phương trình mất đất phổ dụng hiệu chỉnh (RUSLE), bộ chỉ số chất lượng nước (WQI) và chỉ số dễ bị tổn thương tài nguyên nước theo khung UNEP (2008).

So sánh với các nghiên cứu quốc tế điển hình như mô hình JESAT của Wolfgang-Albert Flügel và cs. (2009) tại Nam Phi/Đức hay nghiên cứu lượng hóa dịch vụ hệ sinh thái thủy văn của Aznarez và cs. (2021) tại lưu vực Laguna del Sauce (Uruguay), luận án đã thiết kế một khung công cụ toàn diện và thích ứng cao hơn với bối cảnh dữ liệu quan trắc phân mảnh tại các quốc gia đang phát triển.

Đó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 học thuyết Quản lý Tổng hợp Tài nguyên Nước (IWRM) của GWP và Khung đánh giá tổn thương nguồn nước của UNEP (2008) thông qua việc toán học hóa mối quan hệ giữa biến động cảnh quan bề mặt và tính tổn thương hệ thống sinh thái - xã hội. Luận án bác bỏ giả định truyền thống xem lưu vực sông như một thực thể đồng nhất, thay vào đó chứng minh tính đa tầng và không gian hóa của các ngưỡng dễ bị tổn thương.

Mô hình lý thuyết đề xuất chuẩn hóa cấu trúc đánh giá gồm 4 nhóm chỉ số thành phần:

  1. Sức ép nguồn nước ($RSt$): Biến thiên tỷ lệ khan hiếm ($RS_s$) và biến động dòng chảy tự nhiên ($RS_v$).
  2. Sức ép phát triển ($DP$): Tỷ lệ khai thác nước mặt ($DP_s$) và mức độ tiếp cận nước sạch của dân cư ($DP_d$).
  3. Sức khỏe hệ sinh thái ($EH$): Mức độ ô nhiễm môi trường nước ($EH_p$) và tỷ lệ suy thoái thảm phủ sinh thái ($EH_e$).
  4. Năng lực quản lý ($MC$): Hiệu quả sử dụng nước ($MC_e$), khả năng tiếp cận hạ tầng vệ sinh ($MC_s$) và năng lực kiểm soát mâu thuẫn nguồn nước ($MC_g$).

$$\text{Chỉ số Tổn thương (VI)} = f(RSt, DP, EH, MC)$$

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

Khung phân tích của luận án tích hợp đồng thời 3 nền tảng phương pháp luận:

  • Lý thuyết phân cấp cảnh quan: Ứng dụng thuật toán phân loại xác suất cực đại (Maximum Likelihood Classification - MLC) kết hợp chỉ số thực vật chuẩn hóa (NDVI) để bóc tách 4 thời kỳ dữ liệu Landsat (1994, 2004, 2014, 2020), xác lập ma trận chuyển dịch lớp phủ mặt đất (LC).
  • Lý thuyết thủy văn lưu vực vật lý: Triển khai mô hình SWAT trên cơ sở phân chia lưu vực thành các đơn vị đáp ứng thủy văn (Hydrologic Response Units - HRUs), liên kết trực tiếp với mô hình số độ cao (DEM) và dữ liệu thổ nhưỡng.
  • Lý thuyết quản trị đa tác nhân: Tích hợp dòng chảy môi trường (DCMT) theo Luật Tài nguyên nước 2012 vào bài toán cân bằng nước tổng thể, triệt tiêu nguyên nhân cốt lõi gây mâu thuẫn giữa thủy điện, nông nghiệp và cấp nước đô thị.

Điều kiện biên (boundary conditions) của khung phân tích được xác định rõ ràng: Áp dụng tối ưu cho các lưu vực sông miền núi nhiệt đới gió mùa có độ dốc cao, phân mùa mưa - khô sâu sắc, có sự cạnh tranh gay gắt về sử dụng nước giữa thủy điện thượng nguồn và an ninh nước hạ lưu.

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 quan điểm thực chứng (positivism) kết hợp hệ thống không gian thực nghiệm (spatial empiricism). Luận án triển khai cấu trúc nghiên cứu đa cấp (multi-level design):

  • Cấp độ toàn lưu vực: Đánh giá biến động khí hậu toàn vùng, cân bằng nước vĩ mô và phân tích động lực dòng chảy tại trạm cửa sổ Tà Lài.
  • Cấp độ tiểu lưu vực (Water Management Units - WMUs): Đánh giá chi tiết xói mòn, cân bằng cung - cầu theo ngành kinh tế và phân hạng rủi ro ô nhiễm.
  • Cấp độ vi mô (HRU và điểm quan trắc): Phân tích các quá trình động học dòng chảy mặt, nồng độ phú dưỡng ($NO_2^-$, $NO_3^-$, $NH_4^+$, $PO_4^{3-}$) và chất rắn lơ lửng (TSS).

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu được xây dựng đồng bộ, tuân thủ nguyên tắc tam giác đạc dữ liệu và phương pháp (data and method triangulation):

  1. Tiền xử lý và phân loại ảnh viễn thám: Sử dụng chuỗi ảnh vệ tinh Landsat đa cảm biến (Landsat 5 TM, Landsat 7 ETM+, Landsat 8 OLI) các năm 1994, 2004, 2014 và 2020. Độ chính xác phân loại được kiểm định bằng ma trận nhầm lẫn (confusion matrix), hệ số Kappa đạt chuẩn tin cậy cao ($> 0.85$).
  2. Mô hình hóa xói mòn đất bằng RUSLE: Tích hợp 5 lớp thông tin không gian: $$A = R \times K \times LS \times C \times P$$ Trong đó:
    • $R$: Hệ số xói mòn do mưa được tính từ chuỗi mưa tháng bình quân tại các trạm khí tượng ($MJ\cdot mm\cdot ha^{-1}\cdot h^{-1}\cdot yr^{-1}$).
    • $K$: Hệ số xói mòn đất dựa trên bản đồ thổ nhưỡng.
    • $LS$: Hệ số chiều dài và độ dốc trích xuất từ DEM.
    • $C$: Hệ số thảm phủ xác định qua biến thiên chỉ số NDVI.
    • $P$: Hệ số biện pháp bảo tồn đất canh tác.
  3. Mô phỏng thủy văn và chất lượng nước bằng SWAT: Hiệu chỉnh (calibration) và kiểm định (validation) mô hình thông qua dữ liệu lưu lượng thực đo tại 3 trạm thủy văn chính. Mô hình hóa tải lượng ô nhiễm nông nghiệp, sinh hoạt và công nghiệp.

Data và phân tích

Phân tích chất lượng nước mặt dựa trên QCVN 08-MT:2015/BTNMT và chỉ số WQI theo hướng dẫn của Tổng cục Môi trường. Các công cụ thống kê nâng cao gồm:

  • Phân tích phương sai (ANOVA) nhằm đánh giá sự khác biệt có ý nghĩa thống kê của WQI theo mùa và không gian.
  • Phân tích tương quan Pearson nhằm bóc tách mối liên hệ giữa các thông số lý hóa ($DO, BOD_5, COD, TSS, N-NH_4^+, P-PO_4^{3-}$) với chất lượng nước tổng thể.
  • Phân tích cụm (Cluster Analysis - CA) để nhóm các vị trí quan trắc có mức độ tương đồng về nguồn thải, xác định các điểm nóng ô nhiễm.

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

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

Thứ nhất, biến động lớp phủ và sự thoái hóa cảnh quan rừng diễn ra nghiêm trọng. Trong giai đoạn 1994–2020, diện tích thực phủ rừng trên toàn lưu vực đã suy giảm từ 69,6% (năm 1994) xuống còn 40,2% (năm 2020). Sự suy giảm gần 30% độ che phủ rừng tự nhiên do chuyển đổi sang đất canh tác nông nghiệp và công trình hạ tầng đã trực tiếp kích hoạt nguy cơ xói mòn đất trên quy mô lớn. Kết quả phân tích RUSLE chỉ ra rằng những khu vực có tiềm năng xói mòn đất cao và rất cao chiếm tới 50,48% tổng diện tích toàn lưu vực thượng nguồn.

Thứ hai, biến dạng chế độ thủy văn và gia tăng hiện tượng bồi lắng. Mô phỏng SWAT cho thấy hệ số tương quan $R^2$ giữa lưu lượng dòng chảy mô phỏng và thực đo tại 3 trạm thủy văn Thanh Bình, Đại Nga và Tà Lài đều đạt mức tin cậy cao ($R^2 > 0.8$). Quá trình mất rừng đã làm gia tăng hệ số dòng chảy mặt trong mùa mưa, đồng thời làm suy giảm nghiêm trọng lượng bổ cập nước ngầm, dẫn đến nguy cơ kiệt dòng trong mùa khô. Lượng bồi lắng và tải lượng trầm tích gia tăng mạnh qua 4 kịch bản LULC (1994, 2004, 2014, 2020), làm suy giảm dung tích hữu ích của các hồ chứa thủy điện lớn trên bậc thang sông Đồng Nai.

Thứ ba, suy thoái chất lượng nước và ô nhiễm chất dinh dưỡng. Phân tích chuỗi dữ liệu tại 49 điểm quan trắc (2012–2020) cho thấy chất lượng nước mặt có xu hướng giảm rõ rệt vào mùa khô. Tỷ lệ mẫu vượt ngưỡng quy chuẩn tập trung vào các chỉ tiêu $TSS, NH_4^+, NO_2^-$ và Coliform. Mô hình SWAT đã bóc tách được sự phân bố không gian của các chất ô nhiễm: nồng độ nitrit, nitrate, amoni và phosphate tăng đột biến tại các WMU có mật độ canh tác nông nghiệp cao và các vùng tập trung dân cư đô thị thiếu hệ thống thu gom, xử lý nước thải.

Thứ tư, gia tăng tính dễ bị tổn thương tài nguyên nước ($VI$). Đánh giá hiện trạng cho thấy phần lớn lưu vực có chỉ số tổn thương ở mức trung bình ($0,2 < VI < 0,4$). Tuy nhiên, dưới tác động kép của phát triển kinh tế - xã hội và kịch bản BĐKH RCP4.5, dự báo đến giai đoạn 2020–2030, chỉ số VI sẽ gia tăng lên mức cao ($0,4 < VI < 0,7$), tập trung tại các tiểu lưu vực trọng điểm: Đắk Nông, Đắk R'Keh và La Ngà. Nguyên nhân chính do áp lực khai thác nước tưới cho cây công nghiệp (cà phê, hồ tiêu) tăng vọt trong mùa khô kết hợp với suy giảm dòng chảy tự nhiên.

Chỉ số / Đại lượng Đo lường Giá trị Hiện trạng (1994 - 2020) Dự báo Giai đoạn 2020 - 2030 (RCP4.5) Mức độ Tác động / Ý nghĩa Thực tiễn
Tỷ lệ che phủ rừng Giảm từ 69,6% (1994) xuống 40,2% (2020) Dự báo tiếp tục chịu áp lực do chuyển đổi đất nông nghiệp Mất thảm phủ thực vật gây gia tăng dòng chảy mặt
Diện tích tiềm năng xói mòn cao 50,48% diện tích toàn lưu vực Nguy cơ mở rộng nếu không có giải pháp canh tác bảo tồn Bồi lắng lòng hồ thủy điện, thoái hóa tầng đất mặt
Chỉ số Dễ bị Tổn thương (VI) Mức trung bình: $0,2 < VI < 0,4$ Tăng lên mức cao: $0,4 < VI < 0,7$ Tập trung tại tiểu lưu vực Đắk Nông, Đắk R'Keh, La Ngà
Áp lực Cân bằng Nước mùa khô Căng thẳng cục bộ tại các vùng chuyên canh Thiếu hụt trầm trọng khi tính toán đủ dòng chảy môi trường Gia tăng xung đột giữa phát điện và tưới tiêu

Implications đa chiều

  • Đóng góp lý luận: Cung cấp cơ sở khoa học khẳng định sự cần thiết phải lồng ghép mô hình thủy văn phân bố vào quy trình tính toán chỉ số tổn thương TNN (bổ khuyết cho phương pháp tĩnh truyền thống của UNEP).
  • Đổi mới phương pháp: Thiết lập quy trình chuẩn hóa kết nối GIS – Viễn thám Landsat – RUSLE – SWAT – WQI, có khả năng chuyển giao và áp dụng trực tiếp cho các lưu vực sông liên tỉnh khác tại Việt Nam như lưu vực sông Ba, Sê San, Srêpôk hay Vu Gia – Thu Bồn.
  • Ứng dụng thực tiễn: Hệ thống bản đồ phân vùng rủi ro ô nhiễm nguồn điểm, nguồn diện và bản đồ tổn thương nguồn nước là công cụ trực quan giúp Sở TN&MT Lâm Đồng, Đắk Nông, Bình Phước, Đồng Nai phân bổ hạn ngạch xả thải và cấp phép khai thác TNN hợp lý.
  • Hàm ý chính sách: Khẳng định nguyên tắc quản lý TNN theo ranh giới lưu vực sông thay vì ranh giới hành chính; cung cấp căn cứ tính toán chi trả dịch vụ môi trường rừng (PFES) dựa trên lượng đất xói mòn được giữ lại và lưu lượng dòng chảy duy trì trong mùa kiệt.

Limitations và Future Research

Luận án thẳng thắn thừa nhận 4 giới hạn nghiên cứu:

  1. Độ phân giải không gian của ảnh vệ tinh: Sử dụng ảnh Landsat độ phân giải 30m phù hợp cho quy mô toàn lưu vực nhưng có thể bỏ sót các biến động vi mô của các mảng rừng phân tán và các công trình nhân tạo nhỏ.
  2. Dữ liệu tầng chứa nước ngầm: SWAT mô phỏng dòng chảy ngầm dựa trên các khối hộp giả định (shallow/deep aquifers); chưa tích hợp mô hình dòng thấm 3D chuyên sâu (như MODFLOW) để đánh giá cạn kiệt tầng ngầm do bơm hút tưới cà phê.
  3. Tác động điều tiết của hệ thống hồ chứa: Các kịch bản xả nước thực tế của các nhà máy thủy điện bậc thang chưa được tích hợp hoàn toàn theo thời gian thực do hạn chế trong chia sẻ dữ liệu vận hành.
  4. Quy mô kịch bản BĐKH: Nghiên cứu tập trung vào kịch bản nồng độ đại diện trung bình RCP4.5; chưa so sánh đối chứng chi tiết với kịch bản phát thải cao RCP8.5.

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

  • Tích hợp dữ liệu viễn thám độ phân giải siêu cao (Sentinel-2, PlanetScope) kết hợp thuật toán học máy (Random Forest, SVM, Deep Learning) để nâng cao độ chính xác phân loại LULC.
  • Phát triển mô hình tích hợp SWAT-MODFLOW nhằm đánh giá tương tác nước mặt - nước ngầm toàn diện.
  • Xây dựng nền tảng WebGIS hỗ trợ ra quyết định theo thời gian thực, kết nối mạng lưới trạm quan trắc IoT tự động phục vụ điều hành liên hồ chứa.

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

Nghiên cứu tạo tiền đề nâng cao chất lượng học thuật với tiềm năng trích dẫn cao trong các tạp chí quốc tế thuộc danh mục ISI/Scopus chuyên ngành Thủy văn và Quản lý Môi trường (Journal of Hydrology, Science of The Total Environment, Water Resources Management).

Về mặt chuyển đổi thực tiễn, công trình đóng góp trực tiếp vào việc tái cấu trúc công tác quản lý tài nguyên nước tại lưu vực sông Đồng Nai – vùng kinh tế đóng góp hơn 40% GDP cả nước. Kết quả nghiên cứu là công cụ hỗ trợ trực tiếp cho Bộ Tài nguyên và Môi trường, Ủy ban Bảo vệ Môi trường Lưu vực hệ thống sông Đồng Nai, Ủy ban nhân dân tỉnh Lâm Đồng và các tỉnh lân cận trong việc lập Quy hoạch Tài nguyên Nước thời kỳ 2021–2030, tầm nhìn đến 2050 theo Luật Quy hoạch 2017.

Về mặt xã hội, việc xác định chính xác các điểm nóng dễ bị tổn thương giúp bảo vệ nguồn cấp nước sinh hoạt an toàn cho hơn 10 triệu dân tại TP. Hồ Chí Minh và các đô thị hạ du.

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

  • Nghiên cứu sinh và Giới học thuật: Tiếp cận phương pháp luận chuẩn mực về tích hợp công nghệ không gian (RS/GIS) và mô hình hóa thủy văn (SWAT, RUSLE) để định lượng các thành phần IWRM.
  • Cơ quan Quản lý Nhà nước (Bộ TN&MT, Cục Quản lý TNN, các Sở TN&MT): Sở hữu cơ sở dữ liệu số GIS đồng bộ, bản đồ rủi ro ô nhiễm và công cụ tính toán cân bằng nước phục vụ thẩm định cấp phép, quy hoạch và thanh tra môi trường.
  • Doanh nghiệp và Đơn vị Khai thác Nước (Tập đoàn Điện lực, Công ty Cấp nước, Nông trường Nông nghiệp): Nắm bắt xu thế biến động nguồn nước đến năm 2030 để tối ưu hóa quy trình vận hành hồ chứa, chủ động kế hoạch sản xuất và giảm thiểu rủi ro thiếu hụt nguồn cung trong mùa khô.
  • Cộng đồng Dân cư trên Lưu vực: Thụ hưởng môi trường nước trong sạch hơn nhờ các biện pháp kiểm soát ô nhiễm nguồn diện và giảm thiểu nguy cơ lũ lụt, hạn hán.

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 khung đánh giá tính dễ bị tổn thương tài nguyên nước của UNEP (2008) bằng cách tích hợp trực tiếp các biến động động học bề mặt (LULC từ ảnh viễn thám Landsat) và các phản hồi thủy văn - xói mòn mô phỏng từ mô hình SWAT/RUSLE vào 4 nhóm thông số thành phần ($RSt, DP, EH, MC$). Điều này chuyển đổi khung đánh giá của UNEP từ dạng bảng biểu thống kê hành chính tĩnh sang mô hình phân bố không gian động lực học.

2. Đổi mới phương pháp luận của luận án khi so sánh với các nghiên cứu trước đây? So với các nghiên cứu trong nước (chỉ phân tích đơn lẻ chất lượng nước bằng WQI hoặc đánh giá biến động rừng tĩnh) và các nghiên cứu quốc tế như Mango và cs. (2011) hay Can và cs. (2015), luận án đã thiết lập một chuỗi xử lý khép kín hoàn chỉnh: Viễn thám Landsat (1994–2020) $\rightarrow$ Mô hình mất đất RUSLE $\rightarrow$ Mô hình thủy văn SWAT phân bố $\rightarrow$ Phân tích cụm chất lượng nước $\rightarrow$ Đánh giá cân bằng nước có xét đến Dòng chảy môi trường (DCMT) $\rightarrow$ Bản đồ hóa chỉ số tổn thương VI theo kịch bản BĐKH RCP4.5.

3. Phát hiện bất ngờ nhất từ dữ liệu thực nghiệm là gì? Mặc dù lượng mưa bình quân năm của vùng Tây Nguyên có xu hướng tăng nhẹ tại một số trạm khí tượng (như Bảo Lộc), nhưng chỉ số tổn thương tài nguyên nước ($VI$) tại các tiểu lưu vực Đắk Nông, Đắk R'Keh và La Ngà trong giai đoạn 2020–2030 lại gia tăng đột biến lên mức cao ($0,4 < VI < 0,7$). Dữ liệu chứng minh rằng áp lực gia tăng nhu cầu nước tưới nông nghiệp mùa khô cùng với việc mất 29,4% diện tích rừng tự nhiên (làm triệt tiêu khả năng giữ nước ngầm mùa kiệt) có sức tàn phá an ninh nước vượt trội hơn tác động thuần túy của biến đổi lượng mưa.

4. Luận án có cung cấp quy trình tái lặp (replication protocol) đầy đủ không? Có. Luận án cung cấp chi tiết toàn bộ các bước tính toán hệ số xói mòn ($R, K, LS, C, P$), cấu hình thiết lập mô hình SWAT (phân chia lưu vực, ngưỡng HRU, nguồn tham số thổ nhưỡng/khí tượng), phương pháp nội suy GIS (Kriging/IDW) và quy trình tính toán trọng số cho 4 nhóm thông số của chỉ số VI, cho phép tái lập quy trình trên bất kỳ lưu vực sông nào có điều kiện tương đồng.

5. Chương trình nghiên cứu 10 năm được phác thảo như thế nào? Chương trình 10 năm định hướng: (1) Phát triển mô hình tích hợp SWAT - MODFLOW - 2D Hydrodynamic cho toàn bộ lưu vực sông Đồng Nai từ thượng nguồn đến hạ lưu cửa sông; (2) Tích hợp công nghệ học sâu và viễn thám radar/quang học độ phân giải cao để giám sát biến động mặt phủ theo thời gian thực; (3) Xây dựng hệ thống bản sao số (Digital Twin) lưu vực sông phục vụ vận hành thông minh và thích ứng BĐKH.

Kết luận

  1. Luận án đã hệ thống hóa và phát triển thành công cơ sở khoa học, phương pháp luận tích hợp GIS, Viễn thám và mô hình SWAT nhằm xây dựng Hệ thống Hỗ trợ Ra quyết định hoàn chỉnh cho Quản lý Tổng hợp Tài nguyên Nước (IWRM) tại lưu vực thượng nguồn sông Đồng Nai.
  2. Thiết lập bộ Cơ sở Dữ liệu GIS không gian chuyên đề đa thời gian, tích hợp chuỗi số liệu khí tượng - thủy văn 30 năm (1984–2014) và chất lượng nước 8 năm (2012–2020) tại 49 trạm quan trắc, giải quyết triệt để vấn đề phân tán thông tin trong quản lý TNN.
  3. Định lượng chính xác mức độ suy giảm rừng từ 69,6% (1994) xuống 40,2% (2020) và chứng minh mối quan hệ nhân quả đối với 50,48% diện tích lưu vực có tiềm năng xói mòn cao, làm biến dạng chế độ dòng chảy mặt và gia tăng bồi lắng lòng hồ.
  4. Mô phỏng thành công động lực học chất lượng nước và lưu lượng dòng chảy bằng mô hình SWAT, bóc tách chính xác các nguồn thải ô nhiễm dinh dưỡng ($N, P$) và bùn cát lơ lửng ($TSS$).
  5. Lập bản đồ cân bằng nước chi tiết có tích hợp Dòng chảy môi trường (DCMT) và dự báo chỉ số dễ bị tổn thương TNN ($VI$) đến năm 2030 theo kịch bản RCP4.5, nhận diện chính xác các "điểm nóng" cần ưu tiên hành động tại các tiểu lưu vực Đắk Nông, Đắk R'Keh và La Ngà.
  6. Đề xuất hệ thống giải pháp IWRM khả thi, kết hợp hài hòa giữa giải pháp công trình và phi công trình, tạo luận cứ khoa học vững chắc phục vụ phát triển bền vững lưu vực sông Đồng Nai và cung cấp mô hình mẫu có khả năng nhân rộng trên toàn quốc.