Tổng quan về luận án

Biến đổi khí hậu toàn cầu đang làm gia tăng tần suất, cường độ và tính bất định của các hiểm họa thủy văn cực đoan, đặt ra những thách thức sống còn đối với hệ thống sinh thái nông nghiệp tại các vùng duyên hải nhiệt đới gió mùa. Luận án tiến sĩ ngành Quản lý đất đai (Mã số: 9850103) của tác giả Nguyễn Bích Ngọc (2021) với đề tài: "Nghiên cứu ảnh hưởng của ngập lụt và hạn hán đến sử dụng đất nông nghiệp ở huyện Quảng Điền, tỉnh Thừa Thiên Huế trong bối cảnh biến đổi khí hậu", dưới sự hướng dẫn khoa học của PGS.TS Nguyễn Hữu Ngữ và PGS.TS Trần Thanh Đức tại Trường Đại học Nông Lâm – Đại học Huế, là một công trình tiên phong thiết lập mô hình tích hợp không gian – khí hậu nhằm giải quyết bài toán quản trị tài nguyên đất trước tác động kép của ngập lụt và hạn hán.

Khoảng trống nghiên cứu (research gap) được xác định rõ rệt: Các công trình trước đây phần lớn tiếp cận đơn lẻ hiểm họa lũ lụt (như Karlsen et al., 2008; Islam et al., 2009) hoặc hạn hán (như Niko Wanders et al., 2010) trên quy mô vùng rộng lớn, thiếu sự liên kết không gian động với sự chuyển dịch cơ cấu cây trồng vi mô và chưa lượng hóa tường minh mức độ phơi nhiễm (exposure level) của từng loại hình sử dụng đất (Land Use Types - LUTs) cụ thể tại vùng trũng duyên hải chịu ảnh hưởng mạnh của đầm phá Tam Giang – Cầu Hai.

Nghiên cứu tập trung giải quyết ba câu hỏi nghiên cứu (Research Questions - RQ) và kiểm định ba giả thuyết tương ứng (Hypotheses - H):

  • RQ1: Hiện trạng và các động lực kinh tế - xã hội - tự nhiên nào chi phối sự biến động đất nông nghiệp giai đoạn 2005–2019 tại huyện Quảng Điền?
    • H1: Biến đổi khí hậu (ngập lụt và hạn hán) là một trong những nhân tố tác động có ý nghĩa thống kê cao làm suy giảm diện tích đất canh tác.
  • RQ2: Mức độ tổn thương và phân bố không gian của đất trồng lúa, đất cây hàng năm, đất lâm nghiệp và đất nuôi trồng thủy sản trước các đợt thiên tai lịch sử diễn ra như thế nào?
    • H2: Thuật toán viễn thám đa chỉ số kết hợp mô hình số hóa độ cao DEM và chỉ số chuẩn hóa giáng thủy (SPI) phản ánh chính xác ranh giới ngập và hạn thực tế.
  • RQ3: Xu thế ngập lụt và hạn hán sẽ tác động đến quỹ đất nông nghiệp đến năm 2030 và 2035 theo các kịch bản khí hậu ra sao, và những giải pháp quy hoạch thích ứng nào mang tính khả thi cao nhất?
    • H3: Tái cấu trúc không gian sử dụng đất theo cụm sinh thái nông thái sẽ giảm thiểu tối đa rủi ro thiệt hại mùa màng.

Khung lý thuyết (theoretical framework) của luận án dựa trên Lý thuyết Sử dụng đất bền vững, Thuyết Địa mạo thủy văn ứng dụng và Khung phân tích tổn thương sinh thái nông nghiệp của IPCC. Đóng góp đột phá của công trình được lượng hóa cụ thể: Xác lập bản đồ phân vùng ngập lụt từ chuỗi dữ liệu ảnh vệ tinh Landsat 7/8 với độ phân giải 30m trong các trận lũ cực đoan năm 2017 và 2019; xây dựng bản đồ phân vùng hạn hán vụ Hè Thu năm 2015 bằng chỉ số SPI kết hợp trạm quan trắc mặt đất và mưa vệ tinh TRMM; đồng thời dự báo chi tiết hóa phạm vi ngập lụt đến năm 2030 ứng với kịch bản nước biển dâng thêm 13 cm và hạn hán vụ Hè Thu đến năm 2035 theo kịch bản RCP 4.5. Phạm vi nghiên cứu bao trùm 8.194,13 ha đất nông nghiệp (chiếm 50,26% tổng diện tích tự nhiên) trên toàn bộ 11 đơn vị hành chính cấp xã/thị trấn của huyện Quảng Điền trong chuỗi dữ liệu khí tượng thủy văn 1999–2019.

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 của ba dòng nghiên cứu chính:

