Tổng quan nghiên cứu

Biến đổi khí hậu toàn cầu đang diễn biến phức tạp, trong đó khí mê-tan ($CH_4$) là một trong những tác nhân chính khi đóng góp khoảng 18% vào hiệu ứng nhà kính và sở hữu tiềm năng làm nóng toàn cầu cao gấp 25 lần so với $CO_2$ trong chu kỳ 100 năm. Theo báo cáo kiểm kê khí nhà kính quốc gia, ngành nông nghiệp Việt Nam năm 2010 phát thải khoảng 88,6 triệu tấn $CO_2$ tương đương (chiếm khoảng 33,2% tổng lượng phát thải không tính lĩnh vực sử dụng đất và lâm nghiệp). Trong đó, hoạt động canh tác lúa nước chiếm hơn 57,5% tổng lượng phát thải khí nhà kính toàn ngành nông nghiệp do quá trình phân hủy yếm khí của chất hữu cơ trong môi trường ngập nước kéo dài.

Đồng bằng sông Hồng là vựa lúa lớn thứ hai cả nước với hơn 700.000 ha đất nông nghiệp chuyên trồng lúa và mức độ thâm canh cao nhất Việt Nam. Tuy nhiên, việc tính toán phát thải khí nhà kính nông nghiệp trước đây chủ yếu dựa trên các hệ số phát thải trung bình tĩnh của IPCC (Bậc 1), chưa phản ánh được sự biến động không gian của các tiểu vùng đất đai, khí hậu, tập quán bón phân và chế độ tưới tiêu.

