FracLEACH được sử dụng để làm gì trong kiểm kê khí nhà kính?

FracLEACH được sử dụng để ước tính tỷ lệ nitơ đưa vào đất nông nghiệp bị thất thoát qua rửa trôi và dòng chảy mặt vào các thủy vực. Thông số này là căn cứ trực tiếp để tính toán lượng khí nitơ oxit (N2O) gián tiếp phát thải từ hoạt động bón phân và quản lý đất.

Giá trị mặc định của FracLEACH theo IPCC là bao nhiêu?

Theo Hướng dẫn IPCC 2006, giá trị mặc định bậc 1 của FracLEACH là 0.30 kg N trên mỗi kg N bổ sung. Trong Bản hoàn thiện 2019, giá trị này được cập nhật thành 0.24 cho các vùng ẩm ướt hoặc đất có tưới tiêu, và bằng 0 cho các vùng khô hạn không tưới tiêu.

Tại sao nhiều quốc gia xây dựng giá trị FracLEACH đặc thù thay vì dùng giá trị mặc định?

Giá trị mặc định toàn cầu có khoảng bất định rất lớn và không phản ánh chính xác các đặc trưng thổ nhưỡng, lượng mưa và tập quán canh tác địa phương. Nhiều quốc gia như New Zealand và Đức đã xác định giá trị đặc thù bậc 2 hoặc bậc 3 thấp hơn đáng kể nhằm nâng cao độ chính xác của báo cáo kiểm kê quốc gia.

Nitơ bị rửa trôi chuyển hóa thành N2O như thế nào?

Nitơ rửa trôi khỏi đất chủ yếu dưới dạng ion nitrat hòa tan xâm nhập vào nước ngầm và hệ thống sông ngòi. Trong môi trường nước thiếu oxy, các vi sinh vật thực hiện quá trình khử nitrat kỵ khí, chuyển hóa một phần nitrat thành khí nitơ oxit và phát tán vào khí quyển.

FracLEACH (nitrogen leaching fraction) là gì?

FracLEACH (hoặc FracLEACH-(H)) là thông số kỹ thuật không thứ nguyên trong các mô hình kiểm kê khí nhà kính nông nghiệp, biểu thị tỷ lệ lượng nitơ bổ sung vào đất bị thất thoát qua quá trình rửa trôi vào tầng nước ngầm và chảy tràn bề mặt vào các thủy vực tự nhiên. Khái niệm này được phát triển bởi Ban Liên chính phủ về Biến đổi Khí hậu (IPCC) nhằm định lượng nguồn phát thải gián tiếp của khí nitơ oxit (N2O) từ các hoạt động quản lý đất trồng trọt và chăn nuôi. Bài viết cung cấp phân tích toàn diện về định nghĩa, cơ chế sinh địa hóa, công thức tính toán, diễn biến cập nhật hệ số qua các giai đoạn, thực tiễn áp dụng tại Việt Nam cũng như các hạn chế khoa học cần hoàn thiện.

Định nghĩa và bản chất thông số

Trong quản lý dinh dưỡng cây trồng và kiểm kê phát thải khí nhà kính, nitơ đưa vào đất nông nghiệp qua phân bón khoáng, phân chuồng, chất thải gia súc và phụ phẩm cây trồng không được thực vật hấp thu toàn bộ. Một phần nitơ hòa tan sẽ di chuyển xuống các tầng đất sâu hơn vùng rễ hoạt động hoặc bị cuốn trôi theo dòng chảy mặt khi có mưa lớn và tưới tiêu dư thừa. Thông số FracLEACH xác định tỷ lệ lượng nitơ bị mất này so với tổng lượng nitơ đưa vào hoặc khoáng hóa trong đất.

Bản chất của FracLEACH là một hệ số tỷ lệ được chuẩn hóa, phản ánh nguy cơ tổn thất nitơ của hệ thống canh tác vào môi trường thủy sinh. Đại lượng này đóng vai trò mắt xích then chốt trong chuỗi tính toán kiểm kê khí nhà kính nông nghiệp: từ lượng phân bón sử dụng, mô hình xác định lượng nitơ rửa trôi, sau đó tiếp tục xác định lượng khí nitơ oxit sinh ra ở hạ lưu các lưu vực sông ngòi, tầng chứa nước và vùng cửa sông ven biển.

Cơ chế sinh địa hóa của quá trình rửa trôi nitơ

Quá trình rửa trôi nitơ từ đất nông nghiệp gắn liền với các phản ứng chuyển hóa sinh địa hóa của chu trình nitơ trong đất. Khi phân bón hoặc chất hữu cơ được đưa vào đất, các hợp chất hữu cơ chứa nitơ và ion amoni trải qua quá trình khoáng hóa và nitrat hóa dưới tác động của vi sinh vật đất. Theo phân tích của Butterbach-Bahl và cộng sự (2013), nitrat hóa hiếu khí và khử nitrat kỵ khí là hai con đường sinh học trung tâm chi phối quá trình tạo thành và phát thải khí nitơ oxit.

