Tổng quan về luận án

Quá trình đô thị hóa tại vùng kinh tế trọng điểm phía Nam nói chung và tỉnh Bình Dương nói riêng diễn ra với tốc độ nhanh chóng, trở thành động lực then chốt thúc đẩy chuyển dịch cơ cấu kinh tế từ nông nghiệp sang công nghiệp và dịch vụ. Nằm ở vị trí cửa ngõ chiến lược kết nối trực tiếp với Thành phố Hồ Chí Minh qua trục Quốc lộ 13 và sông Sài Gòn, thị xã Thuận An (nay là thành phố Thuận An) là tâm điểm phát triển công nghiệp - thương mại - dịch vụ đô thị năng động bậc nhất của tỉnh Bình Dương. Tuy nhiên, tốc độ đô thị hóa nhanh chóng và sự gia tăng dân số cơ học đã tạo ra sức ép to lớn lên quỹ đất, dẫn đến sự biến động mạnh mẽ về cơ cấu sử dụng đất. Luận án tiến sĩ chuyên ngành Quản lý đất đai (Mã số: 9.03) của nghiên cứu sinh Đào Đức Hưởng, dưới sự hướng dẫn khoa học của PGS. Nguyễn Hữu Ngữ và PGS. Huỳnh Văn Chương tại Trường Đại học Nông Lâm – Đại học Huế (2021), với đề tài "Ảnh hưởng của đô thị hóa đến biến động sử dụng đất trên địa bàn thị xã Thuận An, tỉnh Bình Dương", là một công trình nghiên cứu mang tính tiên phong, giải quyết trực tiếp bài toán xung đột không gian và tối ưu hóa sử dụng tài nguyên đất trong tiến trình chuyển đổi đô thị.

Khoảng trống nghiên cứu (research gap) cốt lõi được xác định rõ: phần lớn các nghiên cứu trước đây (như Đào Thị Thanh Lam, 2014; Nguyễn Hữu Ngữ và cs., 2016) chỉ dừng lại ở việc thống kê diện tích biến động đơn thuần bằng công nghệ GIS hoặc đánh giá định tính các tác động kinh tế - xã hội. Luận án này đã khắc phục triệt để khoảng trống đó bằng cách thiết lập một khung phân tích tích hợp định lượng đa chiều: kết hợp phương pháp phân tích thành phần chính (PCA) để lượng hóa mức độ đô thị hóa tại cấp cơ sở (56 ấp/khu phố), mô hình hồi quy đa biến ($P\text{-value} \le 0{,}05$) nhằm xác định chính xác các nhân tố động lực tác động, và ứng dụng chuỗi Markov (Markov Chain) để mô phỏng, dự báo biến động không gian sử dụng đất đến năm 2025 và 2030.

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

  • RQ1: Thực trạng, tốc độ và bản chất chuyển dịch cơ cấu sử dụng đất tại thị xã Thuận An diễn ra như thế nào qua các giai đoạn 2005–2010, 2010–2015 và 2015–2020?
    • H1: Tốc độ đô thị hóa tỷ lệ nghịch với quy mô diện tích đất nông nghiệp (NNP) và tỷ lệ thuận với sự mở rộng của đất ở đô thị (ODT), đất sản xuất kinh doanh (CSK) và đất giao thông (GT).
  • RQ2: Những nhân tố kinh tế - xã hội, hạ tầng và quản lý nào của đô thị hóa đóng vai trò chi phối và quyết định đến mức độ chuyển đổi mục đích sử dụng đất tại 56 ấp/khu phố?
    • H2: Sự tập trung công nghiệp, biến động dân số cơ học và hạ tầng giao thông là các thành phần chính (Principal Components) tạo ra sự phân hóa không gian về mức độ đô thị hóa và thúc đẩy chuyển dịch đất đai mạnh nhất.
  • RQ3: Xu hướng biến động cơ cấu sử dụng đất đến năm 2025 và 2030 sẽ diễn biến ra sao dưới tác động của các kịch bản đô thị hóa tiếp diễn?
    • H3: Ma trận xác suất chuyển dịch Markov phản ánh quy luật chuyển đổi ổn định và có thể dự báo chính xác xu thế suy giảm tiệm cận của quỹ đất nông nghiệp phục vụ cho công tác quy hoạch phát triển bền vững.

Khung lý thuyết của luận án được xây dựng dựa trên sự giao thoa giữa Lý thuyết Cực tăng trưởng (Growth Pole Theory - François Perroux, 1955), Lý thuyết Vị thế thuê đất (Bid Rent Theory - William Alonso, 1964) và Lý thuyết Chuyển đổi sử dụng đất (Land Use Transition Theory - Grainger, 1995; Lambin et al., 2001). Phạm vi nghiên cứu bao quát toàn bộ diện tích tự nhiên của thị xã Thuận An, phân tích chi tiết dữ liệu không gian và thuộc tính của 56 ấp/khu phố qua chuỗi thời gian 15 năm (2005–2020), dự báo định lượng cho giai đoạn 2020–2030. Kết quả nghiên cứu không chỉ đóng góp luận cứ khoa học cho ngành quản lý đất đai mà còn cung cấp công cụ dự báo chuẩn xác phục vụ trực tiếp cho công tác quy hoạch, kế hoạch sử dụng đất của tỉnh Bình Dương.

