張 釗， 辛?xí)云?/p>

(中國(guó)農(nóng)業(yè)科學(xué)院農(nóng)業(yè)資源與農(nóng)業(yè)區(qū)劃研究所， 北京 100081)

碳循環(huán)是近年來生態(tài)系統(tǒng)生態(tài)學(xué)的研究熱點(diǎn)，300多年的工業(yè)社會(huì)發(fā)展，極大的改變了人類的生活質(zhì)量和生存方式，同時(shí)也極大的改變了人類和自然的相處模式。化石燃料的利用，礦產(chǎn)資源的開發(fā)，高效農(nóng)業(yè)的普及成為現(xiàn)代社會(huì)的典型特征。在人類文明高速發(fā)展與進(jìn)步的同時(shí)，地球生態(tài)系統(tǒng)也隨之發(fā)生著改變。1958年夏威夷冒納羅亞火山觀測(cè)站最早記錄的CO2濃度數(shù)據(jù)為316 ppm，隨著工業(yè)社會(huì)的發(fā)展,到2014年已經(jīng)出現(xiàn)持續(xù)超過一個(gè)月的400 ppm以上數(shù)據(jù)。近年來全球變化生態(tài)學(xué)家開始重視全球CO2濃度上升對(duì)地球生態(tài)系統(tǒng)的影響，尤其是溫室效應(yīng)和各種生態(tài)系統(tǒng)對(duì)氣候變化的響應(yīng)機(jī)理、機(jī)制。

學(xué)術(shù)界普遍認(rèn)為提高生態(tài)系統(tǒng)碳固定能力和碳儲(chǔ)存能力，是減緩大氣CO2含量增長(zhǎng)，抑制全球變暖的主要手段。因此，定量的測(cè)量、估計(jì)和預(yù)測(cè)生態(tài)系統(tǒng)碳儲(chǔ)量與碳固定量是生態(tài)系統(tǒng)研究的熱點(diǎn)和重點(diǎn)。由于生態(tài)系統(tǒng)具有高度的復(fù)雜性和不確定性，無論是人類活動(dòng)、土地利用變化、氣候變化還是生態(tài)系統(tǒng)自身的變化與演替都會(huì)改變生態(tài)系統(tǒng)中碳循環(huán)的路徑和通量，僅根據(jù)點(diǎn)尺度，離散時(shí)間段的監(jiān)測(cè)數(shù)據(jù)不能滿足當(dāng)前我們對(duì)生態(tài)系統(tǒng)碳循環(huán)在長(zhǎng)時(shí)間尺度，大空間尺度上動(dòng)態(tài)變化研究的需求。建立受氣候、土壤、生物和人類活動(dòng)綜合影響的生態(tài)系統(tǒng)模型不僅有助于客觀認(rèn)識(shí)我國(guó)陸地生態(tài)系統(tǒng)碳循環(huán)的過去、現(xiàn)在和未來，而且可以通過碳循環(huán)模型與氣候模式的雙向耦合對(duì)未來氣候變化做出更為客觀的估計(jì)，以幫助人類制定適應(yīng)氣候變化的措施[1]。

生物地球化學(xué)模型是研究物質(zhì)循環(huán)和追蹤元素運(yùn)動(dòng)軌跡很實(shí)用的模型，生物地化模型可以模擬生態(tài)系統(tǒng)中碳、氮、水、熱量在不同界面中的交換和變化過程。通過建立土壤、植被、大氣等不同界面的元素與水熱庫(kù)，用生理和生態(tài)過程公式來表達(dá)不同庫(kù)之間的動(dòng)態(tài)模擬諸如：光合、呼吸、生長(zhǎng)、凋落、降水、下滲、蒸發(fā)、固氮、分解和礦化等生物化學(xué)和地球化學(xué)過程。生物地球化學(xué)模型一般以日或者月為步進(jìn)，對(duì)生態(tài)系統(tǒng)的所有過程參量進(jìn)行全面模擬，然后對(duì)土壤有機(jī)碳、含水量、植物生物量、CO2濃度等狀態(tài)量進(jìn)行逐日積分得到長(zhǎng)時(shí)間尺度的模擬結(jié)果。同時(shí)結(jié)合大面積高頻采樣數(shù)據(jù)可以在空間上對(duì)生態(tài)系統(tǒng)多指標(biāo)模擬，通過積分得到區(qū)域尺度的生態(tài)系統(tǒng)狀態(tài)估計(jì)。相比較其他生態(tài)系統(tǒng)模型，生物地球化學(xué)模型比統(tǒng)計(jì)模型可以考慮更多的生態(tài)學(xué)過程，生物地球化學(xué)循環(huán)模型要求輸入氣候變量(降水、氣溫等)，植被類型(或者以植被模型作為基礎(chǔ))，有的還需要輸入遙感信息(如光合有效輻射等)，具有輸出生態(tài)系統(tǒng)的碳收支動(dòng)態(tài)的能力。生物地球化學(xué)模型重點(diǎn)考慮植被通過光合作用從大氣中固定碳后，碳在植被-土壤中的傳輸過程，包括自養(yǎng)呼吸、異養(yǎng)呼吸、光合同化產(chǎn)物在各植物器官間的分配、植物死亡、凋落、凋落物分解呼吸、土壤有機(jī)物呼吸和氮礦化等，并模擬土壤有效氮、土壤水分和大氣CO2濃度等對(duì)其的影響。因此是進(jìn)行生態(tài)系統(tǒng)水平碳收支研究的比較適合的一類模型。

1 DNDC模型簡(jiǎn)介