Trong môi trường đất thoáng khí, vi khuẩn nitrat hóa chuyển đổi ion amoni thành ion nitrit và tiếp tục thành ion nitrat. Do keo đất thường mang điện tích âm, các ion nitrat cũng mang điện tích âm nên không bị hấp phụ tĩnh điện bởi phức hệ hấp thu của đất. Do đó, nitrat tồn tại tự do trong dung dịch đất và dễ dàng hòa tan trong nước ngấm. Khi lượng nước ngấm xuống vượt quá sức giữ ẩm đồng ruộng, dòng nước thấm lọc mang theo nitrat di chuyển xuống tầng nước ngầm hoặc xuất lộ ra mương tiêu, suối và sông.

Bouwman và cộng sự (2002) qua phân tích tổng hợp dữ liệu toàn cầu đã chỉ ra rằng mức độ rửa trôi nitơ không chỉ phụ thuộc vào lượng nitơ bón vào đất mà còn chịu sự điều tiết chặt chẽ của kết cấu đất, khả năng giữ nước, chế độ tưới tiêu và cường độ hấp thu của thảm thực vật. Sau khi bị rửa trôi ra khỏi tầng đất mặt, ion nitrat xâm nhập vào môi trường nước ngầm và hệ thống nước mặt. Tại đây, trong các vi môi trường thiếu oxy, vi khuẩn khử nitrat sử dụng nitrat làm chất nhận điện tử, khử dần nitrat thành nitrit, oxit nitric, nitơ oxit và cuối cùng là khí nitơ phân tử. Khí nitơ oxit thoát ra khí quyển từ các thủy vực chính là thành phần phát thải gián tiếp được lượng hóa thông qua FracLEACH.

Phương pháp luận tính toán phát thải theo IPCC

Theo Hướng dẫn kiểm kê khí nhà kính quốc gia của IPCC (2006), phân cấp tính toán phát thải nông nghiệp được chia thành ba bậc:

Bậc 1 (Tier 1): Áp dụng các hệ số phát thải và thông số tỷ lệ mặc định toàn cầu do IPCC cung cấp, không phân biệt chi tiết điều kiện thổ nhưỡng và khí hậu vi mô.
Bậc 2 (Tier 2): Sử dụng các hệ số đặc thù quốc gia hoặc vùng sinh thái dựa trên số liệu nghiên cứu thực nghiệm và điều tra thực địa trong nước.
Bậc 3 (Tier 3): Ứng dụng các mô hình mô phỏng quá trình sinh địa hóa học chi tiết kết hợp dữ liệu không gian có độ phân giải cao và chuỗi thời gian liên tục.

Trong phương pháp luận bậc 1 và bậc 2, lượng nitơ thất thoát qua rửa trôi và dòng chảy mặt, ký hiệu là $N_{LEACH}$ , được xác định theo công thức sau:

N_{LEACH} = N_{inputs} \times Frac_{LEACH-(H)}

Trong đó, biến số $N_{LEACH}$ biểu thị tổng lượng nitơ bị rửa trôi tính bằng kilôgam nitơ; biến số $N_{inputs}$ biểu thị tổng lượng nitơ đầu vào từ phân bón khoáng, phân bón hữu cơ, chất thải chăn nuôi và tàn dư cây trồng bổ sung vào đất; và thông số $Frac_{LEACH-(H)}$ là tỷ lệ nitơ thất thoát qua rửa trôi và chảy tràn ở các khu vực xảy ra quá trình này.

Sau khi xác định được lượng nitơ rửa trôi, lượng phát thải nitơ oxit gián tiếp, ký hiệu là $N_2O_{(L)}-N$ , được tính bằng cách nhân lượng nitơ rửa trôi với hệ số phát thải đặc thù:

N_2O_{(L)}-N = N_{LEACH} \times EF_5

Trong công thức trên, biến số $N_2O_{(L)}-N$ là lượng phát thải nitơ oxit gián tiếp tính theo đơn vị kilôgam nitơ; biến số $N_{LEACH}$ là lượng nitơ bị rửa trôi; và thông số $EF_5$ là hệ số phát thải nitơ oxit từ nitơ rửa trôi và dòng chảy mặt.

Diễn biến cập nhật giá trị FracLEACH qua các giai đoạn

