Tổng quan về luận án

Nghiên cứu "Ảnh hưởng của hạn hán đến quản lý và sử dụng đất trồng lúa tại huyện Hòa Vang, thành phố Đà Nẵng" do nghiên cứu sinh Trần Thị Phượng thực hiện dưới sự hướng dẫn khoa học của PGS. TS Huỳnh Văn Chương (Chuyên ngành Quản lý đất đai, Trường Đại học Nông Lâm – Đại học Huế, 2019) là công trình tiên phong tại miền Trung Việt Nam tiếp cận đa ngành, tích hợp khoa học không gian, mô hình khí hậu và khoa học quản trị đất đai. Bối cảnh nghiên cứu đặt trong giai đoạn đô thị hóa bùng nổ từ khi Đà Nẵng trở thành thành phố trực thuộc Trung ương (1997–2016), kéo theo sự suy giảm mạnh mẽ quỹ đất nông nghiệp ven đô cùng áp lực khốc liệt từ biến đổi khí hậu toàn cầu theo các cảnh báo của Kreft (2013, 2017) và Trần Thục, Koos Neefjes (2015).

Khoảng trống nghiên cứu (research gap) trọng yếu mà luận án định vị là: Các nghiên cứu trước đây (như Nguyễn Lập Dân, 2010; Đào Xuân Học, 2001; Nguyễn Hữu Ngữ và cs, 2014) chủ yếu khảo sát hạn hán ở giác độ khí tượng thuần túy hoặc quy hoạch thủy lợi diện rộng, thiếu vắng sự kết nối không gian vi mô giữa hiện tượng hạn hán với 15 nội dung quản lý nhà nước về đất đai quy định tại Điều 22 Luật Đất đai 2013 và hành vi thích ứng của nông hộ. Luận án đặt ra 4 câu hỏi nghiên cứu và 4 giả thuyết khoa học cốt lõi:

  • Câu hỏi 1 (Q1): Đất trồng lúa tại huyện Hòa Vang biến động theo quy luật không gian - thời gian nào trong 20 năm (1997–2016)? Giả thuyết 1 (H1): Diện tích đất lúa suy giảm mạnh chịu sự chi phối đồng thời của đô thị hóa và tần suất hạn hán gia tăng.
  • Câu hỏi 2 (Q2): Mức độ và phân bố không gian của hạn hán biểu hiện ra sao qua chỉ số SPI và dữ liệu viễn thám? Giả thuyết 2 (H2): Hạn hán phân hóa rõ nét theo mùa vụ (tập trung vụ Hè Thu) và theo địa hình từ miền núi đến đồng bằng.
  • Câu hỏi 3 (Q3): Hạn hán tác động định lượng đến các khâu quản lý nhà nước và sử dụng đất trồng lúa ở các cấp như thế nào? Giả thuyết 3 (H3): Hạn hán làm sai lệch chỉ tiêu quy hoạch, gia tăng tranh chấp nguồn nước và làm suy giảm hiệu quả sử dụng đất lúa cấp nông hộ.
  • Câu hỏi 4 (Q4): Giải pháp thích ứng nào khả thi nhất cho từng tiểu vùng sinh thái đến năm 2035? Giả thuyết 4 (H4): Các gói giải pháp chuyển đổi cơ cấu cây trồng theo Nghị định 35/2015/NĐ-CP chỉ đạt hiệu quả cao khi thỏa mãn đồng thời 5 tiêu chí năng lực địa phương.

Khung lý thuyết của công trình tích hợp Lý thuyết Khoa học Thay đổi Đất (Land Change Science - Lambin và cs, 2003), Lý thuyết Quản lý Hạn hán Tích hợp (Wilhite, 2000) và Khung Thích ứng Sinh kế Nông thôn. Phạm vi nghiên cứu bao quát toàn bộ 11 xã thuộc huyện Hòa Vang trong chuỗi thời gian 20 năm (1997–2016), khảo sát thực địa 347 hộ nông dân canh tác lúa và tham vấn chuyên sâu 35 cán bộ quản lý, dự báo không gian hạn hán theo các kịch bản phát thải khí nhà kính (RCP4.5 và RCP8.5) đến năm 2035. Ý nghĩa khoa học và thực tiễn của công trình thể hiện ở việc lượng hóa tác động khí hậu vào các quyết định quản lý đất đai địa phương, cung cấp luận cứ bảo vệ diện tích 3.004,1 ha đất lúa theo phương án quy hoạch của thành phố Đà Nẵng.

Literature Review và Positioning

Tổng quan y văn quốc tế và trong nước ghi nhận ba dòng nghiên cứu chính. Dòng thứ nhất tập trung vào cơ chế phát sinh và đánh giá định lượng hạn hán. Wilhite (2000) và WMO (2006, 2012) phân loại hạn thành 4 nhóm: hạn khí tượng, hạn nông nghiệp, hạn thủy văn và hạn kinh tế - xã hội. Trong việc đo lường, McKee và cs (1993) phát triển Chỉ số chuẩn hóa giáng thủy (SPI), trong khi Palmer (1965) đề xuất Chỉ số mức độ khắc nghiệt hạn Palmer (PDSI), và Tsakiris và cs (2007) đưa ra Chỉ số nhận diện hạn hán chuẩn hóa (RDIstd). Tranh luận khoa học lớn diễn ra giữa quan điểm sử dụng chỉ số đơn biến dựa trên lượng mưa (McKee và cs, 1993) và chỉ số đa biến tích hợp bốc thoát hơi tiềm năng PET (Palmer, 1965; Tsakiris và cs, 2007). Nghiên cứu của Trần Thị Phượng lựa chọn chỉ số SPI đa thời đoạn nhờ tính linh hoạt thống kê và mức độ phù hợp cao khi mạng lưới trạm quan trắc mặt đất tại địa bàn nghiên cứu còn thưa thớt.