Mục tiêu trọng tâm của nghiên cứu là xây dựng phương pháp tích hợp tư liệu viễn thám MODIS đa thời gian với mô hình sinh địa hóa DNDC (DeNitrification-DeComposition) nhằm chiết tách không gian phân bố mùa vụ lúa và tính toán định lượng phát thải $CH_4$ tại 10 tỉnh thuộc Đồng bằng sông Hồng trong năm 2010. Kết quả nghiên cứu mang ý nghĩa thực tiễn to lớn trong việc số hóa kiểm kê khí nhà kính, phục vụ hoạch định chính sách giảm nhẹ phát thải quốc gia và tối ưu hóa chi phí sản xuất lúa 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 trên nền tảng lý thuyết sinh địa hóa đất ngập nước và nguyên lý viễn thám quang học đa thời gian:

  • Lý thuyết sinh địa hóa và mô hình DNDC: Quá trình sinh khí $CH_4$ diễn ra trong điều kiện yếm khí cực hạn khi thế oxy hóa - khử ($Eh$) của đất giảm xuống dưới mức -150 mV đến -200 mV. Mô hình DNDC mô phỏng chi tiết các quá trình tương tác lý - hóa - sinh trong đất gồm: chu trình các-bon, chu trình ni-tơ, sự phân hủy chất hữu cơ, sự phát triển của hệ vi sinh vật sinh mê-tan (Methanogens) và vi sinh vật oxy hóa mê-tan (Methanotrophs).
  • Cơ chế vận chuyển khí $CH_4$ trên ruộng lúa: Khoảng 90% lượng $CH_4$ thoát vào khí quyển thông qua hệ thống mô khí (aerenchyma) của thân cây lúa, 10% phát tán qua các bọt khí nổi lên mặt nước và chỉ khoảng 1% khuếch tán trực tiếp qua lớp nước bề mặt.
  • Nguyên lý nhận dạng mùa vụ lúa bằng chỉ số thực vật: Cây lúa nước có đặc tính quang phổ độc nhất khi trải qua thời kỳ làm đất và cấy mạ, bề mặt ruộng là hỗn hợp nước và mạ non. Thuật toán phân loại mùa vụ của Xiangming Xiao dựa trên việc phân tích sự giao cắt giữa chỉ số nước bề mặt (LSWI) và các chỉ số thực vật (NDVI, EVI) để nhận diện chính xác thời điểm sạ cấy và thu hoạch.

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

  • Nguồn dữ liệu đầu vào:
    • Dữ liệu viễn thám: Sử dụng chuỗi 46 cảnh ảnh MODIS MOD09A1 (phản xạ bề mặt tổ hợp 8 ngày, phân giải 500m, mảnh h27v06) bao phủ trọn vẹn năm 2010 trên toàn bộ Đồng bằng sông Hồng.
    • Dữ liệu thổ nhưỡng: Bản đồ đất tỷ lệ 1:250.000 của Viện Thổ nhưỡng Nông hóa gồm 35 tiểu đơn vị đất với các thuộc tính lý hóa chi tiết: hàm lượng các-bon hữu cơ (SOC từ 1,19% đến 8,99%), độ pH (4,0 đến 8,5), dung trọng (1,00 đến 1,62 g/cm³) và thành phần cấp hạt sét (5% đến 78%).
    • Dữ liệu khí tượng: Dữ liệu nhiệt độ tối cao, nhiệt độ tối thấp và lượng mưa theo ngày tại 10 trạm khí tượng đại diện (Láng, Phủ Liễn, Hải Dương, Nam Định, Thái Bình, Vĩnh Yên,...).
    • Dữ liệu nông học: Khảo sát nông lịch, lượng phân bón (N dao động từ 90,6 đến 131,1 kg/ha) và chế độ tưới tiêu tại 10 tỉnh thành.
  • Quy mô và phương pháp chọn mẫu: Nghiên cứu bao quát toàn bộ 10 tỉnh Đồng bằng sông Hồng. Việc kiểm chứng thực địa được thực hiện thông qua phương pháp lấy mẫu ngẫu nhiên phân tầng theo các tuyến khảo sát cắt ngang qua các tiểu vùng sinh thái tại Hà Nội, Hải Dương, Thái Bình và Nam Định.
  • Lý do lựa chọn phương pháp: Sự kết hợp giữa MODIS và GIS cho phép khắc phục hạn chế mây mù của ảnh quang học nhiệt đới nhờ chu kỳ lặp 8 ngày, cung cấp bản đồ phân bố diện tích lúa theo từng mùa vụ để tích hợp làm đầu vào không gian cho mô hình DNDC quy mô vùng.

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

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

  • Độ chính xác cao trong phân lập không gian mùa vụ lúa: Việc ứng dụng thuật toán chỉ số kết hợp LSWI, NDVI và EVI đã tách chiết thành công bản đồ diện tích lúa vụ Đông Xuân và vụ Hè Thu năm 2010 trên toàn bộ 10 tỉnh Đồng bằng sông Hồng. Kết quả đối chiếu với số liệu thống kê nông nghiệp cho thấy mức độ tương thích trên 85%, hạn chế tối đa sai số do mây mù che phủ.
  • Sự biến thiên phát thải $CH_4$ rõ rệt giữa hai mùa vụ: Lượng phát thải khí $CH_4$ trong vụ Hè Thu (vụ Mùa) cao hơn đáng kể so với vụ Đông Xuân. Cụ thể, cường độ phát thải trung bình vụ Hè Thu đạt khoảng 150 - 280 kg $CH_4$/ha/vụ, trong khi vụ Đông Xuân dao động từ 80 - 160 kg $CH_4$/ha/vụ (tăng trung bình khoảng 40% đến 60% tùy thuộc vào từng khu vực).
  • Phân hóa không gian phát thải theo đơn vị hành chính và loại đất: Các tỉnh có mật độ phát thải và tổng lượng $CH_4$ cao nhất là Thái Bình, Nam Định và Hải Dương. Nguyên nhân do diện tích thâm canh lúa 2 vụ lớn kết hợp với lượng bón đạm cao (vụ Xuân tại Thái Bình đạt 131,1 kg N/ha, vụ Mùa đạt 116,2 kg N/ha). Những vùng đất phù sa không được bồi đắp hàng năm có hàm lượng sét từ 35% trở lên và hàm lượng các-bon hữu cơ tầng mặt vượt mức 2,5% luôn ghi nhận lượng phát thải cao nhất.
  • Kiểm chứng mô hình bằng đo đạc thực nghiệm: Kết quả mô phỏng từ mô hình DNDC được so sánh với số liệu đo đạc thực tế tại trạm thí nghiệm Hoài Đức (Hà Nội) và các nghiên cứu của Viện Cây lương thực và Cây thực phẩm. Sai số tương đối giữa mô hình và thực nghiệm chỉ dao động trong khoảng 10% đến 15%, khẳng định tính khả thi của mô hình sinh địa hóa tại miền Bắc Việt Nam.

Thảo luận kết quả

Sự chênh lệch phát thải giữa hai vụ lúa bắt nguồn trực tiếp từ điều kiện nhiệt độ và động thái $Eh$ trong đất. Trong vụ Đông Xuân, nhiệt độ trung bình từ tháng 2 đến tháng 4 duy trì ở mức thấp (18°C đến 22°C), làm chậm tốc độ phân giải chất hữu cơ của vi khuẩn yếm khí. Ngược lại, vụ Hè Thu diễn ra trong các tháng 6 đến tháng 9 với nền nhiệt độ cao (28°C đến 33°C) kết hợp lượng mưa lớn chiếm hơn 85% tổng lượng mưa cả năm, kích hoạt mạnh mẽ hoạt động của vi sinh vật sinh mê-tan.

