Tổng quan về luận án

Sự biến đổi sâu sắc của bề mặt Trái đất dưới tác động đồng thời của hoạt động nhân sinh và biến đổi khí hậu (BĐKH) đang đặt ra thách thức sống còn đối với an ninh tài nguyên và phát triển bền vững toàn cầu. Luận án tiến sĩ chuyên ngành Quản lý Tài nguyên và Môi trường (Mã số: 9.01) của nghiên cứu sinh Lê Ngọc Lãm, dưới sự hướng dẫn khoa học của PGS.TS Lê Văn Trung tại Trường Đại học Nông Lâm TP. Hồ Chí Minh (2024), mang tiêu đề "Giám sát thay đổi sử dụng đất trong bối cảnh biến đổi khí hậu tỉnh Bến Tre", đại diện cho một công trình nghiên cứu tiên phong trong việc lượng hóa tương tác phức hợp giữa tai biến khí hậu và động thái chuyển dịch không gian đất đai tại một trong những vùng đồng bằng nhạy cảm sinh thái bậc nhất thế giới.

Khoảng trống nghiên cứu (research gap) cốt lõi xuất phát từ thực trạng: các phương pháp quản lý đất đai truyền thống tại Việt Nam chủ yếu dựa vào công tác thống kê, kiểm kê định kỳ 5 năm một lần theo đơn vị hành chính. Phương thức này bộc lộ độ trễ thông tin lớn, mang tính tĩnh và không phản ánh kịp thời các biến động tức thời do thiên tai cực đoan. Về mặt lý luận học thuật, mặc dù mối quan hệ giữa biến đổi sử dụng đất/lớp phủ mặt đất (LULC) và BĐKH đã được ghi nhận rộng rãi (Foley et al., 2005; Turner et al., 2007), phần lớn các nghiên cứu trong nước chỉ dừng lại ở việc liệt kê định tính các tác nhân hoặc phân tích biến động đơn lẻ mà thiếu khung phương pháp lượng hóa mức độ đóng góp của từng thông số khí hậu cực đoan (nhiệt độ, lượng mưa, xâm nhập mặn, ngập lụt, khô hạn) vào xác suất chuyển đổi không gian của từng loại hình sử dụng đất cụ thể.

Luận án thiết lập hệ thống 4 câu hỏi nghiên cứu và 4 giả thuyết khoa học tương ứng:

  • Câu hỏi nghiên cứu 1 (RQ1): Quy luật chu chuyển và xu thế biến động không gian của các loại hình sử dụng đất chính tại vùng ven biển tỉnh Bến Tre giai đoạn 2009–2019 diễn ra như thế nào?
    Giả thuyết 1 (H1): Diện tích đất trồng lúa truyền thống suy giảm mạnh do xu thế chuyển dịch tất yếu sang đất nuôi trồng thủy sản và đất trồng cây lâu năm dưới áp lực của các yếu tố sinh thái mới.
  • Câu hỏi nghiên cứu 2 (RQ2): Mức độ tương quan và trọng số tác động của các biến số khí hậu cực đoan đến tiềm năng chuyển đổi từng loại đất được lượng hóa ra sao?
    Giả thuyết 2 (H2): Độ mặn trong đất (đo bằng độ dẫn điện EC) và chỉ số khô hạn (TVDI) là hai biến độc lập giữ vai trò chi phối có ý nghĩa thống kê cao nhất ($p < 0.01$) đối với sự suy giảm diện tích canh tác lúa.
  • Câu hỏi nghiên cứu 3 (RQ3): Mô hình tích hợp không gian - toán học nào có khả năng mô phỏng và dự báo chính xác cấu trúc LULC tương lai dưới các kịch bản BĐKH?
    Giả thuyết 3 (H3): Cấu trúc tích hợp giữa chuỗi Markov, mạng nơ-ron nhân tạo (ANN/MLP) và mạng tế bào tự động (Cellular Automata - CA) trong môi trường LCM và MOLUSCE cung cấp độ chính xác tổng thể và chỉ số Kappa vượt trội so với các mô hình ngoại suy tuyến tính cổ điển.
  • Câu hỏi nghiên cứu 4 (RQ4): Cần thiết lập giải pháp công nghệ và quy hoạch nào để kiểm soát không gian đất đai thích ứng với các kịch bản BĐKH dài hạn đến năm 2050?
    Giả thuyết 4 (H4): Hệ thống giám sát dựa trên nền tảng WebGIS mã nguồn mở tích hợp dữ liệu viễn thám đa thời gian cho phép phát hiện sớm biến động và hỗ trợ tái cấu trúc mùa vụ thích ứng hiệu quả với kịch bản phát thải trung bình (RCP 4.5).

