Tổng quan nghiên cứu

Trong cơ cấu phát thải khí nhà kính tại Việt Nam, lĩnh vực nông nghiệp đóng góp khoảng 88,35 triệu tấn CO2 tương đương theo báo cáo kiểm kê quốc gia năm 2010. Đáng chú ý, hoạt động canh tác lúa nước là nguồn phát sinh khí nhà kính lớn nhất trong ngành, chiếm tỷ trọng tới 50,5%, tương đương 44,60 triệu tấn CO2 tương đương. Tuy nhiên, các phương pháp kiểm kê truyền thống trước đây chủ yếu dựa trên hệ số phát thải mặc định bậc 1 của Tổ chức Liên chính phủ về Biến đổi Khí hậu năm 1996, vốn chưa phản ánh được sự biến thiên phức tạp về thổ nhưỡng, thời tiết và tập quán canh tác giữa các tiểu vùng sinh thái.

Nghiên cứu của luận văn tập trung giải quyết bài toán định lượng chi tiết lượng phát thải hai loại khí nhà kính chủ yếu là mêtan (CH4) và oxit nitơ (N2O) từ các hệ thống canh tác lúa nước. Mục tiêu cụ thể là mô phỏng quá trình sinh địa hóa trên từng đơn vị đất và tiểu vùng khí hậu, từ đó quy đổi sang tổng lượng phát thải khí nhà kính tương đương (CO2e) trên phạm vi toàn tỉnh Nam Định.

Phạm vi nghiên cứu bao quát toàn bộ diện tích tự nhiên 1.652,29 km2 của tỉnh Nam Định, gồm 1 thành phố và 9 huyện, kết hợp chuỗi số liệu quan trắc khí tượng và khảo sát thực địa giai đoạn 2013-2015. Kết quả tính toán cung cấp cơ sở khoa học định lượng có độ phân giải cao theo không gian, trực tiếp phục vụ mục tiêu quốc gia về giảm 8% lượng phát thải khí nhà kính vào năm 2030 bằng nguồn lực trong nước (và có thể đạt 25% khi có hỗ trợ quốc tế) theo Quyết định số 2053/QĐ-TTg và Quyết định số 2359/QĐ-CP của Thủ tướng Chính phủ.

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 nền tảng lý thuyết sinh địa hóa đất và khung hướng dẫn kiểm kê khí nhà kính quốc tế của IPCC phiên bản 1996 và 2006, kết hợp hướng dẫn thực hành tốt đối với lĩnh vực sử dụng đất và thay đổi sử dụng đất năm 2003. Khung phân tích vận dụng mô hình sinh địa hóa DNDC (DeNitrification - DeComposition) để mô phỏng động thái chuyển hóa các-bon và nitơ trong hệ sinh thái đất nông nghiệp ngập nước.

Các khái niệm then chốt được chuẩn hóa bao gồm:

  • Quá trình sinh mêtan (Methanogenesis): Quá trình vi sinh vật kỵ khí chuyển hóa cơ chất hữu cơ trong điều kiện khử nghiêm ngặt với thế oxy hóa khử dao động từ -250 mV đến -350 mV, phát sinh khí CH4 thoát vào khí quyển qua hệ thống gian bào cây lúa và khuếch tán bề mặt.
  • Quá trình nitrate hóa và phản nitrate hóa (Nitrification - Denitrification): Chuỗi phản ứng sinh học chuyển đổi ion amoni (NH4+) sang nitrate (NO3-) và khử tiếp thành khí N2O, NO, N2 trong các tầng đất có chế độ oxy biến động.
  • Tiềm năng làm nóng toàn cầu (Global Warming Potential - GWP): Hệ số quy đổi tác động bức xạ của khí CH4 và N2O sang đơn vị CO2 tương đương (CO2e) trong chu kỳ 100 năm.
  • Đơn vị tổ hợp tự nhiên - canh tác: Phân vùng đồng nhất về khí hậu, loại đất và quy trình quản lý phân bón, thủy lợi làm đơn vị cơ sở để tính toán phát thải.

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