Literature Review và Positioning

Tổng quan y văn quốc tế và trong nước về mối quan hệ tương tác giữa đô thị hóa và biến động sử dụng đất cho thấy ba dòng nghiên cứu (research streams) chủ đạo:

  1. Dòng nghiên cứu lý thuyết hình thái và động lực đô thị: Khởi nguồn từ trường phái Chicago với Lý thuyết Sinh thái đô thị (Park & Burgess, 1925), phát triển qua các mô hình cấu trúc đô thị của Hoyt (1939), Harris & Ullman (1945), và được chuẩn hóa trong nghiên cứu của Đàm Trung Phường (2005), Đào Hoàng Tuấn (2008). Dòng nghiên cứu này phân tích quá trình chuyển dịch dân cư, lao động từ nông nghiệp sang phi nông nghiệp và sự bành trướng không gian của các trung tâm đô thị.
  2. Dòng nghiên cứu ứng dụng viễn thám và GIS trong đánh giá biến động lớp phủ/sử dụng đất (LULC): Fu & Ye (1995), Weng (2002), Huỳnh Văn Chương và cs. (2017), Nguyễn Hữu Ngữ và cs. (2016) đã khẳng định tính ưu việt của kỹ thuật chồng ghép bản đồ và viễn thám trong việc trích xuất ma trận biến động không gian, xác định sự biến mất của thảm phủ thực vật và diện tích đất nông nghiệp trước áp lực mở rộng bề mặt không thấm (impervious surface).
  3. Dòng nghiên cứu mô hình hóa và dự báo chuyển đổi đất đai: Clark (1965), Turner (1990), McMillen & McDonald (1991), Levinson & Chen (2005) đã tiên phong ứng dụng chuỗi Markov kết hợp với các mô hình kinh tế lượng không gian để mô phỏng xác suất chuyển dịch các loại đất theo thời gian.

Các cuộc tranh luận học thuật nổi bật tập trung vào hai luồng quan điểm đối lập:

  • Quan điểm thứ nhất (Đô thị hóa tập trung - Compact City): Cho rằng cần tập trung cao độ công nghiệp và dịch vụ vào các vùng lõi nhằm tối ưu hóa hiệu quả kinh tế và hạ tầng, nhưng chấp nhận việc chuyển đổi triệt để đất nông nghiệp vùng ven đô.
  • Quan điểm thứ hai (Đô thị hóa phân tán và đa trung tâm - Polycentric & Sustainable Development): Nhấn mạnh việc bảo tồn không gian xanh, duy trì nông nghiệp đô thị đa chức năng để bảo đảm cân bằng sinh thái và an ninh lương thực.

So sánh với các nghiên cứu quốc tế điển hình:

  • Nghiên cứu của Joanna Wilbers tại Amsterdam (Hà Lan) chứng minh rằng việc áp dụng "Chính sách hiệp ước" (Covenant Policy) giúp duy trì hơn 600 khu vườn sinh thái với 300 ha trên tổng số 21.907 ha diện tích thành phố, chứng minh nông nghiệp vẫn có thể tồn tại hài hòa trong lòng đô thị mật độ cao trên 20.000 người/$\text{km}^2$.
  • Nghiên cứu của Weng (2002) tại Đồng bằng sông Châu Giang (Trung Quốc) chỉ ra rằng quá trình công nghiệp hóa thần tốc đã làm biến dạng cảnh quan nông nghiệp truyền thống, đòi hỏi phải tích hợp GIS và mô hình Markov để kiểm soát tốc độ xâm lấn của đất công nghiệp và dân cư.
  • Nghiên cứu của Hoshino (2001) tại Nhật Bản đã áp dụng thành công phân tích thành phần chính (PCA) để xử lý hiện tượng đa cộng tuyến giữa các biến điều kiện tự nhiên và kinh tế - xã hội trước khi đưa vào mô hình dự báo sử dụng đất.

Luận án của Đào Đức Hưởng đã định vị chính xác vị trí học thuật bằng việc kết hợp mô hình phân tích đa biến PCA với kỹ thuật chồng ghép không gian GIS và chuỗi Markov. Nghiên cứu giải quyết bài toán phức tạp tại một vùng cận đô thị hóa có tốc độ công nghiệp hóa hàng đầu Việt Nam, thu hẹp khoảng cách giữa lý thuyết quy hoạch không gian và thực tiễn quản lý địa chính tại địa phương.

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

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