全球目前有30多種不同的生物地球化學(xué)過程模型[2]。DNDC是目前開發(fā)歷史較長(zhǎng)、使用較多，擴(kuò)展較為豐富，尤其是在我國(guó)大量驗(yàn)正并且廣泛應(yīng)用的一種生物地球化學(xué)過程模型。DNDC模型創(chuàng)建于1992年[3]，最初用于模擬美國(guó)農(nóng)田生態(tài)系統(tǒng)，模型以降水為主要驅(qū)動(dòng)力，主要模擬N2O、CO2、N2在農(nóng)田中的釋放。DNDC整合了土壤環(huán)境中碳氮元素的氧化還原過程，通過大量的實(shí)驗(yàn)數(shù)據(jù)和前人研究構(gòu)建了碳氮元素與土壤氣候環(huán)境的關(guān)系，模型主要追蹤碳素和氮素在土壤中的生物地球化學(xué)過程。過去二十多年中DNDC模型被全球各地的生態(tài)學(xué)者廣泛的應(yīng)用，模型有了很大的改進(jìn)和發(fā)展，加入了許多新的子模塊和更加詳實(shí)的生態(tài)過程參數(shù)，以適應(yīng)不同生態(tài)系統(tǒng)和不同關(guān)注點(diǎn)的研究。

2 DNDC模型發(fā)展與國(guó)內(nèi)學(xué)者在碳循環(huán)領(lǐng)域的應(yīng)用

v7.0之前早期版本的DNDC模型包含3個(gè)基礎(chǔ)模塊：土壤氣候/水熱通量模塊、土壤碳素分解模塊以及土壤脫氮模塊。模型并沒有包含植物生長(zhǎng)模塊，硝化作用使用單一公式近似模擬。

1994年DNDC v7.1 版本加入了一個(gè)以經(jīng)驗(yàn)公式為主的植物生長(zhǎng)模塊，這個(gè)模塊包含了常見的土壤耕作模式和田間管理措施，從而使得DNDC模型可以用于農(nóng)田土壤碳元素的生物地球化學(xué)過程模擬[4]。

2000年基于森林模擬的PnET-N-DNDC模型開發(fā)成功，這個(gè)模型的主要特色是包含一個(gè)光合-蒸散模塊和一個(gè)硝化作用模塊，后者可以用于精確的模擬硝化細(xì)菌的繁殖/死亡速率、亞硝化細(xì)菌的繁殖/死亡速率和整個(gè)硝化作用導(dǎo)致的NO、N2O排放量[5]。這些新的模塊和特性后來被整合到了DNDC v8.0中，PnET-N-DNDC中的厭氧氣球概念，將土壤環(huán)境區(qū)分為有氧環(huán)境和厭氧環(huán)境，加上土壤空隙度的引入，使得模型可以同時(shí)模擬硝化作用和反硝化作用。

2002年，加拿大遙感中心的Zhang等人開發(fā)出crop-DNDC，在DNDC v8.0的基礎(chǔ)上加入了3個(gè)獨(dú)立的植物生長(zhǎng)相關(guān)模塊，3個(gè)模塊包括：土壤氣候與水熱通量模塊、植物生理參數(shù)模塊、分解與硝化模塊。新的模塊使得模型具有追蹤植物生長(zhǎng)過程中碳、氮和水分在系統(tǒng)中的循環(huán)和流動(dòng)的能力[6]。從此模型可以模擬植物生理指數(shù)、葉面積指數(shù)、光合強(qiáng)度、呼吸強(qiáng)度、產(chǎn)物分配、根系生長(zhǎng)和氮元素的吸收。隨后Crop-DNDC的新模塊都被合并入了DNDC v8.2版本，從此DNDC模型的功能性進(jìn)入了一個(gè)新的領(lǐng)域：預(yù)測(cè)作物產(chǎn)量、固碳能力、溫室氣體排放和模擬氣候變化對(duì)農(nóng)業(yè)生產(chǎn)影響。在國(guó)內(nèi)邱建軍、唐華俊[7]最早引入了DNDCv8.2版本用于估算北方農(nóng)牧交錯(cuò)帶內(nèi)蒙古自治區(qū)耕地土壤碳儲(chǔ)量并研究其平衡狀況，結(jié)果表明：該區(qū)耕地0～30 cm土層土壤有機(jī)碳儲(chǔ)量約為46 622.6萬噸，目前土壤有機(jī)碳庫(kù)處于嚴(yán)重負(fù)平衡，年凈損失有機(jī)碳1 335萬噸。指出干旱半干旱農(nóng)牧交錯(cuò)區(qū)應(yīng)推行免耕、休耕等保護(hù)性耕作措施以助于土壤有機(jī)碳的積累。王立剛、邱建軍[8]以高產(chǎn)糧區(qū)河北省曲周縣為例，運(yùn)用DNDC模型，在對(duì)實(shí)測(cè)田間試驗(yàn)數(shù)據(jù)進(jìn)行模型驗(yàn)證的基礎(chǔ)上，對(duì)該地區(qū)農(nóng)業(yè)生態(tài)系統(tǒng)土壤碳氮平衡進(jìn)行了模擬研究。結(jié)果表明：1998年曲周縣47 821 hm2農(nóng)業(yè)耕地土壤的總有機(jī)碳儲(chǔ)量為742.94×106kg，平均每公頃耕地土壤有機(jī)碳儲(chǔ)量為15 536.05 kg。1998年耕地土壤有機(jī)碳為正平衡，1990年為負(fù)平衡，2年土壤有機(jī)碳顯著差異的主要原因在于秸稈還田比例的多少。