Khung lý thuyết của công trình được xây dựng trên sự giao thoa giữa Khoa học Thay đổi Đất đai (Land Change Science), Lý thuyết Vị trí Không gian (Spatial Location Theory) và Lý thuyết Hệ thống Thích ứng Phức hợp Sinh thái - Xã hội (Social-Ecological Complex Adaptive Systems). Phạm vi không gian của luận án tập trung vào 3 huyện duyên hải xung yếu của tỉnh Bến Tre: Ba Tri, Bình Đại và Thạnh Phú. Khung thời gian thực nghiệm bao quát chuỗi số liệu viễn thám và khí tượng 1999–2009–2019, đồng thời thiết lập mô phỏng động thái đất đai cho các chu kỳ tương lai 2029, 2039, 2049 và định hướng chiến lược đến năm 2050.

Literature Review và Positioning

Nghiên cứu sử dụng đất và động thái thay đổi lớp phủ mặt đất (LULC) đã hình thành ba trường phái học thuật chính trên trường quốc tế:

  1. Trường phái mô hình hóa thống kê - thực nghiệm (Empirical-Statistical Modeling): Tiêu biểu là các nghiên cứu của Verburg et al. với mô hình CLUE-S và phiên bản nâng cấp CLUMondo, tập trung vào việc cân bằng giữa nhu cầu sử dụng đất phi không gian (non-spatial demands) và việc phân bổ không gian (spatial allocation) dựa trên hàm hồi quy Logistic đa biến.
  2. Trường phái mô hình hóa toán học rời rạc và ngẫu nhiên (Discrete & Stochastic Modeling): Điển hình là lý thuyết chuỗi Markov kết hợp Mạng tế bào tự động (Markov-CA) của Banos (2001) và Pontius et al. (2004), cho phép mô hình hóa sự tương tác không gian cục bộ (neighborhood interactions) và xác suất chuyển dịch trạng thái theo thời gian.
  3. Trường phái mô hình hóa dựa trên tác tử (Agent-Based Modeling - ABM): Điển hình với các công trình của Castella (2005), Parker (2012), Taillandier et al. (2014) trên nền tảng GAMA, mô phỏng hành vi ra quyết định vi mô của nông hộ và chủ thể quản lý tác động lên cảnh quan vĩ mô.

Trong bức tranh tổng quan đó, tồn tại một cuộc tranh luận học thuật sâu sắc giữa hai luồng quan điểm đối lập:

  • Luồng quan điểm thứ nhất (Virginia, 1997) cho rằng: các động lực kinh tế - xã hội, thể chế và thị trường giữ vai trò áp đảo trong việc định hình LULC, trong khi biến đổi khí hậu chỉ đóng vai trò thứ yếu hoặc là bối cảnh nền.
  • Luồng quan điểm thứ hai (Nyatuame & Agodzo, 2017; IPCC AR4, AR5) lập luận rằng: tại các vùng đồng bằng ven biển và cửa sông trũng thấp, các hiện tượng khí hậu cực đoan (nước biển dâng, hạn hán, xâm nhập mặn) đã vượt qua ngưỡng chống chịu sinh thái, trở thành động lực trực tiếp và mang tính quyết định cưỡng bức sự biến đổi của hệ số sử dụng đất.

Luận án này định vị học thuật chính xác vào điểm giao cắt của cuộc tranh luận trên: tại vùng cửa sông ven biển Bến Tre thuộc Đồng bằng sông Cửu Long (ĐBSCL), các cú sốc khí hậu đóng vai trò "chất xúc tác cưỡng bức", thúc đẩy quá trình chuyển đổi sinh kế và cơ cấu đất đai của con người.

Khi đặt cạnh các công trình quốc tế cùng lĩnh vực:

  • Nghiên cứu của Varun et al. (2014) tại thành phố Muzaffarpur (Ấn Độ) áp dụng mô hình Land Change Modeler (LCM) dựa trên chuỗi số liệu 1988–2010 để dự báo không gian đô thị đến năm 2025 và 2035, song chỉ thuần túy khai thác các biến nhân sinh và khoảng cách hạ tầng giao thông.
  • Nghiên cứu của P. et al. tại Banjarmasin (Indonesia) sử dụng IDRISI LCM và chuỗi Markov (2014–2018) chứng minh 40% đất nông nghiệp chuyển thành đất xây dựng, nhưng chưa tích hợp biến động độ mặn đất (EC) và nhiệt độ bề mặt đất (LST).
  • Nghiên cứu của Noortje (2013) tại Tomohon (Indonesia) tích hợp Quy trình Phân tích Thứ bậc (AHP) phân tích chuyển đổi đất nông nghiệp nhưng dừng lại ở mô hình định tính tĩnh.
  • Hệ thống hỗ trợ ra quyết định CLUES tại New Zealand (Sandy, 2008) liên kết mô hình rửa trôi nitrat với LULC cấp quốc gia nhưng đòi hỏi hạ tầng dữ liệu quan trắc vi mô rất tốn kém.

