Tổng quan nghiên cứu

Đất nông nghiệp canh tác nhờ nước trời hiện chiếm khoảng 80% tổng diện tích đất trồng trọt trên toàn cầu và cung cấp tới 70% sản lượng lương thực cho nhân loại. Trong bối cảnh dân số thế giới được dự báo sẽ tăng từ 7,3 tỷ người lên khoảng 9,7 tỷ người vào năm 2050 và chạm ngưỡng 11,2 tỷ người vào năm 2100, áp lực duy trì an ninh lương thực cùng với biến đổi khí hậu đang đe dọa nghiêm trọng tính bền vững của tài nguyên đất và nước. Khả năng thẩm thấu nước của đất (infiltrability - ký hiệu là Ks) là đại lượng thủy văn then chốt điều tiết tỷ lệ phân chia giữa lượng nước ngấm vào tầng rễ cây và lượng nước chảy tràn trên mặt đất. Khi tốc độ thấm bị suy giảm, nguy cơ xói mòn rửa trôi bề mặt tăng cao, gây thất thoát dinh dưỡng và sụt giảm năng suất cây trồng.

Mục tiêu cốt lõi của nghiên cứu là xây dựng quy trình công nghệ bản đồ hóa và định lượng các biến trạng thái bề mặt đất nông nghiệp từ chuỗi ảnh vệ tinh quang học Sentinel-2, phục vụ việc xác định chính xác cấp độ thẩm thấu nước ở quy mô từng thửa ruộng. Nghiên cứu được triển khai thực nghiệm tại lưu vực Kamech thuộc bán đảo Cap Bon, Tunisia với diện tích 2,63 km² bao gồm khoảng 390 thửa đất canh tác đặc trưng cho vùng Địa Trung Hải. Dữ liệu viễn thám và thực địa được thu thập đồng bộ qua 5 thời điểm quan trắc liên tiếp từ tháng 8 đến tháng 12 năm 2016. Kết quả nghiên cứu đã phân loại thành công 337 trên tổng số 384 thửa đất khảo sát, mang lại công cụ giám sát bề mặt đất có độ lặp lại cao với chu kỳ 5 ngày, thay thế hiệu quả các đợt khảo sát thủ công tốn kém và hỗ trợ đắc lực cho công tác thủy lợi, bảo tồn đất đai.

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 ba trụ cột lý thuyết chuyên sâu:

Mô hình thủy văn bề mặt xác định khả năng thẩm thấu nước ở trạng thái bão hòa (Ks) thông qua hàm tương quan đa biến: Ks = f(RT, Vgt, L, Cx, Faciès). Trong đó, RT biểu thị độ gồ ghề bề mặt (rugosité), Vgt là tỷ lệ che phủ của thực vật tươi (thực vật sống), L đại diện cho thảm tàn dư thực vật khô hoặc rơm rạ (litière), Cx phản ánh tỷ lệ sỏi đá thô trên bề mặt, và Faciès mô tả mức độ đóng mở cơ học của vỏ đất. Các nghiên cứu thổ nhưỡng thực nghiệm khẳng định trạng thái bề mặt phản ánh trực tiếp cấu trúc xốp và độ rỗng vĩ mô của đất canh tác.

Lý thuyết quang phổ viễn thám (Spectroscopy) giải thích sự tương tác giữa bức xạ điện từ mặt trời với các đối tượng thổ nhưỡng và thảm thực vật trong dải bước sóng từ 400 nm đến 2500 nm. Độ phản xạ phổ của đất biến thiên phụ thuộc vào thành phần khoáng vật (như ion carbonat hấp thụ mạnh ở bước sóng khoảng 2340 nm), hàm lượng chất hữu cơ, độ ẩm và độ mịn hạt đất. Ngược lại, thực vật tươi hấp thụ mạnh bức xạ đỏ và phản xạ cực đại ở vùng cận hồng ngoại (NIR), trong khi tàn dư khô phản xạ rõ nét ở các dải hồng ngoại sóng ngắn (SWIR).

Hệ phương pháp học máy thống kê bao gồm các mô hình phân loại có giám sát (Supervised Learning) như Rừng ngẫu nhiên (Random Forest) với cơ chế lấy mẫu lặp lại kết hợp bỏ phiếu đa số và Máy vector hỗ trợ (Support Vector Machine - SVM) sử dụng hàm nhân xuyên tâm (RBF kernel) để tối ưu hóa siêu phẳng phân tách; kết hợp các mô hình phân loại không giám sát (Unsupervised Learning) như K-means và Phân cụm thứ bậc tăng dần (CAH) dựa trên khoảng cách khoảng cách Euclid và tiêu chuẩn Ward sau khi giảm chiều dữ liệu bằng Phân tích thành phần chính (PCA).

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