同年Zhang等人從PnET-DNDC模型繼續(xù)開發(fā)了wet-land DNDC模塊[9]，通過加入土壤水文參數(shù)、土壤溫度參數(shù)、植物生長(zhǎng)狀態(tài)參數(shù)和土壤碳軌跡實(shí)現(xiàn)了對(duì)生態(tài)系統(tǒng)CH4和CO2排放的模擬。其中最重要的一個(gè)核心意義就是能斯特方程和米氏方程被加入了厭氧氣球中。隨著能斯特方程和米氏方程的加入，DNDC模型具有了對(duì)土壤氧化還原電位和底物酶促反應(yīng)產(chǎn)物模擬的能力。

同年Brown等人經(jīng)過5年的開發(fā)完成了UK-DNDC模型[10]。UK-DNDC模型加入了一個(gè)空間化的有機(jī)肥模塊，可以支持耕作參數(shù)、家畜參數(shù)和采食參數(shù)輸入模型。

2004年，DNDC模型更新到v8.5版本[11]，從wet-land DNDC中吸取了能斯特方程和米氏方程，改進(jìn)了原有的厭氧模擬氣球。張凡等[12]以陜西農(nóng)業(yè)生態(tài)系統(tǒng)為對(duì)象，整合DNDC v8.5與陜西農(nóng)業(yè)地理信息系統(tǒng)數(shù)據(jù)庫(kù)，利用陜西地區(qū)氣象和作物資料，對(duì)陜西省2000年作物生長(zhǎng)發(fā)育和土壤碳循環(huán)進(jìn)行了模型模擬，實(shí)例探討了耕作管理對(duì)土壤碳儲(chǔ)量的影響。分析模擬結(jié)果可以得出3點(diǎn)結(jié)論：①2000年陜西耕作土壤總有機(jī)碳儲(chǔ)量約為103 Tg碳，是一個(gè)大氣CO2源，向大氣釋放0.5 Tg碳；②提高作物秸稈還田率是提高陜西農(nóng)田碳庫(kù)儲(chǔ)量有效可行的措施，將作物秸稈還田率從當(dāng)前的15%提高到50%或90%會(huì)使陜西農(nóng)田土壤從大氣CO2源轉(zhuǎn)變?yōu)閰R，每年分別增加土壤有機(jī)碳庫(kù)儲(chǔ)量0.7 Tg碳或2.1 Tg碳；③施用有機(jī)農(nóng)肥(500 kg·hm-2)也會(huì)增加土壤碳輸入，從而提高土壤碳儲(chǔ)量，使陜西農(nóng)田系統(tǒng)轉(zhuǎn)變?yōu)檩^弱的碳源，每年可多固定0.2 Tg碳。

2004年新西蘭的Saggar等人根據(jù)新西蘭的氣候土壤情況對(duì)DNDC做了大量本地化適配，開發(fā)出NZ-DNDC模型[13]，NZ-DNDC模型根據(jù)新西蘭的高有機(jī)碳、長(zhǎng)期放牧草地生態(tài)環(huán)境對(duì)模型進(jìn)行了機(jī)理的優(yōu)化，使之更加符合新西蘭的草地農(nóng)業(yè)生態(tài)系統(tǒng)。

2005年，李長(zhǎng)生等人開發(fā)了Forest-DNDC模型[14]，F(xiàn)orest-DNDC模型整合了PnET-DNDC模型和wet-land DNDC模型，加入了Aber等人開發(fā)的森林生理模塊，實(shí)現(xiàn)了對(duì)森林生物量、森林碳動(dòng)態(tài)、森林生態(tài)系統(tǒng)溫室氣體排放的模擬。模型中加入了大量的森林管理參數(shù)，并且可以模擬森林動(dòng)態(tài)例如：森林采伐、造林、間伐、火燒、灌溉、濕地恢復(fù)、施肥等過程。

同年Kiese等人從PnET-DNDC中進(jìn)一步開發(fā)了Forest-DNDC Tropica[15]，主要用于模擬熱帶雨林生態(tài)系統(tǒng)的N2O排放量。模型針對(duì)熱帶雨林的氣候條件和土壤環(huán)境，優(yōu)化了一些生態(tài)過程參數(shù)。在保留原有物質(zhì)循環(huán)框架的基礎(chǔ)上整合了BIOME-BGC模型的常綠闊葉林模塊。

德國(guó)的Neufeldt等人將EFEM模型和DNDC模型進(jìn)行了耦合[16]，為DNDC模型帶來了一個(gè)以GIS為基礎(chǔ)的區(qū)域模擬模塊，他們用EFEM-DNDC模型對(duì)德國(guó)巴登-符騰堡州的作物和家畜溫室氣體排放進(jìn)行了模擬，獲得了很好的精度。這個(gè)模型和其他DNDC的區(qū)別在于加入了總效益、影子價(jià)格和成本控制等經(jīng)濟(jì)學(xué)參數(shù)。

2006年，比利時(shí)的Beheydt等人根據(jù)DNDC v8.3版繼續(xù)開發(fā)了BE-DNDC模型[17]。BE-DNDC主要模擬大范圍農(nóng)田的N2O排放，其中區(qū)域化的草地農(nóng)田N2O排放量預(yù)測(cè)來源于當(dāng)?shù)卮罅渴占瘮?shù)據(jù)的回歸公式?；貧w公式根據(jù)比利時(shí)不同地區(qū)的22個(gè)長(zhǎng)期試驗(yàn)點(diǎn)數(shù)據(jù)進(jìn)行了優(yōu)化。BE-DNDC同時(shí)能根據(jù)結(jié)果不確定性來輸出土壤有機(jī)碳含量最大最小值。

