Tổng quan nghiên cứu

Quá trình đô thị hóa nhanh chóng tại thành phố Hải Phòng đã tạo ra áp lực chuyển dịch cơ cấu sử dụng đất vô cùng mạnh mẽ, đặc biệt là tại các khu vực ven đô tiếp giáp giữa đô thị và nông thôn. Huyện Kiến Thụy có tổng diện tích tự nhiên 10.751,89 ha cùng quy mô dân số trên 12,5 vạn người, phân bổ tại 17 xã và 1 thị trấn. Sau khi điều chỉnh địa giới hành chính để thành lập các quận Dương Kinh và Đồ Sơn theo Nghị định số 145/2007/NĐ-CP của Chính phủ, Kiến Thụy đóng vai trò là cầu nối kinh tế - sinh thái quan trọng ở phía Đông Nam thành phố Hải Phòng.

Vấn đề nghiên cứu trọng tâm đặt ra là các phương pháp thống kê và kiểm kê đất đai truyền thống theo chu kỳ 5 năm thường tốn kém chi phí, mất nhiều thời gian và không thể hiện trực quan không gian biến động trạng thái lớp phủ. Giai đoạn 2005 - 2010 chứng kiến diện tích đất sản xuất nông nghiệp của huyện sụt giảm 284,3 ha (tương ứng giảm 2,64% tổng diện tích tự nhiên), trong khi diện tích đất phi nông nghiệp tăng nhanh thêm 155,1 ha nhằm phục vụ mở rộng hạ tầng và khu dân cư. Sự thay đổi nhanh chóng này đòi hỏi một công cụ giám sát và dự báo động lực học không gian có độ chính xác cao.

Mục tiêu cụ thể của nghiên cứu là xây dựng quy trình ứng dụng công nghệ viễn thám đa thời gian kết hợp hệ thống thông tin địa lý (GIS) và mô hình chuỗi Markov kết hợp mạng tự động (Cellular Automata - CA) nhằm phân tích biến động lớp phủ mặt đất giai đoạn 2000 - 2005 - 2010, đồng thời mô hình hóa dự báo xu thế biến đổi không gian đến năm 2015 và 2020. Nghiên cứu mang ý nghĩa thực tiễn to lớn khi cung cấp bộ dữ liệu dự báo không gian đạt độ tin cậy trên 85%, hỗ trợ đắc lực cho các nhà hoạch định chính sách tối ưu hóa quy hoạch sử dụng đất trên toàn bộ 107,53 km² lãnh thổ huyện Kiến Thụy theo hướng phát triển bền vững.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Nghiên cứu được xây dựng dựa trên sự giao thoa của ba nền tảng lý thuyết và mô hình tính toán không gian hiện đại:

  1. Lý thuyết phân loại lớp phủ mặt đất (FAO): Dựa trên định nghĩa chuẩn của Tổ chức Lương thực và Nông nghiệp Liên Hợp Quốc, lớp phủ mặt đất là những đối tượng vật chất quan sát được trên bề mặt Trái Đất. Nghiên cứu phân định rõ sự khác biệt giữa lớp phủ mặt đất (trạng thái sinh lý - vật lý thực tế) và hiện trạng sử dụng đất (mục đích kinh tế - xã hội của con người) để thiết lập hệ thống 5 lớp chuyên đề: đất sản xuất nông nghiệp, đất lâm nghiệp, đất ở dân cư, mặt nước và đất chuyên dùng.
  2. Lý thuyết quá trình ngẫu nhiên chuỗi Markov (Markov Chain): Mô hình hóa quá trình chuyển dịch trạng thái không gian rời rạc qua thời gian. Xác suất chuyển đổi của một pixel đất đai sang trạng thái khác tại thời điểm tương lai chỉ phụ thuộc duy nhất vào trạng thái hiện tại của nó mà không phụ thuộc vào quá khứ, cho phép tính toán ma trận xác suất chuyển dịch giữa các thời kỳ 2000 - 2005 và 2005 - 2010.
  3. Mô hình Mạng tự động tế bào (Cellular Automata - CA): Khắc phục nhược điểm phi không gian của chuỗi Markov bằng cách bổ sung quy luật chuyển tiếp cục bộ theo cấu trúc lưới raster 2 chiều. Trạng thái của mỗi ô tế bào (cell) được cập nhật đồng thời dựa trên trạng thái của 4 ô lân cận (quy tắc Von Neumann) hoặc 8 ô lân cận (quy tắc Moore).
  4. Phương pháp Đánh giá đa chỉ tiêu (Multi-Criteria Evaluation - MCE): Đóng vai trò thiết lập các ngưỡng không gian thích nghi, tích hợp các nhân tố tự nhiên, vị trí địa lý ven sông Văn Úc (14,75 km), sông Đa Độ (29 km), mạng lưới giao thông và hành lang quy hoạch phát triển kinh tế - xã hội đến năm 2020.

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