Luận án của Lê Ngọc Lãm tạo ra bước tiến vượt bậc khi tiên phong giải quyết bài toán thiếu hụt trạm quan trắc mặt đất tại các nước đang phát triển bằng cách "vật lý hóa" dữ liệu viễn thám: chuyển đổi ảnh vệ tinh Landsat đa thời gian thành các chỉ số hạn hán TVDI và bản đồ độ mặn EC thông qua các thuật toán hồi quy quang phổ, sau đó nạp trực tiếp vào cấu trúc mô hình LCM và MOLUSCE để điều khiển quá trình mô phỏng không gian đa chiề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 Vị trí Sử dụng Đất cổ điển (Von Thünen Land Rent Theory) và Khoa học Thay đổi Đất đai (Land Change Science) bằng cách đưa các yếu tố rủi ro sinh thái khí hậu cực đoan thành các biến số nội sinh trong hàm quyết định sử dụng đất.

Văn bản luận án khẳng định nền tảng pháp lý và khoa học:

"Đất đai là tài nguyên đặc biệt của quốc gia, nguồn lực quan trọng phát triển đất nước, được quản lý theo pháp luật." (Hiến pháp Việt Nam năm 2013)

Đồng thời chuẩn hóa định nghĩa nền tảng:

"Đất đai bao gồm môi trường vật chất, bao gồm khí hậu, đất đai, thủy văn và thảm thực vật, ở mức độ ảnh hưởng đến tiềm năng sử dụng đất." (FAO, 1976)

và:

"Biến động sử dụng đất được biết như biến động đất đai, đây là một thuật ngữ chung chỉ những thay đổi bề mặt lãnh thổ trái đất xảy ra do tác động của con người." (Từ điển Khoa học Trái đất)

Các mệnh đề lý thuyết (Theoretical Propositions) được xác lập:

  • Mệnh đề 1 (P1): Xác suất chuyển đổi ($P_{ij}$) của một ô không gian từ loại hình sử dụng đất $i$ (ví dụ: đất lúa) sang loại hình $j$ (ví dụ: đất nuôi trồng thủy sản) tỷ lệ thuận với mức độ gia tăng tích lũy của độ dẫn điện đất (EC) và chỉ số khô hạn (TVDI), vượt qua ngưỡng thích nghi sinh học của cây trồng.
  • Mệnh đề 2 (P2): Động thái LULC tại dải ven biển châu thổ không diễn ra theo trật tự tuyến tính mà tuân theo mô hình biến động phi tuyến cục bộ, được dẫn dắt bởi sự tương tác lân cận (neighborhood effect) trong mạng tế bào tự động (CA) dưới sức ép của kịch bản biến đổi khí hậu.
  • Mệnh đề 3 (P3): Quá trình chuyển dịch cơ cấu cây trồng là phản ứng sinh kế thích ứng có điều kiện; khi biên độ xâm nhập mặn vượt ngưỡng chịu đựng của hệ thống ngọt hóa, quá trình chuyển đổi sang đất chuyên tôm hoặc lúa - tôm mang tính chất không thể đảo nghịch trong ngắn hạn.

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

Khung phân tích của luận án tích hợp liền mạch ba trụ cột lý thuyết: (1) Lý thuyết Chuyển tiếp Markov (Markovian Transition Theory); (2) Lý thuyết Tương tác Không gian Tế bào Tự động (Cellular Automata Dynamics); và (3) Lý thuyết Học máy Nhận dạng Mẫu Cảnh quan thông qua Mạng Nơ-ron Nhân tạo Đa lớp (Multilayer Perceptron Neural Network - MLP-NN).

Phương pháp phân tích mang tính đột phá thể hiện ở quy trình 5 giai đoạn:

  1. Không gian hóa các chỉ số khí hậu: Sử dụng dữ liệu phổ phản xạ từ các kênh cận hồng ngoại (NIR), hồng ngoại sóng ngắn (SWIR) và hồng ngoại nhiệt (TIR) của vệ tinh Landsat để chiết xuất các chỉ số thực vật chuẩn hóa (NDVI), chỉ số khác biệt độ mặn (NDSI), nhiệt độ bề mặt đất (LST), từ đó thiết lập ma trận không gian chỉ số khô hạn theo quan hệ nhiệt độ - thực vật (TVDI) và bản đồ độ dẫn điện mặn (EC).
  2. Lượng hóa tiềm năng chuyển đổi: Ứng dụng chỉ số Cramer’s V và mô hình Hồi quy Logistic nhị phân đa biến để sàng lọc, loại bỏ hiện tượng đa cộng tuyến và đánh giá sức mạnh giải thích của từng biến điều kiện tự nhiên đối với biến động đất đai.
  3. Tính toán xác suất chuyển dịch: Thiết lập ma trận chu chuyển đất đai và xác suất quá độ Markov cho giai đoạn 2009–2019.
  4. Phân bổ không gian bằng mạng nơ-ron: Ứng dụng thuật toán MLP trong mô hình MOLUSCE và LCM để huấn luyện dữ liệu, tạo lập bản đồ bề mặt tiềm năng chuyển đổi (Transition Potential Maps).
  5. Mô phỏng động và kiểm định: Khởi tạo mô phỏng CA cho các mốc thời gian 2029, 2039, 2049 và đối chiếu với kịch bản biến đổi khí hậu RCP 4.5 định hướng 2050.

