Tổng quan nghiên cứu

Quá trình đô thị hóa và công nghiệp hóa diễn ra mạnh mẽ tại khu vực ven đô thành phố Hà Nội đã tạo ra sức ép chuyển dịch cơ cấu sử dụng đất vô cùng sâu sắc. Huyện Đông Anh là địa bàn cửa ngõ phía Bắc thủ đô với tổng diện tích tự nhiên đạt 18.230 ha, sở hữu tỷ lệ đất nông nghiệp chiếm tới 23% tổng quỹ đất ngoại thành, đứng thứ hai toàn thành phố sau huyện Sóc Sơn. Đây vốn là vành đai xanh cung cấp nguồn lương thực, thực phẩm và rau màu chủ lực cho toàn bộ khu vực nội đô. Tuy nhiên, số liệu thống kê địa chính cho thấy trong giai đoạn từ năm 1990 đến năm 2011, diện tích đất nông nghiệp của huyện đã bị suy giảm 991,87 ha, tương đương mức giảm bình quân 47,23 ha/năm. Đáng chú ý, trong giai đoạn 2000 đến 2010, huyện ghi nhận tới 39,11% diện tích đất canh tác bị chuyển đổi sang mục đích phi nông nghiệp để xây dựng hạ tầng, khu công nghiệp và khu đô thị mới.

Các phương pháp quản lý đất đai truyền thống dựa trên số liệu kiểm kê định kỳ 5 năm một lần thường bộc lộ độ trễ lớn và chưa phản ánh kịp thời cấu trúc không gian cảnh quan. Nhằm giải quyết bài toán này, đề tài tập trung đo đạc trắc lượng hình thái lớp phủ bằng công nghệ viễn thám đa thời gian trong giai đoạn 1993 đến 2011. Mục tiêu trọng tâm là lượng hóa biến động diện tích, phân tích mức độ phân mảnh, cô lập của các thửa đất canh tác thông qua phương pháp phân loại định hướng đối tượng và các chỉ số cảnh quan. Kết quả nghiên cứu không chỉ thiết lập bộ dữ liệu không gian có độ chính xác phân loại tổng thể đạt trên 91% mà còn cung cấp cơ sở khoa học định lượng phục vụ công tác quy hoạch sử dụng đất bền vững theo định hướng mở rộng đô thị trung tâm Hà Nội đến năm 2030.

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 trên nền tảng lý thuyết sinh thái cảnh quan và lý thuyết hình học fractal, kết hợp với các nguyên lý phản xạ phổ trong viễn thám quang học. Cấu trúc không gian cảnh quan được nhìn nhận dưới dạng mô hình khảm cảnh quan, trong đó các mảng đất nông nghiệp liên tục bị chia cắt bởi quá trình mở rộng hạ tầng nhân tạo. Bốn khái niệm cốt lõi được vận dụng xuyên suốt bao gồm:

  • Lớp phủ mặt đất: Trạng thái sinh học và vật lý hiện hữu trên bề mặt đất, được nhận biết thông qua đặc trưng phản xạ phổ tại các kênh sóng nhìn thấy, cận hồng ngoại và hồng ngoại sóng ngắn.
  • Đối tượng ảnh: Đơn vị không gian cơ bản được hình thành từ quá trình gom cụm các điểm ảnh lân cận có tính đồng nhất về phổ, kết cấu, hình dạng và mối liên hệ ngữ cảnh không gian.
  • Phân loại định hướng đối tượng: Phương pháp phân tích ảnh dựa trên tri thức chuyên gia kết hợp thuật toán phân đoạn đa tỷ lệ, khắc phục triệt để hiện tượng nhiễu muối tiêu của phương pháp phân loại theo điểm ảnh truyền thống.
  • Trắc lượng hình thái cảnh quan: Hệ thống các chỉ số định lượng đo đạc hình học, kích thước, mật độ, độ phức tạp đường biên, mức độ kết nối và độ cô lập của các mảng không gian.

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

Nguồn dữ liệu viễn thám đa thời gian được chuẩn hóa bao gồm ảnh vệ tinh Landsat TM chụp ngày 27 tháng 12 năm 1993, Landsat TM ngày 20 tháng 12 năm 1999, Landsat ETM+ ngày 09 tháng 10 năm 2005 (độ phân giải không gian 30m) và ảnh SPOT 5 ngày 22 tháng 12 năm 2011 (độ phân giải không gian 10m). Toàn bộ dữ liệu ảnh và bản đồ địa hình tỷ lệ 1:10.000, bản đồ hiện trạng sử dụng đất năm 2005 tỷ lệ 1:25.000 đều được nắn chỉnh hình học quy chiếu thống nhất về Hệ tọa độ Quốc gia VN-2000 với sai số trung phương nhỏ hơn 0,5 pixel.