Dòng thứ hai nghiên cứu biến động sử dụng đất (LUCC) và áp dụng công nghệ viễn thám, GIS. Brandt Tso và Mather (2009), Yan Gao và Mas (2008), Mohd Hasmadi và cs (2017) đã chỉ rõ sự dịch chuyển từ phương pháp phân loại điểm ảnh truyền thống (Maximum Likelihood, K-means) sang phân loại hướng đối tượng (OBIA) và tích hợp chuỗi ảnh đa thời gian để triệt tiêu hiện tượng "hạt muối tiêu". Tại Việt Nam, các tác giả Hồ Việt Hoàng và cs (2017), Nguyễn Hữu Ngữ và cs (2016), Trịnh Lê Hùng và cs (2017) đã ứng dụng GIS và chuỗi Markov-CA trong phân tích đất đô thị và phi nông nghiệp, nhưng chưa gắn kết chặt chẽ với động thái thời tiết cực đoan.

                  ┌──────────────────────────────────────────────────────────┐
                  │                 Y VĂN QUỐC TẾ & TRONG NƯỚC               │
                  └─────────────────────────────┬────────────────────────────┘
                                                │
         ┌──────────────────────────────────────┼──────────────────────────────────────┐
         ▼                                      ▼                                      ▼
┌─────────────────────────┐        ┌─────────────────────────┐        ┌─────────────────────────┐
│   Đo lường hạn hán      │        │  Viễn thám & GIS (LUCC) │        │  Quản trị & Thích ứng   │
│  - McKee (1993): SPI    │        │  - Brandt Tso (2009)    │        │  - Lambin (2003)        │
│  - Palmer (1965): PDSI  │        │  - Liu (2017), Wu (2017)│        │  - Wilhite (2000)       │
│  - Tsakiris (2007): RDI │        │  - Mansaray (2017)      │        │  - Luật Đất đai 2013    │
└────────────┬────────────┘        └────────────┬────────────┘        └────────────┬────────────┘
             │                                  │                                  │
             └──────────────────────────────────┼──────────────────────────────────┘
                                                ▼
                  ┌──────────────────────────────────────────────────────────┐
                  │           ĐỊNH VỊ LUẬN ÁN (RESEARCH POSITIONING)          │
                  │ - Tích hợp SPI + Viễn thám độ phân giải cao + GIS 3D     │
                  │ - Lượng hóa tác động hạn hán lên 15 nội dung QLNN đất đai │
                  │ - Đánh giá thích ứng 5 chiều theo 3 vùng địa hình        │
                  └──────────────────────────────────────────────────────────┘

Dòng thứ ba khảo sát tác động của hạn hán đến nông nghiệp và quản lý đất đai. So sánh với hai nghiên cứu quốc tế tiêu biểu:

  1. Nghiên cứu của Batar và cs (2017) tại Garhwal Himalayan (Ấn Độ): Sử dụng chuỗi ảnh Landsat (MSS, TM, OLI) và công cụ LFT v2.0 để phân tích phân mảnh rừng và suy giảm thảm phủ trong 38 năm. Luận án của Trần Thị Phượng kế thừa phương pháp phân tích ma trận chuyển dịch nhưng nâng cao độ mịn không gian bằng cách kết hợp ảnh Landsat TM5 (1997) với ảnh vệ tinh độ phân giải cao RapidEye (2016), bổ sung dữ liệu mưa mô phỏng vệ tinh TRMM (Tropical Rainfall Measuring Mission).
  2. Nghiên cứu của Mansaray và cs (2017) tại vùng Thượng Hải (Trung Quốc): Kết hợp radar khẩu độ tổng hợp Sentinel-1A và quang học Landsat 8 để phân tầng canh tác lúa ven đô. Luận án đã phát triển hướng đi này bằng cách đưa biến số hạn hán SPI nội suy không gian (IDW) lồng ghép trực tiếp vào bản đồ hiện trạng đất lúa, giải quyết bài toán mâu thuẫn giữa quy hoạch sử dụng đất và khả năng cấp nước mặt trong điều kiện biến đổi khí hậu.

Đó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 Lý thuyết Động thái Sử dụng Đất (Land Change Theory) của Eric F. Lambin và cs (2003) và Beatriz da Silveira Pinheiro (2003) bằng cách bổ sung biến số "sốc cực đoan khí hậu cục bộ" (localized climate extreme shocks) vào mô hình tương tác kinh tế - xã hội - môi trường. Trước đây, lý thuyết của Lambin nhấn mạnh các động lực thể chế, nhân khẩu và hạ tầng thị trường; công trình này chứng minh rằng hạn hán không chỉ là yếu tố ngoại sinh thụ động mà là động lực trực tiếp tái cấu trúc không gian đất lúa thông qua cơ chế phá vỡ cân bằng chi phí - lợi nhuận nông hộ.

┌─────────────────────────────────────────────────────────────────────────────────────────────────┐
│                           MÔ HÌNH LÝ THUYẾT MỞ RỘNG CỦA LUẬN ÁN                                 │
├─────────────────────────────────────────────────────────────────────────────────────────────────┤
│                                                                                                 │
│   ┌───────────────────────────┐         ┌──────────────────────────┐                            │
│   │   Động lực Khí hậu        │         │   Động lực Kinh tế-Đô thị│                            │
│   │   - Chỉ số SPI vụ Hè Thu  │         │   - Đô thị hóa Đà Nẵng   │                            │
│   │   - Thiếu hụt nước mặt    │         │   - Chuyển dịch việc làm │                            │
│   └─────────────┬─────────────┘         └────────────┬─────────────┘                            │
│                 │ (P1: Tác động trực tiếp)           │ (P2: Tác động biên)                      │
│                 ▼                                    ▼                                          │
│   ┌────────────────────────────────────────────────────────────────┐                            │
│   │       QUYẾT ĐỊNH QUẢN LÝ VÀ SỬ DỤNG ĐẤT TRỒNG LÚA CỦA NÔNG HỘ   │                            │
│   │       - Bỏ hoang vụ Hè Thu / Thu hẹp quy mô diện tích          │                            │
│   │       - Chuyển đổi sang cây trồng cạn / Đề nghị thu hồi đất    │                            │
│   └────────────────────────────────┬───────────────────────────────┘                            │
│                                    │ (P3: Phản hồi chính sách)                                  │
│                                    ▼                                                            │
│   ┌────────────────────────────────────────────────────────────────┐                            │
│   │       ĐIỀU CHỈNH QUẢN LÝ NHÀ NƯỚC VỀ ĐẤT ĐAI (ĐIỀU 22 LUẬT 2013)│                            │
│   │       - Quy hoạch, kế hoạch sử dụng đất (Nghị định 35/2015)    │                            │
│   │       - Cấp GCNQSDĐ, giải quyết tranh chấp nước tưới, tài chính│                            │
│   └────────────────────────────────────────────────────────────────┘                            │
└─────────────────────────────────────────────────────────────────────────────────────────────────┘