Dữ liệu kết quả được thể hiện trực quan thông qua hai bản đồ chuyên đề GIS về phân bố thông lượng $CH_4$ toàn vùng và các biểu đồ thanh so sánh lượng phát thải giữa 10 tỉnh thành. Phân tích hồi quy trên các bảng số liệu đất cho thấy mối tương quan thuận rõ nét giữa lượng phát thải $CH_4$ với hàm lượng chất hữu cơ (hệ số $R^2 \approx 0,65$) và tương quan nghịch với chế độ rút nước giữa vụ. Kết quả này tương đồng với các nghiên cứu của Viện Nghiên cứu Lúa Quốc tế (IRRI) tại Philippines và Trung Quốc, khẳng định rằng việc ngập nước thường xuyên cả vụ (Continuous Flooding) chính là nguyên nhân đẩy mức phát thải lên đỉnh điểm.

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

  1. Chuyển đổi kỹ thuật tưới sang quy trình Nông - Lộ - Phơi xen kẽ (AWD): Cục Trồng trọt cùng Sở Nông nghiệp & PTNT các tỉnh cần chỉ đạo áp dụng kỹ thuật tưới ngập - khô xen kẽ trên tối thiểu 40% diện tích lúa 2 vụ vào năm 2028. Giải pháp này giúp rút nước mặt ruộng 2 - 3 lần vào giai đoạn đẻ nhánh rộ và sau khi trổ, hướng đến mục tiêu giảm từ 25% đến 35% lượng $CH_4$ phát thải mà vẫn duy trì năng suất lúa ổn định.
  2. Chuẩn hóa quy trình quản lý rơm rạ và phân bón hữu cơ: Trung tâm Khuyến nông các địa phương cần tuyên truyền chấm dứt triệt để việc đốt rơm rạ lộ thiên hoặc vùi tươi rơm rạ vào đất ngập nước trước vụ cấy dưới 20 ngày. Khuyến khích sử dụng chế phẩm vi sinh ủ rơm rạ hoặc thu gom làm nguyên liệu nấm/năng lượng sinh học, đặt mục tiêu giảm 20% lượng phát tán khí nhà kính hữu cơ trước năm 2027.
  3. Điều chỉnh liều lượng phân đạm hóa học: Các hợp tác xã nông nghiệp tại Thái Bình, Hải Dương, Nam Định cần áp dụng bảng so màu lá lúa (LCC) để tối ưu hóa lượng phân đạm bón vào đất, duy trì dưới ngưỡng 100 kg N/ha/vụ, giúp ngăn chặn sự ức chế vi khuẩn oxy hóa mê-tan và giảm 15% lượng phát thải liên quan đến chu trình ni-tơ.
  4. Ứng dụng công nghệ GIS - Viễn thám trong giám sát phát thải thường niên: Bộ Tài nguyên & Môi trường cần phối hợp cùng các viện nghiên cứu xây dựng hệ thống tự động cập nhật bản đồ phát thải khí nhà kính dựa trên dữ liệu ảnh vệ tinh quang học và radar định kỳ hàng năm, hoàn thiện cơ sở dữ liệu số phục vụ báo cáo đóng góp do quốc gia tự quyết định (NDC) trước năm 2030.

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

  • Cơ quan quản lý nhà nước về môi trường và nông nghiệp: Cán bộ Cục Biến đổi khí hậu (Bộ TN&MT), Vụ Khoa học Công nghệ và Môi trường (Bộ Nông nghiệp & PTNT) sử dụng kết quả nghiên cứu làm căn cứ khoa học để hoàn thiện hệ số phát thải $CH_4$ Bậc 2 và Bậc 3 trong kiểm kê khí nhà kính quốc gia.
  • Chuyên gia và kỹ sư công nghệ địa không gian (GIS & Viễn thám): Các nhà nghiên cứu ứng dụng chuỗi ảnh quang học đa thời gian (MODIS, Sentinel, Landsat) và mô hình hóa không gian để giám sát mùa màng, lập bản đồ thảm phủ thực vật và đánh giá tác động môi trường nông nghiệp.
  • Nhà khoa học, giảng viên và sinh viên ngành Môi trường - Thổ nhưỡng: Các trường đại học khối tài nguyên, địa lý và nông nghiệp sử dụng tài liệu làm giáo trình tham khảo chuyên sâu về mô hình hóa sinh địa hóa đất ngập nước bằng phần mềm DNDC.
  • Các tổ chức phi chính phủ và đơn vị phát triển dự án tín chỉ carbon: Các nhà phát triển dự án carbon nông nghiệp tham khảo phương pháp luận đo đếm, báo cáo và thẩm tra (MRV) để xây dựng phương án giảm phát thải trong sản xuất lúa theo tiêu chuẩn quốc tế.

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