Nghiên cứu sử dụng phương pháp mô hình hóa kết hợp với Hệ thống Thông tin Địa lý (GIS) và thực nghiệm đồng ruộng nhằm chuyển đổi các thông số điểm thành bản đồ phát thải liên tục.

Về nguồn dữ liệu và cỡ mẫu:

  • Dữ liệu khí tượng đầu vào gồm chuỗi số liệu ngày đo đạc liên tục trong 3 năm (2013-2015) tại 4 trạm khí tượng đại diện là Nam Định, Văn Lý, Ninh Bình và Thái Bình, tương ứng hơn 1.095 bản ghi dữ liệu cho mỗi biến số nhiệt độ tối cao, nhiệt độ tối thấp, độ ẩm, lượng mưa, bức xạ mặt trời và tốc độ gió.
  • Dữ liệu không gian bao gồm bản đồ hiện trạng sử dụng đất và bản đồ thổ nhưỡng tỉnh Nam Định tỷ lệ 1:50.000 với 8 nhóm đất trồng lúa chính được chuẩn hóa từ cơ sở dữ liệu của Tổng cục Quản lý đất đai và Bảo tàng Đất Việt Nam.
  • Dữ liệu kiểm chứng thực địa được thu thập từ 2 điểm quan trắc đo đạc phát thải khí nhà kính tại Thịnh Long và Rạng Đông thuộc khuôn khổ dự án nghiên cứu biến đổi khí hậu nông nghiệp phối hợp cùng Đại sứ quán Na Uy năm 2014.

Phương pháp chọn mẫu và phân vùng không gian sử dụng thuật toán đa giác Thiessen để xây dựng bản đồ phân vùng khí hậu, sau đó tiến hành kỹ thuật chồng xếp không gian (Spatial Overlay) giữa 3 lớp thông tin bản đồ khí hậu, bản đồ đất và bản đồ chế độ canh tác để hình thành các đơn vị mô phỏng đồng nhất.

Lý do lựa chọn mô hình DNDC kết hợp GIS là khả năng mô phỏng chi tiết các quá trình sinh địa hóa theo bước thời gian hàng ngày, khắc phục nhược điểm áp dụng hệ số tĩnh của phương pháp IPCC bậc 1. Mô hình được hiệu chỉnh và đánh giá độ tin cậy thông qua hệ số xác định R2 và chỉ số hiệu quả Nash - Sutcliffe (NSI), đảm bảo cả hai chỉ số đều đạt ngưỡng lớn hơn 0,5 trước khi áp dụng mô phỏng diện rộng trên toàn tỉnh.

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

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

Mô phỏng DNDC kết hợp phân tích không gian đã mang lại các phát hiện định lượng then chốt về cường độ phát thải khí nhà kính trong canh tác lúa tại Nam Định:

Thứ nhất, lượng phát thải khí CH4 và N2O có sự phân hóa rõ nét theo từng loại đất và tiểu vùng sinh thái. Các diện tích đất phù sa được bồi và đất phù sa gley có mức phát thải mêtan trung bình năm cao hơn rõ rệt so với đất cát ven biển và đất mặn nhiều, do hàm lượng các-bon hữu cơ trong đất phù sa cao hơn tạo nguồn cơ chất dồi dào cho vi khuẩn kỵ khí.

Thứ hai, chế độ quản lý nước tác động mang tính quyết định đến tổng lượng khí phát sinh. So sánh kết quả mô hình hóa cho thấy chế độ tưới ngập nước thường xuyên sinh ra lượng CH4 chiếm tới hơn 90% tổng phát thải khí nhà kính từ ruộng lúa, tương đương với bức tranh kiểm kê toàn quốc khi lúa tưới ngập chiếm khoảng 2,0 triệu tấn CH4 trong khi lúa nhờ nước trời chỉ phát sinh 0,16 triệu tấn CH4.