Dòng nghiên cứu Tác giả tiêu biểu Trọng tâm học thuật Hạn chế tồn tại
1. Viễn thám & Nhận diện ngập lụt Rouse et al. (1973), Karlsen et al. (2008), Islam et al. (2009), Dong et al. (2014) Ứng dụng chỉ số phản xạ phổ NDVI, EVI, LSWI để tách bề mặt nước tự do và nước ngập dưới tán thực vật Thiếu kiểm chứng đồng thời với mô hình số hóa độ cao DEM vi địa hình và số liệu thủy văn trạm mặt đất
2. Chỉ số khí tượng & Đánh giá hạn hán Palmer (1965), McKee et al. (1993), Niko Wanders et al. (2010), Tsakiris et al. (2007) Phát triển chỉ số PDSI, SPI, RDI để định lượng thiếu hụt ẩm tích lũy theo các thang thời gian Thường nội suy trạm thưa, chưa tích hợp dữ liệu mưa viễn thám đa nguồn để giải quyết sai số địa hình
3. Chuyển dịch sử dụng đất nông nghiệp Lê Quang Trí (2005), Lê Thị Giang (2012), Nguyễn Thị Hồng Hạnh (2014) Phân tích động lực chuyển đổi đất nông nghiệp dưới tác động của kinh tế - xã hội và thị trường Tiếp cận nhân tố khí hậu ở mức định tính, chưa mô hình hóa thành biến số độc lập trong phân tích hồi quy

Cuộc tranh luận học thuật (scholarly debate) nảy sinh giữa hai trường phái: Một trường phái cho rằng diễn biến thị trường và quy hoạch hạ tầng đô thị hóa giữ vai trò tuyệt đối trong chuyển đổi cơ cấu đất (quan điểm kinh tế học thuần túy); trường phái thứ hai nhấn mạnh các cú sốc khí hậu cực đoan mới là lực đẩy quyết định làm suy thoái độ phì nhiêu và biến đổi quỹ đất sản xuất (quan điểm sinh thái học thích ứng). Luận án của Nguyễn Bích Ngọc định vị vững chắc tại điểm giao thoa của hai trường phái này, chứng minh rằng mặc dù yếu tố kinh tế - xã hội dẫn dắt quy hoạch, nhưng chính các xung đột thủy văn (ngập lụt và khô hạn) là tác nhân kích hoạt các đợt chuyển đổi đất đột ngột.

Công trình tạo bước tiến vượt bậc khi so sánh với các nghiên cứu quốc tế:

  • So với nghiên cứu của White & Lewis (2011) tại lưu vực Great Artesian (Úc) vốn sử dụng ảnh vệ tinh độ phân giải siêu cao để quan sát phản ứng của thực vật ngập nước nhưng chi phí dữ liệu lớn, luận án đã tối ưu hóa chuỗi ảnh Landsat TM/ETM+/OLI hoàn toàn miễn phí của USGS với thuật toán ngưỡng kép kết hợp LSWI và NDVI, đạt độ chính xác tương đương trên diện tích canh tác phức tạp.
  • So với các nghiên cứu ngưỡng hạn tại đồng bằng Bắc Trung Quốc (lượng mưa <350 mm/năm gây giảm năng suất và <250 mm/năm gây hạn hán trầm trọng), nghiên cứu này đã bản địa hóa chỉ số SPI theo các tháng sinh trưởng then chốt của cây lúa vụ Hè Thu miền Trung Việt Nam, chỉ ra rằng giá trị SPI < -1,5 vào tháng 6 là giới hạn sinh thái gây tổn thương nghiêm trọng cho toàn bộ hệ thống canh tá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 sâu sắc thêm Lý thuyết Địa tô và Giá trị Đất đai của Karl Marx và William Petty khi tích hợp các biến số rủi ro khí hậu cực đoan vào hàm giá trị sử dụng đất. Trong điều kiện biến đổi khí hậu, độ phì tự nhiên và độ phì kinh tế của đất không còn là hàm số tuyến tính của chi phí lao động và phân bón, mà bị biến dạng nghiêm trọng bởi tần suất rửa trôi (do ngập lụt) và tích lũy muối/khô kiệt nước mao dẫn (do hạn hán).