Mô hình lý thuyết thiết lập 3 mệnh đề khoa học (Propositions):

  • Mệnh đề P1: Mức độ hạn hán khí tượng (tính bằng SPI) tương tác phi tuyến tính với phân vùng địa hình, tạo ra áp lực đào thải diện tích đất lúa không đồng đều giữa các tiểu vùng.
  • Mệnh đề P2: Tác động tiêu cực của hạn hán làm khuếch đại tốc độ chuyển đổi đất trồng lúa sang đất phi nông nghiệp tại các vùng chịu áp lực đô thị hóa mạnh.
  • Mệnh đề P3: Sự đứt gãy giữa chỉ tiêu kế hoạch sử dụng đất hành chính và điều kiện tự nhiên tạo ra độ trễ quản lý, dẫn đến hiện tượng đất lúa bị bỏ hoang cục bộ nhưng không thể chuyển đổi chính thức.

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

Khung phân tích của đề tài tích hợp liên ngành 3 trụ cột lý thuyết: (1) Thuyết Khí tượng - Thủy văn ứng dụng (đánh giá cán cân nước qua SPI); (2) Thuyết Cảm thám không gian và Địa tin học (OBIA và GIS đa tầng); và (3) Thuyết Thể chế Quản lý Đất đai (State Land Management Framework).

                               ┌─────────────────────────────────────────┐
                               │   KHUNG PHÂN TÍCH ĐỘC ĐÁO CỦA ĐỀ TÀI    │
                               └────────────────────┬────────────────────┘
                                                    │
                 ┌──────────────────────────────────┼──────────────────────────────────┐
                 ▼                                  ▼                                  ▼
   ┌───────────────────────────┐      ┌───────────────────────────┐      ┌───────────────────────────┐
   │    TRỤ CỘT KHÍ TƯỢNG      │      │    TRỤ CỘT ĐỊA TIN HỌC    │      │     TRỤ CỘT THỂ CHẾ       │
   │  - Số liệu 4 trạm mặt đất │      │  - Landsat TM5 (1997)     │      │  - Điều 22 Luật Đất đai   │
   │  - 4 trạm mô phỏng TRMM   │      │  - RapidEye (2016)        │      │  - Nghị định 35/2015/NĐ-CP│
   │  - Tính toán chỉ số SPI   │      │  - Giải đoán & Kiểm định  │      │  - Đánh giá 15 nội dung   │
   └─────────────┬─────────────┘      └─────────────┬─────────────┘      └─────────────┬─────────────┘
                 │                                  │                                  │
                 └──────────────────────────────────┼──────────────────────────────────┘
                                                    ▼
                               ┌─────────────────────────────────────────┐
                               │  MÔ HÌNH CHỒNG GHÉP KHÔNG GIAN - THỜI   │
                               │  GIAN & HỒI QUY TUYẾN TÍNH ĐA BIẾN      │
                               └────────────────────┬────────────────────┘
                                                    ▼
                               ┌─────────────────────────────────────────┐
                               │  ĐÁNH GIÁ TÍNH KHẢ THI THÍCH ỨNG 5 CHIỀU │
                               │  (Tài chính, Kỹ thuật, Lao động,        │
                               │   Quản lý, Hưởng lợi) theo 3 tiểu vùng  │
                               └─────────────────────────────────────────┘

Điểm đột phá phương pháp luận là việc bổ sung 4 trạm mưa ảo TRMM vào mạng lưới 4 trạm quan trắc mặt đất hiện có (Đà Nẵng, Ái Nghĩa, Tam Kỳ, Thượng Nhật), cho phép nội suy không gian IDW đạt hệ số tương quan $R^2 > 0.85$, khắc phục triệt để khoảng trống thiếu trạm khí tượng nội huyện. Điều kiện biên xác lập cho mô hình là vùng nông nghiệp ven đô nhiệt đới gió mùa có chế độ tưới phụ thuộc 100% vào nguồn nước mặt hồi lưu và hệ thống hồ chứa cục bộ.

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

Thiết kế nghiên cứu

Nghiên cứu vận dụng triết lý Thực chứng Phản biện (Critical Realism), kết hợp chặt chẽ giữa phân tích định lượng không gian khách quan và khảo sát thực chứng xã hội học. Thiết kế nghiên cứu đa tầng (multi-level design) vận hành đồng bộ trên 3 cấp độ phân tích: cấp vĩ mô (toàn huyện Hòa Vang và quan hệ liên vùng), cấp trung mô (11 xã phân chia thành 3 vùng địa hình: Miền núi, Trung du, Đồng bằng) và cấp vi mô (nông hộ sản xuất lúa).

Quy mô mẫu điều tra nông hộ được xác định chuẩn xác thông qua công thức Slovin (1984): $$n = \frac{N}{1 + N \cdot e^2}$$ Với tổng thể $N = 2.650$ hộ có đất trồng lúa trực tiếp bị hạn hán tại 11 xã, biên độ sai số chấp nhận $e = 0{,}05$ (độ tin cậy 95%), quy mô mẫu tối thiểu tính toán là 347 hộ. Tác giả áp dụng kỹ thuật chọn mẫu ngẫu nhiên phân tầng (Stratified Random Sampling) phân bổ cụ thể theo 3 vùng địa lý:

  • Nhóm xã Miền núi (Hòa Ninh, Hòa Phú, Hòa Liên): $n_1 = 173$ hộ (chiếm 49,86%).
  • Nhóm xã Trung du (Hòa Nhơn, Hòa Sơn): $n_2 = 111$ hộ (chiếm 31,99%).
  • Nhóm xã Đồng bằng (Hòa Châu, Hòa Tiến): $n_3 = 63$ hộ (chiếm 18,15%).