Điều kiện biên (Boundary conditions): Khung phân tích được tối ưu hóa cho các vùng đồng bằng ngập lũ và ven biển nhiệt đới gió mùa, nơi có địa hình bằng phẳng, cao độ thấp (< 2m so với mực nước biển), chịu ảnh hưởng mạnh của chế độ bán nhật triều và xâm nhập mặn theo mùa.

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 hành trên nền tảng triết học Thực chứng Luận (Positivism) kết hợp với Thuyết Hiện thực Không gian Định lượng (Spatial Quantitative Realism), sử dụng phương pháp luận hỗn hợp (mixed-methods) lấy phân tích không gian số lượng lớn làm trung tâm và hiệu chỉnh bằng điều tra thực địa xã hội học sinh kế.

Thiết kế nghiên cứu đa cấp độ (multi-level design) được triển khai trên 3 tầng bậc:

  • Cấp độ vi mô (Pixel level): Ô lưới không gian độ phân giải $30\text{m} \times 30\text{m}$ (đồng bộ giữa dữ liệu ảnh Landsat, DEM và các biến giải thích).
  • Cấp độ trung mô (Landscape/District level): Phạm vi ranh giới 3 huyện duyên hải Ba Tri, Bình Đại, Thạnh Phú với đặc trưng sinh thái ngọt - lợ - mặn phân hóa rõ rệt.
  • Cấp độ vĩ mô (Regional level): Toàn bộ cấu trúc địa bàn tỉnh Bến Tre trong mối liên kết vùng hạ lưu sông Tiền và sông Hàm Luông thuộc ĐBSCL.

Quy trình nghiên cứu rigorous

Quy trình thu thập và xử lý dữ liệu được kiểm soát nghiêm ngặt theo các tiêu chuẩn viễn thám và trắc địa quốc tế:

  • Nguồn dữ liệu vệ tinh: Ảnh Landsat 5 TM (thời điểm năm 1999, 2009) và Landsat 8 OLI/TIRS (thời điểm năm 2019) được thu nhận trong cùng mùa khô (nhằm hạn chế tối đa độ che phủ của mây và bắt trọn đỉnh điểm xâm nhập mặn/khô hạn). Ảnh được nắn chỉnh hình học, hiệu chỉnh bức xạ khí quyển (Atmospheric Correction) và chuyển đổi giá trị bức xạ số (DN) sang nhiệt độ bức xạ bề mặt (Brightness Temperature) và độ phản xạ bề mặt (Surface Reflectance).
  • Thu thập dữ liệu thực địa: Thiết lập bộ dữ liệu mẫu huấn luyện (Region of Interest - ROI) và mẫu kiểm chứng độc lập phân bố đều khắp các loại hình sử dụng đất tại 3 huyện nghiên cứu. Đo đạc mẫu độ mặn tầng đất mặt ($0 - 20\text{cm}$) để xác định giá trị độ dẫn điện $EC\ (\text{dS/m})$ bằng thiết bị chuyên dụng, phục vụ xây dựng phương trình tương quan phản xạ phổ.
  • Kỹ thuật tam giác đạc dữ liệu (Data Triangulation): Kết hợp chéo giữa số liệu thống kê kiểm kê đất đai định kỳ của Sở Tài nguyên và Môi trường tỉnh Bến Tre, số liệu quan trắc khí tượng thủy văn (nhiệt độ, lượng mưa, mực nước trạm Ba Tri) và dữ liệu giải đoán ảnh viễn thám.
  • Đánh giá độ tin cậy và giá trị phân loại: Sử dụng thuật toán Phân loại Cực đại Gần nhất (Maximum Likelihood Classification - MLC) kết hợp giải đoán bằng mắt có sự hỗ trợ của chuyên gia. Đánh giá độ chính xác thông qua Ma trận Sai số (Confusion Matrix). Kết quả giải đoán ảnh năm 2009 và 2019 đều đạt Độ chính xác Tổng thể (Overall Accuracy) $> 85%$ và Hệ số Kappa Coefficient $> 0.82$, đáp ứng tiêu chuẩn khắt khe trong mô hình hóa LULC.

Data và phân tích

