Tổng quan về luận án

Biến động sử dụng đất (Land Use and Land Cover Change - LUCC) tại các dải đồng bằng ven biển châu thổ nhiệt đới là tâm điểm phản ánh sự tương tác phức tạp giữa các quá trình sinh - địa vật lý và áp lực nhân sinh. Luận án tiến sĩ Địa lý của nghiên cứu sinh Phạm Thị Làn, chuyên ngành Bản đồ, viễn thám và Hệ thông tin địa lý (Mã số: 62 44 02 14), bảo vệ tại Trường Đại học Khoa học Tự nhiên – Đại học Quốc gia Hà Nội (2017) dưới sự hướng dẫn của GS. Võ Chí Mỹ và PGS. Phạm Văn Cự, đã tạo dựng một dấu ấn khoa học tiên phong. Công trình giải quyết câu hỏi then chốt: Làm thế nào để định lượng hóa chính xác mối quan hệ hữu cơ giữa các yếu tố tự nhiên nền tảng với động lực biến động lớp phủ/sử dụng đất trong bối cảnh chịu sự chi phối mạnh mẽ của các chính sách phát triển kinh tế - xã hội?

Khoảng trống nghiên cứu (research gap) cốt lõi xuất phát từ thực tiễn học thuật: các nghiên cứu trước đây thường phân tách ranh giới chuyên ngành, hoặc thuần túy phân tích định tính mối quan hệ tự nhiên – nhân tạo, hoặc xem xét biến động sử dụng đất chủ yếu dưới lăng kính kinh tế - xã hội (Lambin & Geist, 2006; Eric Lambin, 2007). Một số công trình quốc tế như Schneider & Pontius (2001) hay Hietel (2004) cho rằng biến động sử dụng đất chỉ có mối quan hệ lỏng lẻo với điều kiện tự nhiên; trong khi các nghiên cứu của Pan et al. (1999), Chen (2001) và Fu (2006) lại khẳng định mối liên hệ vô cùng chặt chẽ. Sự mâu thuẫn lý thuyết này đặt ra yêu cầu phải kiểm chứng thực nghiệm tại các hệ sinh thái ven biển có độ nhạy cảm cao.

Luận án xác lập hệ thống 03 câu hỏi và giả thuyết nghiên cứu trọng tâm:

  1. Câu hỏi 1 (Q1): Phương pháp phân loại viễn thám nào cho phép trích xuất chính xác chuỗi dữ liệu hiện trạng sử dụng đất đa thời gian tại vùng ven biển có độ phân mảnh cao?
    Giả thuyết 1 (H1): Tích hợp phân loại định hướng đối tượng (OBIA) và phân vùng mẫu (per-field classification) từ ảnh vệ tinh Landsat đa thời gian sẽ nâng cao độ chính xác giải đoán vượt trội so với các thuật toán dựa trên điểm ảnh (pixel-based) truyền thống.
  2. Câu hỏi 2 (Q2): Các nhân tố tự nhiên đóng vai trò kiểm soát hay chỉ là điều kiện nền tảng thích ứng đối với sự chuyển dịch cơ cấu sử dụng đất?
    Giả thuyết 2 (H2): Địa mạo (cồn cát cổ, bãi bồi ven biển) và thổ nhưỡng (đất mặn, đất phù sa, đất phèn, đất cát) có mối tương quan phi tuyến chặt chẽ, chi phối không gian phân bố và xác suất chuyển đổi của các loại hình sử dụng đất.
  3. Câu hỏi 3 (Q3): Việc tích hợp các thông số tự nhiên vào mô hình học máy phi tuyến có nâng cao độ chính xác dự báo không gian sử dụng đất tương lai?
    Giả thuyết 3 (H3): Mô hình tích hợp Mạng nơ-ron nhân tạo đa lớp (Artificial Neural Network / Multi-Layer Perceptron) và Tế bào tự hành kết hợp chuỗi Markov (CA-Markov) trong môi trường Land Change Modeler (LCM) sẽ mô phỏng chuẩn xác động lực biến động sử dụng đất đến năm 2021.

Khung lý thuyết của công trình kế thừa hệ thống phân loại chức năng lớp phủ của Anderson et al. (1976), Jensen et al. (2003), lý thuyết hệ thống cảnh quan của Nguyễn Cao Huần (2005), Phạm Hoàng Hải (1997) và khung tương tác nhân tố ngoại sinh - nội sinh của Satoshi Hoshino (1996).