┌─────────────────────────────────────────────────────────────────────────────────────────────────┐
│                                THIẾT KẾ MẪU VÀ QUY TRÌNH THU THẬP                               │
├─────────────────────────────────────────────────────────────────────────────────────────────────┤
│                                                                                                 │
│   ┌─────────────────────────────────────────────────────────────────────────────────────────┐   │
│   │ TỔNG THỂ KHẢO SÁT: N = 2.650 hộ nông dân canh tác lúa bị hạn hán tại 11 xã               │   │
│   └────────────────────────────────────────────┬────────────────────────────────────────────┘   │
│                                                │ Công thức Slovin (e = 0.05, Độ tin cậy 95%)    │
│                                                ▼                                                │
│   ┌─────────────────────────────────────────────────────────────────────────────────────────┐   │
│   │ MẪU ĐIỀU TRA ĐẠT CHUẨN: n = 347 hộ (Chọn mẫu ngẫu nhiên phân tầng)                      │   │
│   └────────────────────────────────────────────┬────────────────────────────────────────────┘   │
│                                                │                                                │
│         ┌──────────────────────────────────────┼──────────────────────────────────────┐         │
│         ▼                                      ▼                                      ▼         │
│   ┌───────────────────────────┐  ┌───────────────────────────┐  ┌───────────────────────────┐   │
│   │ XÃ MIỀN NÚI (173 hộ)      │  │ XÃ TRUNG DU (111 hộ)      │  │ XÃ ĐỒNG BẰNG (63 hộ)      │   │
│   │ Hòa Ninh, Hòa Phú, HòaLiên│  │ Hòa Nhơn, Hòa Sơn         │  │ Hòa Châu, Hòa Tiến        │   │
│   └───────────────────────────┘  └───────────────────────────┘  └───────────────────────────┘   │
│                                                │                                                │
│                                                ▼                                                │
│   ┌─────────────────────────────────────────────────────────────────────────────────────────┐   │
│   │ THAM VẤN CHUYÊN GIA & CÁN BỘ: 35 phiếu hợp lệ (42 phát ra) + 4 cuộc thảo luận nhóm FGDs │   │
│   │ - Cấp huyện: Phòng TN&MT (2), Phòng NN&PTNT (2), VP Đăng ký đất đai, Thủy lợi           │   │
│   │ - Cấp xã: 31 cán bộ phụ trách địa chính, nông nghiệp, lãnh đạo 11 xã                    │   │
│   └─────────────────────────────────────────────────────────────────────────────────────────┘   │
└─────────────────────────────────────────────────────────────────────────────────────────────────┘

Song song đó, nghiên cứu tiến hành 4 cuộc thảo luận nhóm tập trung (FGD) gồm 1 cuộc cấp huyện và 3 cuộc cấp cụm xã; cùng 35 phiếu tham vấn chuyên sâu đối với cán bộ quản lý (2 lãnh đạo/chuyên viên Phòng TN&MT, 2 cán bộ Phòng NN&PTNT, 31 cán bộ địa chính và lãnh đạo 11 xã).

Quy trình nghiên cứu rigorous

Quy trình giải đoán ảnh viễn thám và xử lý địa không gian được kiểm soát nghiêm ngặt theo các bước:

  1. Tiền xử lý ảnh: Hiệu chỉnh hình học, nắn chỉnh tọa độ về hệ quy chiếu VN-2000 kinh tuyến trục $107^\circ 45'$, hiệu chỉnh bức xạ và khí quyển cho ảnh Landsat TM5 (năm 1997, độ phân giải 30m) và ảnh RapidEye (năm 2016, độ phân giải 5m).
  2. Phân loại ảnh: Áp dụng thuật toán phân loại có kiểm định Maximum Likelihood kết hợp kiểm chứng thực địa bằng thiết bị GPS cầm tay Garmin.
  3. Đánh giá độ chính xác phân loại: Sử dụng ma trận nhầm lẫn (Confusion Matrix). Kết quả giải đoán ảnh năm 1997 đạt độ chính xác tổng thể $88{,}5%$ (Hệ số $\text{Kappa} = 0{,}85$); ảnh năm 2016 đạt độ chính xác tổng thể $92{,}3%$ (Hệ số $\text{Kappa} = 0{,}90$), vượt xa ngưỡng tiêu chuẩn tối thiểu ($80%$) của USGS.
  4. Tính toán chỉ số SPI: Xử lý chuỗi số liệu lượng mưa tháng trong 20 năm bằng thuật toán phân phối xác suất Gamma theo chuẩn WMO, xuất ra các chuỗi giá trị SPI 3 tháng (vụ Đông Xuân và vụ Hè Thu) và SPI 12 tháng.
  5. Nội suy không gian và mô phỏng tương lai: Ứng dụng kỹ thuật IDW trong môi trường phần mềm ArcGIS 10.5 để chồng ghép lớp bản đồ SPI lên bản đồ phân bố các công trình thủy lợi (hồ chứa, đập dâng, trạm bơm) và bản đồ hiện trạng sử dụng đất trồng lúa. Tích hợp kịch bản biến đổi khí hậu phiên bản 2016 của Bộ TN&MT để mô phỏng mức độ rủi ro hạn hán giai đoạn 2016–2035 theo hai kịch bản RCP4.5 và RCP8.5.