Luận án mở rộng và làm phong phú thêm Lý thuyết Vị thế thuê đất (Bid Rent Theory - Alonso, 1964) và Lý thuyết Cực tăng trưởng (Growth Pole Theory - Perroux, 1955) trong bối cảnh các nền kinh tế chuyển đổi tại Đông Nam Á. Nghiên cứu chứng minh rằng giá trị địa tô chênh lệch (differential rent) không chỉ đơn thuần phụ thuộc vào khoảng cách hình học tới trung tâm hành chính truyền thống, mà bị tái cấu trúc mạnh mẽ bởi mạng lưới giao thông kết nối liên vùng (Quốc lộ 13, đường Mỹ Phước – Tân Vạn) và sự hình thành các cực tăng trưởng công nghiệp tập trung (KCN VSIP 1, KCN Việt Hương).

                        ┌────────────────────────────────────────────────────────┐
                        │              TIẾN TRÌNH ĐÔ THỊ HÓA THUẬN AN             │
                        │ (Dân số cơ học, KCN/CCN, Hạ tầng kỹ thuật, TMDV-CNXD)   │
                        └───────────────────────────┬────────────────────────────┘
                                                    │
                                                    ▼
                        ┌────────────────────────────────────────────────────────┐
                        │          CƠ CHẾ TÁI CẤU TRÚC ĐỊA TÔ & KHÔNG GIAN       │
                        │  (Bid Rent Theory & Growth Pole: Cực công nghiệp/Giao  │
                        │   thông đẩy giá trị sinh lợi đất phi nông nghiệp tăng) │
                        └───────────────────────────┬────────────────────────────┘
                                                    │
                                                    ▼
                        ┌────────────────────────────────────────────────────────┐
                        │             BIẾN ĐỘNG CƠ CẤU SỬ DỤNG ĐẤT               │
                        │   - Đất Nông nghiệp (NNP) giảm mạnh                    │
                        │   - Đất Phi nông nghiệp (ODT, CSK, CDG, GT) tăng vọt   │
                        │   - Đất Chưa sử dụng (CSD) tiệm cận triệt tiêu         │
                        └───────────────────────────┬────────────────────────────┘
                                                    │
                                                    ▼
                        ┌────────────────────────────────────────────────────────┐
                        │       KHUNG ĐỊNH LƯỢNG TÍCH HỢP ĐÁNH GIÁ & DỰ BÁO      │
                        │   1. PCA: Lượng hóa mức độ ĐTH tại 56 ấp/khu phố       │
                        │   2. Multivariate Regression (p ≤ 0.05): Xác định lực  │
                        │   3. Markov Chain: Dự báo ma trận chuyển dịch đến 2030 │
                        └────────────────────────────────────────────────────────┘

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

  • Proposition 1 (P1): Dưới tác động lan tỏa từ cực tăng trưởng công nghiệp, đường cong vị thế trả giá thuê đất (bid-rent curve) của đất sản xuất kinh doanh phi nông nghiệp (CSK) và đất ở đô thị (ODT) vượt trội hoàn toàn so với đất nông nghiệp (NNP), dẫn đến quá trình chuyển đổi mục đích sử dụng đất mang tính tất yếu kinh tế.
  • Proposition 2 (P2): Sự biến động của không gian sử dụng đất diễn ra theo các giai đoạn phân kỳ: giai đoạn chuyển dịch diện rộng theo chiều ngang (horizontal expansion) được thay thế dần bởi giai đoạn tái cấu trúc sử dụng đất theo chiều sâu (vertical intensification) khi quỹ đất tự nhiên cạn kiệt.
  • Proposition 3 (P3): Quá trình chuyển dịch cơ cấu sử dụng đất mang tính chất ngẫu nhiên có điều kiện, tuân theo quy luật chuỗi Markov với ma trận xác suất chuyển dịch $[\gamma_{ij}]$ phản ánh hành vi chuyển đổi trạng thái của các mảnh thửa đất đai qua từng chu kỳ kinh tế - xã hội.

Sự chuyển dịch mô thức (paradigm shift) mà luận án mang lại là việc thay thế tư duy quản lý đất đai tĩnh (dựa trên các chỉ tiêu hành chính phân bổ cơ học) bằng mô thức Quản trị thích ứng động (Dynamic Adaptive Spatial Governance), tích hợp không gian hóa các chỉ số kinh tế - xã hội và xác suất toán học vào việc lập quy hoạch sử dụng đất.

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 lý thuyết: Lý thuyết Chuyển dịch sử dụng đất (Lambin & Meyfroidt), Lý thuyết Hệ thống phức hợp thích ứng (Complex Adaptive Systems), và Lý thuyết Quy hoạch không gian chức năng.