Mô hình lý thuyết đề xuất các mệnh đề khoa học (propositions):

  • Mệnh đề P1: Độ bền vững của hệ thống sử dụng đất nông nghiệp ven biển đồng bằng tỷ lệ nghịch với biên độ dao động giữa chỉ số ẩm bề mặt cực đại ($LSWI_{max}$) trong mùa mưa và chỉ số chuẩn hóa giáng thủy cực tiểu ($SPI_{min}$) trong mùa khô.
  • Mệnh đề P2: Điểm chuyển dịch hình thái canh tác (agricultural land-use regime shift) xảy ra khi tần suất ngập lụt vượt quá ngưỡng chịu đựng sinh lý của giống cây trồng chủ lực trong 3 vụ liên tiếp hoặc giá trị SPI rơi vào ngưỡng khô hạn nặng liên tục quá 45 ngày.

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

Khung phân tích của luận án tích hợp liên ngành ba trụ cột: Khoa học Thay đổi Sử dụng Đất (Land Change Science), Viễn thám Thủy văn Ứng dụng và Kinh tế lượng Không gian. Khung phân tích thiết lập các bước phân tách đa tầng:

  1. Tầng vật lý - khí hậu: Xác định trường không gian ngập lụt (bằng tổ hợp ảnh viễn thám phản xạ hồng ngoại cận và sóng ngắn kết hợp DEM) và trường không gian hạn hán (bằng nội suy địa thống kê Kriging từ chỉ số SPI tính toán từ chuỗi trạm quan trắc kết hợp dữ liệu vệ tinh TRMM).
  2. Tầng hiện trạng sử dụng đất: Chồng xếp các lớp bản đồ chuyên đề để xác định ma trận chuyển dịch của 4 loại đất chính: Đất trồng lúa (LUA), Đất trồng cây hàng năm khác (BHK), Đất lâm nghiệp (LNP) và Đất nuôi trồng thủy sản (NTS).
  3. Tầng thích ứng thể chế và cộng đồng: Định hình không gian quy hoạch dựa trên điều kiện biên (boundary conditions) là vùng trũng ngập hạ lưu sông Bồ và cồn cát ven phá Tam Giang.

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

+---------------------------------------------------------------------------------------+
|                              QUY TRÌNH NGHIÊN CỨU TÍCH HỢP                            |
+---------------------------------------------------------------------------------------+
|                                                                                       |
|   [DỮ LIỆU ĐẦU VÀO]                                                                   |
|   - Viễn thám: Landsat 7/8 (30m), Mưa vệ tinh TRMM 3B43                               |
|   - Địa không gian: SRTM DEM, Bản đồ hiện trạng & quy hoạch đất                       |
|   - Khí tượng thủy văn: Chuỗi số liệu mưa trạm Phú Ốc, Huế (1999-2019)                |
|   - Điều tra xã hội học: 11 xã, 381 phiếu khảo sát hộ + 45 chuyên gia cán bộ          |
|   [XỬ LÝ DỮ LIỆU CHUYÊN SÂU]                                                          |
|   - ENVI 5.x: Giải đoán chỉ số NDVI, LSWI, trích xuất mặt nước ngập lũ               |
|   - SPIGenerator v1.0: Tính chuỗi chỉ số SPI (thời đoạn 1, 3, 6, 12 tháng)            |
|   - ArcGIS 10.x: Nội suy Kriging không gian hạn hán, mô phỏng DEM ngập lũ             |
|   - SPSS: Phân tích nhân tố, tương quan Pearson, Hồi quy tuyến tính đa biến           |
|   [KẾT QUẢ ĐẦU RA]                                                                    |
|   - Bản đồ phân vùng ngập lũ thực tế (2017, 2019) & Dự báo 2030 (+13cm NBD)          |
|   - Bản đồ phân vùng hạn hán vụ Hè Thu (2015) & Dự báo 2035 (Kịch bản RCP 4.5)        |
|   - Hệ thống giải pháp phân vùng thích ứng cho 3 tiểu vùng sinh thái nông nghiệp      |
|                                                                                       |
+---------------------------------------------------------------------------------------+

Thiết kế nghiên cứu

Luận án tuân thủ thế giới quan thực chứng (positivism) kết hợp chủ nghĩa hiện thực phản tư (critical realism). Thiết kế nghiên cứu hỗn hợp (mixed-methods design) kết hợp chặt chẽ giữa định lượng không gian địa lý và phân tích thống kê kinh tế lượng với khảo sát định tính có sự tham gia của các bên liên quan (Stakeholder Participatory Approach). Thiết kế đa cấp độ (multi-level design) được cấu trúc: Cấp huyện (vĩ mô), Cấp xã/cụm xã (trung mô) và Cấp nông hộ sản xuất (vi mô).