Nghiên cứu sử dụng nguồn dữ liệu viễn thám đa thời gian độ phân giải cao và trung bình gồm ảnh vệ tinh Landsat TM/ETM+ (độ phân giải không gian 30 m) và ảnh vệ tinh SPOT 3/5 (độ phân giải không gian từ 15 m đến 20 m) thu chụp tại 3 mốc thời gian: năm 2000, 2005 và 2010. Dữ liệu bổ trợ bao gồm bản đồ địa hình tỷ lệ 1:25.000, bản đồ hiện trạng sử dụng đất năm 2005 và 2010, cùng số liệu niên giám thống kê huyện Kiến Thụy.

Cỡ mẫu nghiên cứu được thiết lập với 250 điểm kiểm chứng thực địa trên toàn bộ 18 xã, thị trấn. Phương pháp chọn mẫu là chọn mẫu phân tầng ngẫu nhiên (stratified random sampling) theo từng loại hình lớp phủ. Phương pháp xử lý ảnh số sử dụng thuật toán phân loại có kiểm định Maximum Likelihood (xác suất cực đại) trên phần mềm IDRISI Andes 15. Lý do lựa chọn phương pháp này là khả năng tối ưu hóa hàm mật độ xác suất của các vector phổ phản xạ, giúp phân tách hiệu quả các lớp phủ có đặc tính phổ gần nhau. Sau đó, phương pháp so sánh sau phân loại (post-classification comparison) được ứng dụng để lập ma trận chéo biến động diện tích, kết hợp module Land Change Modeler (LCM) để chạy mô hình dự báo không gian CA-Markov cho giai đoạn 2010 - 2020.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Quá trình phân tích dữ liệu không gian và thống kê biến động lớp phủ giai đoạn 2000 - 2010 trên địa bàn huyện Kiến Thụy mang lại 4 phát hiện cốt lõi:

  • Sự suy giảm liên tục của quỹ đất sản xuất nông nghiệp: Năm 2005, đất nông nghiệp chiếm 5.623,6 ha (52,30% tổng diện tích tự nhiên). Đến năm 2010, diện tích này giảm xuống còn 5.339,3 ha (chiếm 49,66%), tức giảm ròng 284,3 ha chỉ trong 5 năm. Xu hướng chuyển dịch chủ yếu diễn ra tại các xã ven đô như Thuận Thiên, Hữu Bằng và Thanh Sơn.
  • Sự phục hồi và gia tăng ấn tượng của diện tích rừng ngập mặn: Đất lâm nghiệp tăng trưởng vượt bậc từ 337,5 ha năm 2005 lên 646,2 ha vào năm 2010, tương ứng mức tăng 308,7 ha (tăng 91,47%). Diện tích tăng thêm tập trung gần như toàn bộ tại khu vực bãi bồi ven biển xã Đại Hợp (chiếm 632,2 ha năm 2010) nhờ các chương trình trồng mới rừng ngập mặn phòng hộ.
  • Tốc độ mở rộng nhanh của đất phi nông nghiệp và đất chuyên dùng: Tổng diện tích đất phi nông nghiệp tăng từ 3.934,5 ha (năm 2005) lên 4.089,5 ha (năm 2010), tăng thêm 155,1 ha. Trong đó, đất chuyên dùng tăng mạnh từ 1.659,4 ha lên 1.784,1 ha (tăng 124,7 ha) và đất ở tăng từ 1.173,4 ha lên 1.187,1 ha (tăng 13,7 ha) để đáp ứng nhu cầu mở rộng các trục lộ giao thông và công trình công cộng.
  • Chuyển đổi mạnh mẽ nhóm đất chưa sử dụng sang nuôi trồng thủy sản: Đất chưa sử dụng và mặt nước ven biển giảm mạnh 366,1 ha (từ 450,0 ha năm 2005 xuống 84,0 ha năm 2010). Phần lớn diện tích này đã được chuyển hóa thành các đầm nuôi trồng thủy sản nước lợ và mặn, nâng tổng diện tích nuôi trồng thủy sản của huyện từ 404,9 ha lên 591,8 ha (tăng 186,9 ha).