Quy trình phân loại định hướng đối tượng được thực hiện trên phần mềm eCognition 8.64 với thuật toán phân đoạn đa độ phân giải. Thông số phân đoạn đối với ảnh SPOT 5 được tối ưu hóa với tham số tỷ lệ là 15, trọng số hình dạng 0,7 và độ chặt 0,3; đối với ảnh Landsat, tham số tỷ lệ được thiết lập ở mức 5, trọng số hình dạng 0,5 và độ chặt 0,2. Quá trình chiết tách thông tin tích hợp các chỉ số phổ chuyên biệt gồm NDVI, VI nhằm tách lớp thực vật, NDBI và UI để nhận diện đất xây dựng, cùng chỉ số LSWI để trích xuất mặt nước.

Công tác kiểm chứng độ chính xác trong phòng được tiến hành với cỡ mẫu 51 ô mẫu kiểm chứng được chọn theo phương pháp ngẫu nhiên phân tầng. Cỡ mẫu này được xác định dựa trên công thức thống kê với mức độ tin cậy 95% (hệ số Z = 1,96), sai số cho phép 10% và độ chính xác kỳ vọng 85%, chiếm quy mô tương đương 2% tổng diện tích tự nhiên của huyện. Kết quả ma trận sai số năm 2005 đạt độ chính xác tổng quát 91,6% và hệ số Kappa đạt 0,844, khẳng định độ tin cậy cao của phương pháp phân loại. Ngoài ra, nghiên cứu đã tiến hành 27 điểm kiểm chứng thực địa vào tháng 12 năm 2012 kết hợp điều tra phiếu nông lịch nông hộ nhằm xác định lịch canh tác 3 vụ/năm (2 vụ lúa và 1 vụ rau màu).

Để tính toán trắc lượng hình thái cảnh quan, lớp đất nông nghiệp sau khi chiết tách được đưa vào công cụ Patch Analysis 4 trên nền tảng ArcGIS 9.3 để xử lý 31 chỉ số hình thái. Lý do áp dụng phương pháp phân tích thành phần chính PCA là nhằm loại bỏ các biến số dư thừa, chọn lọc ra 6 chỉ số đại diện có tương quan cao nhất gồm: LPI (tỷ lệ mảng lớn nhất), AWMSI (độ phức tạp hình dạng), MPS (kích thước mảng trung bình), MPI (mức độ tiếp giáp), MNN (khoảng cách lân cận gần nhất) và TCA (tổng diện tích vùng lõi).

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 không gian và thống kê trắc lượng cho thấy sự chuyển đổi sâu sắc của đất nông nghiệp huyện Đông Anh trong 18 năm nghiên cứu:

  • Suy giảm mạnh diện tích canh tác truyền thống: Diện tích đất sản xuất nông nghiệp từ mức 9.989,02 ha năm 1995 đã giảm xuống còn 8.600,25 ha vào năm 2011, trung bình mỗi năm mất 86,80 ha. Trong đó, đất trồng lúa chịu tác động nặng nề nhất khi giảm từ 9.062,90 ha năm 1995 xuống còn 7.506,00 ha năm 2011, tương ứng tốc độ suy giảm 97,30 ha/năm. Đáng chú ý, giai đoạn 2006 đến 2011 ghi nhận diện tích lúa giảm tới 591,69 ha, đạt tốc độ suy giảm kỷ lục 118,34 ha/năm.
  • Tái cơ cấu cơ cấu nội bộ đất nông nghiệp: Trong khi đất lúa sụt giảm, diện tích đất trồng cây lâu năm tăng mạnh gấp 4,7 lần từ 43,04 ha năm 1995 lên 203,51 ha năm 2011. Diện tích mặt nước phục vụ nuôi trồng thủy sản cũng tăng gần gấp đôi từ 316,33 ha lên 613,34 ha trong cùng giai đoạn, phản ánh xu hướng chuyển dịch kinh tế nông nghiệp ven đô từ trồng trọt thuần túy sang các mô hình trang trại giá trị cao.
  • Phân mảnh cảnh quan và biến đổi hình thái nghiêm trọng: Chỉ số kích thước mảng trung bình MPS và tổng diện tích vùng lõi TCA sụt giảm liên tục, trong khi chỉ số hình thái có trọng số diện tích AWMSI tăng rõ rệt từ năm 1993 đến 2011. Điều này chứng minh các cánh đồng canh tác tập trung đã bị chia cắt thành các mảnh nhỏ với đường biên méo mó, phức tạp. Chỉ số khoảng cách lân cận trung bình MNN gia tăng đáng kể phản ánh mức độ cô lập ngày càng cao giữa các dải đất nông nghiệp còn lại.
  • Tình trạng bỏ hoang và suy thoái cục bộ: Toàn huyện ghi nhận 84,27 ha đất nông nghiệp canh tác kém hiệu quả hoặc bị bỏ hoang, nằm phân tán rải rác trên địa bàn 12 xã trọng điểm gồm Dục Tú, Liên Hà, Vân Hà, Thụy Lâm, Xuân Nộn, Việt Hùng, Uy Nỗ, Vân Nội, Xuân Canh, Nguyên Khê, Võng La và Hải Bối.