┌─────────────────────────────────────────────────────────────────────────────────────────────────┐
│                       QUY TRÌNH XỬ LÝ VIỄN THÁM & ĐỊA TIN HỌC (GIS)                             │
├─────────────────────────────────────────────────────────────────────────────────────────────────┤
│                                                                                                 │
│   ┌─────────────────────────────────────┐         ┌─────────────────────────────────────┐       │
│   │ Dữ liệu ảnh Landsat TM5 (1997, 30m) │         │ Dữ liệu ảnh RapidEye (2016, 5m)     │       │
│   └──────────────────┬──────────────────┘         └──────────────────┬──────────────────┘       │
│                      │                                               │                          │
│                      ▼                                               ▼                          │
│   ┌─────────────────────────────────────────────────────────────────────────────────────┐       │
│   │ Tiền xử lý: Hiệu chỉnh khí quyển, bức xạ, nắn hình học về hệ quy chiếu VN-2000      │       │
│   └──────────────────────────────────────────┬──────────────────────────────────────────┘       │
│                                              ▼                                                  │
│   ┌─────────────────────────────────────────────────────────────────────────────────────┐       │
│   │ Phân loại có kiểm định (Maximum Likelihood) & Đánh giá ma trận nhầm lẫn             │       │
│   │ - Ảnh 1997: Độ chính xác tổng thể = 88.5%, Hệ số Kappa = 0.85                       │       │
│   │ - Ảnh 2016: Độ chính xác tổng thể = 92.3%, Hệ số Kappa = 0.90                       │       │
│   └──────────────────────────────────────────┬──────────────────────────────────────────┘       │
│                                              ▼                                                  │
│   ┌─────────────────────────────────────────────────────────────────────────────────────┐       │
│   │ Tích hợp dữ liệu SPI (4 trạm mặt đất + 4 trạm TRMM) qua thuật toán nội suy IDW      │       │
│   └──────────────────────────────────────────┬──────────────────────────────────────────┘       │
│                                              ▼                                                  │
│   ┌─────────────────────────────────────────────────────────────────────────────────────┐       │
│   │ Chồng ghép không gian (Overlay) trong ArcGIS 10.5: Lớp đất lúa x Lớp SPI x Thủy lợi │       │
│   └──────────────────────────────────────────┬──────────────────────────────────────────┘       │
│                                              ▼                                                  │
│   ┌─────────────────────────────────────────────────────────────────────────────────────┐       │
│   │ Dự báo nguy cơ hạn hán đến năm 2035 theo kịch bản BĐKH RCP4.5 và RCP8.5             │       │
│   └─────────────────────────────────────────────────────────────────────────────────────┘       │
└─────────────────────────────────────────────────────────────────────────────────────────────────┘

Data và phân tích

Dữ liệu xã hội học được kiểm định độ tin cậy thang đo thông qua hệ số Cronbach's Alpha ($\alpha > 0{,}78$ cho tất cả các nhóm chỉ báo đánh giá năng lực thích ứng). Mô hình hồi quy tuyến tính đa biến (Multivariate Linear Regression) được thiết lập trên phần mềm SPSS 22.0 nhằm xác định mức độ tác động của các nhân tố đến sự suy giảm diện tích đất trồng lúa: $$Y = \beta_0 + \beta_1 X_1 + \beta_2 X_2 + \beta_3 X_3 + \beta_4 X_4 + \epsilon$$ Trong đó:

  • $Y$: Biến động diện tích đất trồng lúa;
  • $X_1$: Tần suất và cường độ hạn hán (giá trị SPI);
  • $X_2$: Tốc độ đô thị hóa và phát triển công nghiệp;
  • $X_3$: Chênh lệch thu nhập giữa sản xuất lúa và ngành nghề phi nông nghiệp;
  • $X_4$: Hiệu lực thực thi chính sách và quy hoạch đất lúa địa phương.

Kết quả kiểm định F đạt giá trị thống kê có ý nghĩa ($p < 0{,}001$), hệ số phóng đại phương sai $\text{VIF} < 2{,}1$ khẳng định không xảy ra hiện tượng đa cộng tuyến. Các kiểm định Robustness thông qua việc thay đổi điểm cắt phân vùng hạn hán chứng minh mô hình có độ vững chắc rất cao.

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

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

Thứ nhất, lượng hóa chính xác sự suy giảm nghiêm trọng của quỹ đất trồng lúa. Dữ liệu chiết tách viễn thám và thống kê đất đai chứng minh diện tích đất trồng lúa huyện Hòa Vang giảm mạnh từ 5.124,6 ha (năm 1997) xuống còn 3.910,2 ha (năm 2016), tức mất đi 1.214,4 ha (giảm $23{,}7%$). Đất lúa chuyển đổi phần lớn sang đất phi nông nghiệp phục vụ các dự án hạ tầng đô thị, khu công nghiệp và một phần chuyển thành đất trồng cây hàng năm hoặc bị bỏ hoang do thiếu nước.

   Diện tích đất trồng lúa huyện Hòa Vang (1997 - 2016)
   ┌───────────────────────────────────────────────────────────────────────────────────────┐
   │ 1997: ██████████████████████████████████████████████████ 5.124,6 ha                  │
   │ 2016: ██─────────────────────────────────────────────── 3.910,2 ha (Giảm 23,7%)      │
   │ Chuyển đổi: 1.214,4 ha (Chủ yếu sang đất phi nông nghiệp & bỏ hoang do thiếu nước)    │
   └───────────────────────────────────────────────────────────────────────────────────────┘

Thứ hai, xác lập quy luật không gian và mùa vụ của hạn hán. Phân tích chỉ số SPI trong giai đoạn 1997–2016 chỉ ra rằng hạn hán tại Hòa Vang có tính quy luật mùa vụ cực kỳ nghiêm trọng. Trong khi vụ Đông Xuân hầu như không xuất hiện hạn nặng (SPI dao động từ $-0{,}5$ đến $1{,}2$), thì vụ Hè Thu ghi nhận các đợt hạn khốc liệt với tần suất xuất hiện ngày càng dày đặc, đặc biệt vào các năm 1998, 2002, 2005, 2014 và 2015 với giá trị $\text{SPI} \le -1{,}5$ (hạn nặng đến hạn đặc biệt nghiêm trọng).

Thứ ba, phân hóa mức độ tổn thương theo 3 vùng địa hình sinh thái. Kết quả khảo sát 347 nông hộ cho thấy:

  • Vùng núi (Hòa Ninh, Hòa Phú, Hòa Liên): Tỷ lệ diện tích lúa bị hạn chiếm $38{,}2%$ diện tích gieo trồng vụ Hè Thu do phụ thuộc hoàn toàn vào các khe suối tự nhiên và đập dâng nhỏ có lưu lượng suy kiệt nhanh.
  • Vùng trung du (Hòa Nhơn, Hòa Sơn): Tỷ lệ diện tích bị hạn là $29{,}5%$, chịu tác động kép từ khô hạn nguồn nước và sỏi sạn hóa tầng mặt đất.
  • Vùng đồng bằng (Hòa Châu, Hòa Tiến): Tỷ lệ bị hạn là $18{,}6%$, tuy có hệ thống trạm bơm điện trên sông Yên và sông Cầu Đỏ nhưng thường xuyên bị tê liệt do hiện tượng xâm nhập mặn sâu vào mùa khô.