Quy trình nghiên cứu rigorous

Quy trình thu thập và xử lý dữ liệu được tiến hành nghiêm ngặt qua 4 kênh tam giác đạc (triangulation):

  • Thu thập dữ liệu thứ cấp: Chuỗi số liệu khí tượng thủy văn 1999–2019 tại trạm Phú Ốc trên sông Bồ và các trạm phụ cận; bản đồ kiểm kê đất đai các năm 2005, 2010, 2015, 2019 của huyện Quảng Điền; kịch bản biến đổi khí hậu và nước biển dâng của Bộ Tài nguyên và Môi trường (2016).
  • Khảo sát dữ liệu sơ cấp: Triển khai phỏng vấn sâu cán bộ quản lý (45 cán bộ Phòng TN&MT, Nông nghiệp, UBND 11 xã) và điều tra bảng hỏi bán cấu trúc đối với 381 hộ gia đình phân bổ theo phương pháp chọn mẫu phân tầng ngẫu nhiên (stratified random sampling) tại các vùng sinh thái đặc trưng: vùng ngập lũ ven sông, vùng cồn cát ven biển và vùng ven đầm phá.
  • Độ tin cậy và giá trị: Bảng hỏi Likert 5 mức độ được kiểm định độ tin cậy nhất quán nội tại qua hệ số Cronbach's Alpha ($\alpha > 0,8$). Giá trị hội tụ và giá trị phân biệt được khẳng định thông qua thảo luận nhóm tập trung (FGD) tại các điểm nóng thiên tai.

Data và phân tích

Quy trình trích xuất và xử lý dữ liệu không gian được vận hành trên các phần mềm chuyên dụng:

  • ArcGIS for DesktopENVI: Xử lý chuỗi ảnh Landsat 7 ETM+ và Landsat 8 OLI/TIRS (độ phân giải 30m) vào các thời điểm đỉnh lũ tháng 10 và tháng 11 các năm 2017 và 2019. Sử dụng thuật toán chỉ số phản xạ: $$\text{NDVI} = \frac{\text{NIR} - \text{RED}}{\text{NIR} + \text{RED}}$$ $$\text{LSWI} = \frac{\text{NIR} - \text{SWIR}}{\text{NIR} + \text{SWIR}}$$ Kết hợp quy tắc quyết định không gian: Điểm ảnh ngập lụt được xác định khi thỏa mãn điều kiện $LSWI - NDVI > 0$ hoặc $LSWI > 0$ kết hợp với cao trình địa hình từ mô hình SRTM DEM.
  • SPIGenerator v1.0: Tính toán chỉ số chuẩn hóa giáng thủy SPI từ chuỗi số liệu lượng mưa tháng và mưa trích xuất từ vệ tinh TRMM (Tropical Rainfall Measuring Mission - sản phẩm TRMM 3B43) độ phân giải $0,25^\circ \times 0,25^\circ$. Mối tương quan giữa mưa vệ tinh TRMM và mưa quan trắc trạm thực địa đạt hệ số xác định cao ($R^2 > 0,85$), cho phép hiệu chỉnh và nội suy không gian liên tục bằng kỹ thuật Địa thống kê Kriging.
  • Phân tích hồi quy đa biến: Mô hình xác định mức độ tác động của 5 nhóm yếu tố đến biến động diện tích đất nông nghiệp: $$Y = \beta_0 + \beta_1 X_{\text{TuNhien}} + \beta_2 X_{\text{KinhTe}} + \beta_3 X_{\text{XaHoi}} + \beta_4 X_{\text{HaTang}} + \beta_5 X_{\text{KhiHau}} + \epsilon$$ Kết quả kiểm định F và hệ số phóng đại phương sai (VIF < 2) bảo đảm mô hình không vi phạm hiện tượng đa cộng tuyến.

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

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

  1. Biến động đất đai và vị thế của biến đổi khí hậu: Giai đoạn 2015–2019, tổng diện tích đất nông nghiệp của huyện Quảng Điền giảm ròng 124,79 ha, tập trung cao nhất tại các xã Quảng Công, Quảng Lợi, Quảng Ngạn và Quảng An. Mô hình hồi quy đa biến và xếp hạng thang đo Likert chỉ rõ: Biến đổi khí hậu (ngập lụt và hạn hán) là yếu tố tác động mạnh thứ 3 trong 5 nhóm nhân tố chính, trực tiếp thúc đẩy quá trình hoang hóa tạm thời hoặc buộc người dân phải chuyển đổi cơ cấu cây trồng.