Trong tài liệu Hướng dẫn IPCC 2006, Volume 4, Chương 11, giá trị mặc định bậc 1 của FracLEACH-(H) được quy định là 0.30 kg N trên mỗi kg N bổ sung, đi kèm khoảng bất định rất rộng từ 0.10 đến 0.80. Hệ số phát thải EF5 tương ứng khi đó có giá trị mặc định là 0.0075 kg N2O-N trên mỗi kg N rửa trôi. Mức mặc định 0.30 này giả định rằng một phần ba tổng lượng phân bón nitơ đưa vào đất nông nghiệp đều bị mất vào nguồn nước ở bất kỳ khu vực nào có phát sinh dòng chảy.

Tuy nhiên, các bằng chứng thực nghiệm trong giai đoạn sau đó cho thấy giả định trên không phản ánh đúng thực tế tại nhiều vùng khí hậu. Để khắc phục sai lệch này, trong Bản hoàn thiện 2019 của Hướng dẫn IPCC 2006 tại Bảng 11.3, IPCC đã cập nhật và phân tách giá trị FracLEACH-(H) theo điều kiện lượng mưa và thủy văn:

Đối với vùng khí hậu ẩm ướt hoặc đất canh tác có tưới tiêu (nơi lượng mưa hoặc lượng nước tưới vượt quá lượng bốc thoát hơi nước tiềm năng), giá trị mặc định giảm xuống 0.24, với khoảng bất định dao động từ 0.01 đến 0.45.
Đối với vùng khí hậu khô hạn không có hệ thống tưới tiêu bổ sung, giá trị mặc định được xác định bằng 0, bởi vì tại các vùng này nước ngấm không đủ để vận chuyển nitrat vượt qua tầng rễ cây xuống tầng nước ngầm.
Hệ số phát thải tổng hợp EF5 cho nitơ rửa trôi cũng được cập nhật từ 0.0075 lên 0.011 kg N2O-N trên mỗi kg N rửa trôi.

Bên cạnh các giá trị mặc định của IPCC, nhiều quốc gia đã tiến hành nghiên cứu thực địa chuyên sâu nhằm xác định giá trị FracLEACH bậc 2 đặc thù cho điều kiện trong nước. Bảng dưới đây so sánh các giá trị FracLEACH giữa các hướng dẫn quốc tế và các nghiên cứu điển hình:

Nguồn tài liệu	Phạm vi áp dụng	Giá trị FracLEACH	Khoảng bất định
IPCC (2006)	Mặc định bậc 1 toàn cầu	0.30	0.10 - 0.80
IPCC (2019)	Vùng khí hậu ẩm ướt hoặc đất có tưới tiêu	0.24	0.01 - 0.45
IPCC (2019)	Vùng khí hậu khô hạn không tưới tiêu	0	Không áp dụng
Thomas và cộng sự (2005)	Đồng cỏ chăn thả New Zealand	0.07	Nghiên cứu bậc 2 đặc thù
Eysholdt và cộng sự (2022)	Đất nông nghiệp toàn nước Đức	0.099	Mô hình RAUMIS-mGROWA-DENUZ

Nghiên cứu của Thomas và cộng sự (2005) tại New Zealand đã chứng minh rằng trong điều kiện đồng cỏ chăn thả được quản lý chặt chẽ, tỷ lệ rửa trôi thực tế chỉ đạt mức 0.07, thấp hơn nhiều so với giá trị mặc định 0.30 của IPCC. Tương tự, nghiên cứu của Eysholdt và cộng sự (2022) tại Đức thông qua chuỗi mô hình không gian kết hợp RAUMIS-mGROWA-DENUZ đã xác định giá trị FracLEACH bình quân toàn quốc là 0.099 kg N trên mỗi kg N đầu vào, chứng minh sự phân hóa địa phương sâu sắc của quá trình rửa trôi.

Các yếu tố chi phối mức độ rửa trôi nitơ

Mức độ biến thiên của FracLEACH trong thực tế chịu sự tác động tương hỗ của nhiều nhóm yếu tố tự nhiên và nhân tác:

Đặc tính thổ nhưỡng: Đất có thành phần cơ giới thô như đất cát, đất xốp có khả năng giữ nước kém và tốc độ thấm lọc cao, dẫn đến nguy cơ rửa trôi nitrat lớn hơn nhiều so với đất sét hoặc đất giàu chất hữu cơ.
Chế độ thủy văn và lượng mưa: Lượng mưa tích lũy và lượng nước tưới vượt ngưỡng trữ ẩm của tầng đất mặt là động lực vật lý chính đẩy dung dịch đất chứa nitrat xuống sâu. Các đợt mưa lớn tập trung ngay sau khi bón phân gây ra tổn thất rửa trôi đột biến.
Chủng loại và liều lượng phân bón: Bón phân nitơ khoáng vượt quá nhu cầu hấp thu tức thời của cây trồng tạo ra lượng nitơ dư thừa trong đất, làm tăng mạnh tỷ lệ nitrat hòa tan dễ bị rửa trôi.
Hệ thống cây trồng và kỹ thuật canh tác: Cây trồng có bộ rễ ăn sâu và sinh khối lớn hấp thu nitrat hiệu quả hơn, làm giảm nồng độ nitrat trong nước ngấm. Việc trồng cây che phủ vụ đông giúp giữ lại lượng nitơ dư thừa sau vụ thu hoạch chính, hạn chế rửa trôi trong mùa mưa.