Về phạm vi không gian và thời gian, luận án khảo sát thực nghiệm trên hai huyện ven biển điển hình của Đồng bằng sông Hồng (ĐBSH): huyện Thái Thụy (tỉnh Thái Bình) – đặc trưng bởi thùy châu thổ Diêm Điền cổ đang trong giai đoạn lụi tàn, chịu ảnh hưởng thủy văn sông Thái Bình với hướng vận chuyển trầm tích chủ đạo Đông Nam; và huyện Giao Thủy (tỉnh Nam Định) – đặc trưng bởi thùy châu thổ Ba Lạt phát triển mạnh mẽ, bồi tụ tích cực bởi hệ thống sông Hồng với dòng trầm tích vận chuyển chủ đạo hướng Đông Bắc. Dữ liệu viễn thám bao trùm chuỗi thời gian 24 năm từ 1989 đến 2013 (các mốc 1989, 2001, 2005, 2013) và mở rộng dự báo xu hướng đến năm 2021.

Literature Review và Positioning

Tổng quan y văn quốc tế và trong nước về biến động lớp phủ/sử dụng đất chỉ ra ba trường phái học thuật chủ đạo:

Trường phái thứ nhất nhấn mạnh vai trò quyết định của các động lực kinh tế - xã hội (Socio-economic driving forces). Lambin et al. (2001, 2003, 2007) phân định rạch ròi giữa nguyên nhân trực tiếp (proximate causes) ở cấp vi mô (hộ gia đình, nông trại) như sự mở rộng nông nghiệp, khai thác gỗ, và nguyên nhân gián tiếp (underlying driving forces) ở cấp vĩ mô như chính sách kinh tế, nhân khẩu học, công nghệ và thể chế. A. Khatoon-Abali et al. (2011) khi nghiên cứu vùng Kalardasht (Iran) giai đoạn 1985–2011 đã chỉ ra chính sách đất đai bất cập, tình trạng nghèo đói và áp lực thị trường bất động sản là nguyên nhân trực tiếp tàn phá cảnh quan.

Trường phái thứ hai nhấn mạnh vai trò kiểm soát của các yếu tố sinh thái cảnh quan và địa vật lý (Biophysical determinants). Jan Kabrda (2004, 2007) khi phân tích dữ liệu biến động sử dụng đất tại Cộng hòa Séc đã áp dụng hệ số tương quan Pearson và mô hình hồi quy bội (Multiple Regression Analysis), chứng minh độ cao, độ dốc, đặc tính thạch quyển và khí hậu tạo lập ngưỡng giới hạn sinh thái quyết định tính ổn định của cảnh quan nông nghiệp. Tương tự, Ataollah Kelarestaghi et al. (2011) tại Bắc Iran hay Shu-Wei Huang & Hsiu-I Hsiech (2012) tại dải ven biển Đài Loan đều khẳng định địa mạo, độ dốc, khoảng cách tới bờ biển và mức độ nhạy cảm môi trường chi phối cấu trúc phân bố không gian đất đai.

Trường phái thứ ba tích hợp mô hình hóa không gian và động lực học hệ sinh thái. Satoshi Hoshino (1996) tại vùng Kansai (Nhật Bản) đã xây dựng khung phân tích đột phá: phân tách nhân tố làm biến động sử dụng đất thành bốn nhóm bao gồm lực dẫn dắt ngoại sinh (driving forces), lực thích ứng nội sinh (counter forces), điều kiện tự nhiên nền tảng và hiện trạng sử dụng đất ban đầu. Hoshino sử dụng phân tích thành phần chính (PCA) kết hợp hồi quy Logistic đa thức (Multinomial Logistic Regression) để xử lý tương quan phi tuyến. Cùng hướng tiếp cận này, Walsh et al. (2001) tại huyện Nang Rong (Thái Lan) dùng phân tích tương quan chính tắc (Canonical Correlation Analysis) để lượng hóa quan hệ dân số - môi trường; còn Veldkamp & Fresco (1997) tại Costa Rica dùng phân tích nhân tố trên sáu cấp độ không gian khác nhau để nhận diện cấu trúc phụ thuộc của thảm phủ thực vật vào khả năng tưới tiêu và độ cao.

Tại Việt Nam, các công trình của Nguyễn Cao Huần (2005) với mô hình tích hợp ALES-GIS, Vũ An Tuân (2004) tại lưu vực sông Trà Khúc, Phạm Bách Việt & Hồ Đình Duẩn (2006) tại Hà Nội, Tống Sĩ Sơn (2009) với kỹ thuật phân tích véc-tơ biến động (CVA) dựa trên chỉ số NDVI và SAVI, Nguyễn Thị Thu Hiền et al. (2014) và Nguyễn Hải Hòa (2016) tại Tiên Yên (Quảng Ninh) đã ứng dụng kỹ thuật phân tích sau phân loại (post-classification comparison) để phát hiện biến động rừng ngập mặn và đất đô thị.