"Trong giai đoạn 2015 - 2019, diện tích nông nghiệp của huyện biến động giảm 124,79 ha chủ yếu ở các xã Quảng Công, Quảng Lợi, Quảng Ngạn, xã Quảng An... Yếu tố biến đổi khí hậu, cụ thể là ngập lụt và hạn hán cũng có những tác động đáng kể đến sử dụng đất nông nghiệp trên địa bàn huyện Quảng Điền và yếu tố này ảnh hưởng ở mức độ thứ 3 trong nhóm 5 các yếu tố ảnh hưởng." (Trích tóm tắt luận án).

  1. Cơ chế tác động phân hóa của ngập lụt: Phân tích không gian từ ảnh viễn thám các đợt lũ lịch sử năm 2017 và 2019 cho thấy tính quy luật khốc liệt:
    • Trong mỗi đợt ngập lũ, hơn 80% diện tích bị ngập là đất trồng lúa, tiếp đến là đất trồng cây hàng năm khác (>10%) và đất nuôi trồng thủy sản (2,64%).
    • Không gian ngập tập trung với độ sâu từ 1,0 m đến 2,5 m tại các xã hạ lưu ven sông Bồ và ven phá: Quảng An, Quảng Phước, Quảng Thành và Quảng Thọ.

"Mỗi đợt ngập lụt, hơn 80% diện tích bị ngập là đất trồng lúa, 10% là diện tích đất trồng cây hàng năm khác, 2,64% là diện tích đất nuôi trồng thủy sản. Các vùng trũng thấp ven sông và ven phá chịu mức ngập sâu và thời gian ngập kéo dài nhất." (Trích kết quả nghiên cứu Chương 3).

  1. Cực đoan hạn hán vụ Hè Thu: Đánh giá hạn hán năm 2015 dựa trên phân tích chỉ số SPI chứng minh:
    • Tổng diện tích đất nông nghiệp bị hạn ở mức trung bình chiếm 24,21%, hạn nặng chiếm 12,93% và hạn nhẹ chiếm 11,12%.
    • Trong chuỗi các tháng từ tháng 5 đến tháng 8, tháng 6 được xác định là tháng có chỉ số SPI thấp nhất, tạo đỉnh hạn hán khốc liệt trùng với thời kỳ làm đòng - trổ bông của lúa vụ Hè Thu, gây tổn thất nặng nề cho các xã cát ven biển (Quảng Lợi, Quảng Thái, Quảng Công, Quảng Ngạn).

"Đã xây dựng được bản đồ phân vùng hạn hán năm 2015 dựa vào chỉ số SPI với các mức hạn nhẹ, trung bình và nặng. Trong đó, diện tích đất nông nghiệp bị hạn chủ yếu ở mức trung bình chiếm 24,21%, tiếp đến diện tích bị hạn nặng chiếm 12,93%... từ tháng 5 đến tháng 8 thì tháng 6 được đánh giá là tháng hạn nhất trong vụ Hè Thu." (Trích tóm tắt luận án).

  1. Dự báo không gian ngập lụt đến năm 2030: Dưới kịch bản nước biển dâng thêm 13 cm, diện tích ngập lụt mở rộng thêm hàng trăm hecta. Đất trồng lúa tiếp tục là đối tượng bị nhấn chìm lớn nhất, lan rộng không gian ngập vào sâu các vùng đất chuyên màu và vùng nuôi thủy sản hạ triều tại 4 xã trọng điểm: Quảng An, Quảng Phước, Quảng Thành, Quảng Thọ.

  2. Dự báo không gian hạn hán đến năm 2035: Theo kịch bản nồng độ khí nhà kính trung bình (RCP 4.5), nền nhiệt gia tăng cùng sự sụt giảm lượng mưa mùa khô làm cho cấp độ hạn nặng bao phủ phần lớn diện tích đất lâm nghiệp (rừng phi lao phòng hộ trên cát), đất trồng cây hàng năm và đất lúa tại dải cồn cát ven biển thuộc các xã Quảng Lợi, Quảng Thái, Quảng Công và Quảng Ngạn.