同年，DNDC模型更新到v9.0版[18]，模型內(nèi)部詳細(xì)細(xì)化了自由氨基的動(dòng)態(tài)流程、硝化和亞硝化作用。使得模型對(duì)氮循環(huán)的模擬精度大幅提高。李虎等[19]以黃淮海平原河北省范圍內(nèi)的農(nóng)田土壤為研究對(duì)象，通過與田間實(shí)際觀測(cè)數(shù)據(jù)進(jìn)行比較發(fā)現(xiàn)，DNDC 9.0能夠較好地反映農(nóng)田土壤溫室氣體CO2和N2O的排放通量，可以用來模擬估算農(nóng)田土壤CO2和N2O的排放通量.根據(jù)模型估算，2003年河北省111個(gè)縣市農(nóng)業(yè)土壤CO2排放量約3.758×106t碳，各縣市總的N2O排放量40.345×106kg氮.全省釋放的CO2和N2O中有40%左右來自冬小麥/夏玉米地。邱建軍等[20]將DNDCv9.0與GIS數(shù)據(jù)庫(kù)進(jìn)行結(jié)合，估測(cè)了全國(guó)六大農(nóng)業(yè)區(qū)域土壤有機(jī)碳含量變化對(duì)農(nóng)作物產(chǎn)量的影響。研究認(rèn)為初始土壤有機(jī)碳每增加1 g碳每千克，東北地區(qū)的玉米每公頃產(chǎn)量將增加176 kg，北方地區(qū)的玉米小麥輪作田每公頃產(chǎn)量將增加454 kg，西北地區(qū)的玉米每公頃產(chǎn)量增加328 kg，中南地區(qū)的單播水稻每公頃將增產(chǎn)185 kg，東部地區(qū)的兩季稻每公頃將增產(chǎn)266 kg，西南地區(qū)的麥稻輪作田每公頃將增產(chǎn)229 kg。王淑芳等[21]運(yùn)用DNDCv9.0對(duì)密云水庫(kù)上游地區(qū)農(nóng)田土壤有機(jī)碳儲(chǔ)量及其變化進(jìn)行模擬研究，結(jié)果表明：2006年密云水庫(kù)上游地區(qū)214 920 hm2農(nóng)田土壤(0～25 cm)的總有機(jī)碳儲(chǔ)量為7 646×106kg，其中位于河北省境內(nèi)的該地區(qū) 63.1%的農(nóng)田儲(chǔ)存了全區(qū) 68.1%的 SOC；平均每公頃農(nóng)田SOC儲(chǔ)量為 35 576.1 kg，低于全國(guó)平均水平；由于化肥和有機(jī)肥投入的增加，經(jīng)一年耕種后，2006年該地區(qū)農(nóng)田土壤有機(jī)碳儲(chǔ)量增加 142.5×106kg，整個(gè)地區(qū)及各區(qū)縣農(nóng)田土壤碳收支均為正，是大氣 CO2的一個(gè)匯。馬玉芳等[22]利用甘肅農(nóng)業(yè)大學(xué)定西旱農(nóng)綜合試驗(yàn)站的2001-2008年長(zhǎng)期定位試驗(yàn)數(shù)據(jù)對(duì)DNDC9.0模型進(jìn)行了驗(yàn)證，結(jié)果表明：秸稈還田或覆蓋可提高作物秸稈與根系殘留的外源碳攜入量，也能提高土壤異氧呼吸對(duì)內(nèi)源碳的消耗。該區(qū)域2008年表層農(nóng)田土壤有機(jī)碳總儲(chǔ)量為2.8×109kg，平均土壤有機(jī)碳密度為2.33 kg·m-2。預(yù)測(cè)2008-2037年土壤有機(jī)碳儲(chǔ)量，得出增加有機(jī)肥能夠顯著增加土壤有機(jī)碳的積累，其次為免耕同時(shí)增加秸稈還田率50%，而單獨(dú)實(shí)施這兩種措施，土壤有機(jī)碳的積累速度較慢。陳晨等[23]在利用中科院安塞水土保持綜合實(shí)驗(yàn)站長(zhǎng)期定位數(shù)據(jù)，應(yīng)用DNDC9.0模型研究黃土丘陵溝壑區(qū)坡耕地農(nóng)田土壤碳庫(kù)儲(chǔ)量動(dòng)態(tài)變化及碳循環(huán)特征，結(jié)果表明：施肥可提高作物殘?bào)w與根系分泌物的外源碳攜入量，也能提高土壤異氧呼吸對(duì)內(nèi)源碳的消耗，且施用有機(jī)肥后提高效果顯著?？傮w上，單施有機(jī)肥、有機(jī)肥配施氮肥40年后各處理土壤有機(jī)碳分別比初始值提高了90.29%，86.46%。不施肥和單施氮肥 40年后各處理土壤有機(jī)碳分別比初始值降低了3.52%，0.38%。王立剛等[24]利用DNDC9.0，選擇6個(gè)典型的種植模式點(diǎn)，考察改變田間管理措施(秸稈還田、增施有機(jī)肥、改變化肥施用等)對(duì)土壤有機(jī)碳含量影響的長(zhǎng)期效應(yīng)，研究結(jié)果表明：增加秸稈還田比例的固碳效應(yīng)在眾多的農(nóng)田管理措施中體現(xiàn)最突出，以基本管理措施20年后有機(jī)碳含量為對(duì)照，當(dāng)秸稈還田率增加到80%時(shí)，齊齊哈爾、平?jīng)龊蛙平貐^(qū)增加的土壤有機(jī)碳含量都在3 000 kgC·hm-2以上，北京郊區(qū)、江寧和鹽亭地區(qū)每公頃可以增加碳10 000 kg以上?；适┯昧繛楸镜字档?0%配合秸稈還田率50%和化肥施用量為本底值的50%增施有機(jī)肥1 000 kgC·hm-2的兩項(xiàng)措施，在各個(gè)地區(qū)都表現(xiàn)出比基本管理措施能存儲(chǔ)更多的碳素，因此，在我國(guó)現(xiàn)在的管理措施下，減少化肥的同時(shí)增加秸稈還田比例和增施有機(jī)肥是既增加土壤有機(jī)碳含量又減少化肥污染的“雙贏”措施。陳長(zhǎng)青[25]等為探討長(zhǎng)期不同施肥條件下土壤有機(jī)碳的變化規(guī)律及DNDC模型的應(yīng)用，利用封丘石灰性潮土不同施肥措施下的長(zhǎng)期定位試驗(yàn)數(shù)據(jù)，選取CK、NPK、1/2NPK+ 1/2OM、OM 4個(gè)處理，分析了15年來作物產(chǎn)量及土壤有機(jī)碳的變化特征，并用DNDC模型預(yù)測(cè)了試驗(yàn)地近100年(2000～2099)的土壤有機(jī)碳變化趨勢(shì)，結(jié)果表明：有機(jī)肥和化肥配施能獲得較高的作物產(chǎn)量，并能在百年尺度上有效地增加土壤有機(jī)碳含量，從而提高土壤的可持續(xù)利用能力。