Thực tiễn áp dụng và bối cảnh tại Việt Nam

Tại Việt Nam, nông nghiệp đóng vai trò quan trọng trong cơ cấu phát thải khí nhà kính quốc gia, trong đó phát thải khí nitơ oxit từ đất canh tác là một hợp phần trọng yếu. Theo Báo cáo Cập nhật Hai năm một lần lần thứ ba (BUR3) năm 2020 của Bộ Tài nguyên và Môi trường nộp cho Công ước Khung của Liên Hợp Quốc về Biến đổi Khí hậu (UNFCCC), công tác kiểm kê khí nhà kính nông nghiệp quốc gia áp dụng phương pháp luận bậc 1 của IPCC 2006.

Trong báo cáo BUR3 (2020), lượng nitơ rửa trôi từ đất nông nghiệp được tính toán dựa trên số liệu thống kê lượng phân đạm hóa học và nguồn phân bón hữu cơ sử dụng trên quy mô các tỉnh, kết hợp thông số FracLEACH mặc định. Do điều kiện khí hậu nhiệt đới gió mùa tại Việt Nam có lượng mưa hàng năm lớn nhưng phân bố không đều giữa mùa mưa và mùa khô, việc áp dụng hệ số mặc định toàn cầu tĩnh có thể chưa phản ánh chính xác lượng nitơ rửa trôi thực tế giữa các hệ canh tác khác nhau, chẳng hạn như giữa ruộng lúa nước ngập triều và đất trồng cây công nghiệp vùng cao nguyên.

Các nhà nghiên cứu trong nước và các cơ quan chuyên môn đang hướng tới mục tiêu xây dựng các hệ số đặc thù quốc gia bậc 2 cho các vùng sinh thái nông nghiệp trọng điểm. Việc xác định giá trị FracLEACH phù hợp cho từng vùng đất sẽ giúp nâng cao độ chính xác của báo cáo kiểm kê quốc gia, đồng thời cung cấp căn cứ khoa học để xây dựng các biện pháp giảm thiểu phát thải thông qua kỹ thuật bón phân hợp lý và quản lý nước tưới tiết kiệm.

Hạn chế khoa học và hướng nghiên cứu mở

Mặc dù FracLEACH là công cụ không thể thiếu trong các khuôn khổ kiểm kê khí nhà kính tiêu chuẩn, thông số này tồn tại một số hạn chế nội tại đã được cộng đồng khoa học quốc tế thảo luận sâu rộng:

Trước hết, phê bình của Nevison (2000) đã chỉ ra rằng việc sử dụng một hệ số FracLEACH cố định tuyến tính bỏ qua mối quan hệ phi tuyến tính phức tạp giữa lượng phân bón dư thừa và tỷ lệ nitơ rửa trôi. Khi lượng phân bón vượt xa năng suất hấp thu tối đa của cây trồng, tỷ lệ rửa trôi thực tế tăng theo cấp số mũ chứ không duy trì ở một tỷ lệ phần trăm cố định.

Thứ hai, khoảng bất định của FracLEACH trong các hướng dẫn bậc 1 vẫn ở mức rất cao, làm giảm độ tin cậy của các ước tính phát thải nitơ oxit gián tiếp trên phạm vi quốc gia và lưu vực sông. Sự khác biệt giữa giá trị mặc định của IPCC với kết quả mô hình hóa chi tiết tại New Zealand theo Thomas và cộng sự (2005) hay tại Đức theo Eysholdt và cộng sự (2022) cho thấy tính tất yếu của việc chuyển dịch từ phương pháp hệ số cố định sang phương pháp bậc 3 ứng dụng các mô hình động lực học đất và thủy văn.

Hướng nghiên cứu mở tập trung vào việc tích hợp dữ liệu viễn thám, bản đồ thổ nhưỡng số và các mô hình sinh thái nông nghiệp để ước tính FracLEACH theo thời gian thực. Cách tiếp cận này cho phép định lượng chính xác hiệu quả của các biện pháp quản lý dinh dưỡng tiên tiến, đóng góp thiết thực vào mục tiêu giảm phát thải khí nhà kính nông nghiệp toàn cầu.