Nghiên cứu sử dụng nguồn dữ liệu kép gồm ảnh vệ tinh và khảo sát thực địa:

Dữ liệu viễn thám bao gồm 5 cảnh ảnh vệ tinh Sentinel-2 cấp độ 2A (phản xạ đáy khí quyển - BoA) chụp vào các ngày 04/08/2016, 03/10/2016, 02/11/2016, 22/11/2016 và 02/12/2016. Mỗi ảnh trích xuất 10 dải phổ hữu ích với độ phân giải không gian từ 10m đến 20m, bao quát toàn bộ 15.127 pixel trên lưu vực nghiên cứu.

Dữ liệu kiểm chứng mặt đất gồm 34 thửa đất chuẩn được quan trắc định kỳ hàng tháng từ năm 2006 trong khuôn khổ đài quan sát ORE OMERE, tương ứng với tập dữ liệu 1.264 pixel thực địa. Diện tích trung bình mỗi thửa đạt 5.129 m² (dao động từ 654 m² đến 12.272 m²). Các chuyên gia tiến hành đánh giá 5 biến số trạng thái: Thực vật tươi (6 cấp độ từ 0% đến 100%), Thực vật khô (4 cấp độ từ 5% đến 100%), Hoạt động làm đất (3 lớp: cày mới TR, cày cũ TA, không cày NT), Độ gồ ghề (3 lớp: RT0 dưới 2cm, RT2 từ 2-5cm, RT5 trên 5cm) và Cấu trúc Faciès (6 mức độ từ F0 mở hoàn toàn đến F2 đóng kín vón cục).

Quy trình chọn mẫu và phân tách dữ liệu được thực hiện theo nguyên tắc phân tầng ngẫu nhiên: Với cách tiếp cận theo pixel, 70% số pixel được đưa vào tập huấn luyện và 30% dành cho tập kiểm thử độc lập. Với cách tiếp cận theo thửa, 2/3 số thửa được dùng để huấn luyện và 1/3 số thửa dùng để kiểm thử. Kỹ thuật kiểm định chéo 10 lần (10-fold cross-validation với 90% dữ liệu huấn luyện và 10% dữ liệu thẩm định nội bộ) được áp dụng để triệt tiêu hiện tượng quá khớp (overfitting). Các chỉ số đánh giá độ chính xác bao gồm Độ chính xác tổng thể (Overall Accuracy), Hệ số Cohen's Kappa và Chỉ số Rand hiệu chỉnh (Adjusted Rand Index - ARI).

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

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

Mô hình học máy giám sát đạt hiệu suất phân loại vượt trội so với phân cụm tự động. Khi áp dụng thuật toán Random Forest theo cách tiếp cận đơn thời điểm (mono-date) và đa thời điểm không thứ tự (multi-date non-ordonnée) ở cấp độ pixel, mô hình đã nhận diện chính xác 337 trên tổng số 384 thửa đất nông nghiệp toàn lưu vực, đạt tỷ lệ tin cậy xấp xỉ 87,8%.

Phân loại thảm phủ thực vật đạt độ khớp cao nhất trong các biến khảo sát. Thảm thực vật tươi (Vgtv) và thực vật khô (Vgts) được phân biệt rõ rệt trên ảnh Sentinel-2 nhờ sự tương phản phổ mạnh mẽ giữa kênh cận hồng ngoại B8, B8a và các kênh khả kiến. Ngược lại, việc phân loại độ gồ ghề (RT) và cấu trúc bề mặt (Faciès) đòi hỏi phải kết hợp các kênh hồng ngoại sóng ngắn SWIR (B11 và B12) mới đạt độ chính xác trên 75%.

Phân tích mức độ suy giảm độ chính xác trung bình (Mean Decrease in Accuracy - MDA) từ thuật toán Random Forest chỉ ra rằng các dải phổ B11, B12 và các dải biên đỏ Red Edge (B5, B6, B7) giữ vai trò quan trọng nhất trong việc bóc tách trạng thái cơ học của bề mặt đất, vượt qua tầm ảnh hưởng của các kênh màu quang học truyền thống.

Thuật toán không giám sát K-means và CAH (sau khi giảm chiều qua 3 trục PCA giữ lại 95% phương sai) chỉ đạt chỉ số tương đồng Adjusted Rand Index (ARI) khoảng 50% đối với thực vật và 40% đối với trạng thái làm đất. Lý do là sự pha trộn phổ tự nhiên giữa đất ẩm, tàn dư hữu cơ và sỏi đá trên cùng một pixel khiến các cụm phổ không phản ánh đơn nhất một biến số thủy văn độc lập.

Thảo luận kết quả