Quy trình xử lý số liệu thống kê không gian và học máy được triển khai trên nền tảng phần mềm mã nguồn mở QGIS kết hợp plugin MOLUSCE và phần mềm chuyên dụng IDRISI TerrSet (LCM):

  1. Mô hình Hồi quy Logistic Nhị phân Đa biến (Binary Logistic Regression): Thiết lập phương trình hàm logit: $$\ln\left(\frac{P}{1 - P}\right) = \beta_0 + \beta_1 X_1 + \beta_2 X_2 + \dots + \beta_k X_k$$ trong đó biến phụ thuộc là sự xuất hiện (1) hoặc không xuất hiện (0) của biến động chuyển loại đất; các biến độc lập giải thích ($X_i$) gồm: Nhiệt độ trung bình năm, Lượng mưa trung bình tháng, Khoảng cách đến mạng lưới sông ngòi/bờ biển, Mức độ ngập lụt theo mô hình DEM, Chỉ số khô hạn TVDI, và Chỉ số độ mặn EC.
  2. Kiểm tra hiện tượng Đa cộng tuyến (Multicollinearity Test): Tính toán ma trận hệ số tương quan Pearson và Hệ số Phóng đại Biến thiên (Variance Inflation Factor - VIF). Toàn bộ các biến được lựa chọn đưa vào mô hình cuối cùng đều có hệ số $\text{VIF} < 2.5$, loại trừ hoàn toàn nguy cơ sai lệch ước lượng do đa cộng tuyến.
  3. Đánh giá tiềm năng chuyển đổi qua Cramer’s V: Xác định chỉ số Cramer’s V cho từng lớp dữ liệu nhân tố; các biến có giá trị $\text{V} > 0.15$ được xác nhận có lực phân biệt tốt và được nạp vào mạng nơ-ron.
  4. Cấu trúc Mạng Nơ-ron Nhân tạo (ANN/MLP): Mạng nơ-ron truyền thẳng đa lớp trong LCM/MOLUSCE được cấu hình với: Tầng vào (Input layer) gồm các nơ-ron tương ứng với các biến khí hậu và không gian; Tầng ẩn (Hidden layer) sử dụng hàm kích hoạt Sigmoidal phi tuyến; Tầng ra (Output layer) xuất bản đồ xác suất tiềm năng chuyển đổi. Quá trình huấn luyện đạt chỉ số Root Mean Square (RMS) error $< 0.01$ và độ chính xác xác thực chéo $> 80%$.
  5. Mô phỏng Không gian Markov-CA: Xác suất chuyển tiếp được lặp qua các chu kỳ thời gian với bộ lọc không gian ma trận kề $3 \times 3$ và $5 \times 5$, cho phép tái tạo chân thực hình thái phân mảnh hoặc tập trung của cảnh quan.

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

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

  1. Sự sụp đổ không gian của đất chuyên trồng lúa và sự bành trướng của đất nuôi trồng thủy sản: Giai đoạn 2009–2019 ghi nhận làn sóng chuyển dịch đất đai chưa từng có tại 3 huyện duyên hải. Hàng ngàn hecta đất trồng lúa năng suất thấp tại Ba Tri, Bình Đại và Thạnh Phú đã chuyển đổi mục đích sang nuôi trồng thủy sản nước lợ/mặn (chuyên tôm, cua) và đất trồng cây lâu năm chịu hạn/mặn (như cây dừa). Ma trận chu chuyển đất đai phản ánh tỷ lệ thất thoát đất lúa diễn ra với tốc độ trên $3.5%/\text{năm}$.
  2. Xâm nhập mặn (EC) và Khô hạn (TVDI) là hai biến số cưỡng bức có ý nghĩa thống kê quyết định: Kết quả phân tích hồi quy Logistic xác nhận độ mặn trong đất ($EC > 4\ \text{dS/m}$) và chỉ số khô hạn TVDI đạt mức cực đoan có hệ số tương quan dương vượt trội ($p < 0.001$, $\text{Odds Ratio} > 2.8$) đối với việc từ bỏ canh tác lúa chuyển sang nuôi tôm. Điều này bác bỏ giả thuyết cho rằng thay đổi đất đai chỉ thuần túy do quy hoạch hành chính.
  3. Phân hóa không gian khô hạn theo không - thời gian: Bản đồ TVDI chiết xuất từ mối quan hệ phân tán tán xạ giữa nhiệt độ bề mặt đất (LST) và chỉ số thực vật (NDVI) cho thấy diện tích khô hạn nghiêm trọng tập trung chủ yếu vào các dải giồng cát ven biển thuộc huyện Ba Tri và Thạnh Phú trong các đợt El Niño đạt đỉnh.
  4. Dự báo biến động LULC đến 2029, 2039 và 2049: Kết quả mô phỏng Markov-CA-ANN chỉ ra rằng nếu giữ nguyên xu thế hiện tại kết hợp kịch bản nước biển dâng và triều cường RCP 4.5 đến năm 2050, diện tích đất lúa nước ngọt tại khu vực ven biển Bến Tre sẽ co cụm cục bộ chỉ còn duy trì tại các tiểu vùng có đê bao ngọt hóa hoàn chỉnh; đất nuôi trồng thủy sản nước mặn và lợ sẽ chiếm ưu thế tuyệt đối trong cơ cấu đất nông nghiệp.
  5. Độ chính xác kiểm định mô phỏng vượt chuẩn: Bản đồ mô phỏng hiện trạng năm 2019 (được sinh ra từ mô hình huấn luyện trên dữ liệu 1999–2009) khi đối chiếu với bản đồ hiện trạng thực tế năm 2019 đạt chỉ số Kappa không gian $> 0.85$, chứng minh độ tin cậy khoa học tuyệt đối của thuật toán tích hợp.