Tính độc đáo trong phương pháp tiếp cận thể hiện ở quy trình 3 giai đoạn:

  1. Giai đoạn 1 - Lượng hóa mức độ đô thị hóa bằng PCA: Thu thập hệ thống các chỉ tiêu phản ánh tỷ trọng sử dụng đất, cơ cấu lao động, thu nhập và mật độ dân số tại 56 ấp/khu phố. Ứng dụng PCA để rút gọn biến, loại trừ đa cộng tuyến và tính toán điểm chỉ số đô thị hóa chuẩn hóa cho từng đơn vị hành chính vi mô.
  2. Giai đoạn 2 - Lượng hóa tương quan không gian và nhân tố ảnh hưởng: Sử dụng mô hình hồi quy đa biến với kiểm định chặt chẽ ($P\text{-value} \le 0{,}05$) để xác lập mối quan hệ hàm số giữa mức độ đô thị hóa và diện tích các loại đất biến động.
  3. Giai đoạn 3 - Mô phỏng và dự báo xác suất bằng Chuỗi Markov: Thiết lập ma trận xác suất chuyển dịch $P = [\gamma_{ij}]$ từ dữ liệu chồng ghép bản đồ GIS qua các thời kỳ thực tế (2005–2010, 2010–2015, 2015–2020), kiểm định độ chính xác của mô hình dự báo năm 2020 trước khi ngoại suy dự báo cho năm 2025 và 2030.

Điều kiện biên (boundary conditions) của khung phân tích được xác lập rõ: áp dụng tối ưu cho các đô thị công nghiệp - dịch vụ chuyển đổi nhanh tại các vùng kinh tế động lực, nơi có hệ thống dữ liệu bản đồ địa chính chính quy và số liệu thống kê kinh tế - xã hội cấp cơ sở đầy đủ, đáng tin cậy.

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

Thiết kế nghiên cứu

Nghiên cứu đứng trên lập trường triết học Thực chứng luận (Positivism) kết hợp với Hiện thực phê phán (Critical Realism). Nghiên cứu thừa nhận sự tồn tại của các quy luật khách quan chi phối biến động không gian đất đai, đồng thời giải mã các cơ chế cấu trúc xã hội, thể chế và hành vi con người tác động đằng sau các biến số thống kê.

Thiết kế nghiên cứu hỗn hợp định lượng không gian (Spatial Quantitative Mixed-Methods) được triển khai qua thiết kế đa cấp độ (multi-level design):

  • Cấp độ vĩ mô (Macro-level): Toàn bộ địa bàn thị xã Thuận An và mối liên kết vùng với tỉnh Bình Dương và TP.HCM.
  • Cấp độ trung mô (Meso-level): 10 đơn vị hành chính cấp xã/phường (gồm 9 phường: An Phú, An Thạnh, Bình Chuẩn, Bình Hòa, Bình Nhâm, Hưng Định, Lái Thiêu, Thuận Giao, Vĩnh Phú và 1 xã: An Sơn).
  • Cấp độ vi mô (Micro-level): 56 ấp/khu phố – đơn vị phân tích cơ sở cho phép mô hình PCA bóc tách chi tiết sự phân hóa không gian đô thị hóa.

Quy trình nghiên cứu rigorous

Chiến lược thu thập dữ liệu bảo đảm tính chuẩn xác và toàn diện từ các nguồn dữ liệu chính thống:

  • Dữ liệu bản đồ và địa chính: Bản đồ địa chính chính quy tỷ lệ 1:500, 1:1.000, 1:2.000; bản đồ hiện trạng sử dụng đất các năm 2005, 2010, 2015 và 2020 của thị xã Thuận An do Sở Tài nguyên và Môi trường tỉnh Bình Dương thẩm định; dữ liệu mốc địa chính cơ sở.
  • Dữ liệu kinh tế - xã hội: Niên giám thống kê thị xã Thuận An, báo cáo tình hình giải quyết thủ tục hành chính về đất đai, cấp Giấy chứng nhận quyền sử dụng đất (GCNQSD), cơ cấu các ngành kinh tế, biến động lao động và thu nhập giai đoạn 2005–2020.
  • Quy trình kiểm chuẩn (Triangulation): Thực hiện tam giác giác hóa dữ liệu (data triangulation) giữa ảnh viễn thám/bản đồ số GIS, hồ sơ địa chính lưu trữ và kết quả điều tra, khảo sát thực địa tại các tổ dân phố.
  • Kiểm định độ tin cậy và giá trị đo lường: Hệ số tải nhân tố (Factor Loading) của các chỉ tiêu trong mô hình PCA đều đạt ngưỡng tiêu chuẩn $> 0{,}5$. Kiểm định tính phù hợp của dữ liệu thông qua chỉ số KMO (Kaiser-Meyer-Olkin $> 0{,}6$) và kiểm định Bartlett ($p < 0{,}001$). Độ tin cậy thang đo đạt chuẩn với hệ số Cronbach’s Alpha $> 0{,}7$.

Data và phân tích