Thảo luận kết quả

Nguyên nhân chính dẫn đến sự biến đổi mạnh mẽ của lớp phủ mặt đất tại Kiến Thụy xuất phát từ sự tương tác phức tạp giữa các yếu tố quy hoạch kinh tế - xã hội và tác động tự nhiên. Việc triển khai các dự án giao thông trọng điểm như đường cao tốc Hà Nội - Hải Phòng, nâng cấp đường tỉnh cùng định hướng chuyển đổi huyện thuần nông sang công nghiệp - dịch vụ đã thúc đẩy quá trình chuyển đổi đất lúa sang đất phi nông nghiệp. Ngược lại, điều kiện tự nhiên bờ biển dài khoảng 4 km tại xã Đại Hợp cùng chính sách bảo vệ đê điều cửa sông Văn Úc đã tạo điều kiện mở rộng thành công vành đai rừng ngập mặn ven biển.

Khi so sánh với các nghiên cứu biến động đô thị tại thành phố Đà Nẵng hay thành phố Vinh, xu thế chuyển đổi đất nông nghiệp ven đô tại Kiến Thụy có nét tương đồng về tốc độ nhưng khác biệt rõ rệt về tính bền vững sinh thái thông qua việc chủ động gia tăng diện tích rừng phòng hộ. Toàn bộ dữ liệu chuyển dịch không gian này có thể được trình bày và trực quan hóa sinh động thông qua bảng ma trận chuyển dịch trạng thái 5x5 cùng biểu đồ cột kép so sánh tỷ trọng cơ cấu đất đai 2005 - 2010. Kết quả mô phỏng trên bản đồ raster 30m x 30m của mô hình CA-Markov năm 2010 khi kiểm chứng với ảnh phân loại thực tế cho chỉ số tương đồng không gian cao, khẳng định tính khả thi vượt trội của phương pháp trong dự báo lớp phủ đến năm 2020.

Đề xuất và khuyến nghị