Tại sao ảnh viễn thám MODIS được lựa chọn thay vì các ảnh có độ phân giải cao hơn như Landsat hay SPOT? Ảnh MODIS sở hữu tần suất chụp lặp lại rất cao (1 - 2 ngày/lần) và cung cấp sản phẩm chuẩn tổ hợp 8 ngày (MOD09A1), giúp loại bỏ hiệu quả hiện tượng mây mù dày đặc tại khu vực nhiệt đới gió mùa của Đồng bằng sông Hồng. Trong khi đó, ảnh Landsat có chu kỳ chụp 16 ngày nên thường xuyên bị mây che phủ trong giai đoạn gieo cấy lúa.

Mô hình DNDC cần những nhóm tham số đầu vào cơ bản nào để chạy mô phỏng phát thải? Mô hình DNDC yêu cầu 4 nhóm thông tin chính: Dữ liệu khí tượng hàng ngày (nhiệt độ tối cao, tối thấp, lượng mưa); Thuộc tính đất tầng mặt (pH, dung trọng, hàm lượng các-bon hữu cơ SOC, tỷ lệ sét); Nông lịch và đặc tính sinh trưởng của giống lúa; Chế độ canh tác (lịch cày xới, lượng phân bón N-P-K, thời gian ngập và rút nước).

Tại sao phần lớn khí $CH_4$ từ ruộng lúa lại phát thải qua thân cây lúa thay vì bốc hơi qua mặt nước? Cây lúa nước sở hữu hệ thống mô khí (aerenchyma) phát triển từ rễ lên đến bẹ lá và phiến lá để cung cấp oxy cho rễ ngập nước. Cấu trúc rỗng này vô tình tạo thành một "ống dẫn khí" trực tiếp đưa $CH_4$ từ vùng rễ yếm khí thoát thẳng ra ngoài khí quyển, chiếm tới 90% tổng lượng phát thải của ruộng lúa.

Lượng phát thải $CH_4$ ở vụ Hè Thu khác biệt như thế nào so với vụ Đông Xuân tại miền Bắc? Phát thải $CH_4$ vụ Hè Thu cao hơn vụ Đông Xuân khoảng 40% đến 60%. Nguyên nhân là do nhiệt độ vụ Hè Thu cao (trên 30°C) thúc đẩy tốc độ phân hủy chất hữu cơ của vi sinh vật, đồng thời lượng mưa mùa hè chiếm hơn 85% tổng lượng mưa cả năm tạo môi trường ngập nước yếm khí liên tục.

Làm thế nào để kiểm chứng độ tin cậy của mô hình DNDC khi tính toán ở quy mô vùng rộng lớn? Độ tin cậy được kiểm chứng bằng cách so sánh đầu ra mô phỏng với dữ liệu đo đạc thực tế tại các trạm thực nghiệm nông nghiệp (sử dụng phương pháp buồng đo kín) và đối chiếu diện tích lúa chiết tách từ viễn thám với số liệu thống kê chính thức của Tổng cục Thống kê, duy trì sai số dưới 15%.

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 ảnh viễn thám MODIS đa thời gian và mô hình sinh địa hóa DNDC để giám sát không gian và định lượng phát thải khí $CH_4$ từ lúa nước trên toàn vùng Đồng bằng sông Hồng.
  • Thuật toán phân loại mùa vụ dựa trên chỉ số kết hợp LSWI/NDVI/EVI đạt độ chính xác cao trên 85%, phản ánh chính xác lịch gieo cấy và thu hoạch của 2 vụ lúa chính trong năm 2010.
  • Phát hiện sự phân hóa phát thải sâu sắc theo mùa vụ và không gian, trong đó vụ Hè Thu phát thải cao hơn vụ Đông Xuân từ 40% đến 60%; các tỉnh Thái Bình, Nam Định và Hải Dương là những tâm điểm phát thải lớn nhất.
  • Mô hình DNDC cho thấy độ thích ứng cao với điều kiện nông học miền Bắc Việt Nam, có sai số so với đo đạc thực nghiệm chỉ từ 10% đến 15%, đủ điều kiện áp dụng cho kiểm kê khí nhà kính Bậc 3.
  • Trong giai đoạn 2026 - 2030, việc mở rộng mô hình kết hợp với dữ liệu viễn thám độ phân giải cao (Sentinel-1 SAR, Sentinel-2) sẽ là hướng đi then chốt để xây dựng hệ thống bản đồ số phát thải khí nhà kính tự động phục vụ nông nghiệp phát thải thấp. Hãy áp dụng ngay quy trình tưới Nông - Lộ - Phơi và quản lý phân bón thông minh để chung tay giảm nhẹ phát thải khí nhà kính quốc gia.