Thảo luận kết quả

Nguyên nhân căn bản dẫn đến sự suy giảm và manh mún đất nông nghiệp tại Đông Anh xuất phát từ việc triển khai hàng loạt dự án hạ tầng giao thông động lực và các khu công nghiệp tập trung. Việc mở rộng Quốc lộ 3, tuyến đường cao tốc Thăng Long - Nội Bài, xây dựng cầu Nhật Tân, cầu Đông Trù cùng dự án đường 5 kéo dài đã tạo nên các trục chia cắt không gian sinh thái nông nghiệp. Song song đó, sự hình thành Khu công nghiệp Bắc Thăng Long, các cụm công nghiệp vừa và nhỏ tại xã Nguyên Khê, cụm làng nghề Vân Hà và các khu đô thị mới tại xã Liên Hà, Kim Chung đã trực tiếp chuyển đổi hàng trăm hecta đất bờ xôi ruộng mật thành đất phi nông nghiệp.

Khi so sánh với các nghiên cứu đô thị hóa tại vùng đồng bằng Châu thổ sông Châu Giang của Trung Quốc hay vùng ven đô bang Arizona tại Hoa Kỳ, biến đổi hình thái đất đai ở Đông Anh thể hiện đặc trưng rõ nét của mô hình đô thị hóa dạng dải và nhảy cóc. Quá trình phát triển các khu dân cư tự phát ven các trục lộ làm gia tăng chỉ số mật độ cạnh biên và đẩy nhanh tốc độ phân mảnh đất nông nghiệp nhanh hơn so với mô hình đô thị hóa tập trung. Tình trạng 84,27 ha đất nông nghiệp bị bỏ hoang cũng phản ánh tâm lý người nông dân chờ đợi đền bù giải tỏa hoặc chuyển đổi sinh kế sang dịch vụ thương mại khi diện tích canh tác bình quân chỉ còn 0,051 ha/lao động nông nghiệp (theo thống kê năm 2009).

Các phát hiện này có thể được trực quan hóa tối ưu thông qua hệ thống bản đồ ma trận chuyển dịch lớp phủ giai đoạn 1993 - 2011, kết hợp với các biểu đồ cột thể hiện tốc độ suy giảm diện tích lúa theo từng giai đoạn 5 năm và biểu đồ radar so sánh 6 chỉ số hình thái không gian giữa 23 xã và 1 thị trấn. Cách trình bày này giúp các nhà quản lý nhận diện trực tiếp những điểm nóng phân mảnh không gian cần ưu tiên can thiệp.

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