Nhằm quản lý và sử dụng bền vững tài nguyên đất đai huyện Kiến Thụy giai đoạn 2015 - 2020, nghiên cứu đề xuất 4 nhóm giải pháp hành động cụ thể:

  • Thiết lập ranh giới bảo vệ nghiêm ngặt quỹ đất nông nghiệp chất lượng cao: Ủy ban nhân dân huyện Kiến Thụy phối hợp với Phòng Tài nguyên và Môi trường cần khoanh định và giữ ổn định tối thiểu 4.800 ha đất trồng 2 vụ lúa và rau màu chuyên canh đến năm 2020; ban hành quy định nghiêm cấm việc chuyển mục đích sử dụng đất nông nghiệp màu mỡ sang đất ở tự phát dọc các tuyến đường liên xã.
  • Mở rộng và nâng cao chức năng phòng hộ của hệ sinh thái rừng ngập mặn: Sở Nông nghiệp và Phát triển nông thôn thành phố Hải Phòng kết hợp cùng chính quyền xã Đại Hợp triển khai kế hoạch khoanh nuôi, bảo vệ và trồng mới rừng ngập mặn bãi bồi, hướng tới mục tiêu duy trì ổn định quy mô trên 650 ha rừng phòng hộ ven biển giai đoạn 2015 - 2020 nhằm ứng phó với nước biển dâng và triều cường.
  • Kiểm soát chặt chẽ chỉ tiêu mở rộng đất phi nông nghiệp theo quy hoạch: Ủy ban nhân dân huyện Kiến Thụy cần giới hạn tổng diện tích đất phi nông nghiệp dưới ngưỡng 4.500 ha vào năm 2020; ưu tiên bố trí các khu công nghiệp, cụm công nghiệp tập trung bám sát hành lang kinh tế sông Đa Độ và các trục cao tốc đã được phê duyệt trong Nghị quyết Đại hội Đảng bộ huyện.
  • Xây dựng hệ thống cơ sở dữ liệu GIS và giám sát lớp phủ mặt đất định kỳ: Phòng Tài nguyên và Môi trường huyện Kiến Thụy cần chuyển giao và ứng dụng công nghệ viễn thám kết hợp phần mềm GIS mã nguồn mở hoặc thương mại để cập nhật biến động đất đai hàng năm (thay vì 5 năm một lần), duy trì độ chính xác của cơ sở dữ liệu địa chính số hóa trên toàn bộ 10.751,89 ha diện tích tự nhiên.

Đối tượng nên tham khảo luận văn

Nội dung nghiên cứu và khung phương pháp luận của luận văn mang lại giá trị thực tiễn cho 4 nhóm đối tượng chính:

  • Cơ quan quản lý nhà nước về đất đai và quy hoạch: Ủy ban nhân dân huyện Kiến Thụy, Sở Tài nguyên và Môi trường, Sở Xây dựng thành phố Hải Phòng có thể sử dụng kết quả dự báo không gian năm 2020 làm căn cứ định lượng để lập, điều chỉnh quy hoạch sử dụng đất và kiểm soát biến động đô thị ven đô.
  • Chuyên viên kỹ thuật GIS và Viễn thám: Các kỹ sư trắc địa, địa chính và cán bộ kỹ thuật tại các trung tâm viễn thám có thể khai thác quy trình tích hợp giữa thuật toán Maximum Likelihood, phân tích chuỗi Markov và module CA-Markov trong IDRISI Andes 15 để áp dụng cho các bài toán phân tích không gian tương tự.
  • Giảng viên, nhà nghiên cứu và học viên cao học: Các cá nhân nghiên cứu thuộc chuyên ngành Địa chính, Quản lý Đất đai, Địa lý Tài nguyên tại các trường đại học có thể sử dụng luận văn như tài liệu tham khảo chuẩn mực về mô hình hóa biến đổi không gian cảnh quan và đánh giá đa chỉ tiêu MCE.
  • Các tổ chức quản lý môi trường và thích ứng biến đổi khí hậu: Các ban quản lý rừng phòng hộ và tổ chức bảo tồn sinh thái ven biển có thể ứng dụng số liệu diễn thế rừng ngập mặn Đại Hợp để xây dựng các dự án bảo vệ dải ven bờ cửa sông Văn Úc.