Luận án của Phạm Thị Làn định vị chính xác khoảng trống học thuật: Trong khi các công trình đi trước tại Việt Nam chủ yếu dừng lại ở việc mô tả thống kê biến động diện tích, chồng xếp bản đồ định tính hoặc đánh giá thích nghi đất đai tĩnh, luận án đã tiến một bước dài khi thiết lập mô hình định lượng hóa xác suất chuyển dịch sử dụng đất dưới tác động đồng thời của động lực trầm tích ven bờ, quá trình xói lở/bồi tụ, phân hóa địa mạo thùy châu thổ và tính chất thổ nhưỡng mặn - phèn.

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

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

Luận án đóng góp quan trọng vào lý thuyết Địa lý tổng hợp và Sinh thái học cảnh quan thông qua việc làm sáng tỏ cơ chế phản hồi tương hỗ giữa hệ thống tự nhiên và can thiệp nhân sinh tại vùng bờ biển (Coastal zone).

Thứ nhất, luận án kiểm chứng và bổ sung cho mô hình lý thuyết của Lambin & Geist (2006) và Hoshino (1996): Điều kiện tự nhiên không thuần túy là bối cảnh tĩnh (static background) mà là bộ biến số kiểm soát động (dynamic controlling variables) thiết lập nên các ngưỡng sinh thái (ecological thresholds). Tác giả đã bảo vệ một quan điểm học thuật sâu sắc:

"Đổi Mới cũng không thể đưa thực vật ngập mặn vào đất phù sa để trồng. Đổi Mới hay bất kỳ chính sách phát triển kinh tế - xã hội nào muốn chuyển đổi sử dụng đất đều phải tuân thủ quy luật tự nhiên nhằm hướng tới mục tiêu phát triển bền vững."

Thứ hai, công trình hệ thống hóa và lượng hóa quy luật diễn thế cảnh quan nhân sinh vùng ven biển Bắc Bộ:

"Sử dụng đất biến đổi theo lộ trình “lúa lấn cói, cói lấn vẹt, vẹt lấn biển”. Sự chuyển đổi sử dụng đất mạnh mẽ trong suốt gần 30 năm qua là kết quả của sự những lần thay đổi chính sách phát triển kinh tế - xã hội, mốc chuyển đổi chính sách đặc biệt quan trọng và làm thay đổi lớn các loại hình sử dụng đất là chính sách Đổi Mới và các dự án về trồng rừng ngập mặn."

Mô hình lý thuyết của luận án xác lập các mệnh đề khoa học (propositions):

  • Mệnh đề 1 (P1): Xác suất chuyển đổi đất nông nghiệp sang đất nuôi trồng thủy sản tỷ lệ thuận với mức độ nhiễm mặn của thổ nhưỡng và khoảng cách tiếp cận nguồn nước lợ tại các bãi triều thấp.
  • Mệnh đề 2 (P2): Động lực bồi tụ tại các cửa sông lớn (thùy Ba Lạt) tạo điều kiện tiên quyết mở rộng đai rừng ngập mặn tiên phong (sú, vẹt), đóng vai trò bệ đỡ sinh thái cho quá trình quai đê lấn biển nhân tạo.
  • Mệnh đề 3 (P3): Động lực xói lở và suy thoái trầm tích (thùy Diêm Điền) gia tăng tính dễ tổn thương của các mô hình canh tác lúa nước, thúc đẩy chuyển đổi thích ứng sang đất phi nông nghiệp hoặc nuôi trồng thủy sản chuyên canh.
flowchart TD
    subgraph Natural_Drivers ["HỆ THỐNG ĐIỀU KIỆN TỰ NHIÊN (NỘI SINH)"]
        A1["Địa mạo (Cồn cát cổ, Bãi bồi ven biển)"]
        A2["Thổ nhưỡng (Đất mặn, Phù sa, Phèn, Cát)"]
        A3["Động lực bờ biển (Bồi tụ thùy Ba Lạt, Xói lở thùy Diêm Điền)"]
        A4["Khoảng cách tới đường bờ & nguồn nước"]
    end

    subgraph Socioeconomic_Drivers ["CHÍNH SÁCH & KINH TẾ (NGOẠI SINH)"]
        B1["Chính sách Đổi mới kinh tế (1986)"]
        B2["Dự án Trồng rừng ngập mặn & Quai đê lấn biển"]
        B3["Chuyển dịch cơ cấu Thủy sản - Công nghiệp"]
    end

    subgraph Integration_Model ["MÔ HÌNH HÓA ĐỊNH LƯỢNG (GIS - VIỄN THÁM)"]
        C1["Hồi quy Logistic Đa biến (SPSS/R)"]
        C2["Mạng Nơ-ron Nhân tạo Đa lớp (MLP)"]
        C3["Tế bào Tự hành - Chuỗi Markov (CA-Markov)"]
    end

    subgraph LUCC_Transitions ["MA TRẬN CHUYỂN DỊCH SỬ DỤNG ĐẤT"]
        D1["Đất trồng lúa & hoa màu"]
        D2["Đất nuôi trồng thủy sản"]
        D3["Đất rừng ngập mặn"]
        D4["Đất trống / Bãi triều"]
        D5["Đất dân cư & xây dựng"]
    end

    Natural_Drivers --> Integration_Model
    Socioeconomic_Drivers --> Integration_Model
    Integration_Model --> LUCC_Transitions
    LUCC_Transitions -.->|"Phản hồi sinh thái (Xâm nhập mặn, Xói lở)"| Natural_Drivers

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