Quy trình phân tích dữ liệu ứng dụng đồng bộ các công cụ phần mềm chuyên dụng: ArcGIS 10.2 (xử lý không gian, chồng ghép bản đồ, biên tập bản đồ biến động), SPSS (chạy PCA, phân tích tương quan và hồi quy đa biến), và MS Excel / Ngôn ngữ lập trình mô hình hóa (xây dựng ma trận Markov).

Phân tích PCA xử lý nhóm chỉ tiêu kinh tế - xã hội, lao động và cơ cấu đất đai:

  • Rút gọn các biến gốc về các trục thành phần chính có giá trị Eigenvalues $> 1$.
  • Tổng phương sai trích giải thích đạt từ $68%$ đến trên $80%$, phản ánh đầy đủ thông tin của tập dữ liệu ban đầu.

Mô hình chuỗi Markov được thiết lập dựa trên phương trình ma trận xác suất chuyển đổi: $$\mathbf{V}{t_2} = \mathbf{V}{t_1} \times \mathbf{P}$$ Trong đó:

  • $\mathbf{V}{t_1} = [V_1, V_2, \dots, V_k]{t_1}$ là vector diện tích các loại đất tại thời điểm đầu kỳ $t_1$.
  • $\mathbf{V}{t_2} = [V_1, V_2, \dots, V_k]{t_2}$ là vector diện tích các loại đất tại thời điểm cuối kỳ $t_2$.
  • $\mathbf{P} = [\gamma_{ij}]$ là ma trận xác suất chuyển dịch trạng thái với $\sum_{j=1}^k \gamma_{ij} = 1$ và $0 \le \gamma_{ij} \le 1$.

Các kiểm tra độ bền vững (robustness checks) được thực hiện bằng cách đối chiếu kết quả dự báo diện tích các loại đất năm 2020 từ chuỗi dữ liệu 2010–2015 với số liệu hiện trạng thực tế năm 2020. Kết quả kiểm chứng cho thấy sai số tương đối rất nhỏ ($< 3%$), khẳng định mô hình có độ tin cậy cao để tiến hành dự báo cho các mốc 2025 và 2030.

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

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

Thứ nhất, sự đảo chiều mang tính bước ngoặt trong cơ cấu 3 nhóm đất chính giai đoạn 2005–2020: Quá trình đô thị hóa thần tốc đã dẫn đến sự chuyển dịch mạnh mẽ từ nhóm đất nông nghiệp sang nhóm đất phi nông nghiệp. Đất nông nghiệp (NNP) sụt giảm nghiêm trọng từ diện tích chiếm ưu thế tuyệt đối trước năm 2005 xuống mức rất thấp vào năm 2020. Ngược lại, đất phi nông nghiệp tăng vọt, trong đó đất ở tại đô thị (ODT), đất sản xuất kinh doanh phi nông nghiệp (CSK), đất chuyên dùng (CDG) và đất giao thông (GT) tăng trưởng đột biến. Đất chưa sử dụng (CSD) bị khai thác triệt để và tiệm cận về mức 0 ha.

Thứ hai, phân hóa không gian mức độ đô thị hóa tại 56 ấp/khu phố: Kết quả chạy mô hình PCA giai đoạn 2010, 2015 và 2020 bóc tách rõ sự phân hóa không gian đô thị hóa thành 3 nhóm rõ rệt:

  • Nhóm đô thị hóa cao: Tập trung tại các phường trung tâm công nghiệp - dịch vụ như Lái Thiêu, An Phú, Thuận Giao, Bình Hòa với tỷ lệ lao động phi nông nghiệp đạt trên $90%$, mật độ dân số vượt $4.000\text{ người/km}^2$, và diện tích đất phi nông nghiệp chiếm đa số.
  • Nhóm đô thị hóa trung bình: Các phường An Thạnh, Bình Chuẩn, Vĩnh Phú, Hưng Định đang trong giai đoạn chuyển đổi mạnh mẽ.
  • Nhóm đô thị hóa thấp/ven đô sinh thái: Điển hình là xã An Sơn và phường Bình Nhâm, nơi vẫn duy trì tỷ trọng đất vườn cây ăn trái và nông nghiệp sinh thái tương đối lớn do vị trí nằm ven sông Sài Gòn.

Thứ ba, xác định phương trình hồi quy và các nhân tố chi phối có ý nghĩa thống kê: Mô hình hồi quy đa biến xác định các biến độc lập có hệ số hồi quy mang ý nghĩa thống kê cao ($P\text{-value} \le 0{,}05$) chi phối biến động đất đai:

  • Tốc độ tăng trưởng dân số cơ học và mật độ dân số (hệ số tác động dương mạnh nhất đến tăng diện tích đất ở ODT).
  • Diện tích khu, cụm công nghiệp và cơ sở sản xuất kinh doanh (thúc đẩy trực tiếp sự mở rộng của đất CSK và đất giao thông GT).
  • Tỷ lệ chuyển dịch lao động từ nông nghiệp sang công nghiệp - xây dựng (CNXD) và thương mại - dịch vụ (TMDV).