Thứ ba, kỹ thuật bón phân nitơ khoáng và quản lý phụ phẩm hữu cơ quyết định trực tiếp cường độ phát thải N2O. Việc bón vãi phân đạm trên mặt nước làm tăng tỷ lệ mất nitơ qua bay hơi NH3 và phản nitrate hóa, trong khi các công thức bón vùi sâu kết hợp phân đạm dạng amoni sunfat giúp giảm phát thải khí nhà kính quy đổi xuống khoảng 30% đến 50% so với bón phân hữu cơ tươi không ủ hoai.

Thứ tư, phân bố phát thải theo địa giới hành chính chỉ ra rằng các huyện có quy mô thâm canh lúa lớn như Hải Hậu, Nghĩa Hưng, Ý Yên và Trực Ninh đóng góp tỷ trọng lớn nhất vào tổng lượng CO2e toàn tỉnh, phản ánh sự tương quan thuận giữa diện tích gieo cấy lúa 2 vụ và mật độ phát thải khí nhà kính.

Thảo luận kết quả

Cơ chế sinh địa hóa giải thích sự phát thải vượt trội của CH4 tại các chân ruộng trũng là do việc ngập nước kéo dài làm triệt tiêu nguồn oxy hòa tan, hạ thấp thế oxy hóa khử xuống dưới -250 mV. Khi đó, vi khuẩn sinh mêtan chiếm ưu thế tuyệt đối trong việc khử CO2 và phân giải axit axetic. Tại các vùng đất mặn ven 72 km bờ biển của Nam Định, sự hiện diện của ion sunfat (SO42-) trong nước biển xâm nhập đã kìm hãm một phần hoạt động của vi sinh vật sinh mêtan do sự cạnh tranh cơ chất của vi khuẩn khử sunfat.

So sánh với các nghiên cứu quốc tế tại đồng bằng sông Nile (Ai Cập) và các vùng lúa nước tại Trung Quốc hay Đức, kết quả mô phỏng tại Nam Định hoàn toàn tương đồng về động thái phát thải theo mùa vụ. Cụ thể, lượng phát thải CH4 đạt đỉnh vào giai đoạn lúa đẻ nhánh rộ và làm đòng khi sinh khối rễ phát triển mạnh nhất, cung cấp dịch tiết rễ làm thức ăn cho vi sinh vật.

Toàn bộ dữ liệu phát thải trong nghiên cứu được thể hiện trực quan thông qua hệ thống bản đồ chuyên đề GIS, bao gồm bản đồ đơn vị tổ hợp điều kiện tự nhiên, bản đồ phân bố phát thải CH4 (kg/ha/năm), bản đồ phát thải N2O (kg/ha/năm) và bản đồ tổng lượng phát thải khí nhà kính quy đổi CO2e. Các dữ liệu này được tổng hợp thành các bảng ma trận tra cứu chéo theo 10 đơn vị hành chính và 8 loại đất chính, giúp cơ quan quản lý dễ dàng nhận diện các điểm nóng phát thải để ưu tiên can thiệp kỹ thuật.

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