DNDC-Europe模型開發(fā)于2008年[26]，主要用途是模擬農(nóng)業(yè)生產(chǎn)、農(nóng)業(yè)經(jīng)濟(jì)效益與溫室氣體排放之間的關(guān)系。在整合了農(nóng)業(yè)經(jīng)濟(jì)模型CAPRI之后，Leip等人用DNDC-Europe模型進(jìn)行農(nóng)業(yè)市場(chǎng)、環(huán)境政策、交易模式、農(nóng)戶、經(jīng)銷商、消費(fèi)者行為和農(nóng)業(yè)生產(chǎn)及農(nóng)業(yè)生態(tài)系統(tǒng)變化的細(xì)微關(guān)系，并且提供了一個(gè)全球尺度政策對(duì)市場(chǎng)分配影響的預(yù)測(cè)。

DNDC模型對(duì)溫室氣體的模擬最早源于水稻田相關(guān)研究，2008年Fumoto等人基于李長(zhǎng)生等在2004年的工作進(jìn)一步發(fā)展出了DNDC-Rice模型[27]，主要用于模擬不同氣候、土壤和管理措施條件下水稻田CO2、無機(jī)碳和甲烷的排放， 這個(gè)模型在2010年進(jìn)一步改進(jìn)支持模擬水田不同水淹制度下的水稻田甲烷排放。

2009年，Grote等人將DNDC模型連接到MoBiLE (Modular Biosphere Simulation Environment)模型研究框架中，MoBiLE框架可以在研究生態(tài)系統(tǒng)某一參數(shù)或者某一方面時(shí)為模型提供高效的參數(shù)選擇、初始化和校正驗(yàn)證，并且DNDC在框架內(nèi)將被視為一個(gè)更大模型的子模塊，從而可以和更多的模型連接進(jìn)行多種不同用途的模擬[28]。2011年Grote等人在MoBiLE框架內(nèi)將PSIM模型連接入DNDC，細(xì)化了DNDC模型的植物生長(zhǎng)模塊，對(duì)森林生態(tài)系統(tǒng)碳通量和碳庫(kù)的細(xì)節(jié)進(jìn)行了詳細(xì)的研究[29]。Wolf等人將PnEF-N-DNDC連接入MoBiLE框架用于模擬溫帶生態(tài)系統(tǒng)融雪期的N2O排放，并且制作了最大雪深和生長(zhǎng)季末生物量之間的數(shù)學(xué)關(guān)系[30]。