Nhằm cân bằng giữa mục tiêu phát triển đô thị công nghiệp và bảo vệ không gian sinh thái nông nghiệp bền vững, nghiên cứu đề xuất 4 nhóm giải pháp chiến lược:

  • Thiết lập ranh giới bảo tồn và vành đai nông nghiệp sinh thái: Ủy ban nhân dân huyện Đông Anh cần phối hợp với Sở Quy hoạch - Kiến trúc Hà Nội khoanh vùng bảo vệ nghiêm ngặt tối thiểu 7.000 ha đất chuyên canh lúa chất lượng cao và vùng sản xuất rau màu an toàn tại các xã có điều kiện đất đai màu mỡ như Xuân Nộn, Thụy Lâm, Vân Nội. Cần áp dụng chỉ tiêu mật độ mảng tối thiểu nhằm ngăn chặn việc chia cắt cơ học các cánh đồng lớn trước năm 2025.
  • Ứng dụng công nghệ viễn thám và GIS trong giám sát biến động định kỳ: Phòng Tài nguyên và Môi trường huyện Đông Anh cần thiết lập quy trình quan trắc lớp phủ tự động sử dụng ảnh vệ tinh SPOT, Sentinel định kỳ 6 tháng một lần. Công cụ này cho phép phát hiện sớm các vi phạm lấn chiếm đất nông nghiệp, theo dõi chính xác biến động chuyển đổi mục đích sử dụng đất và kiểm soát triệt để diện tích 84,27 ha đất nông nghiệp đang bị bỏ hoang để có biện pháp phục hồi sản xuất kịp thời trong vòng 12 đến 24 tháng.
  • Chuyển đổi mô hình canh tác sang nông nghiệp đô thị công nghệ cao: Khuyến khích các hợp tác xã nông nghiệp chuyển đổi ít nhất 30% diện tích đất trồng lúa kém hiệu quả sang mô hình trồng hoa, cây cảnh, rau an toàn ứng dụng nhà màng và nuôi trồng thủy sản sinh thái. Mục tiêu nâng cao giá trị sản xuất trên một hecta đất canh tác từ 150 triệu đồng/ha lên 350 - 500 triệu đồng/ha trước năm 2030, giải quyết bài toán việc làm cho hơn 108.000 lao động nông nghiệp địa phương.
  • Hoàn thiện khung thể chế và quy hoạch tích hợp không gian: Ủy ban nhân dân thành phố Hà Nội cần ban hành cơ chế bồi thường, tích tụ ruộng đất linh hoạt, đồng thời tích hợp các chỉ số trắc lượng hình thái cảnh quan vào quy chuẩn lập quy hoạch sử dụng đất cấp huyện. Điều này giúp kiểm soát ngưỡng phân mảnh không gian của đất công cộng và đất cây xanh theo Đồ án Quy hoạch chung xây dựng Thủ đô Hà Nội đến năm 2030 và tầm nhìn đến năm 2050.

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

Nghiên cứu mang giá trị thực tiễn và học thuật cao, là tài liệu tham khảo hữu ích cho các nhóm đối tượng:

  • Cơ quan quản lý nhà nước về tài nguyên và quy hoạch: Sở Tài nguyên và Môi trường, Sở Nông nghiệp và Phát triển nông thôn Hà Nội, Phòng Tài nguyên và Môi trường huyện Đông Anh có thể sử dụng trực tiếp kết quả phân loại và ma trận chuyển đổi lớp phủ làm dữ liệu đầu vào cho công tác lập quy hoạch, kế hoạch sử dụng đất giai đoạn 2021 - 2030.
  • Viện nghiên cứu, Giảng viên và Học viên cao học: Các nhà khoa học thuộc chuyên ngành Bản đồ, Viễn thám, Hệ thông tin địa lý (GIS), Quản lý đất đai và Khoa học môi trường có thể ứng dụng quy trình công nghệ phân loại định hướng đối tượng eCognition và phương pháp phân tích chỉ số FRAGSTATS vào các nghiên cứu biến đổi cảnh quan tại các vùng sinh thái khác.
  • Đơn vị tư vấn quy hoạch đô thị và kiến trúc cảnh quan: Các chuyên gia quy hoạch sử dụng phương pháp trắc lượng hình thái để đánh giá mức độ gắn kết không gian, từ đó thiết kế các hành lang xanh, đệm sinh thái giữa các khu đô thị mới và vùng nông nghiệp truyền thống.
  • Các tổ chức phát triển nông nghiệp và dự án quốc tế: Các ban quản lý dự án về thích ứng biến đổi khí hậu, sinh kế cộng đồng và an ninh lương thực đô thị (như các dự án của Danida, Nafosted) có thể tham khảo phương pháp tiếp cận không gian định lượng để đánh giá tính tổn thương của tài nguyên đất đai.

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

Phương pháp phân loại định hướng đối tượng (OBIA) có ưu thế gì vượt trội so với phương pháp phân loại theo pixel truyền thống? Phương pháp định hướng đối tượng giải quyết triệt để hạn chế của phân loại theo điểm ảnh đối với ảnh độ phân giải cao nhờ việc gộp các pixel có cùng đặc tính thành các mảng đối tượng. Quá trình này tích hợp không chỉ giá trị phản xạ phổ mà cả hình dạng, kết cấu và ngữ cảnh không gian, giúp phân biệt rõ ràng giữa đất trống ven sông và đất xây dựng, đạt độ chính xác tổng quát lên tới 91,6% và hệ số Kappa 0,844.