Implications đa chiều

  • Về mặt lý luận: Khẳng định tính khả thi của việc tích hợp các mô hình vật lý viễn thám (LST, TVDI, Salinity Index) vào các mô hình học máy phân bổ không gian (Markov-CA-ANN), mở ra hướng đi mới cho Khoa học Biến động Đất đai tại các vùng chịu tổn thương sinh thái.
  • Về mặt phương pháp luận: Cung cấp một khung quy trình hoàn chỉnh, kinh tế và có thể nhân rộng dựa trên phần mềm mã nguồn mở QGIS và dữ liệu vệ tinh miễn phí Landsat/Sentinel, giải phóng sự phụ thuộc vào các phần mềm thương mại đắt tiền cho các viện nghiên cứu và cơ quan quản lý tại các nước đang phát triển.
  • Về mặt quản lý và ứng dụng thực tiễn: Đề xuất 4 nhóm giải pháp chiến lược khả thi cao:
    1. Giám sát theo phương án quy hoạch: Đối soát liên tục giữa bản đồ mô phỏng không gian và quy hoạch sử dụng đất đã được phê duyệt nhằm phát hiện sớm các điểm nóng vi phạm quy hoạch tự phát.
    2. Hệ thống WebGIS mã nguồn mở: Xây dựng cổng thông tin địa lý trực tuyến chia sẻ dữ liệu biến động đất đai thời gian thực phục vụ liên ngành.
    3. Hệ thống giám sát viễn thám cảnh báo sớm mặn - hạn: Tự động cập nhật chỉ số TVDI và EC định kỳ theo mùa khô để cảnh báo nguy cơ suy thoái đất.
    4. Tái cấu trúc mùa vụ và không gian nông nghiệp: Chuyển đổi mô hình sản xuất từ độc canh lúa sang các mô hình thích ứng thông minh với khí hậu (lúa - tôm, chuyên tôm rừng, trồng dừa xen cây ăn trái chịu mặn).
  • Về mặt chính sách: Cung cấp cơ sở khoa học định lượng vững chắc giúp Ủy ban Nhân dân tỉnh Bến Tre và Bộ Tài nguyên & Môi trường điều chỉnh quy hoạch sử dụng đất tầm nhìn 2030–2050 phù hợp với Nghị quyết 120/NQ-CP về phát triển bền vững ĐBSCL thích ứng với BĐKH.

Limitations và Future Research