2010年DNDC模型更新到v9.3版，改進(jìn)了土壤蒸發(fā)的計(jì)算方式和精度。馬玉芳等[31]應(yīng)用DNDC9.3對(duì)甘肅農(nóng)業(yè)大學(xué)定西旱農(nóng)綜合試驗(yàn)站的2001—2008年長(zhǎng)期定位試驗(yàn)的4個(gè)處理的長(zhǎng)期(100 年)模擬表明：免耕覆蓋秸稈(NTS)和傳統(tǒng)耕作+秸稈還田(TS)能夠較大幅度地提高土壤有機(jī)碳含量。而傳統(tǒng)耕作不覆蓋(T)處理下土壤有機(jī)碳呈下降趨勢(shì)。土壤有機(jī)碳含量對(duì)土壤性狀以及耕作管理措施變化的響應(yīng)模擬研究表明，秸稈還田和施用有機(jī)肥是最有效的提高土壤有機(jī)碳含量的耕作措施，而土壤性狀，尤其是初始有機(jī)碳含量是影響土壤有機(jī)碳變化的最主要敏感因素。孫園園等[32]利用DNDC9.3模型模擬川中丘陵區(qū)不同輪作制度下稻田CO2排放情況，結(jié)果表明：PF、RR和RW 3種輪作制度下RW和RR處理CO2排放量均顯著高于PF，但3種輪作制度下單季作物CO2排放量差異不大，RR處理的單季作物的CO2排放量最小，其年總作物產(chǎn)量居中，RW處理產(chǎn)量最高。根際呼吸是土壤呼吸的主要部分，整個(gè)生長(zhǎng)期PF，RR，RW的根呼吸貢獻(xiàn)率平均為59.14%～62.96%。徐勝祥等[33]以江蘇省稻田為對(duì)象，整合DNDC9.3模型和100∶1土壤數(shù)據(jù)庫(kù)，定量估算少耕、免耕以及綜合措施少耕+桔稈還田下江蘇省稻田土壤的固碳潛力，模擬結(jié)果表明相對(duì)于傳統(tǒng)耕作，采用少耕、免耕和少耕+秸稈還田均可明顯地增加稻田的積累，長(zhǎng)期模擬結(jié)果表明：免耕和少耕+秸稈還田下固碳能力是單純少耕的2倍左右。張明圓等[34]研究在DNDC9.3模型調(diào)整的基礎(chǔ)上選取了土壤有機(jī)碳(SOC)、土壤碳密度(SCD)、土壤呼吸(SR)以及生物量碳(BC)4個(gè)指標(biāo)對(duì)模型在華北麥-玉兩熟農(nóng)田的適用性進(jìn)行驗(yàn)證，并用該模型模擬當(dāng)?shù)赝寥捞純?chǔ)量(SCS)動(dòng)態(tài)變化以及溫室氣體排放特征，結(jié)果表明：2001-2010年SOC和SCS逐年遞增；對(duì)未來100年模擬發(fā)現(xiàn)，15 年旋耕(RT)和翻耕(CT)處理SOC增長(zhǎng)迅速，而免耕(NT)SOC的劇烈增長(zhǎng)趨勢(shì)要持續(xù)近40 年；對(duì)比各處理100年碳儲(chǔ)量變化可知，前20 年CT處理SCS最大，20年后NT處理SCS最大；各處理土壤全球變暖潛勢(shì)(GWP)大小為CT>RT>NT。Kang等[35]運(yùn)用DNDC9.3模擬分析了羊草生態(tài)系統(tǒng)2004-2005年凈CO2通量(NEE)及其組分(GPP、Re)的動(dòng)態(tài)變化過程。通過與渦相關(guān)實(shí)測(cè)數(shù)據(jù)對(duì)比發(fā)現(xiàn)羊草草原由2004年生長(zhǎng)季(濕潤(rùn)年)的碳匯轉(zhuǎn)變?yōu)?005年(極端干早年)碳源的過程。2005年的極端干旱氣候嚴(yán)重影響草原生產(chǎn)力，進(jìn)而減少了土壤的凋落物輸入量，最終導(dǎo)致生態(tài)系統(tǒng)碳的凈損失。進(jìn)一步的28年圍封和放牧情景模擬發(fā)現(xiàn)：與圍欄封育相比，中度放牧一方面提高了植物的光合作用速率，顯著增加了系統(tǒng)的固碳能力，與2006-2007年放牧樣地實(shí)測(cè)結(jié)果一致；另一方面主要影響土壤異養(yǎng)呼吸，對(duì)植物自養(yǎng)呼吸影響不大。氣候變化和放牧管理措施將對(duì)內(nèi)蒙古羊草草原碳收支產(chǎn)生顯著的影響。王德營(yíng)等[36]探討黑土有機(jī)碳的長(zhǎng)期變化規(guī)律及DNDC模型在土壤有機(jī)碳預(yù)測(cè)方面的適用性，結(jié)果表明：有機(jī)肥的施用對(duì)土壤有機(jī)碳含量的影響最顯著，且這種影響具有持久性。模擬4種施肥情境下未來100年的土壤有機(jī)碳變化情況表明：不施肥處理土壤有機(jī)碳含量略有下降;單施化肥處理(NPK)土壤有機(jī)碳含量較為穩(wěn)定，并未出現(xiàn)土壤有機(jī)碳含量下降;配施有機(jī)肥(NPKM)和增施有機(jī)肥(M2+NPK)處理土壤有機(jī)碳含量增加明顯。

Kr?bel等[37]于2011年給DNDC模型外掛了一個(gè)CSW子模塊，用于高精度模擬加拿大春小麥生長(zhǎng)和小麥農(nóng)田生態(tài)系統(tǒng)氮固定。目前這一子模塊在國(guó)內(nèi)尚無人引進(jìn)。

Haas等人于2013年將DNDC和Forest-DNDC在MoBiLE框架中合并到一個(gè)綜合的土壤生物地化模型中，稱為L(zhǎng)andscape DNDC[38]，主要用于模擬生態(tài)系統(tǒng)中的碳氮元素周轉(zhuǎn)以及不同土地利用方式和不同土地利用變化下土壤碳氮儲(chǔ)量的變化。Landscape DNDC主要加入了一個(gè)基于DNDC結(jié)構(gòu)的植物生長(zhǎng)生理模塊，允許土地的利用方式和作物、森林和草原在短時(shí)間內(nèi)切換。該模塊基于物理和化學(xué)基本公式驅(qū)動(dòng)，描述了土壤環(huán)境中土壤化學(xué)驅(qū)動(dòng)下的微生物碳氮周轉(zhuǎn)過程和植被生長(zhǎng)動(dòng)態(tài)。模型可以同時(shí)支持三維區(qū)域尺度和點(diǎn)尺度的模擬。

為了研究永久凍土中的CH4排放，Zhang等人于2012年開發(fā)了NEST-DNDC模型[39]，模型模擬了植物群落尺度中的生物物理和生物化學(xué)過程，并且可以根據(jù)地域特征擴(kuò)展到生態(tài)系統(tǒng)尺度。NEST-DNDC模型整合了永凍土模型Northern Ecosystem Soil Temperature，使得模型可以詳細(xì)模擬永凍土的水熱條件和生物地球化學(xué)過程。借助預(yù)設(shè)的多種土壤質(zhì)地和砂石分層數(shù)據(jù)，模型也能模擬非永凍土生態(tài)系統(tǒng)例如森林、苔原等生境。