Khung phân tích của luận án tích hợp ba trụ cột phương pháp luận:

  1. Lý thuyết tiếp cận lớp phủ suy luận sử dụng đất (Jensen et al., 2003): Do tư liệu Landsat không trực tiếp phân biệt được các thực thể pháp lý hoặc kinh tế chi tiết (như đất chuyên dùng, đất an ninh quốc phòng), luận án thiết lập "chìa khóa giải đoán" phân tầng: từ các thuộc tính vật lý của lớp phủ (nước, thực vật, bề mặt nhân tạo, đất trống) kết hợp với các dữ liệu bổ trợ GIS (bản đồ địa hình, địa mạo, thổ nhưỡng) để ngoại suy ra chức năng sử dụng đất.
  2. Khung phân tích nhân quả Hoshino - Lambin: Tách biệt rõ ràng biến phụ thuộc (biến động sử dụng đất dạng nhị phân/rời rạc) và biến độc lập (các lớp chuyên đề tự nhiên định lượng).
  3. Mô hình hóa động lực học không gian phi tuyến trong nền tảng LCM (Land Change Modeler): Kết hợp trí tuệ nhân tạo (Mạng nơ-ron MLP) để tạo lập bản đồ tiềm năng chuyển đổi (Transition Potential Maps) và thuật toán CA-Markov để phân bổ không gian tối ưu.

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

Thiết kế nghiên cứu

Nghiên cứu được thiết kế theo lập trường bản thể luận thực tế khách quan (Objectivism/Positivism) kết hợp chủ nghĩa hiện thực phê phán (Critical Realism) trong phân tích tương tác không gian. Thiết kế phương pháp luận hỗn hợp định lượng (Quantitative Mixed-Method Design) đa tầng bậc được triển khai chặt chẽ, kết hợp công nghệ xử lý ảnh viễn thám số, hệ thống thông tin địa lý (GIS) và mô hình hóa toán xác suất.