Thứ tư, kết quả dự báo chuỗi Markov giai đoạn 2020–2025 và 2025–2030: Dự báo cho thấy xu hướng thu hẹp đất nông nghiệp vẫn tiếp tục diễn ra nhưng với tốc độ chậm lại do quỹ đất nông nghiệp có thể chuyển đổi đã cạn dần. Đất ở đô thị (ODT) và đất phát triển hạ tầng giao thông (GT) tiếp tục chiếm tỷ trọng chi phối trong cơ cấu sử dụng đất năm 2025 và hướng tới năm 2030, tiệm cận ngưỡng đô thị hóa hoàn toàn của một thành phố công nghiệp - dịch vụ văn minh, hiện đại.

Thứ năm, phát hiện hiện tượng phân mảnh không gian "da báo" mang tính nghịch lý: Trong giai đoạn 2005–2015, việc quy hoạch xây dựng chưa bám sát nhu cầu thực tế và thiếu kiểm soát đã tạo ra tình trạng chuyển mục đích sử dụng đất phân tán, manh mún, hình thành các khu đất xen cài giữa nông nghiệp và nhà ở tự phát, làm suy giảm hiệu quả sử dụng đất và gia tăng chi phí đầu tư hạ tầng kỹ thuật đồng bộ.

Implications đa chiều

  • Về mặt lý thuyết: Luận án cung cấp bằng chứng thực nghiệm củng cố lý thuyết chuyển đổi sử dụng đất vùng ven đô tại các quốc gia đang phát triển; hoàn thiện phương pháp luận lượng hóa đô thị hóa cấp vi mô thông qua tích hợp PCA - GIS - Markov.
  • Về mặt phương pháp luận: Thiết lập quy trình chuẩn tắc xử lý dữ liệu viễn thám kết hợp với thống kê đa biến có thể chuyển giao và áp dụng trực tiếp cho các địa phương khác có điều kiện tương đồng tại vùng Đông Nam Bộ và đồng bằng sông Hồng.
  • Về mặt thực tiễn và chính sách:
    • Chính sách quy hoạch không gian: Đề xuất chuyển từ phát triển đô thị dàn trải theo chiều ngang sang phát triển đô thị nén, tập trung, tăng cường hệ số sử dụng đất, khuyến khích nhà ở cao tầng tại các phường lõi.
    • Chính sách bảo tồn sinh thái: Cần khoanh vùng bảo vệ nghiêm ngặt diện tích đất vườn cây ăn trái đặc sản tại xã An Sơn và phường Bình Nhâm gắn với phát triển du lịch sinh thái ven sông Sài Gòn, tạo "lá phổi xanh" điều hòa khí hậu đô thị.
    • Chính sách hạ tầng kỹ thuật và xã hội: Định hướng quỹ đất hợp lý cho giao thông tĩnh, công viên cây xanh, trường học, trạm y tế nhằm giải quyết tình trạng quá tải hạ tầng do gia tăng dân số cơ học tại các phường tập trung đông công nhân như An Phú, Bình Hòa.

Limitations và Future Research

Nghiên cứu thẳng thắn thừa nhận một số hạn chế mang tính khách quan:

  1. Giới hạn về dữ liệu chuỗi thời gian liên tục: Nghiên cứu dựa trên các mốc kiểm kê đất đai 5 năm/lần (2005, 2010, 2015, 2020), do đó chưa phản ánh được các biến động mang tính tức thời giữa các năm đơn lẻ trong chu kỳ.
  2. Giới hạn về các biến phi cấu trúc: Chưa lượng hóa được đầy đủ tác động của tâm lý đầu cơ đất đai, biến động giá đất thị trường phi chính thức và sự thay đổi đột xuất của các dự án hạ tầng giao thông liên vùng lớn vào mô hình toán học.
  3. Giới hạn về yếu tố sinh thái - khí hậu: Mô hình Markov chưa tích hợp biến động thủy văn, ngập lụt do triều cường sông Sài Gòn và biến đổi khí hậu tác động đến sự thích nghi sử dụng đất.

Chương trình nghiên cứu tiếp theo (Future Research Agenda):

  • Tích hợp mô hình tự động tế bào (Cellular Automata - CA) và Mô hình đa tác tử (Agent-Based Model - ABM): Kết hợp CA-Markov để mô phỏng chính xác vị trí không gian (spatial allocation) của từng ô lưới biến động đất đai thay vì chỉ dự báo quy mô diện tích.
  • Tích hợp yếu tố rủi ro ngập lụt và biến đổi khí hậu: Đánh giá tổn thương và khả năng thích ứng của quỹ đất đô thị ven sông trong bối cảnh nước biển dâng và triều cường.
  • Mở rộng phạm vi nghiên cứu liên vùng: Nghiên cứu so sánh chuỗi biến động đất đai dọc hành lang kinh tế Đông Nam Bộ (TP.HCM – Thuận An – Dĩ An – Thủ Dầu Một – Biên Hòa).

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