+-----------------------------------------------------------------------------------------+
|                    TỔNG HỢP MA TRẬN TỔN THƯƠNG THEO TIỂU VÙNG ĐỊA LÝ                    |
+------------------------------------+-----------------------------+----------------------+
| Tiểu vùng sinh thái nông nghiệp    | Hiểm họa chính & Xã đại diện| Loại đất tổn thương  |
+------------------------------------+-----------------------------+----------------------+
| 1. Vùng hạ lưu trũng thấp ven sông | Ngập lụt sâu & ngập kéo dài | Đất lúa (LUA >80%),  |
|    và đầm phá nội đồng             | (Quảng An, Quảng Phước,     | Nuôi trồng thủy sản  |
|                                    |  Quảng Thành, Quảng Thọ)    | (NTS), Đất hoa màu   |
+------------------------------------+-----------------------------+----------------------+
| 2. Vùng cồn cát & đồng bằng ven    | Hạn hán gay gắt (Tháng 6),  | Đất lâm nghiệp (LNP),|
|    biển chắn sóng gió              | Thiếu nước tưới vụ Hè Thu   | Đất cây hàng năm     |
|                                    | (Quảng Lợi, Quảng Thái,     | khác (BHK), Đất cát  |
|                                    |  Quảng Công, Quảng Ngạn)    | hoang hóa            |
+------------------------------------+-----------------------------+----------------------+
| 3. Vùng trung gian chuyển tiếp     | Ngập úng cục bộ vụ Đông Xuân| Đất chuyên màu,      |
|    và ven đê bao nội đồng          | & Hạn nhẹ xen kẽ Hè Thu     | Đất lúa tái cơ cấu   |
|                                    | (Quảng Phú, Quảng Vinh, TT) |                      |
+------------------------------------+-----------------------------+----------------------+

Implications đa chiều

  • Về mặt lý thuyết: Bổ sung cơ sở dữ liệu thực nghiệm kiểm chứng mô hình tương tác giữa chỉ số viễn thám bề mặt nước (LSWI/NDVI) và chỉ số khí tượng thủy văn (SPI) trong đánh giá sức tải sinh thái đất.
  • Về mặt phương pháp luận: Thiết lập quy trình công nghệ chuẩn hóa, chi phí thấp, cho phép các địa phương sử dụng ảnh Landsat và dữ liệu TRMM mã nguồn mở để giám sát rủi ro thiên tai định kỳ mà không phụ thuộc vào hệ thống thiết bị đo đạc đắt tiền.
  • Về mặt thực tiễn quy hoạch: Cung cấp luận cứ trực tiếp để UBND huyện Quảng Điền xây dựng Quy hoạch sử dụng đất giai đoạn 2021–2030; phân định ranh giới cách ly an toàn sinh thái; chuyển đổi linh hoạt cơ cấu mùa vụ (tránh tháng 6 cho vụ Hè Thu và thu hoạch sớm trước tháng 10 cho vụ Mùa).
  • Về mặt chính sách: Đề xuất danh mục ưu tiên đầu tư công trình (nâng cấp đê bao nội đồng ven phá Tam Giang, nạo vét trục tiêu thoát lũ sông Bồ, xây dựng hồ trữ nước ngọt vùng cát) kết hợp phi công trình (bảo hiểm nông nghiệp chỉ số thời tiết, khuyến nông chuyển đổi giống chịu mặn/chịu hạn).

Limitations và Future Research

Mặc dù đạt được những kết quả xuất sắc, luận án vẫn tồn tại một số hạn chế mang tính khách quan:

  1. Độ phân giải thời gian và không gian của ảnh quang học Landsat (30m, chu kỳ quay lại 16 ngày) dễ bị gián đoạn bởi mây mù dày đặc trong mùa mưa bão miền Trung, dẫn đến việc nắm bắt thời điểm đỉnh lũ tuyệt đối có độ trễ nhất định.
  2. Dữ liệu lượng mưa vệ tinh TRMM có độ phân giải thô ($0,25^\circ \approx 27\text{ km}$), đòi hỏi các thuật toán nội suy làm mịn có thể chưa phản ánh hết các vi khí hậu cục bộ tại vùng giáp ranh giữa biển và đầm phá.
  3. Luận án chưa lượng hóa đầy đủ tổn thất kinh tế lượng vi mô (micro-econometric damage function) trên từng mét vuông đất bị suy thoái chất lượng đất (độ phì, độ mặn, hóa chua phèn) sau các đợt ngập hạn.