graph LR
    subgraph Data_Acquisition ["THU THẬP & TIỀN XỬ LÝ DỮ LIỆU"]
        T1["Chuỗi ảnh Landsat TM/ETM+/OLI (1989, 2001, 2005, 2013)"]
        T2["Bản đồ Địa hình 1:50.000 (VN-2000)"]
        T3["Bản đồ Thổ nhưỡng & Địa mạo 1:50.000"]
        T4["Niên giám thống kê huyện (Thái Thụy, Giao Thủy)"]
    end

    subgraph Image_Processing ["CHIẾT XUẤT THÔNG TIN SỬ DỤNG ĐẤT"]
        P1["Nắn chỉnh hình học tọa độ VN-2000"]
        P2["Phân đoạn ảnh đa tỷ lệ (Multiresolution Segmentation)"]
        P3["Phân loại Định hướng Đối tượng (OBIA)"]
        P4["Hiển chỉnh Phân vùng mẫu (Per-field) + Dữ liệu Vector"]
    end

    subgraph Spatial_Analysis ["PHÂN TÍCH QUAN HỆ & MÔ HÌNH HÓA"]
        M1["Chồng xếp ma trận biến động sau phân loại"]
        M2["Phân tích Hồi quy Logistic Đa biến"]
        M3["Huấn luyện Mạng Nơ-ron Nhân tạo (MLP) trong LCM"]
        M4["Mô phỏng Động lực Không gian CA-Markov (2013 -> 2021)"]
    end

    subgraph Validation ["KIỂM ĐỊNH ĐỘ CHÍNH XÁC"]
        V1["Hệ số Kappa (Kstandard, Kno, Klocation)"]
        V2["Chỉ số ROC (Area Under Curve)"]
        V3["So sánh tương quan Niên giám thống kê"]
    end

    Data_Acquisition --> Image_Processing
    Image_Processing --> Spatial_Analysis
    Spatial_Analysis --> Validation

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu bao gồm 4 giai đoạn logic khép kín:

  1. Chuẩn hóa và tiền xử lý tư liệu viễn thám:

    • Sử dụng 04 cảnh ảnh Landsat (Path 126 / Row 046) bao gồm: Landsat TM (ngày 23/11/1989), Landsat ETM+ (ngày 16/11/2001), Landsat ETM+ có kênh toàn sắc Pan 15m (ngày 10/10/2005), và Landsat 8 OLI (ngày 08/10/2013).
    • Nắn chỉnh hình học và chuyển đổi hệ quy chiếu toàn bộ các cảnh ảnh về Hệ tọa độ Quốc gia VN-2000 (kinh tuyến trục phù hợp với tỉnh Thái Bình và Nam Định) dựa trên bản đồ địa hình 1:50.000 xuất bản năm 2000 của Cục Đo đạc Bản đồ Quân đội, đảm bảo sai số vị trí điểm kiểm tra (RMSE) nhỏ hơn 0.5 pixel (< 15m).

  2. Chiết xuất thông tin lớp phủ/sử dụng đất:

    • Áp dụng phương pháp phân loại định hướng đối tượng (OBIA) trên phần mềm chuyên dụng eCognition/Definiens. Quá trình bắt đầu bằng thuật toán phân mảnh đa độ phân giải (multiresolution segmentation) để gộp nhóm các pixel đồng nhất thành các đối tượng ảnh (image objects) dựa trên trọng số phổ (color/spectral) và hình học (shape, compactness, smoothness).
    • Thiết lập cây quyết định mờ (fuzzy logic rules) kết hợp các chỉ số phổ: NDVI (Normalized Difference Vegetation Index), MNDWI (Modified Normalized Difference Water Index), SAVI (Soil-Adjusted Vegetation Index), cùng các thông số quan hệ không gian đối tượng lân cận (relations to neighbor objects).
    • Tích hợp kỹ thuật phân vùng mẫu (per-field classification) bằng cách liên kết ranh giới tự nhiên của các đối tượng ảnh với ranh giới bản đồ số dạng vector nhằm loại bỏ hiệu ứng "hạt tiêu" (salt-and-pepper noise), nâng cao độ chính xác giải đoán đối với các thửa đất nuôi trồng thủy sản và thổ cư xen kẹt.

  3. Phân tích tương quan thống kê không gian:

    • Thiết lập mô hình hồi quy Logistic đa biến (Multivariate Logistic Regression) để lượng hóa mức độ tác động của 4 biến độc lập tự nhiên: (1) Loại hình địa mạo, (2) Loại đất thổ nhưỡng, (3) Khoảng cách tới đường bờ biển hiện tại, và (4) Biên độ biến động đường bờ (xói lở/bồi tụ giai đoạn 1989–2008).
    • Phương trình hồi quy Logistic tổng quát: $$\ln\left(\frac{P_i}{1 - P_i}\right) = \beta_0 + \beta_1 X_{1i} + \beta_2 X_{2i} + \dots + \beta_k X_{ki}$$ Trong đó $P_i$ là xác suất biến đổi của loại hình sử dụng đất thứ $i$, $X_{ki}$ là các biến độc lập tự nhiên, $\beta_k$ là các hệ số hồi quy được ước lượng bằng phương pháp hợp lý cực đại (Maximum Likelihood Estimation).

  4. Mô hình hóa dự báo động lực sử dụng đất bằng LCM (Land Change Modeler):

    • Ứng dụng cấu trúc mạng nơ-ron nhân tạo lan truyền ngược nhiều lớp (Multi-Layer Perceptron - MLP) trong module LCM của phần mềm IDRISI để huấn luyện mẫu chuyển dịch giữa các giai đoạn 1989–2001 và 2001–2013.
    • Sử dụng thuật toán Tế bào tự hành kết hợp chuỗi Markov (CA-Markov) để mô phỏng ma trận xác suất chuyển dịch và dự báo cấu trúc không gian sử dụng đất đến năm 2021.

Data và phân tích