Những chỉ số trắc lượng hình thái nào phản ánh rõ nét nhất tình trạng phân mảnh đất nông nghiệp tại Đông Anh? Thông qua phân tích thành phần chính PCA từ 31 biến số, 6 chỉ số nhạy cảm nhất được xác định là LPI (tỷ lệ mảng lớn nhất), AWMSI (độ phức tạp hình học), MPS (kích thước mảng trung bình), MPI (mức độ tiếp giáp), MNN (khoảng cách lân cận) và TCA (diện tích lõi). Sự sụt giảm của MPS, TCA kết hợp với sự gia tăng của AWMSI và MNN phản ánh trực tiếp tình trạng các cánh đồng bị chia nhỏ và cô lập.

Giai đoạn nào ghi nhận tốc độ suy giảm diện tích đất nông nghiệp cao nhất tại huyện Đông Anh? Giai đoạn 2007 đến 2011 chứng kiến tốc độ mất đất nông nghiệp nhanh nhất với 312,66 ha bị chuyển đổi mục đích, cao gấp hơn 20 lần so với giai đoạn 2000 - 2005. Riêng diện tích đất lúa giai đoạn 2006 - 2011 đã sụt giảm 591,69 ha với tốc độ bình quân lên tới 118,34 ha/năm do hàng loạt dự án hạ tầng lớn như cầu Nhật Tân và các khu công nghiệp được đẩy mạnh giải phóng mặt bằng.

Nguyên nhân chính dẫn đến hiện tượng hơn 84 ha đất nông nghiệp bị bỏ hoang tại huyện Đông Anh là gì? Tình trạng 84,27 ha đất nông nghiệp canh tác kém hiệu quả tại 12 xã bắt nguồn từ việc các thửa đất nằm xen kẹt giữa các công trình xây dựng, bị cắt đứt hệ thống thủy lợi tưới tiêu và ô nhiễm cục bộ. Thêm vào đó, diện tích canh tác quá nhỏ (chỉ khoảng 0,051 ha/lao động) mang lại hiệu quả kinh tế thấp khiến người dân có xu hướng bỏ ruộng để chuyển sang làm công nhân hoặc kinh doanh dịch vụ.

Làm thế nào để đảm bảo tính đồng nhất khi sử dụng nhiều loại ảnh vệ tinh có độ phân giải khác nhau (Landsat và SPOT 5)? Nghiên cứu áp dụng quy trình nắn chỉnh hình học chuẩn xác theo bản đồ địa hình 1:10.000 về hệ tọa độ VN-2000 với sai số nhỏ hơn 0,5 pixel. Quá trình phân đoạn ảnh trên eCognition được điều chỉnh tham số tỷ lệ linh hoạt phù hợp với từng độ phân giải (tham số 15 cho SPOT 5 10m và tham số 5 cho Landsat 30m), kết hợp chuẩn hóa các kênh chỉ số thực vật và thổ nhưỡng để đảm bảo khả năng so sánh chuỗi thời gian liên tục.

Kết luận

  • Nghiên cứu đã xây dựng thành công quy trình công nghệ kết hợp phân loại định hướng đối tượng (OBIA) và phân tích trắc lượng hình thái cảnh quan trên chuỗi ảnh vệ tinh Landsat TM, ETM+ và SPOT 5 giai đoạn 1993 - 2011.
  • Kết quả kiểm chứng đạt độ chính xác phân loại tổng thể 91,6% và hệ số Kappa 0,844, chứng minh tính ưu việt của phương pháp đối với các cảnh quan ven đô phức tạp và manh mún.
  • Lượng hóa chi tiết sự suy giảm của 991,87 ha đất nông nghiệp trong 21 năm, chỉ ra sự chuyển dịch cơ cấu nội bộ với xu hướng giảm mạnh đất lúa (97,30 ha/năm) và gia tăng đất cây lâu năm cùng diện tích nuôi trồng thủy sản.
  • Xác định quy luật biến đổi cấu trúc cảnh quan thông qua 6 chỉ số cốt lõi (LPI, AWMSI, MPS, MPI, MNN, TCA), làm rõ tác động tiêu cực của các dự án công nghiệp và hạ tầng giao thông gây phân mảnh và cô lập không gian canh tác.
  • Định hình khung giải pháp khoa học gồm bảo vệ 7.000 ha vùng sản xuất trọng điểm, ứng dụng GIS giám sát định kỳ và chuyển đổi nông nghiệp công nghệ cao, đóng góp trực tiếp vào mục tiêu quy hoạch thủ đô bền vững.

Các nhà quản lý, nhà quy hoạch và các nhóm nghiên cứu quan tâm có thể tiếp tục mở rộng phương pháp này kết hợp với các nguồn dữ liệu viễn thám độ phân giải siêu cao và mô hình mô phỏng không gian để dự báo biến động sử dụng đất cho các giai đoạn tiếp theo.