Luận án tạo ra tác động sâu rộng trên cả 5 phương diện:

  • Ảnh hưởng học thuật (Academic Impact): Đóng góp công trình nghiên cứu chuẩn mực cho ngành Quản lý đất đai và Địa lý kinh tế; mở ra hướng nghiên cứu kết hợp viễn thám, GIS và mô hình toán kinh tế lượng cho các nghiên cứu sinh và học viên cao học.
  • Chuyển đổi thực tiễn quản lý địa phương (Industry & Planning Transformation): Cung cấp cơ sở dữ liệu số hóa và ma trận biến động đất đai chuẩn xác, hỗ trợ trực tiếp cho Phòng Tài nguyên và Môi trường thị xã Thuận An và Sở Tài nguyên và Môi trường tỉnh Bình Dương trong việc rà soát quy hoạch.
  • Ảnh hưởng chính sách công (Policy Influence): Là tài liệu tham khảo khoa học quan trọng cho UBND tỉnh Bình Dương trong quá trình xây dựng Quy hoạch tỉnh Bình Dương thời kỳ 2021–2030, tầm nhìn đến năm 2050, đặc biệt trong phương án phân bổ chỉ tiêu đất khu công nghiệp, đất ở đô thị và bảo vệ hành lang xanh.
  • Lợi ích xã hội (Societal Benefits): Góp phần hạn chế tình trạng lãng phí đất đai, giảm thiểu xung đột thu hồi đất nông nghiệp, đề xuất giải pháp tạo việc làm và chuyển đổi nghề nghiệp bền vững cho bộ phận nông dân bị thu hồi đất.
  • Ý nghĩa quốc tế (International Relevance): Cung cấp bức tranh thực nghiệm điển hình từ một "vùng lõi công nghiệp" của Việt Nam vào mạng lưới nghiên cứu đô thị hóa toàn cầu của các tổ chức quốc tế như WB, ADB, UN-Habitat.

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

  • Nghiên cứu sinh và Giới học thuật: Tiếp cận phương pháp luận tích hợp PCA - GIS - Markov hoàn chỉnh; kế thừa hệ thống cơ sở lý luận về chuyển dịch đất đai đô thị.
  • Cơ quan Quản lý Nhà nước và Hoạch định Chính sách (Sở TN&MT, Sở Xây dựng, UBND cấp huyện/tỉnh): Nhận được luận cứ khoa học và công cụ dự báo định lượng chính xác để ban hành các quyết định giao đất, cho thuê đất, chuyển mục đích sử dụng đất và lập quy hoạch, kế hoạch sử dụng đất có tính khả thi cao.
  • Các Đơn vị Tư vấn Quy hoạch và Doanh nghiệp Phát triển Đô thị: Nắm bắt được quy luật và xu hướng chuyển đổi không gian, xác định các phân vùng tiềm năng để phát triển dự án nhà ở, khu dịch vụ thương mại đồng bộ.
  • Cộng đồng Dân cư Địa phương: Hưởng lợi gián tiếp từ việc nâng cao chất lượng môi trường sống, hạ tầng kỹ thuật và xã hội được đầu tư hợp lý, hạn chế ngập úng, ô nhiễm môi trường và ách tắc giao thô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ào?

Đóng góp độc đáo nhất là việc mở rộng Lý thuyết Vị thế thuê đất (Bid Rent Theory - Alonso, 1964) và Lý thuyết Chuyển đổi sử dụng đất (Lambin et al., 2001) trong bối cảnh đô thị công nghiệp vệ tinh Đông Nam Á. Luận án chứng minh rằng giá trị địa tô không chỉ bị chi phối bởi khoảng cách tới trung tâm đô thị truyền thống mà bị tái định hình bởi các cực tăng trưởng hạ tầng kỹ thuật và khu công nghiệp tập trung. Nghiên cứu đã toán học hóa mối quan hệ giữa các chỉ số đô thị hóa đa biến với tốc độ biến động không gian đất đai thông qua khung tích hợp PCA và chuỗi Markov.

2. Điểm mới về mặt phương pháp nghiên cứu so với các công trình trước đây là gì?

So với các nghiên cứu trước đây (như Đào Thị Thanh Lam, 2014 chỉ dùng thống kê GIS thuần túy; Huỳnh Văn Chương và cs., 2017 chỉ chạy Markov cấp huyện), điểm đột phá của luận án là:

  • Triển khai phân tích PCA chi tiết đến cấp độ 56 ấp/khu phố, cho phép bóc tách bản chất đô thị hóa vi mô mà các nghiên cứu cấp huyện/tỉnh bị làm phẳng (spatial aggregation bias).
  • Kết hợp mô hình hồi quy đa biến với điều kiện kiểm định nghiêm ngặt ($P\text{-value} \le 0{,}05$) nhằm loại bỏ biến nhiễu trước khi nạp dữ liệu vào ma trận xác suất chuyển đổi Markov, nâng cao độ chính xác dự báo lên mức sai số $< 3%$.