Dựa trên kết quả định lượng bằng mô hình DNDC, bốn nhóm giải pháp hành động cụ thể được khuyến nghị nhằm giảm thiểu phát thải khí nhà kính trong nông nghiệp tỉnh Nam Định:

  1. Chuyển đổi quy trình thủy lợi sang kỹ thuật tưới ngập khô xen kẽ (AWD): Sở Nông nghiệp và Phát triển nông thôn Nam Định chỉ đạo các công ty khai thác công trình thủy lợi phối hợp cùng hợp tác xã nông nghiệp áp dụng kỹ thuật tưới rút nước giữa vụ cho ít nhất 60% diện tích lúa 2 vụ. Giải pháp này hướng tới mục tiêu giảm 30% đến 40% lượng phát thải khí CH4 trên mỗi hecta canh tác, thực hiện đồng bộ trong giai đoạn 2026-2030.

  2. Cải tiến kỹ thuật bón phân và dinh dưỡng cây trồng: Trung tâm Khuyến nông tỉnh phối hợp Chi cục Trồng trọt và Bảo vệ thực vật triển khai kỹ thuật bón phân đạm sâu, sử dụng phân bón tan chậm và bổ sung phân amoni sunfat thay thế một phần urê thông thường, khống chế lượng đạm nguyên chất dưới 100 kg N/ha/vụ. Mục tiêu đặt ra là cắt giảm 20% đến 25% lượng phát thải N2O trực tiếp từ đất ruộng, hoàn thành mở rộng kỹ thuật trước năm 2028.

  3. Thu gom và xử lý sinh học phụ phẩm rơm rạ: Hội Nông dân tỉnh phát động phong trào chấm dứt hoàn toàn tập quán đốt rơm rạ lộ thiên ngoài đồng, chuyển sang sử dụng chế phẩm vi sinh để ủ phân hữu cơ hiếu khí hoặc cày vùi rơm rạ kết hợp bón vôi xử lý sớm trước khi gieo cấy ít nhất 20 ngày. Mục tiêu làm giảm 15% tổng lượng CO2e phát sinh từ phân hủy hữu cơ kỵ khí, triển khai ngay trong giai đoạn 2026-2027.

  4. Số hóa hệ thống đo đạc, báo cáo và thẩm tra (MRV) kiểm kê khí nhà kính: Sở Tài nguyên và Môi trường tỉnh Nam Định chủ trì xây dựng cơ sở dữ liệu GIS tích hợp mô hình DNDC để định kỳ 2 năm một lần cập nhật bản đồ phát thải khí nhà kính cấp huyện và xã. Hoạt động này tạo lập dữ liệu cơ sở cho 10 đơn vị hành chính, hoàn thiện hệ thống quan trắc phát thải vào năm 2030 phục vụ thị trường tín chỉ carbon.

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

Nội dung và kết quả của công trình nghiên cứu mang giá trị ứng dụng cao cho 4 nhóm đối tượng trọng tâm:

  • Cơ quan quản lý nhà nước về tài nguyên, môi trường và nông nghiệp: Cán bộ chuyên trách tại Sở Tài nguyên và Môi trường, Sở Nông nghiệp và Phát triển nông thôn có thể sử dụng phương pháp luận và dữ liệu bản đồ để xây dựng Kế hoạch hành động ứng phó biến đổi khí hậu cấp tỉnh, tích hợp chỉ tiêu giảm phát thải vào quy hoạch sử dụng đất giai đoạn 2026-2030.

  • Các nhà khoa học, giảng viên và học viên cao học: Các chuyên gia nghiên cứu trong lĩnh vực Biến đổi khí hậu, Khoa học Môi trường, Nông học và Quản lý Đất đai có thể tham khảo quy trình tích hợp mô hình sinh địa hóa DNDC với GIS, phương pháp hiệu chỉnh mô hình bằng chỉ số NSI và R2 để mở rộng ứng dụng cho các lưu vực sông khác.

  • Tổ chức phát triển quốc tế và doanh nghiệp phát triển tín chỉ carbon: Các tổ chức phi chính phủ và định chế tài chính khí hậu có thể khai thác bộ số liệu đường phát thải cơ sở (Baseline) để thiết kế các dự án tín chỉ carbon nông nghiệp theo tiêu chuẩn quốc tế và Điều 6 của Thỏa thuận Paris.

  • Cán bộ khuyến nông và quản lý hợp tác xã nông nghiệp: Kỹ thuật viên nông nghiệp tại các địa phương thuộc vùng đồng bằng sông Hồng có thể sử dụng các kết quả so sánh kỹ thuật canh tác để trực tiếp xây dựng tài liệu hướng dẫn nông dân điều tiết nước và bón phân giảm phát thải.

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

Mô hình DNDC mang lại ưu thế gì vượt trội so với phương pháp tính hệ số mặc định của IPCC? Mô hình DNDC mô phỏng các chu trình sinh hóa theo ngày dựa trên 4 biến số động là khí tượng, đặc tính đất, sinh lý cây trồng và biện pháp quản lý nông học. Nhờ đó, mô hình phản ánh chính xác sự sai khác về lượng phát thải giữa từng thửa ruộng và từng mùa vụ, khắc phục tính cào bằng của hệ số tĩnh IPCC bậc 1.