Thứ tư, phát hiện nghịch lý đứt gãy trong công tác quản lý nhà nước về đất đai. Khảo sát 35 cán bộ quản lý chỉ ra rằng hạn hán tác động tiêu cực mạnh nhất đến 4/15 nội dung quản lý nhà nước theo Điều 22 Luật Đất đai: (1) Lập và quản lý quy hoạch, kế hoạch sử dụng đất (điểm đánh giá $4{,}45/5$); (2) Đăng ký đất đai và cấp Giấy chứng nhận quyền sử dụng đất ($4{,}12/5$); (3) Quản lý việc thu hồi, chuyển mục đích sử dụng đất ($4{,}08/5$); và (4) Giải quyết tranh chấp đất đai và nguồn nước tưới ($3{,}95/5$). Hạn hán dẫn đến tình trạng nông dân tự ý bỏ vụ hoặc tự phát chuyển đổi cây trồng, làm vô hiệu hóa các chỉ tiêu kế hoạch sử dụng đất cấp huyện được phê duyệt.

┌─────────────────────────────────────────────────────────────────────────────────────────────────┐
│        MỨC ĐỘ ẢNH HƯỞNG CỦA HẠN HÁN ĐẾN 4 NỘI DUNG QUẢN LÝ NHÀ NƯỚC VỀ ĐẤT ĐAI (THANG 5 ĐIỂM)   │
├─────────────────────────────────────────────────────────────────────────────────────────────────┤
│                                                                                                 │
│  1. Lập & quản lý quy hoạch, kế hoạch SDĐ:    ███████████████████████████████████████▌ 4.45/5   │
│  2. Đăng ký đất đai & cấp GCNQSDĐ:            ██████████████████████████████████     4.12/5   │
│  3. Thu hồi & chuyển mục đích SDĐ:            ████████████████████████████████       4.08/5   │
│  4. Giải quyết tranh chấp đất đai & nguồn nước:████████────────────────────────────── 3.95/5   │
│                                                                                                 │
└─────────────────────────────────────────────────────────────────────────────────────────────────┘

Thứ năm, mô phỏng nguy cơ hạn hán mở rộng đến năm 2035. Theo kịch bản RCP4.5, diện tích đất trồng lúa có nguy cơ hạn hán vụ Hè Thu tại Hòa Vang sẽ tăng lên $1.420{,}5$ ha ($36{,}3%$ diện tích đất lúa còn lại); theo kịch bản phát thải cao RCP8.5, diện tích bị hạn sẽ lên tới $1.895{,}8$ ha ($48{,}5%$), đe dọa trực tiếp mục tiêu giữ ổn định 3.004,1 ha đất lúa đến năm 2020 và định hướng 2030 của thành phố Đà Nẵng.

Implications đa chiều

  • Về mặt lý luận: Khẳng định tính tất yếu của việc tích hợp biến động khí hậu cực đoan vào mô hình kinh tế đất đai, bổ sung cơ sở khoa học cho lý thuyết quy hoạch sử dụng đất thích ứng (adaptive land use planning).
  • Về phương pháp luận: Thiết lập quy trình chuẩn mực kết hợp giữa chỉ số SPI, dữ liệu viễn thám đa thời gian và công cụ GIS 3D, có thể chuyển giao và áp dụng hiệu quả cho toàn bộ dải đồng bằng duyên hải miền Trung.
  • Về thực tiễn quản lý và chính sách: Đề xuất UBND huyện Hòa Vang và Sở TN&MT Đà Nẵng điều chỉnh linh hoạt cơ chế chuyển đổi cơ cấu cây trồng theo Nghị định 35/2015/NĐ-CP:
    1. Phân vùng miền núi: Chuyển đổi $100%$ diện tích lúa 1 vụ thường xuyên thiếu nước sang mô hình nông lâm kết hợp, trồng cây ăn quả có múi hoặc ngô sinh khối.
    2. Phân vùng trung du: Áp dụng công thức luân canh 1 vụ lúa (Đông Xuân) - 1 vụ hoa màu chịu hạn (lạc, mè, đậu đỗ) trong vụ Hè Thu.
    3. Phân vùng đồng bằng: Giữ vững 2 vụ lúa chuyên canh thông qua đầu tư nâng cấp hệ thống cống ngăn mặn An Trạch và đập dâng giữ ngọt trên sông Yên.

Limitations và Future Research

Nghiên cứu tồn tại một số hạn chế nhất định cần được hoàn thiện trong các giai đoạn tiếp theo:

  1. Hạn chế về mạng lưới quan trắc khí tượng mặt đất: Do trên địa bàn huyện Hòa Vang không có trạm khí tượng chuyên dụng, tác giả phải thu thập dữ liệu từ 4 trạm lân cận và sử dụng số liệu mưa vệ tinh TRMM. Mặc dù hệ số kiểm định tương quan cao, việc thiếu số liệu vi khí hậu cục bộ tại từng thung lũng miền núi có thể tạo ra sai số biên nhất định trong các vi vùng hẹp.
  2. Độ phân giải không gian của ảnh vệ tinh thời kỳ 1997: Ảnh Landsat TM5 với độ phân giải 30m chưa thể hiện được các thửa ruộng lúa bậc thang có diện tích nhỏ hơn 0,09 ha tại các xã Hòa Bắc và Hòa Phú, trong khi ảnh RapidEye 2016 có độ phân giải cao 5m.
  3. Chưa lượng hóa toàn diện tổn thất kinh tế gián tiếp: Nghiên cứu tập trung vào thiệt hại năng suất lúa và biến động diện tích, chưa tính toán sâu giá trị suy thoái dịch vụ hệ sinh thái đất ngập nước và chi phí cơ hội chuyển đổi nghề nghiệp của lao động nông thôn.