3. Phát hiện nào trong dữ liệu gây bất ngờ và có tính phản trực giác nhất?

Phát hiện bất ngờ nhất là hiện tượng suy thoái và phân mảnh đất nông nghiệp diễn ra trước khi đất thực sự được xây dựng (Pre-development Land Fragmentation). Dữ liệu chồng ghép bản đồ cho thấy hàng trăm hecta đất nông nghiệp ven đô bị bỏ hoang hoặc phân lô manh mún do tâm lý chờ đợi quy hoạch và chuyển mục đích ngầm, dẫn đến việc đất đai bị suy giảm nghiêm trọng chức năng sinh thái và năng suất cây trồng từ nhiều năm trước khi có quyết định thu hồi đất chính thức của cơ quan nhà nước.

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ó. Luận án cung cấp đầy đủ:

  • Bộ công thức toán học tính toán tỷ lệ đô thị hóa ($T$), tốc độ đô thị hóa ($M$), chỉ số đô thị - nông thôn ($URR_t$), hệ phương trình ma trận xác suất chuyển đổi chuỗi Markov.
  • Danh mục biến số chuẩn hóa đầu vào cho mô hình PCA (tỷ trọng đất đai, cơ cấu lao động, thu nhập, mật độ dân số) và tiêu chí chọn hệ số tải Factor Loading $> 0{,}5$.
  • Quy trình chuẩn hóa dữ liệu bản đồ địa chính trên ArcGIS 10.2 và quy trình kiểm định hồi quy đa biến trên SPSS với $P\text{-value} \le 0{,}05$.

5. Luận án đã phác thảo chương trình nghiên cứu 10 năm tiếp theo như thế nào?

Chương trình 10 năm tiếp theo tập trung vào 3 trọng tâm:

  1. Phát triển mô hình dự báo không gian thời gian thực (Real-time Spatial Simulation) kết hợp viễn thám độ phân giải cao và trí tuệ nhân tạo (AI/Machine Learning) để tự động nhận diện biến động đất đai.
  2. Xây dựng khung đánh giá dịch vụ hệ sinh thái đất (Soil Ecosystem Services Valuation) lượng hóa giá trị mất mát khi chuyển đổi đất nông nghiệp sang đất công nghiệp.
  3. Nghiên cứu cơ chế tài chính đất đai và công cụ thu hồi giá trị gia tăng từ đất (Land Value Capture) phục vụ tái đầu tư hạ tầng xanh thích ứng với biến đổi khí hậu tại lưu vực sông Sài Gòn.

Kết luận

Luận án tiến sĩ của Đào Đức Hưởng đã hoàn thành xuất sắc các mục tiêu nghiên cứu đề ra với 6 đóng góp cốt lõi:

  1. Lượng hóa toàn diện thực trạng đô thị hóa và chuyển dịch cơ cấu kinh tế - sử dụng đất tại thị xã Thuận An giai đoạn 2005–2020 qua chuỗi số liệu thực chứng đáng tin cậy.
  2. Xây dựng thành công bản đồ biến động sử dụng đất chi tiết qua 3 thời kỳ (2005–2010, 2010–2015, 2015–2020) bằng công nghệ GIS, bóc tách chính xác dòng chuyển dịch giữa các nhóm đất.
  3. Ứng dụng thành công phương pháp PCA để phân hạng và lập bản đồ mức độ đô thị hóa tại 56 ấp/khu phố, khắc phục triệt để hạn chế cào bằng dữ liệu không gian.
  4. Xác lập mô hình hồi quy đa biến ($P\text{-value} \le 0{,}05$) chỉ rõ các động lực cốt lõi chi phối biến động đất đai, trong đó dân số cơ học và hạ tầng KCN đóng vai trò quyết định.
  5. Ứng dụng chuỗi Markov dự báo chính xác cơ cấu sử dụng đất năm 2025 và 2030, làm cơ sở khoa học định hình quy hoạch sử dụng đất dài hạn.
  6. Đề xuất hệ thống giải pháp quản lý đất đai và quy hoạch đồng bộ, hài hòa giữa mục tiêu phát triển đô thị công nghiệp nén và bảo tồn dải sinh thái nông nghiệp ven sông Sài Gòn.

Công trình tạo tiền đề khai mở 3 hướng nghiên cứu mới về mô phỏng tế bào không gian CA-Markov, kinh tế học địa tô ven đô, và quản trị đất đai thích ứng biến đổi khí hậu; khẳng định vị thế học thuật vững chắc và giá trị thực tiễn lâu dài cho sự phát triển bền vững của hệ thống đô thị Việt Nam.