為了滿足對(duì)家畜生產(chǎn)中的溫室氣體和氨氣排放的模擬，李長(zhǎng)生等[40]于2012年開發(fā)了Manure-DNDC模型。Manure-DNDC將一系列的生物地球化學(xué)過程整合到一個(gè)可計(jì)算框架中，來反映有機(jī)肥在農(nóng)場(chǎng)中的周轉(zhuǎn)周期，以便預(yù)測(cè)家畜有機(jī)肥處理系統(tǒng)中的溫室氣體和氨氣釋放。鑒于有機(jī)肥中有機(jī)物的特性與土壤有機(jī)質(zhì)類似，Manure-DNDC直接借鑒了DNDC中的土壤有機(jī)質(zhì)動(dòng)態(tài)模塊。環(huán)境因子和生物地化反應(yīng)用于判斷生產(chǎn)中各環(huán)節(jié)的CO2、N2O、CH4和NH3釋放量。飼養(yǎng)場(chǎng)、堆肥場(chǎng)、污水池、厭氧消化池以及整個(gè)養(yǎng)殖場(chǎng)都直接和有機(jī)肥動(dòng)態(tài)模塊關(guān)聯(lián)，每一個(gè)環(huán)節(jié)都可以獨(dú)立或者整體計(jì)算有機(jī)肥的存儲(chǔ)發(fā)酵以及氣體排放，正如DNDC模型可以模擬多種耕作土壤一樣，Manure-DNDC可以模擬多種不同架構(gòu)的養(yǎng)殖場(chǎng)。

2012年，DNDC更新到v9.4版，將Manure-DNDC中的土壤銨根離子動(dòng)態(tài)模塊整合入模型[40]。陳杰華等[41]分析了15年來土壤有機(jī)碳的變化特征，并用DNDC模型模擬了長(zhǎng)期試驗(yàn)點(diǎn)15年和100年有機(jī)碳變化趨勢(shì)，結(jié)果表明：經(jīng)過15年的耕種，4種耕作制度的土壤有機(jī)碳含量都較原始有機(jī)碳含量有所上升。T2處理的土壤有機(jī)碳增幅最大，分別為CK，T1和T3處理的131%，142%，153%。T3處理的有機(jī)碳含量最低。100年尺度下四種耕作制度土壤有機(jī)碳含量都會(huì)上升。薛靜等[42]基于河北省曲周試驗(yàn)站的試驗(yàn)數(shù)據(jù)，對(duì)DNDC模型的參數(shù)進(jìn)行詳盡的靈敏度分析和標(biāo)定，經(jīng)過詳細(xì)參數(shù)標(biāo)定后的DNDC模型將為深入探討不同農(nóng)田管理措施下，實(shí)現(xiàn)作物增產(chǎn)和農(nóng)業(yè)減排雙贏局面提供科學(xué)有效的途徑。倪元龍等[43]在1982 年水稻田土壤有機(jī)碳庫(kù)的基礎(chǔ)上，利用 DNDC9.4模型結(jié)合3種比例尺矢量圖斑單元和不同分辨率的柵格單元模擬了2000年太湖地區(qū)水稻田表層土壤有機(jī)碳庫(kù)，對(duì)比分析了不同分辨率柵格單元與3種比例尺矢量圖斑單元的水稻土類型數(shù)量、面積、表層土壤有機(jī)碳儲(chǔ)量以及有機(jī)碳密度的變化特征，并以矢量圖斑單元獲得的這4個(gè)指標(biāo)結(jié)果為基準(zhǔn)，用相對(duì)變異百分?jǐn)?shù)(VIV)來判別基于3種比例尺DNDC模型模擬的最佳柵格單元分辨率，結(jié)果表明：最佳柵格模擬單元分辨率分別為0.2 km×0.2 km，2 km×2 km，17 km×17 km，既能保證模型模擬過程中的精度要求，又可以避免數(shù)據(jù)冗余，提高模擬效率。建立的土壤碳庫(kù)模擬研究的比例尺與其最佳柵格單元分辨率對(duì)應(yīng)轉(zhuǎn)換關(guān)系，對(duì)區(qū)域土壤碳模擬研究具有重要參考價(jià)值。