Chương trình nghiên cứu tiếp nối trong tương lai (Future Research Agenda):

  • Tích hợp cảm biến Radar khẩu độ tổng hợp (SAR): Ứng dụng chuỗi ảnh vệ tinh Sentinel-1A/1B (sóng C-band) có khả năng xuyên mây và chụp đêm để theo dõi ngập lụt theo thời gian thực (real-time flood hydrodynamics).
  • Ứng dụng Trí tuệ Nhân tạo & Học máy (Machine Learning): Triển khai các thuật toán Random Forest, XGBoost hoặc Mạng nơ-ron tích chập (CNN) để dự báo không gian hạn hán SPI đa bước thời gian kết hợp dữ liệu độ ẩm đất SMAP.
  • Mô hình hóa kinh tế sinh thái động: Tích hợp mô hình SWAT (Soil and Water Assessment Tool) với mô hình kinh tế nông hộ để đánh giá tác động dây chuyền của chính sách hỗ trợ tín dụng xanh và chuyển dịch cơ cấu cây trồng.

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

  • Tác động học thuật: Thiết lập chuẩn mực nghiên cứu mới trong ngành Quản lý đất đai tại Việt Nam khi kết hợp viễn thám không gian với quản trị tài nguyên; đóng góp tiềm năng vào các chỉ số trích dẫn khoa học quốc tế về quản lý rủi ro thiên tai vùng ven biển Đông Nam Á.
  • Chuyển đổi kinh tế - ngành nông nghiệp: Giúp ngành nông nghiệp huyện Quảng Điền và tỉnh Thừa Thiên Huế giảm thiểu tỷ lệ thất thoát sản lượng lúa hàng năm (ước tính bảo vệ an toàn trên 4.000 ha lúa trũng); chuyển đổi hiệu quả hàng trăm hecta đất cát khô hạn sang mô hình nông nghiệp công nghệ cao trồng rau an toàn, sen cao sản kết hợp nuôi cá xen canh.
  • Tác động chính sách công: Trở thành tài liệu định hướng nền tảng để Sở Tài nguyên và Môi trường tỉnh Thừa Thiên Huế nhân rộng mô hình phân vùng thích ứng khí hậu cho các huyện ven biển lân cận như Phong Điền, Phú Vang, Phú Lộc.
  • Lợi ích xã hội và cộng đồng: Trực tiếp nâng cao năng lực chống chịu (resilience capacity) cho hơn 13.000 hộ nông dân huyện Quảng Điền, bảo vệ sinh kế bền vững trước các dị thường thời tiết ngày càng cực đoan.

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

  • Nghiên cứu sinh và Giới học giả trẻ: Tiếp cận một quy trình phương pháp luận mẫu mực kết hợp GIS, viễn thám, chỉ số khí tượng và phân tích thống kê xã hội học trong quản lý đất đai.
  • Các nhà khoa học cao cấp và Chuyên gia biến đổi khí hậu: Nguồn tư liệu kiểm chứng sinh động về diễn biến ngập hạn vùng hạ lưu các con sông ngắn dốc miền Trung Việt Nam.
  • Cơ quan Quản lý Nhà nước (Sở TN&MT, Sở NN&PTNT, UBND các cấp): Sở hữu bản đồ số hóa phân vùng rủi ro chi tiết đến cấp xã làm công cụ ra quyết định quy hoạch đất và phân bổ ngân sách đầu tư thủy lợi.
  • Hợp tác xã nông nghiệp và Nông dân địa phương: Nắm bắt lịch thời vụ tối ưu và định hướng lựa chọn cơ cấu cây trồng/con nuôi né tránh đỉnh ngập lũ tháng 10–11 và đỉnh hạn hán tháng 6.

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 Sử dụng đất bền vững bằng việc tích hợp cấu trúc đa rủi ro khí hậu (ngập lụt và hạn hán) vào cơ chế định hình giá trị không gian đất nông nghiệp. Luận án chỉ ra rằng sự suy giảm diện tích đất canh tác ở vùng duyên hải không đơn thuần do động lực đô thị hóa mà là quá trình đào thải sinh thái tự nhiên khi các chỉ số cực đoan ($LSWI > 0$ mùa mưa lũ và $SPI < -1,5$ mùa khô) vượt ngưỡng chống chịu của hệ sinh thái nông nghiệp truyền thống.