Luận án thẳng thắn thừa nhận các giới hạn học thuật nội tại:

  • Độ phân giải không gian và thời gian của dữ liệu ảnh: Việc sử dụng ảnh Landsat ($30\text{m} \times 30\text{m}$) có thể bỏ sót các biến động vi mô của các thửa đất manh mún dưới $900\text{m}^2$. Chu kỳ lặp 16 ngày của Landsat dễ bị gián đoạn bởi mây mù trong mùa mưa.
  • Biến số kinh tế - xã hội dạng không gian: Mặc dù các yếu tố tự nhiên được lượng hóa rất chi tiết, các dữ liệu về giá đất, biến động thu nhập nông hộ, chính sách tín dụng vi mô và tập quán canh tác cục bộ chưa được số hóa thành các lớp bản đồ vector/raster có cùng độ phân giải.
  • Giả định quán tính của chuỗi Markov: Ma trận xác suất Markov mặc định quy luật chuyển dịch trong quá khứ tiếp tục chi phối tương lai, chưa tính hết các cú sốc chính sách đột ngột (ví dụ: công trình siêu đê biển ngăn mặn hoàn thành đóng cửa toàn bộ).

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

  1. Tích hợp dữ liệu viễn thám radar đa thời gian (Sentinel-1 SAR) kết hợp ảnh quang học độ phân giải cao (Sentinel-2 MSI $10\text{m}$, PlanetScope $3\text{m}$) nhằm khắc phục triệt để bài toán mây mù và nâng cao độ phân giải không gian.
  2. Mở rộng khung mô hình bằng cách lai ghép giữa mô hình Markov-CA-ANN với Mô hình Đa tác tử (GAMA/ABM) nhằm mô phỏng đồng thời hành vi tâm lý thích ứng của người nông dân và phản ứng giá cả nông sản thị trường.
  3. Đánh giá chu chuyển trữ lượng carbon trong đất và sinh khối (Biomass & Soil Organic Carbon) dưới tác động của sự chuyển dịch từ đất lúa sang đất nuôi tôm và đất rừng ngập mặn.
  4. Xây dựng nền tảng học sâu (Deep Learning) sử dụng mạng tích chập không - thời gian (Spatial-Temporal CNN-LSTM) để tự động hóa toàn bộ chu trình giám sát LULC từ luồng dữ liệu lớn (Big Data) điện toán đám mây Google Earth Engine (GEE).

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

  • Tác động học thuật: Công trình tạo tiền lệ phương pháp luận chuẩn mực cho các nghiên cứu sinh và nhà khoa học trong nước và quốc tế nghiên cứu về đồng bằng châu thổ sông Hồng, sông Mekong, sông Chao Phraya (Thái Lan) và sông Ayeyarwady (Myanmar). Dự kiến công trình sẽ thu hút trích dẫn học thuật cao trong các tạp chí chuyên ngành ISI/Scopus (Q1/Q2) về Land Use Policy, Applied Geography, và Science of The Total Environment.
  • Tác động kinh tế ngành: Giúp ngành nông nghiệp và thủy sản tỉnh Bến Tre cắt giảm hàng trăm tỷ đồng thiệt hại hàng năm do canh tác lúa cố thủ tại các vùng đã bị nhiễm mặn nghiêm trọng, đồng thời định hướng phân vùng đầu tư hạ tầng thủy lợi cấp - thoát nước mặn/ngọt tách biệt.
  • Tác động chính sách: Hỗ trợ trực tiếp Sở Tài nguyên và Môi trường, Sở Nông nghiệp và Phát triển Nông thôn tỉnh Bến Tre xây dựng phương án Quy hoạch tỉnh Bến Tre thời kỳ 2021–2030, tầm nhìn đến năm 2050 theo Quyết định số 1399/QĐ-TTg của Thủ tướng Chính phủ.
  • Tác động xã hội và môi trường: Góp phần bảo vệ sinh kế bền vững cho hơn 500.000 cư dân thuộc 3 huyện ven biển Ba Tri, Bình Đại, Thạnh Phú, hạn chế tình trạng di cư khí hậu cực đoan và bảo tồn đa dạng sinh học hệ sinh thái rừng ngập mặn ven biển.

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

  • Nghiên cứu sinh và Giảng viên ngành Tài nguyên, Môi trường, Trắc địa - Viễn thám, GIS: Tiếp cận một khung phương pháp luận lai ghép hoàn chỉnh, quy trình xử lý ảnh và kịch bản thực thi mô hình hóa toán học chặt chẽ.
  • Các nhà Hoạch định Chính sách và Cơ quan Quản lý Nhà nước (Bộ TN&MT, UBND các tỉnh ĐBSCL): Sở hữu bộ công cụ khoa học có độ tin cậy cao để thẩm định, phân kỳ và điều chỉnh quy hoạch, kế hoạch sử dụng đất dựa trên bằng chứng dữ liệu số.
  • Doanh nghiệp Nông nghiệp, Thủy sản và Chế biến Xuất khẩu: Định vị chính xác các vùng nguyên liệu ổn định (vùng chuyên canh dừa, vùng thâm canh tôm công nghệ cao, vùng lúa sinh thái) để lập dự án đầu tư chuỗi giá trị dài hạn.
  • Cộng đồng Nông dân địa phương: Hưởng lợi từ các thông tin dự báo sớm về xâm nhập mặn, hạn hán và định hướng mùa vụ, chuyển đổi mô hình kinh tế đúng thời điểm, giảm thiểu rủi ro mất trắng mùa màng.

Câu hỏi chuyên sâu

1. Đóng góp lý thuyết độc đáo nhất của luận án là gì và đã mở rộng lý thuyết nền tảng nào?
Trả lời: Luận án đã mở rộng Lý thuyết Khoa học Thay đổi Đất đai (Turner et al., 2007) bằng việc lượng hóa cụ thể các chỉ số phản xạ quang phổ đại diện cho BĐKH (TVDI và EC mặn) thành các biến độc lập liên tục điều khiển xác suất chuyển tiếp của mô hình Markov-CA-ANN. Nghiên cứu đã chứng minh rằng trong điều kiện sinh thái dễ bị tổn thương, các biến cố khí hậu đóng vai trò là "biến cưỡng bức trạng thái" (state-forcing variables) làm thay đổi căn bản cấu trúc bề mặt tiềm năng chuyển đổi đất đai.

2. Điểm cải tiến phương pháp luận cốt lõi của nghiên cứu khi so sánh với các nghiên cứu tiền nhiệm?
Trả lời: So với các nghiên cứu trước đây tại Việt Nam vốn chỉ áp dụng mô hình Markov-CA đơn thuần (như Huỳnh Văn Chương et al., 2017 tại Nha Trang; Trịnh Lê Hùng et al., 2017 tại Hà Nội) chỉ dựa vào khoảng cách không gian đô thị, luận án của NCS. Lê Ngọc Lãm đã thiết lập liên kết thành công giữa viễn thám vật lý (Physical Remote Sensing - chiết xuất LST, NDVI, NDSI để tính toán TVDI và EC) với mô hình học máy phân bổ không gian MOLUSCE và LCM. Đây là bước đột phá kỹ thuật giải quyết trọn vẹn bài toán thiếu dữ liệu quan trắc trạm bề mặt.