Câu hỏi thường gặp

  • Mô hình CA-Markov mang lại ưu điểm vượt trội nào so với các phương pháp dự báo truyền thống? Mô hình kết hợp hoàn hảo khả năng tính toán xác suất biến đổi theo thời gian của chuỗi Markov với năng lực mô phỏng động lực học không gian theo cấu trúc lân cận của Cellular Automata. Nhờ đó, mô hình không chỉ dự báo chính xác quy mô diện tích mà còn chỉ rõ vị trí không gian cụ thể của từng pixel đất biến đổi trên bản đồ số.
  • Tại sao cần tích hợp phương pháp đánh giá đa chỉ tiêu MCE vào mô hình dự báo? Thuật toán Markov thuần túy chỉ thực hiện nội suy tuyến tính thống kê mà không tính đến các yếu tố cản trở hay thúc đẩy thực tế. Phương pháp MCE cho phép tích hợp các yếu tố tự nhiên, khoảng cách đến sông ngòi (như sông Đa Độ 29 km), mạng lưới đường giao thông và chính sách quy hoạch nhằm thiết lập các ngưỡng thích nghi không gian chuẩn xác.
  • Nguyên nhân chính khiến đất nông nghiệp huyện Kiến Thụy giảm 284,3 ha giai đoạn 2005 - 2010 là gì? Sự sụt giảm này chủ yếu do tác động của quá trình đô thị hóa nhanh chóng tại Hải Phòng và việc điều chỉnh địa giới hành chính theo Nghị định số 145/2007/NĐ-CP. Đất nông nghiệp được thu hồi để phát triển cơ sở hạ tầng giao thông, công trình công cộng và mở rộng diện tích đất ở dân cư từ 1.173,4 ha lên 1.187,1 ha.
  • Phương pháp kiểm định độ chính xác của bản đồ phân loại ảnh vệ tinh được thực hiện ra sao? Nghiên cứu sử dụng 250 điểm kiểm chứng thực tế tại địa phương kết hợp phân tích ma trận sai số (confusion matrix) trên ảnh SPOT và Landsat đã nắn chỉnh hình học. Thuật toán phân loại xác suất cực đại Maximum Likelihood kết hợp giải đoán chuyên gia bằng mắt giúp loại bỏ triệt để hiện tượng nhầm lẫn phổ, đảm bảo độ chính xác trên 85%.
  • Kết quả dự báo biến động lớp phủ đến năm 2020 giúp ích gì cho quy hoạch của huyện Kiến Thụy? Bản đồ dự báo cung cấp kịch bản trực quan về sự thu hẹp của đất nông nghiệp và hướng mở rộng của đất phi nông nghiệp, giúp chính quyền chủ động phân bổ chỉ tiêu sử dụng đất, tránh quy hoạch chồng chéo và bảo vệ vững chắc 646,2 ha rừng ngập mặn phòng hộ ven biển.

Kết luận

  • Nghiên cứu đã xây dựng thành công quy trình công nghệ tích hợp Viễn thám đa thời gian, GIS và mô hình CA-Markov trong việc đánh giá và dự báo biến động lớp phủ mặt đất tại huyện Kiến Thụy.
  • Định lượng chính xác bức tranh biến động giai đoạn 2005 - 2010 trên tổng diện tích 10.751,89 ha, nổi bật với sự suy giảm 284,3 ha đất sản xuất nông nghiệp và sự gia tăng vượt bậc 308,7 ha đất rừng ngập mặn ven biển.
  • Ứng dụng thành công phương pháp đánh giá đa chỉ tiêu MCE để thiết lập các ngưỡng không gian, tạo ra bản đồ dự báo biến đổi lớp phủ mặt đất huyện Kiến Thụy năm 2020 có độ tin cậy khoa học cao.
  • Cung cấp cơ sở khoa học và dữ liệu số hóa trực quan phục vụ công tác điều chỉnh quy hoạch tổng thể phát triển kinh tế - xã hội huyện Kiến Thụy tầm nhìn 2020 theo hướng công nghiệp hóa và bảo vệ môi trường sinh thái.
  • Thiết lập tiền đề phương pháp luận vững chắc để mở rộng ứng dụng mô hình CA-Markov cho các quận, huyện ven biển khác thuộc vùng Duyên hải Bắc Bộ.

Về kế hoạch triển khai tiếp theo trong giai đoạn 2015 - 2020, địa phương cần nhanh chóng tích hợp bản đồ dự báo không gian vào hệ thống thông tin địa lý cấp huyện, định kỳ thu nhận và giải đoán ảnh vệ tinh mới để hiệu chỉnh các kịch bản phát triển đất đai. Các nhà quản lý và cơ quan chức năng hãy chủ động ứng dụng công nghệ viễn thám và mô hình không gian CA-Markov ngay hôm nay nhằm nâng cao hiệu quả quản trị tài nguyên đất đai và kiến tạo không gian sống bền vững cho cộng đồng.