Hướng nghiên cứu tương lai:

  • Ứng dụng công nghệ bay không người lái (UAV) tích hợp cảm biến siêu phổ (Hyperspectral) và dữ liệu vệ tinh Radar Sentinel-1/2 thời gian thực để giám sát ẩm độ đất và sức khỏe cây lúa ở cấp độ từng thửa ruộng.
  • Xây dựng mô hình cân bằng nước động (Hydro-economic dynamic model) kết hợp dự báo khí hậu khu vực (CWRF) để tối ưu hóa bài toán phân bổ nguồn nước mặt giữa nông nghiệp, công nghiệp và cấp nước đô thị Đà Nẵng.
  • Mở rộng phân tích kinh tế lượng không gian (Spatial Econometrics) nhằm đánh giá tác động lan tỏa của chính sách dồn điền đổi thửa và tích tụ đất đai trong bối cảnh thích ứng rủi ro khí hậu.

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

Luận án tạo ra tác động học thuật và ứng dụng thực tiễn sâu rộng. Về mặt học thuật, công trình đóng góp bộ cơ sở dữ liệu số hóa hoàn chỉnh về biến động đất lúa và chỉ số hạn hán trong suốt 20 năm, mở ra hướng trích dẫn quan trọng cho các nghiên cứu về địa lý tài nguyên, quy hoạch môi trường và biến đổi khí hậu tại miền Trung.

┌─────────────────────────────────────────────────────────────────────────────────────────────────┐
│                              HỆ THỐNG TÁC ĐỘNG ĐA CHIỀU CỦA LUẬN ÁN                             │
├─────────────────────────────────────────────────────────────────────────────────────────────────┤
│                                                                                                 │
│   ┌───────────────────────────┐         ┌──────────────────────────┐                            │
│   │ GIÁ TRỊ HỌC THUẬT         │         │ TÁC ĐỘNG CHÍNH SÁCH      │                            │
│   │ - Bộ dữ liệu số hóa 20 năm│         │ - Cơ sở điều chỉnh QHHĐ  │                            │
│   │ - Quy trình tích hợp SPI, │         │   3.004,1 ha đất lúa     │                            │
│   │   Viễn thám & GIS 3D      │         │ - Kế hoạch thực hiện NĐ35│                            │
│   └─────────────┬─────────────┘         └────────────┬─────────────┘                            │
│                 │                                    │                                          │
│                 ▼                                    ▼                                          │
│   ┌────────────────────────────────────────────────────────────────┐                            │
│   │       HIỆU QUẢ KINH TẾ - XÃ HỘI & CHUYỂN ĐỔI SẢN XUẤT           │                            │
│   │       - Giảm thiểu 35-40% thiệt hại kinh tế do hạn vụ Hè Thu   │                            │
│   │       - Cải thiện sinh kế cho 2.650 hộ dân nông nghiệp ven đô  │                            │
│   │       - Thúc đẩy các chuỗi giá trị nông nghiệp hữu cơ chịu hạn │                            │
│   └────────────────────────────────────────────────────────────────┘                            │
└─────────────────────────────────────────────────────────────────────────────────────────────────┘

Về mặt chính sách, kết quả nghiên cứu là tài liệu khoa học trực tiếp phục vụ UBND thành phố Đà Nẵng, Sở TN&MT, Sở NN&PTNT trong việc rà soát quy hoạch sử dụng đất giai đoạn 2021–2030, tầm nhìn 2050; cung cấp căn cứ phân bổ nguồn vốn đầu tư công cho các công trình cấp bách phòng chống hạn theo Quyết định số 158/2008/QĐ-TTg và Chỉ thị 04/CT-TTg của Thủ tướng Chính phủ.

Về mặt xã hội, các gói giải pháp chuyển đổi mùa vụ giúp hơn 2.650 hộ dân tại 11 xã chủ động giảm thiểu từ $35%$ đến $40%$ thiệt hại kinh tế trong các mùa khô hạn cực đoan, bảo đảm ổn định đời sống xã hội nông thôn trong quá trình đô thị hóa.

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

  • Nghiên cứu sinh và Giới học thuật: Tiếp cận quy trình kết nối dữ liệu địa không gian đa độ phân giải với phân tích chính sách quản lý đất đai; khai thác phương pháp luận đánh giá tính khả thi thích ứng 5 chiều (tài chính, kỹ thuật, lao động, quản lý, hưởng lợi).
  • Cơ quan Quản lý Nhà nước (Sở TN&MT, Phòng TN&MT): Sở hữu bản đồ số hóa vùng nguy cơ hạn hán tỷ lệ 1/25.000, hỗ trợ ban hành quyết định cho phép chuyển mục đích sử dụng đất hoặc lập kế hoạch điều tiết nước chính xác.
  • Ngành Nông nghiệp và PTNT, Công ty Khai thác Thủy lợi: Xác định chính xác các điểm nghẽn thủy lợi và các vùng trạm bơm bị nhiễm mặn để bố trí lịch gieo cấy vụ Hè Thu linh hoạt.
  • Nông dân và Hợp tác xã Nông nghiệp tại huyện Hòa Vang: Được định hướng các mô hình cây trồng cạn có giá trị kinh tế cao (lạc, rau hoa công nghệ cao), giúp tăng thu nhập trên một đơn vị diện tích gấp 1,5–2 lần so với độc canh cây lúa bấp bênh nguồ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 Động thái Thay đổi Đất (Land Change Theory) của Eric F. Lambin và cs (2003) bằng việc xác lập mô hình phản hồi thể chế - sinh kế dưới tác động của các cú sốc khí hậu cục bộ. Thay vì chỉ xem hạn hán là biến ngoại cảnh đơn thuần, công trình chứng minh hạn hán vụ Hè Thu đóng vai trò là "chất xúc tác gia tốc" (accelerating catalyst) làm bộc lộ các mâu thuẫn thể chế trong công tác quản lý nhà nước về đất đai (Điều 22 Luật Đất đai 2013), thúc đẩy sự dịch chuyển cơ cấu sử dụng đất ven đô nhanh hơn tốc độ quy hoạch dự kiến.

2. Đổ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 tiền nhiệm? So sánh với các nghiên cứu của Nguyễn Hữu Ngữ và cs (2014) tại Quảng Nam (chỉ sử dụng nội suy SPI 4 trạm quan trắc mặt đất) và Batar và cs (2017) tại Ấn Độ (chỉ dùng viễn thám Landsat), luận án tạo đột phá bằng phương pháp tích hợp 3 cấp độ:

  • Bổ sung 4 trạm mưa ảo TRMM vào mạng lưới quan trắc, nâng hệ số xác định không gian $R^2$ từ $0{,}65$ lên $0{,}88$;
  • Kết hợp đa nguồn ảnh viễn thám (Landsat TM5 phân giải 30m cho năm 1997 và RapidEye phân giải 5m cho năm 2016), đạt hệ số tin cậy $\text{Kappa} \ge 0{,}85$;
  • Chồng ghép trực tiếp lớp nguy cơ hạn hán với dữ liệu khảo sát vi mô 347 hộ và 15 nội dung quản lý đất đai.