3. Phát hiện nào trong dữ liệu thực nghiệm mang tính bất ngờ hoặc trái với trực giác thông thường nhất?
Trả lời: Phát hiện trái trực giác là: Quá trình chuyển đổi từ đất lúa sang đất nuôi trồng thủy sản không chỉ diễn ra tại các vùng sát mép biển mà còn xảy ra mạnh mẽ tại các vùng nội đồng sâu ven các nhánh sông lớn (Hàm Luông, Cổ Chiên). Dữ liệu chỉ ra rằng quá trình truyền triều mang nước mặn lấn sâu vào hệ thống kênh rạch vào mùa khô kết hợp với chỉ số khô hạn TVDI cao đã triệt tiêu hoàn toàn hiệu quả của các công trình ngăn mặn cục bộ, buộc nông dân phải tự phát phá vỡ thế độc canh cây lúa để đào ao nuôi tôm nước lợ.

4. Luận án có cung cấp quy trình nhân bản (Replication Protocol) hoàn chỉnh cho các vùng địa lý khác không?
Trả lời: Có. Luận án xây dựng quy trình 5 bước chuẩn hóa hoàn toàn dựa trên dữ liệu viễn thám công khai (Landsat/Sentinel) và phần mềm mã nguồn mở (QGIS/MOLUSCE). Toàn bộ phương trình chuyển đổi quang phổ, công thức tính toán chỉ số độ mặn EC, thuật toán phân cấp TVDI và quy trình huấn luyện mạng nơ-ron đều được mô tả chi tiết, cho phép áp dụng ngay lập tức cho các tỉnh ven biển khác tại ĐBSCL (như Trà Vinh, Sóc Trăng, Bạc Liêu, Cà Mau, Kiên Giang) cũng như các vùng cửa sông trên toàn thế giới.

5. Chương trình nghị sự nghiên cứu 10 năm tới (10-year Research Agenda) được định hình như thế nào từ kết quả luận án?
Trả lời: Chương trình nghị sự tập trung vào 3 trọng tâm: (1) Nâng cấp hệ thống lên nền tảng Điện toán Đám mây Big Data (Google Earth Engine) giám sát thời gian thực chu kỳ 5 ngày; (2) Tích hợp mô hình hóa rủi ro đa thiên tai (Multi-hazard risk modeling) kết hợp kịch bản ngập triều dâng bão và sạt lở bờ sông bờ biển; (3) Tích hợp chỉ số vốn xã hội và hành vi thích ứng kinh tế nông hộ vào mạng nơ-ron sâu để dự báo điểm bùng phát sinh thái (ecological tipping points) của hệ thống cảnh quan.

Kết luận

Luận án tiến sĩ của Lê Ngọc Lãm đã đạt được những cống hiến học thuật và thực tiễn xuất sắc:

  1. Hệ thống hóa và hoàn thiện cơ sở khoa học: Xác lập luận cứ định lượng vững chắc về mối quan hệ nhân quả tương hỗ giữa BĐKH, xâm nhập mặn, khô hạn và biến động sử dụng đất tại vùng đồng bằng ven biển.
  2. Đột phá phương pháp tích hợp RS - GIS - Math Models: Xây dựng thành công quy trình công nghệ kết hợp ảnh Landsat, chỉ số TVDI, chỉ số độ dẫn điện EC với các mô hình toán LCM, MOLUSCE, Markov-CA và mạng nơ-ron MLP.
  3. Đánh giá toàn diện bức tranh biến động 2009–2019: Làm rõ bản chất và quy luật chu chuyển không gian của các loại hình đất nông nghiệp, đất phi nông nghiệp và đất nuôi trồng thủy sản tại 3 huyện ven biển Ba Tri, Bình Đại, Thạnh Phú.
  4. Dự báo chuẩn xác các kịch bản không gian 2029–2049: Cung cấp bức tranh phân bổ đất đai đáng tin cậy phục vụ định hướng phát triển bền vững đến năm 2050 trong bối cảnh kịch bản phát thải RCP 4.5.
  5. Kiến tạo 4 giải pháp công nghệ - quản lý mang tính ứng dụng cao: Đóng góp trực tiếp hệ thống công cụ hỗ trợ ra quyết định, nền tảng WebGIS và chiến lược tái cơ cấu mùa vụ thích ứng BĐKH cho tỉnh Bến Tre.
  6. Mở ra 3 hướng nghiên cứu liên ngành mới: Mở đường cho các nghiên cứu chuyên sâu về viễn thám mặn - hạn thời gian thực, mô hình hóa tác tử sinh thái - xã hội và đánh giá dịch vụ hệ sinh thái đất ngập nước ven biển.

Công trình là một tài liệu tham khảo mẫu mực, đóng góp to lớn vào kho tàng tri thức Khoa học Quản lý Tài nguyên và Môi trường Việt Nam trong kỷ nguyên biến đổi toàn cầu.