Tại sao đất trồng lúa nước lại phát sinh lượng khí mêtan cao hơn các loại cây trồng cạn? Đất lúa ngập nước tạo môi trường yếm khí tuyệt đối, đẩy thế oxy hóa khử xuống mức -250 mV đến -350 mV. Ở điều kiện này, các vi khuẩn sinh mêtan phân hủy chất hữu cơ để tạo ra khí CH4, sau đó khuếch tán qua mô khí của thân cây lúa và thoát ra ngoài không khí.

Độ tin cậy của kết quả mô phỏng trong luận văn được bảo đảm bằng cách nào? Các tham số mô hình DNDC được hiệu chỉnh trực tiếp với số liệu đo đạc thực địa tại 2 trạm Thịnh Long và Rạng Đông. Độ chính xác được kiểm chứng thông qua chỉ số hiệu quả Nash - Sutcliffe (NSI) và hệ số tương quan R2 đều đạt trên 0,5, khẳng định sai số mô phỏng nằm trong giới hạn khoa học cho phép.

Biện pháp tưới ngập khô xen kẽ giúp cắt giảm phát thải khí nhà kính như thế nào? Việc chủ động rút nước để đất nứt chân chim làm gián đoạn trạng thái yếm khí, đưa oxy vào tầng đất mặt và nâng cao thế oxy hóa khử. Quá trình này ức chế vi khuẩn sinh mêtan, giúp giảm từ 30% đến 40% lượng phát thải CH4 mà vẫn đảm bảo năng suất lúa.

Quy trình nghiên cứu trong luận văn có thể áp dụng cho các tỉnh khác tại đồng bằng sông Hồng không? Quy trình tích hợp mô hình DNDC và công nghệ GIS có tính tổng quát cao. Các địa phương như Thái Bình, Hải Phòng, Ninh Bình hoàn toàn có thể tái sử dụng khung phương pháp này bằng cách nạp dữ liệu khí tượng từ 3 đến 4 trạm đo cục bộ và bản đồ thổ nhưỡng tỷ lệ 1:50.000 tương ứng.

Kết luận

  • Tích hợp thành công mô hình sinh địa hóa DNDC và hệ thống GIS để mô phỏng chính xác động thái phát thải khí nhà kính CH4 và N2O trên toàn bộ đất lúa tỉnh Nam Định.
  • Thiết lập bộ bản đồ chuyên đề và dữ liệu định lượng về tổng phát thải CO2e phân theo 10 đơn vị hành chính và 8 nhóm đất canh tác chính.
  • Xác định mối tương quan giữa điều kiện ngập nước, liều lượng phân bón nitơ với cường độ phát thải khí nhà kính, chứng minh tiềm năng giảm phát thải lên tới 40% qua các gói kỹ thuật cải tiến.
  • Cung cấp cơ sở khoa học thực tiễn phục vụ công tác kiểm kê khí nhà kính định kỳ cấp quốc gia theo Quyết định số 2359/QĐ-CP và Quyết định số 2053/QĐ-TTg.
  • Mở ra hướng tiếp cận công nghệ số hóa trong giám sát và xây dựng kịch bản giảm nhẹ biến đổi khí hậu cho ngành trồng trọt.

Đóng góp lớn nhất của luận văn là đã chuyển đổi phương pháp kiểm kê phát thải từ dạng thống kê tĩnh sang mô hình hóa sinh địa hóa có độ phân giải không gian cao tại vùng trọng điểm lúa Nam Định. Trong giai đoạn 2026-2030, các cơ quan quản lý và địa phương cần đẩy mạnh nhân rộng mô hình giám sát số hóa này trên quy mô toàn vùng đồng bằng sông Hồng, đồng thời triển khai rộng rãi các gói kỹ thuật tưới ngập khô xen kẽ để hiện thực hóa cam kết phát thải ròng bằng không trong nông nghiệp xanh bền vững.