2014年，DNDC更新到v9.5版，增強(qiáng)了植物生長(zhǎng)模塊和水文動(dòng)態(tài)模塊，并且針對(duì)溫室氣體排放進(jìn)行了相關(guān)參數(shù)的優(yōu)化[44]。韓東亮等[45]選擇國(guó)家灰漠土長(zhǎng)期肥力監(jiān)測(cè)試驗(yàn)中的5種處理(hNPKM，NPKM，NPKS、NPK、CK)，采用DNDC9.5模型對(duì)土壤表層土壤有機(jī)碳150年(2010-2160年)的長(zhǎng)期變化趨勢(shì)進(jìn)行預(yù)測(cè)，結(jié)果表明：化肥與有機(jī)肥配施土壤有機(jī)碳的固碳速率最高，化肥平衡施肥、化肥與秸稈配施也具有固碳效果，但固碳速率小于化肥有機(jī)肥配施，而不施肥土壤有機(jī)碳逐漸減少，因此合理施用有機(jī)肥和化肥能夠?qū)崿F(xiàn)作物增產(chǎn)和土壤培肥的雙贏。韓東亮等[46]以國(guó)家灰漠土肥力與肥料效益觀測(cè)試驗(yàn)站22年長(zhǎng)期定位試驗(yàn)監(jiān)測(cè)數(shù)據(jù)為基礎(chǔ)，利用DNDC模型對(duì)9個(gè)不同施肥處理的耕層(0～20 cm)土壤有機(jī)碳(SOC)演變進(jìn)行模擬和驗(yàn)證，并預(yù)測(cè)了未來30年(2010-2040)SOC的變化趨勢(shì)；同時(shí)也對(duì)CO2的釋放量和變化速率進(jìn)行了模擬，結(jié)果表明：持續(xù)施用氮磷鉀肥(NPK)不能使SOC含量提高；而氮磷鉀肥配施有機(jī)肥(NPKM)或秸稈(NPKS)可使SOC含量比2010年增加；雖然土壤有機(jī)碳升高越快，CO2的排放量也越多，但有機(jī)碳變化率超過25.1%時(shí)則表現(xiàn)為CO2的固定。楊黎等[47]在對(duì)DNDC模型進(jìn)行驗(yàn)證的基礎(chǔ)上，應(yīng)用模型研究不同施氮和秸稈還田措施對(duì)東北地區(qū)春玉米農(nóng)田固碳和氧化亞氮排放的長(zhǎng)期綜合影響。陳海心等[48]利用DNDC9.5模型模擬了陜西楊凌地區(qū)不施肥(CK)、單施氮磷鉀無機(jī)肥(NPK)、無機(jī)肥配施秸稈(SNPK)、無機(jī)肥配施低量有機(jī)肥(M1NPK)及無機(jī)肥配施高量有機(jī)肥(M2NPK)5個(gè)處理的定位試驗(yàn)對(duì)土壤有機(jī)碳變化的影響。從模型輸出20年來土壤碳含量的變化看，施肥能提高作物地上部、地下部生物量以及土壤異氧呼吸作用，即同時(shí)提高了外源碳攜入量和內(nèi)源碳消耗量，但各施肥處理均在碳儲(chǔ)量上表現(xiàn)為盈余狀態(tài)。

3 結(jié)語(yǔ)

DNDC模型經(jīng)過25年的發(fā)展，從最初的農(nóng)田N2O排放模擬模型，逐步開發(fā)和發(fā)展到涵蓋森林、濕地、草原等多種生態(tài)系統(tǒng)類型模擬，并且可以模擬系統(tǒng)中碳氮元素運(yùn)動(dòng)變化軌跡等多種生態(tài)學(xué)過程，經(jīng)過多點(diǎn)的驗(yàn)證與校正，DNDC模型可以在多種生態(tài)系統(tǒng)中闡明物質(zhì)循環(huán)的過程，同時(shí)可以模擬3種主要溫室氣體的排放量，可以很方便的進(jìn)行溫室氣體效應(yīng)的模擬和估算。在氮元素的追蹤和模擬方面有著大量成熟的模擬案例和溫室氣體估算，在碳元素動(dòng)態(tài)的模擬和分析方面DNDC也出現(xiàn)越來越多的成果。

過去DNDC模型在碳循環(huán)模擬中應(yīng)用的研究在數(shù)量和尺度上都不如在氮素模擬中廣泛，70%的DNDC相關(guān)文獻(xiàn)都是用于模擬氮素循環(huán)或者溫室氣體排放。這一方面是因?yàn)镈NDC最早從農(nóng)田溫室氣體排放模擬起源，因此有著大量的相關(guān)用戶和驗(yàn)證經(jīng)驗(yàn)；另一方面，碳元素運(yùn)動(dòng)和變化在時(shí)間空間上的連續(xù)性較強(qiáng)，模型對(duì)碳循環(huán)過程模擬，對(duì)驗(yàn)證數(shù)據(jù)的采樣頻率，采樣范圍，采樣精度要求較高。傳統(tǒng)的箱式法測(cè)土壤呼吸與生態(tài)系統(tǒng)呼吸，誤差較大而且時(shí)相頻率無法做到較高頻的連續(xù)監(jiān)測(cè)。近年來，隨著Li-co 8150野外連續(xù)自動(dòng)檢測(cè)儀的普及，渦度相關(guān)通量技術(shù)的廣泛應(yīng)用甚至大孔徑閃爍儀的建立，長(zhǎng)時(shí)間較大范圍的碳動(dòng)態(tài)監(jiān)測(cè)和高質(zhì)量低誤差數(shù)據(jù)的獲取成為了現(xiàn)實(shí)。結(jié)合DNDC模型的幾大優(yōu)勢(shì)：多指標(biāo)全面的生態(tài)系統(tǒng)模擬，易于獲得的初始化狀態(tài)參數(shù)，大量根據(jù)多種生態(tài)類型驗(yàn)證的預(yù)設(shè)參數(shù)，詳細(xì)精準(zhǔn)的土壤氧化還原點(diǎn)電位反應(yīng)機(jī)理，基于能斯特方程的土壤化學(xué)平衡，支持模擬多種農(nóng)田管理方式對(duì)系統(tǒng)的影響，日尺度步進(jìn)的模擬過程和方便實(shí)用的蒙特卡洛不確定性模擬，使得DNDC模型成為了未來進(jìn)行生態(tài)系統(tǒng)碳循環(huán)模擬與研究又一個(gè)很好的工具。從模型輸出的生態(tài)系統(tǒng)各組分參數(shù)對(duì)環(huán)境變化的響應(yīng)和反饋中，研究者可以更好地理解土壤化學(xué)過程的機(jī)理與環(huán)境因子的互作，還可以評(píng)估不同生態(tài)系統(tǒng)管理和不同土地利用方式下溫室氣體的排放，土壤有機(jī)碳的流失/積累。DNDC模型對(duì)研究生物地球化學(xué)循環(huán)和生態(tài)系統(tǒng)可持續(xù)發(fā)展具有積極的指導(dǎo)作用。