3. Phát hiện bất ngờ nhất có sự hỗ trợ của dữ liệu thực nghiệm là gì? Nghịch lý tại vùng đồng bằng huyện Hòa Vang (Hòa Châu, Hòa Tiến): Mặc dù là vùng có hệ thống trạm bơm điện kiên cố và gần các trục sông lớn (sông Yên, sông Cầu Đỏ), tỷ lệ thiệt hại sản xuất lúa vụ Hè Thu trong các năm hạn nặng lại không giảm tương ứng so với vùng trung du. Dữ liệu thực nghiệm chứng minh rằng hiện tượng suy giảm lưu lượng nước ngọt thượng nguồn kết hợp thủy triều đẩy mặn sâu vào đất liền làm độ mặn vượt ngưỡng chịu đựng của lúa ($> 1‰$), biến hệ thống trạm bơm thành các công trình "vô hiệu hóa" trong thời điểm đỉnh hạn.

┌─────────────────────────────────────────────────────────────────────────────────────────────────┐
│              CƠ CHẾ NGHỊCH LÝ THIỆT HẠI HẠN - MẶN VÙNG ĐỒNG BẰNG VEN ĐÔ                         │
├─────────────────────────────────────────────────────────────────────────────────────────────────┤
│                                                                                                 │
│   Nắng nóng cực đoan (Vụ Hè Thu) ──► Suy giảm lưu lượng dòng chảy ngọt thượng nguồn             │
│                                                   │                                             │
│                                                   ▼                                             │
│   Thủy triều đẩy nước biển sâu vào sông Yên/Cầu Đỏ ──► Độ mặn nguồn nước vượt ngưỡng (> 1‰)    │
│                                                   │                                             │
│                                                   ▼                                             │
│   Hệ thống trạm bơm điện ngừng vận hành ──► Đất lúa đồng bằng bị khô hạn & nhiễm mặn nghiêm trọng│
│                                                                                                 │
└─────────────────────────────────────────────────────────────────────────────────────────────────┘

4. Giao thức lặp lại nghiên cứu (Replication Protocol) có được cung cấp đầy đủ không? Luận án cung cấp giao thức chuẩn mực, minh bạch và hoàn toàn có khả năng tái lập:

  • Công thức tính toán toán học chi tiết cho chỉ số SPI (phân phối xác suất Gamma);
  • Bảng tham số phân loại ảnh viễn thám, tọa độ các điểm kiểm chứng mặt đất GCPs, ma trận kiểm định sai số nhầm lẫn;
  • Bảng câu hỏi phỏng vấn chuẩn hóa 347 hộ (Phụ lục 07) và khung phỏng vấn sâu 35 chuyên gia;
  • Bộ tiêu chí lượng hóa 5 bậc đánh giá tính khả thi của các giải pháp thích ứng.

5. Chương trình nghị sự nghiên cứu 10 năm (10-year Research Agenda) được vạch ra như thế nào? Tác giả định hình lộ trình nghiên cứu 10 năm gồm 3 trọng tâm:

  • Giai đoạn 2020–2023: Xây dựng hệ thống WebGIS cảnh báo sớm hạn hán và xâm nhập mặn phục vụ điều hành sản xuất lúa thời gian thực cho thành phố Đà Nẵng;
  • Giai đoạn 2024–2027: Mô hình hóa tác động tích hợp của đô thị hóa và kịch bản nước biển dâng lên thị trường chuyển nhượng quyền sử dụng đất nông nghiệp;
  • Giai đoạn 2028–2030: Hoàn thiện khung thể chế kinh tế tuần hoàn nước nông nghiệp và thí điểm chi trả dịch vụ môi trường đối với đất trồng lúa thích ứng khí hậu.

Kết luận

  1. Công trình đã lượng hóa chính xác sự biến động diện tích đất trồng lúa tại huyện Hòa Vang trong 20 năm (1997–2016), ghi nhận mức suy giảm 1.214,4 ha ($23{,}7%$) dưới áp lực kép của đô thị hóa và khô hạn.
  2. Ứng dụng thành công chỉ số chuẩn hóa giáng thủy (SPI) kết hợp chuỗi ảnh viễn thám Landsat TM5 và RapidEye, xác định rõ quy luật hạn hán tập trung khốc liệt vào vụ Hè Thu với mức độ phân hóa sâu sắc theo 3 vùng địa hình sinh thái.
  3. Chỉ ra 4 nhóm yếu tố cốt lõi chi phối biến động đất lúa gồm: hạn hán, thể chế chính sách quản lý, chênh lệch thu nhập và tốc độ đô thị hóa; chứng minh hạn hán tác động tiêu cực trực tiếp đến 4/15 nội dung quản lý nhà nước về đất đai.
  4. Dự báo không gian rủi ro hạn hán đến năm 2035 theo hai kịch bản biến đổi khí hậu RCP4.5 (nguy cơ hạn $36{,}3%$) và RCP8.5 (nguy cơ hạn $48{,}5%$), cung cấp cơ sở khoa học để giữ vững diện tích 3.004,1 ha đất lúa theo quy hoạch thành phố Đà Nẵng.
  5. Xây dựng ma trận giải pháp thích ứng khả thi 5 chiều (tài chính, kỹ thuật, lao động, quản lý, hưởng lợi) đặc thù cho từng tiểu vùng địa hình, tạo bước đột phá trong thực thi chuyển đổi cơ cấu cây trồng theo Nghị định 35/2015/NĐ-CP.
  6. Mở ra 3 hướng nghiên cứu mới: Ứng dụng viễn thám siêu phổ thời gian thực trong quản trị tài nguyên đất; Mô hình hóa kinh tế - thủy văn liên vùng; và Cơ chế bảo hiểm rủi ro khí hậu đối với đất trồng lúa ven đô.