Độ tin cậy của dữ liệu được thẩm định nghiêm ngặt qua quy trình tam giác hóa (triangulation):

  • Đối chiếu số liệu thống kê: So sánh diện tích từng loại hình sử dụng đất chiết xuất từ ảnh vệ tinh với số liệu tổng hợp trong Niên giám thống kê huyện Giao Thủy và huyện Thái Thụy tại các mốc 1989, 2001, 2005 và 2013.
  • Kiểm định thực địa: Sử dụng dữ liệu khảo sát GPS thực địa và các bản đồ hiện trạng sử dụng đất 1:50.000 năm 2000, 2005 của Sở Tài nguyên và Môi trường tỉnh Thái Bình, Nam Định.
  • Kiểm định mô hình dự báo: Đánh giá độ chính xác mô phỏng không gian thông qua các chỉ số thống kê tiêu chuẩn: Chỉ số diện tích dưới đường cong ROC (Receiver Operating Characteristic - AUC) đạt giá trị > 0.82 cho các lớp chuyển đổi chính; hệ số Kappa tổng thể ($K_{standard}$, $K_{no}$, $K_{location}$) đạt trên 85%, khẳng định tính chuẩn xác và khả năng khái quát hóa cao của mô hình.

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

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

  1. Sự bùng nổ của đất nuôi trồng thủy sản và sự suy giảm đất trồng trọt:
    Giai đoạn 1989–2013 chứng kiến cuộc chuyển dịch cơ cấu sử dụng đất ngoạn mục nhất trong lịch sử đới bờ ĐBSH. Tại huyện Thái Thụy và Giao Thủy, diện tích đất nông nghiệp (đặc biệt là đất trồng lúa và hoa màu) liên tục suy giảm để nhường chỗ cho đất phi nông nghiệp, hạ tầng dân sinh và đất nuôi trồng thủy sản. Quá trình chuyển đổi từ đất lúa - cói sang đầm tôm/cua diễn ra mạnh mẽ nhất trong giai đoạn 2001–2005, phù hợp với chính sách chuyển dịch cơ cấu nông nghiệp sau Đổi mới.

  2. Địa mạo và thổ nhưỡng đóng vai trò kiểm soát xác suất chuyển đổi đất đai:
    Kết quả hồi quy Logistic khẳng định mối tương quan có ý nghĩa thống kê cao ($p < 0.01$) giữa các biến tự nhiên với xác suất biến đổi:

    • Các dải cồn cát cổ ven biển có tính ổn định cao, là nền tảng tự nhiên hình thành các điểm dân cư truyền thống và đất hoa màu, hầu như không bị chuyển đổi sang đất ngập nước.
    • Vùng bãi bồi ven biển và đất mặn nhiều, đất mặn trung bình có xác suất chuyển đổi sang đất nuôi trồng thủy sản cao gấp 3.8 lần so với các đơn vị đất phù sa ngọt không mặn.
    • Các diện tích đất phèn tiềm tàng ven sông bị chuyển đổi mạnh từ lúa năng suất thấp sang đầm nuôi thủy sản nước lợ hoặc rừng ngập mặn phòng hộ.

  3. Sự phân hóa sâu sắc theo động lực trầm tích thùy châu thổ:

    • Tại huyện Giao Thủy (thùy Ba Lạt - sông Hồng), lượng phù sa dồi dào với hướng vận chuyển trầm tích Đông Bắc đã bồi tụ bãi triều Cồn Lu, Cồn Ngạn với tốc độ hàng chục mét mỗi năm. Hiện tượng "đất mới nổi cuối Cồn Lu" ghi nhận sự gia tăng đất chưa sử dụng, tạo tiền đề cho diện tích rừng ngập mặn tiên phong mở rộng tự nhiên.
    • Ngược lại, tại huyện Thái Thụy (thùy Diêm Điền - sông Thái Bình), trạng thái lụi tàn của thùy châu thổ dẫn đến hiện tượng xói lở cục bộ đường bờ biển, làm suy giảm diện tích bãi bồi tự nhiên, buộc các can thiệp nhân sinh phải tập trung vào việc gia cố đê kè biển và chuyển dịch mạnh sang nuôi ngao/vạng ngoài bãi triều thấp.

  4. Tác động nhân sinh tái cấu trúc điều kiện tự nhiên:
    Luận án phát hiện quy luật tác động ngược chiều (feedback loop): Việc con người đẩy mạnh quai đê lấn biển, đắp đập ngăn mặn xây hồ chứa nước ngọt, và phá bỏ rừng ngập mặn nguyên sinh để đào ao tôm công nghiệp đã làm thay đổi căn bản chế độ thủy động lực ven bờ, gia tăng hiện tượng xói lở đáy và xâm nhập mặn cục bộ trong tầng nước ngầm.