Sự biến động theo mùa của thảm phủ thực vật giải thích rõ nét độ chính xác của các mô hình viễn thám. Vào tháng 8 và tháng 9, thời điểm nông dân chuẩn bị đất sau thu hoạch, độ che phủ thực vật tươi rất thấp (chủ yếu thuộc lớp dưới 5%), trong khi tàn dư thực vật khô chiếm ưu thế với phổ phản xạ biến thiên mạnh. Đến tháng 11 và tháng 12, khi cây trồng vụ đông bắt đầu nảy mầm và phát triển, phổ phản xạ diệp lục tăng vọt giúp thuật toán phân tách cực kỳ chuẩn xác các lớp che phủ từ 25% đến trên 75%.

So sánh giữa hai cấp độ không gian cho thấy phân loại cấp độ pixel (pixel-level) phản ánh trung thực sự không đồng nhất vi mô trong thửa ruộng, bắt trọn các vệt cày và ổ trũng tích tụ rơm rạ. Trong khi đó, phân loại cấp độ thửa (parcel-level) gộp các pixel theo giá trị trung bình hoặc đa số, tuy làm mượt bản đồ và loại bỏ nhiễu hạt tiêu nhưng lại làm giảm nhẹ độ chính xác cục bộ khoảng 3% đến 5% ở những thửa có diện tích lớn trên 10.000 m².

Để trực quan hóa hiệu quả phân loại, dữ liệu viễn thám được tổng hợp qua ma trận nhầm lẫn (confusion matrix), làm nổi bật sự tách biệt giữa các lớp trạng thái đối lập (như cày xới F0 và đóng vỏ F2). Bản đồ sai khác không gian (difference map) giữa phân loại đơn thời điểm và đa thời điểm cho thấy chuỗi thời gian liên tục giúp hạn chế đáng kể các sai số do phản xạ góc chiếu mặt trời thay đổi giữa các tháng. Biểu đồ phân bố diện tích thửa đất chỉ ra rằng các thửa có quy mô từ 3.000 m² đến 8.000 m² đạt độ chính xác dự báo ổn định nhất.

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

Tích hợp dữ liệu Radar khẩu độ tổng hợp Sentinel-1: Khuyến nghị các cơ quan nghiên cứu viễn thám phối hợp chuỗi ảnh quang học Sentinel-2 với dữ liệu phân cực Radar SAR Sentinel-1 băng C nhằm cải thiện độ chính xác phân loại độ gồ ghề bề mặt (RT) và độ ẩm lớp đất mặt lên trên 92%, hoàn thành tích hợp thử nghiệm trong vòng 12 tháng tới.

Xây dựng hệ thống tự động hóa bản đồ thấm nước theo thời gian thực: Đề xuất ban quản lý lưu vực sông và cơ quan thủy lợi thiết lập quy trình tính toán tự động chỉ số thấm Ks theo chu kỳ cập nhật 5 ngày của vệ tinh, giúp cảnh báo sớm nguy cơ lũ quét và xói mòn đất trước các đợt mưa lớn với thời gian phản hồi dưới 24 giờ.

Bổ sung các chỉ số phổ chuyên biệt cho tàn dư thực vật: Nhóm phát triển thuật toán cần đưa thêm các chỉ số viễn thám nâng cao như NDTI (Normalized Difference Tillage Index), STI (Soil Tillage Index) kết hợp cùng toàn bộ phổ phản xạ nguyên bản (full spectrum) nhằm nâng cao độ nhạy nhận diện lớp thực vật khô (Vgts) thêm 15% đến 20% trong giai đoạn mùa khô.

Mở rộng phạm vi kiểm chứng sang các hệ sinh thái nông nghiệp nhiệt đới: Viện nghiên cứu nông học cần nhân rộng mô hình phân loại từ lưu vực thử nghiệm Kamech (2,63 km²) sang các vùng đồng bằng thâm canh lúa và cây công nghiệp có diện tích trên 100 km² tại Đông Nam Á trước quý 4 năm 2028, hiệu chỉnh lại các ngưỡng phân lớp Faciès để tương thích với kỹ thuật canh tác ngập nước.

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

Nhà nghiên cứu viễn thám và khoa học không gian: Tiếp cận quy trình trích xuất thông tin quang phổ đa thời gian, phương pháp tối ưu hóa kênh ảnh Red Edge/SWIR và kỹ thuật so sánh đối chuẩn hiệu năng giữa các thuật toán học máy kinh điển.

Kỹ sư nông nghiệp và chuyên gia quản lý tài nguyên đất: Ứng dụng mô hình đánh giá biến trạng thái mặt đất để chẩn đoán sức khỏe thổ nhưỡng, theo dõi tàn dư hữu cơ sau thu hoạch và định hướng các biện pháp canh tác bảo tồn đất giảm thiểu xói mòn.