2. Điểm đổi mới phương pháp luậ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 của Karlsen et al. (2008) (chỉ dùng NDVI) và Niko Wanders et al. (2010) (chỉ phân tích chỉ số khí tượng độc lập), luận án đã tạo đột phá khi ghép nối ma trận không gian đa nguồn: Dùng cặp chỉ số phổ kép $LSWI - NDVI$ từ ảnh Landsat 7/8 kết hợp mô hình độ cao DEM để bóc tách chính xác diện tích mặt nước ngập dưới tán lúa, đồng thời chuẩn hóa chuỗi dữ liệu mưa vệ tinh TRMM bằng mạng lưới trạm quan trắc mặt đất (trạm Phú Ốc) để xây dựng bản đồ hạn hán SPI liên tục trong không gian địa lý.

3. Phát hiện nào gây bất ngờ nhất về mặt số liệu thực nghiệm?
Phát hiện bất ngờ nhất là sự đối lập gay gắt về không gian tổn thất: Trong khi hiểm họa ngập lụt tấn công tập trung vào đất trồng lúa (>80% diện tích ngập) tại các xã phù sa trũng thấp (Quảng An, Quảng Phước, Quảng Thành), thì hiểm họa hạn hán lại tàn phá nặng nề nhất đối với đất lâm nghiệp và đất hoa màu trên dải cồn cát ven biển (Quảng Lợi, Quảng Thái, Quảng Công), với đỉnh hạn rơi cục bộ và khốc liệt vào duy nhất tháng 6 của vụ Hè Thu, chứ không dàn trải đều trong suốt mùa khô.

4. Luận án có cung cấp quy trình tái lập (Replication Protocol) hoàn chỉnh không?
Có. Luận án trình bày chi tiết từng bước thuật toán: Từ tiền xử lý ảnh viễn thám (chỉnh sai hình học, bức xạ, khí quyển trên ENVI), công thức tính toán chỉ số LSWI, NDVI, đến quy trình nhập dữ liệu mưa vào phần mềm SPIGenerator v1.0, các thông số bán kính tìm kiếm trong mô hình nội suy Kriging trên ArcGIS và cấu trúc ma trận hồi quy đa biến trong SPSS. Bất kỳ nhà nghiên cứu nào cũng có thể áp dụng nguyên vẹn giao thức này cho các lưu vực sông ven biển khác.

5. Chương trình nghiên cứu 10 năm tới được phác thảo thế nào?
Chương trình nghiên cứu mở rộng trong thập kỷ tới định hướng: (1) Nâng cấp hệ thống giám sát ngập hạn tự động bằng điện toán đám mây Google Earth Engine; (2) Tích hợp dữ liệu viễn thám Radar Sentinel-1 SAR để giám sát ngập lũ trong điều kiện mưa bão 100% mây phủ; (3) Phát triển các kịch bản thích ứng dựa trên giải pháp sinh thái tự nhiên (Nature-based Solutions - NbS) kết hợp định giá chi phí - lợi ích (CBA) đa chiều cho toàn dải ven biển miền Trung.

Kết luận

  1. Luận án đã giải quyết trọn vẹn và xuất sắc bài toán khoa học về đánh giá tác động kép của ngập lụt và hạn hán đến sử dụng đất nông nghiệp tại địa bàn đặc thù huyện Quảng Điền, tỉnh Thừa Thiên Huế.
  2. Thiết lập cơ sở khoa học thực chứng khẳng định biến đổi khí hậu là nhân tố tác động mạnh thứ 3 làm biến động giảm 124,79 ha đất nông nghiệp giai đoạn 2015–2019.
  3. Ứng dụng thành công công nghệ GIS và viễn thám đa phổ (Landsat, TRMM) xây dựng bộ bản đồ phân vùng ngập lụt (2017, 2019) và hạn hán SPI (2015), lượng hóa chính xác tỷ lệ tổn thương của từng loại đất (đất lúa chiếm >80% diện tích ngập; 24,21% đất bị hạn trung bình và 12,93% bị hạn nặng).
  4. Dự báo không gian chuẩn xác nguy cơ ngập lụt đến năm 2030 (kịch bản nước biển dâng +13 cm) và hạn hán đến năm 2035 (kịch bản RCP 4.5), nhận diện chính xác các xã trọng điểm bị đe dọa.
  5. Đề xuất hệ thống giải pháp công trình và phi công trình đồng bộ, phân hóa cụ thể cho 3 tiểu vùng sinh thái nông nghiệp, cung cấp luận cứ thực tiễn trực tiếp cho công tác lập quy hoạch sử dụng đất cấp huyện giai đoạn 2021–2030.
  6. Mở ra các hướng nghiên cứu liên ngành đột phá, nâng cao chuẩn mực học thuật trong quản lý đất đai và quản trị rủi ro thiên tai ven biển nhiệt đới.