Implications đa chiều

  • Đóng góp lý luận: Chứng minh sinh động nguyên lý địa lý tổng hợp: tính phân hóa của cảnh quan tự nhiên chi phối trực tiếp không gian kinh tế nhân sinh; đồng thời cung cấp bằng chứng thực nghiệm giải quyết tranh luận học thuật giữa trường phái "tự nhiên lỏng lẻo" (Schneider & Pontius) và "tự nhiên quyết định" (Pan et al.).
  • Đổi mới phương pháp luận: Chuẩn hóa quy trình tích hợp giữa phân loại viễn thám định hướng đối tượng (OBIA), kỹ thuật phân vùng mẫu (per-field) với mô hình học máy Mạng nơ-ron đa lớp (MLP) trong LCM, có thể chuyển giao áp dụng cho toàn bộ 28 tỉnh, thành ven biển Việt Nam và các vùng châu thổ nhiệt đới Đông Nam Á.
  • Ứng dụng thực tiễn và quản trị: Cung cấp cơ sở khoa học để chính quyền tỉnh Thái Bình và Nam Định rà soát ranh giới 3 loại rừng (đặc dụng, phòng hộ, sản xuất), khoanh định không gian bảo vệ nghiêm ngặt Vườn quốc gia Xuân Thủy, và quy hoạch phân vùng nuôi trồng thủy sản bền vững thích ứng với biến đổi khí hậu và mực nước biển dâng.

Limitations và Future Research

Mặc dù đạt được những kết quả đột phá, luận án thẳng thắn chỉ ra các giới hạn nghiên cứu:

  1. Độ phân giải không gian của tư liệu viễn thám: Việc sử dụng ảnh Landsat độ phân giải 30m gặp khó khăn nhất định trong việc phân tách các thửa đất thổ cư xen kẽ vườn cây ăn quả hoặc các công trình công cộng quy mô nhỏ dưới 900 $m^2$.
  2. Số lượng biến số độc lập: Mô hình mới chỉ tích hợp 4 nhóm biến tự nhiên chủ chốt (địa mạo, thổ nhưỡng, biến động đường bờ, khoảng cách bờ biển), chưa tích hợp định lượng dữ liệu vi khí hậu biến thiên theo mùa và số liệu cao độ số (DEM) độ phân giải siêu cao (LiDAR).
  3. Giả định ngoại suy kinh tế - xã hội: Mô hình Markov giả định tốc độ chuyển đổi quá khứ tiếp tục duy trì trong tương lai mà chưa mô phỏng hết các kịch bản đột biến chính sách (chẳng hạn quy hoạch đột xuất các khu kinh tế ven biển lớn).

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

  • Tích hợp chuỗi dữ liệu vệ tinh đa phổ độ phân giải cao (Sentinel-2A/B 10m, PlanetScope 3m) kết hợp dữ liệu siêu cao tần SAR (Sentinel-1) để giám sát độ ẩm và độ mặn đất theo thời gian thực.
  • Mở rộng mô hình Agent-Based Modeling (ABM) kết hợp CA-Markov để mô phỏng hành vi ra quyết định của từng nhóm cộng đồng ngư dân ven biển.
  • Đánh giá định lượng dịch vụ hệ sinh thái (Ecosystem Services Valuation) gắn với thị trường tín chỉ carbon (Blue Carbon) của rừng ngập mặn ĐBSH.

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