Cán bộ quy hoạch thủy lợi và phòng chống thiên tai: Khai thác bản đồ phân vùng khả năng thẩm thấu nước cấp độ thửa ruộng để tính toán chính xác lưu lượng dòng chảy mặt, phục vụ thiết kế kênh mương và hồ điều hòa giảm thiểu ngập lụt cục bộ.

Học viên cao học và sinh viên chuyên ngành Khoa học dữ liệu / Tin học môi trường: Tham khảo cấu trúc xử lý dữ liệu không gian hoàn chỉnh trên môi trường lập trình R (sử dụng các gói caret, FactoMineR, randomForest) kết hợp phần mềm mã nguồn mở QGIS để thực hiện các đồ án phân loại chuỗi thời gian thực tế.

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

Dữ liệu ảnh vệ tinh Sentinel-2 có những ưu thế kỹ thuật nào trong việc theo dõi trạng thái bề mặt đất? Sentinel-2 sở hữu 13 kênh phổ với độ phân giải không gian cao từ 10m đến 20m và tần suất lặp lại 5 ngày. Đặc biệt, sự hiện diện của 3 kênh Red Edge và 2 kênh SWIR giúp phát hiện chính xác hàm lượng diệp lục, độ ẩm và tàn dư hữu cơ trên bề mặt đất canh tác.

Tại sao phân loại cấp độ pixel lại đạt hiệu quả phân biệt vi địa hình tốt hơn cấp độ thửa ruộng? Phân loại theo pixel giữ nguyên giá trị phản xạ quang phổ độc lập của từng ô lưới 10x10m, giúp nhận diện rõ sự phân hóa không đồng đều giữa các vệt cày xới và vùng tích tụ tàn dư thực vật. Trong khi đó, cấp độ thửa gộp chung dữ liệu làm mất đi các chi tiết vi mô cục bộ này.

Vì sao thuật toán Random Forest lại vượt trội hơn hẳn mô hình SVM và các phương pháp phân cụm K-means? Random Forest xử lý xuất sắc các tập dữ liệu có biến số tương quan cao, chống nhiễu tốt nhờ cơ chế tập hợp nhiều cây quyết định ngẫu nhiên và cung cấp chỉ số đo lường tầm quan trọng của từng dải phổ (MDA), giúp mô hình đạt độ chính xác phân loại thành công 337/384 thửa đất.

Năm biến trạng thái bề mặt trong nghiên cứu ảnh hưởng như thế nào đến khả năng thẩm thấu nước của đất? Thực vật tươi và tàn dư khô giúp cản lực hạt mưa rơi, giảm nén chặt đất và giữ ẩm; hoạt động làm đất cùng độ gồ ghề tạo các rãnh giữ nước tạm thời; trong khi cấu trúc Faciès quyết định độ rỗng bề mặt cho phép nước mưa ngấm sâu vào đất hay chảy tràn.

Phương pháp phân tích đa thời gian có thứ tự (multi-date ordonnée) đóng vai trò gì trong giám sát nông nghiệp? Phương pháp này phân tích đạo hàm biến thiên trạng thái giữa hai thời điểm liên tiếp (t1 và t2), giúp xác định xu hướng tăng trưởng, suy giảm hoặc ổn định của thảm thực vật và độ gồ ghề, từ đó phát hiện nhanh các hoạt động gieo cấy hoặc thu hoạch của nông dân.

Kết luận

  • Luận văn đã thiết lập thành công khung phương pháp luận kết hợp chuỗi ảnh quang học Sentinel-2 với các mô hình học máy để giải mã 5 biến số trạng thái quyết định khả năng thấm của đất.
  • Thuật toán Random Forest chứng minh tính ưu việt vượt trội khi phân loại chính xác 337 trên tổng số 384 thửa đất nông nghiệp trên lưu vực thực nghiệm 2,63 km².
  • Nghiên cứu khẳng định vai trò cốt lõi của các kênh phổ hồng ngoại sóng ngắn SWIR và Red Edge trong việc định lượng độ gồ ghề và cấu trúc vỏ đất thổ nhưỡng.
  • Bốn quy trình phân tách dữ liệu đa thời gian (mono-date, multi-date non-ordonnée, semi-ordonnée và ordonnée) cung cấp bộ giải pháp linh hoạt cho nhiều bài toán giám sát mùa vụ khác nhau.
  • Kết quả mở ra triển vọng to lớn trong việc thay thế các khảo sát thực địa tốn kém bằng hệ thống viễn thám tự động hóa, hỗ trợ quản lý tài nguyên nông nghiệp bền vững.

Quý độc giả, các nhà nghiên cứu và kỹ sư nông nghiệp quan tâm đến phương pháp bản đồ hóa trạng thái đất và mô hình hóa thủy văn viễn thám hãy tiếp tục theo dõi, chia sẻ và ứng dụng bộ giải pháp này vào các đề án quản lý lưu vực thực tế!