Công trình của nghiên cứu sinh Phạm Thị Làn tạo ra tác động học thuật và thực tiễn sâu rộng:

  • Tác động học thuật: Đặt nền móng phương pháp luận chuẩn xác cho các nghiên cứu sinh và học viên cao học ngành Địa lý, Bản đồ học, Viễn thám & GIS; cung cấp tài liệu tham khảo giá trị cho các giáo trình chuyên khảo về Viễn thám ứng dụng và Quản lý dải ven bờ.
  • Tác động chính sách và quy hoạch: Cung cấp luận cứ khoa học trực tiếp cho Sở Tài nguyên và Môi trường, Sở Nông nghiệp và Phát triển nông thôn các tỉnh Thái Bình, Nam Định trong công tác lập Quy hoạch sử dụng đất cấp huyện giai đoạn 2021–2030, tầm nhìn 2050.
  • Bảo tồn và thích ứng biến đổi khí hậu: Hỗ trợ kỹ thuật cho Ban quản lý Khu Dự trữ sinh quyển thế giới liên tỉnh Châu thổ sông Hồng trong việc thiết lập hành lang xanh bảo vệ rừng ngập mặn phòng hộ sóng thần, giảm thiểu rủi ro thiên tai ven bờ.

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

  • Nghiên cứu sinh và Giới nghiên cứu Địa tin học: Tiếp cận quy trình trích xuất thông tin OBIA kết hợp hồi quy Logistic và học máy MLP-Markov chuẩn mực.
  • Nhà hoạch định chính sách cấp tỉnh/huyện: Nhận diện chính xác quy luật chuyển dịch đất đai để ban hành các quy định phân vùng sử dụng đất phù hợp với sức tải sinh thái.
  • Chuyên gia Nông - Thủy sản: Hiểu rõ mối liên hệ giữa độ mặn đất, địa mạo bãi triều với năng suất thủy sản để tối ưu hóa mô hình canh tác tôm - rừng bền vững.
  • Cộng đồng dân cư ven biển: Hưởng lợi gián tiếp từ các chính sách quy hoạch hợp lý, giảm thiểu thiệt hại do sạt lở bờ sông, bờ biển và xâm nhập mặn.

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

  1. Đóng góp lý thuyết độc đáo nhất của luận án là gì?
    Đó là việc chứng minh và lượng hóa thành công vai trò "kiểm soát giới hạn" của địa mạo và thổ nhưỡng đối với việc chuyển đổi sử dụng đất vùng ven biển, bổ sung luận cứ phản biện cho khung lý thuyết tương tác môi trường - kinh tế của Lambin (2007) và Hoshino (1996).

  2. Điểm mới về phương pháp luận so với các nghiên cứu trước đây?
    Khác với các nghiên cứu của Schneider & Pontius (2001) dùng phương pháp điểm ảnh đơn lẻ hay Nguyễn Thị Thu Hiền et al. (2014) chỉ so sánh sau phân loại đơn thuần, luận án đã tích hợp thành công OBIA, phân vùng mẫu (per-field) và học máy Mạng nơ-ron MLP trong LCM, giúp nâng cao độ chính xác giải đoán và dự báo không gian đạt hệ số Kappa trên 85%.

  3. Phát hiện bất ngờ và thú vị nhất rút ra từ dữ liệu thực nghiệm?
    Mặc dù chính sách Đổi mới thúc đẩy mạnh mẽ chuyển đổi sang thủy sản trên toàn vùng, nhưng tốc độ và hình thái chuyển đổi tại Thái Thụy và Giao Thủy lại diễn dịch hoàn toàn khác biệt do sự đối lập về động lực thùy châu thổ (thùy Diêm Điền xói lở - lụi tàn đối lập với thùy Ba Lạt bồi tụ - phát triển).

  4. Luận án có cung cấp quy trình tái lặp (Replication Protocol) hoàn chỉnh không?
    Hoàn toàn có. Toàn bộ các tham số phân mảnh ảnh, bảng quy tắc phân loại mờ, giá trị ngưỡng phân loại phổ Landsat, cấu trúc biến hồi quy Logistic và quy trình huấn luyện mạng nơ-ron MLP đều được mô tả chi tiết, rõ ràng trong các sơ đồ và bảng biểu kỹ thuật.

  5. Chương trình nghiên cứu 10 năm tiếp theo được định hình ra sao?
    Tập trung phát triển mô hình tích hợp Viễn thám radar đa thời gian, Trí tuệ nhân tạo không gian (GeoAI) và Mô hình dựa trên tác tử (Agent-Based Model) nhằm giám sát và dự báo tự động biến động đường bờ và suy thoái đất ngập mặn dưới kịch bản nước biển dâng cực đoan.

Kết luận

Luận án tiến sĩ của Phạm Thị Làn đã hoàn thành xuất sắc các mục tiêu đề ra với 05 đóng góp mang tính đột phá:

  1. Xác lập cơ sở khoa học và quy trình công nghệ chuẩn mực tích hợp phương pháp phân loại định hướng đối tượng (OBIA) và phân vùng mẫu (per-field classification) trên tư liệu ảnh Landsat đa thời gian, chiết xuất chính xác chuỗi dữ liệu hiện trạng sử dụng đất 1989–2013 tại ven biển ĐBSH.
  2. Định lượng hóa thành công mối quan hệ hữu cơ chặt chẽ giữa các yếu tố tự nhiên (địa mạo, thổ nhưỡng, biến động bờ biển, khoảng cách đới bờ) với biến động sử dụng đất thông qua mô hình hồi quy Logistic đa biến.
  3. Làm sáng tỏ sự chi phối của động lực tiến hóa thùy châu thổ (thùy Ba Lạt bồi tụ và thùy Diêm Điền suy thoái) đối với sự phân hóa không gian của các mô hình kinh tế nông nghiệp - thủy sản.
  4. Ứng dụng thành công mô hình tích hợp Mạng nơ-ron nhân tạo (MLP) và CA-Markov trong nền tảng LCM để dự báo đáng tin cậy xu hướng chuyển dịch sử dụng đất đến năm 2021.
  5. Đề xuất hệ thống định hướng quy hoạch và giải pháp sử dụng đất bền vững, hài hòa giữa mục tiêu tăng trưởng kinh tế với bảo tồn sinh thái và thích ứng biến đổi khí hậu cho dải ven biển châu thổ sông Hồng.