王 茜，查元源

（武漢大學(xué)水資源工程與調(diào)度全國重點實驗室，武漢 430072）

0 引 言

農(nóng)田生態(tài)系統(tǒng)與人類生命、生活、生產(chǎn)息息相關(guān)，它為我們提供賴以生存的糧食，同時也消耗著全球淡水抽取量的70%和總耗水量的90%[1]。作為水-糧食-能源三者相互聯(lián)系的復(fù)雜系統(tǒng)的縮影，準(zhǔn)確刻畫農(nóng)田生態(tài)系統(tǒng)中涉及的水、能量及CO2通量變化過程，為區(qū)域水資源配置、糧食生產(chǎn)、生態(tài)環(huán)境保護(hù)提供理論依據(jù)，對保障國家水安全、糧食安全、生態(tài)安全有著重要意義。陸面模式作為模擬陸面上發(fā)生的各種物理、化學(xué)、生物和水文過程，以及這些過程與大氣的相互作用的有效工具，歷經(jīng)簡單模式階段與生物大氣模式階段，進(jìn)入了新一代多模式階段[2]。從簡單地考慮水量平衡，到將陸氣間水分與熱量的交換過程進(jìn)行耦合；從不考慮植被的作用，到顯式引入植被的生物物理作用，再到通過耦合葉片的光合作用與蒸騰作用實現(xiàn)陸氣交換過程的水碳耦合，模式的物理機制逐漸完善，對陸面過程與陸氣交換的刻畫也逐漸細(xì)化。新一代陸面模式中，Noah-MP[3]模型集合多種參數(shù)化方案選項，提供多達(dá)4 608 種參數(shù)化方案組合，方便使用者依據(jù)適用性靈活選擇，具有廣闊的應(yīng)用前景。

隨著陸面模式的不斷完善，參數(shù)不確定性成為了限制模型模擬性能的主要因素。進(jìn)入大數(shù)據(jù)時代，衛(wèi)星、無人機等多種遙感技術(shù)快速發(fā)展，使得基于天-空-地多源觀測數(shù)據(jù)建立高分辨率的土壤屬性與植被屬性等地表參數(shù)庫成為可能[4]。因此，對陸面模式進(jìn)行參數(shù)敏感性分析，識別出影響模型模擬效果的敏感參數(shù)，不僅能夠指導(dǎo)參數(shù)優(yōu)化以提升模型模擬性能，也能為地表參數(shù)庫的構(gòu)建指引方向。

近年來，許多學(xué)者展開了Noah-MP 模型的參數(shù)敏感性分析研究。在密西西比河流域評估Noah-MP 陸面模式模擬水文變量的表現(xiàn)時，CAI 等[5]簡單地對3 個徑流模擬關(guān)鍵參數(shù)進(jìn)行了敏感性測試，并在此基礎(chǔ)上以納什效率系數(shù)（NSE）為性能指標(biāo)進(jìn)行參數(shù)優(yōu)化，從而提升模型徑流模擬能力。ARSENAULT 等[6]利用10 個國際FluxNet 站點的數(shù)據(jù)對動態(tài)植被參數(shù)、靜態(tài)植被參數(shù)及土壤參數(shù)進(jìn)行了全局敏感性分析，找出了12 個敏感參數(shù)，并指出針對特定區(qū)域的具體問題需要對參數(shù)分布范圍進(jìn)行進(jìn)一步約束。而HUO等[7]重點關(guān)注草及常綠闊葉林2種典型植被，提出了一種空間抽樣與參數(shù)抽樣相結(jié)合的新方法，在美國大陸的網(wǎng)格單元中篩選出植被總初級生產(chǎn)力（GPP）和潛熱通量（LE）模擬的敏感參數(shù)。SHU 等[8]在黃土高原的研究顯示，相比于直接同化葉面積指數(shù)（LAI），基于參數(shù)敏感性分析的參數(shù)優(yōu)化能更大幅度地降低葉面積指數(shù)的模擬值與觀測值的均方根誤差（RMSE）。上述研究均表明，受不同的下墊面及氣候區(qū)的影響，需開展針對性的Noah-MP 模型參數(shù)敏感性分析試驗，以助挖掘參數(shù)優(yōu)化對模型模擬性能提升的巨大潛力。目前，已有大量有關(guān)Noah-MP 模型在華北平原農(nóng)田生態(tài)系統(tǒng)的應(yīng)用型研究[9-11]，但在該區(qū)域開展的參數(shù)敏感性分析研究未見報道。華北平原是確保全國糧食安全的重要基地[12]，基于Noah-MP 模型開展參數(shù)敏感性分析，有助于準(zhǔn)確量化這一典型農(nóng)業(yè)生態(tài)系統(tǒng)模式與大氣間的輻射傳輸、物質(zhì)和能量交換，對推動區(qū)域經(jīng)濟(jì)、社會和生態(tài)可持續(xù)性協(xié)調(diào)發(fā)展有著促進(jìn)作用。

1 數(shù)據(jù)與方法

1.1 觀測站點與數(shù)據(jù)

中國國家生態(tài)系統(tǒng)觀測研究網(wǎng)絡(luò)（CERN）的禹城農(nóng)業(yè)綜合試驗站（116.570 2°E，36.829 0°N，平均海拔22.0 m）與欒城農(nóng)業(yè)生態(tài)系統(tǒng)試驗站（114.699 3°E，37.893 0°N，平均海拔50.1 m）均位于華北平原，屬暖溫帶半濕潤季風(fēng)氣候區(qū)，以冬小麥-夏玉米為主要作物，一年2 季?？紤]到華北平原輪作特點，即每年10月至次年6月為冬小麥生長季，6月至10月為夏玉米生長季，結(jié)合數(shù)據(jù)質(zhì)量與可獲取性等原因，本研究以2020年6月10日至2021年6月9日的輪作年為研究時段。獲取的觀測數(shù)據(jù)包括基本氣象數(shù)據(jù)、多層土壤水（0.1、0.2、0.3、0.4和0.5 m處）及水熱碳通量。研究時段內(nèi)的數(shù)據(jù)存在大量缺失，需要進(jìn)行插補處理。插補過程參考R 語言的fluxnetLSM[13]程序包，該程序包主要用于插補遵循Fluxnet 2015 協(xié)議的數(shù)據(jù)。處理后的氣象數(shù)據(jù)仍存在大量連續(xù)缺省值難以作為模型所需輸入數(shù)據(jù)，故不加以利用，而土壤水與通量觀測數(shù)據(jù)缺失率低于15%，可以用于模型敏感性分析檢驗。

1.2 Noah-MP模型

Noah-MP 模型（The Community Noah Land Surface Model with Multi-Parameterization Options）是在Noah 3.0 版本的基礎(chǔ)上進(jìn)行改進(jìn)和完善的新一代陸面模式。Noah 模型將陸面簡單視為0～0.1、0.1～0.4、0.4～1.0 和1.0～2.0 m 的4 層土壤，能量平衡計算在土壤與大氣的交界面展開，Noah-MP 模型在此基礎(chǔ)上調(diào)整模型架構(gòu)，將冠層與地表分開，引入3 層雪層，同時對積雪/融雪、產(chǎn)匯流、植被生長等關(guān)鍵過程進(jìn)行優(yōu)化，并為表1 所示的12 個陸面過程提供多種參數(shù)化方案選項。Noah-MP 模型能很好地模擬地表能量收支平衡及陸地水文循環(huán)過程，已被作為陸面參數(shù)化方案嵌入GFS、CFS 及WRF 等氣候/天氣預(yù)測模式中。為簡化輸入輸出文件讀取過程，修改網(wǎng)格化模型的輸入輸出端口，得到用于站點敏感性分析的單點模型（https://github.com/wangxi1998/NoahMP-Sensitivity-NCP/tree/main/POINT_NOAHMP）。

表1 模型參數(shù)化方案設(shè)置Tab.1 Model parameterization scheme

1.3 Noah-MP模型參數(shù)敏感性分析設(shè)計

1.3.1 模型模擬試驗的設(shè)計

模型所需的氣象數(shù)據(jù)來自GLDAS-2.1[14,15]，空間分辨率為0.25°，時間分辨率為3 h。該數(shù)據(jù)為Noah-MP 等陸面模式的“在線”模擬集合，故常用作Noah-MP 模型的氣象輸入數(shù)據(jù)。利用GEE 平臺提取2 個站點在研究時段的數(shù)據(jù)（2 m 氣溫、2 m 比濕、10 m 風(fēng)速、氣壓、下行短波輻射、下行長波輻射及降水），將其插值得到0.5 h分辨率的數(shù)據(jù)作為輸入。模型初始場參考模型指導(dǎo)手冊確定，分層土壤溫度與地表溫度取全球均值（274.0、278.0、280.0、284.0 和287.0 K），土壤含水率為飽和含水率的40%（取為0.24 m3/m3）。

選用的參數(shù)化方案及參數(shù)參考ARSENAULT 等[6]和HUO等[7]的研究（見表1 和表2）。選取的42 個參數(shù)包含30 個植被相關(guān)參數(shù)與12 個土壤相關(guān)參數(shù)；植被相關(guān)參數(shù)中有12 個與輻射傳輸過程密切相關(guān)，其余參數(shù)為動態(tài)植被模擬模塊涉及參數(shù)。模型模擬過程中，模型將土地利用類型與土壤類型作為索引條件，再在模型自帶的查找表中匹配對應(yīng)參數(shù)。因此，大部分參數(shù)的取值上下限通過查找表確定。少量參數(shù)間具有一定大小關(guān)系，如冠層頂部高度HVT 與冠層底部高度HVB、田間持水量SMCREF 與孔隙度SMCMAX、凋萎系數(shù)SMCWLT與田間持水量SMCREF 及蒸發(fā)停止的土壤含水率SMCDRY，需要對它們的取值范圍進(jìn)行約束。

表2 敏感性分析選用的參數(shù)及其取值范圍Tab.2 Selected parameters and the value ranges for sensitivity analysis

模型運行步長為0.5 h，并利用研究時段內(nèi)的氣象數(shù)據(jù)進(jìn)行多輪預(yù)熱。為了對模型參數(shù)敏感性進(jìn)行系統(tǒng)性的分析，本研究選取潛熱通量（LE）、顯熱通量（H）、凈生態(tài)系統(tǒng)碳交換量（NEE）、第1 層土壤含水率（SM1）及第2 層土壤含水率（SM2）進(jìn)行分析。ARSENAULT 等[6]通過數(shù)值試驗發(fā)現(xiàn)當(dāng)樣本量接近1 500 時，全局敏感度基本收斂。更重要的是，那些在小樣本量下表現(xiàn)出較高靈敏度的參數(shù)在較大樣本量下也表現(xiàn)出較高靈敏度，反之亦然。為了減少計算成本，先利用ARSENAULT 等[6]研究中的算例對42 個參數(shù)抽樣分別抽樣250次、500 次、1 000 次與1 500 次，探究合理的抽樣次數(shù)。有研究表明，動態(tài)植被模塊開啟時，氣孔阻抗的土壤水分限制因子的參數(shù)化方案對模型參數(shù)的取值有較大影響。因此，除對2個站點進(jìn)行參數(shù)敏感性分析外，以CLM 和Noah 2 種氣孔阻抗的土壤水分限制因子方案為代表，探索基于土壤水勢與基于土壤含水率的參數(shù)化方案對參數(shù)敏感性的影響。

1.3.2 參數(shù)敏感性分析評估方法

本研究采用Matlab 工具箱中的拉丁超立方抽樣法（LHS）對參數(shù)抽樣，抽樣得到的參數(shù)代入模型運行后，基于Sobol 方差分析計算參數(shù)的敏感性指標(biāo)（實現(xiàn)代碼獲取網(wǎng)址為https://github.com/wangxi1998/NoahMP-Sensitivity-NCP/tree/main）。全局敏感度Ti（Ti越大代表該模型輸出變量對第i個參數(shù)越敏感）的計算方法如下[6,16]：

式中：x 為進(jìn)行敏感性分析的N 個參數(shù)的向量（N=42）；xi為第i個參數(shù)，x～i為除xi外的N-1個其他參數(shù)的向量。蒙特卡洛抽樣的集合平均計算方法如下：

式中：M 為拉丁超立方抽樣次數(shù)；xm為第m 次抽樣過程得到的參數(shù)向量；是拉丁超立方抽樣的2個平行抽樣與的參數(shù)；f(x)為性能指標(biāo)。

性能指標(biāo)選用衡量預(yù)測值與真實值之間偏差的均方根誤差，計算方法如下：

式中：K 為模型輸出變量的時間維度；Rk為第k 個模型運行步長時刻的輸出結(jié)果；Ok為對應(yīng)時刻的觀測值。

2 結(jié)果與分析

2.1 模型預(yù)熱輪數(shù)與抽樣次數(shù)M的確定

模型預(yù)熱以相鄰2輪的平均相對誤差達(dá)0.001為結(jié)束條件，判斷標(biāo)準(zhǔn)見式(4)。2站點的預(yù)熱過程類似，此處僅展示禹城站結(jié)果（見圖1），圖1中結(jié)果均為各輪預(yù)熱后模型的最終結(jié)果。

圖1 禹城站預(yù)熱過程Fig.1 Spin-up process of Yucheng station

式中：Rk為第k 個模型運行步長時刻的輸出結(jié)果；K 為模型輸出變量的時間維度；n為變量達(dá)到穩(wěn)定所需預(yù)熱輪數(shù)。

不同輸出變量達(dá)到穩(wěn)定所需預(yù)熱時間的結(jié)果如表3所示。結(jié)果表明，土壤含水率（SM1和SM2）的模擬結(jié)果達(dá)到穩(wěn)定所需的時間最短，熱通量（LE、H和土壤熱通量G）相較需要更長的時間，其中顯熱通量的預(yù)熱結(jié)果不穩(wěn)定，預(yù)熱多輪后其相對誤差成周期性循環(huán)，取進(jìn)入循環(huán)前一年作為最終所需預(yù)熱輪數(shù)；NEE 模擬穩(wěn)定所需時間需要100 a 以上，與其值較小（10-4量級），達(dá)到數(shù)值穩(wěn)定更為困難不無關(guān)系。此結(jié)果與密西西比流域的研究結(jié)果存在差異，CAI 等[5]認(rèn)為土壤水預(yù)熱需要8 a，而顯熱與潛熱通量預(yù)熱的時間更短，大約需要4 a。

圖2 為不同抽樣次數(shù)的輸出變量敏感參數(shù)分析結(jié)果。結(jié)果表明：對于不同抽樣次數(shù)而言，隨著抽樣次數(shù)的減少，大部分參數(shù)的全局敏感度大小發(fā)生變化，尤其是植被相關(guān)參數(shù)（如QE25、SLA、LTOVRC 等）變化較為顯著，全局敏感度較高的參數(shù)變化較明顯。盡管如此，無論抽樣次數(shù)如何變化，對輸出變量不敏感的參數(shù)依舊只有非常小的全局敏感度。因此，在不考慮敏感性參數(shù)精確排序的條件下，若只為確定哪些參數(shù)更重要、更為敏感，拉丁超立方抽樣法（LHS）的抽樣次數(shù)M 可以取參數(shù)數(shù)量N 的5 倍左右，這與HUO 等[7]的研究結(jié)論具有一致性。下文中M 取為250 以節(jié)約計算成本。

圖2 不同抽樣次數(shù)的輸出變量參數(shù)敏感參數(shù)分析Fig.2 Parameter sensitivity analysis of output variables for different sampling number

2.2 站點參數(shù)敏感性結(jié)果分析

為了對參數(shù)的敏感性進(jìn)行區(qū)分，簡單地將全局敏感度高于0.10 的參數(shù)劃分為高敏感性參數(shù)，將全局敏感度高于0.05且不高于0.10 的參數(shù)劃分為中敏感性參數(shù)，禹城站與欒城站的敏感性參數(shù)劃分結(jié)果見表4。

表4 禹城站與欒城站敏感性參數(shù)匯總（btr = 2 : CLM）Tab.4 Summary of sensitivity parameters for two stations（btr = 2 : CLM）

綜合2 個站點的結(jié)果，由于LE 與H、SM1 與SM2 涉及的物理過程一致，因此其敏感性參數(shù)也大致相同。對LE 與H 來說，模擬過程不僅包含水量平衡還涉及到冠層能量平衡過程，因此它們對植被動態(tài)生長過程中的植物器官周轉(zhuǎn)率LTOVRC及25 ℃時葉綠素吸收的光量子效率QE25也很敏感。LE、H與土壤有關(guān)的敏感參數(shù)主要為BEXP、SMCMAX、SMCREF、SMCWLT，這些參數(shù)都與植被作用下的土壤水狀態(tài)密切相關(guān)。故而它們也是SM1 的敏感參數(shù)。與SM1 不同的是，由于SM2所處位置更深，土壤水與植被的相互作用更小，因此對凋萎系數(shù)SMCWLT 及田間持水量SMCREF 這類與植物相關(guān)的土壤參數(shù)敏感性較弱，而對土壤孔徑分布指數(shù)BEXP 及土壤孔隙率SMCMAX這類土壤物理性質(zhì)參數(shù)更敏感。

NEE 的敏感性參數(shù)為微生物呼吸參數(shù)MRP 和葉組織轉(zhuǎn)化速率LTOVRC，但NEE 對MRP 的敏感性顯著高于其他敏感參數(shù)，這一方面可以證實農(nóng)田生態(tài)系統(tǒng)的固碳潛力主要集中在農(nóng)田土壤[17]，另一方面與模型對農(nóng)田生態(tài)系統(tǒng)的作物碳同化過程模擬能力較弱有關(guān)，Noah-MP 模型的模擬結(jié)果低估了冬小麥與夏玉米生長旺季（4 月和8 月）時農(nóng)田生態(tài)系統(tǒng)的碳匯作用（見圖3）。同時，LTOVRC 是一個與葉片物質(zhì)量更新周轉(zhuǎn)進(jìn)入土壤有機碳有關(guān)的參數(shù)，這同樣能說明在模型模擬農(nóng)田生態(tài)系統(tǒng)NEE 過程中，農(nóng)田土壤碳庫的變化過程占主導(dǎo)，而與光合作用碳同化過程有關(guān)的作物碳庫變化過程較弱。與光合作用相關(guān)的VCMX25、QE25 及SLA，雖然也是重要參數(shù)，但NEE 對它們敏感性較低。

圖3 禹城站與欒城站NEE模擬值與真實值的比較Fig.3 Comparison of simulated and observed NEE for two stations

整體上，2個站點的敏感參數(shù)基本相同（見圖4）。結(jié)合所有輸出變量來看，2 個站點的區(qū)別體現(xiàn)在土壤孔隙率SMCMAX、田間持水量SMCREF、凋萎系數(shù)SMCWLT 和土壤孔徑分布指數(shù)BEXP 這4 個參數(shù)上。相比禹城站，欒城站模擬的LE、H 對SMCMAX 的變化更不敏感，而SM1、SM2 受SMCREF 和SMCWLT 的影響更大。SMCMAX 用于確定土壤水達(dá)到靜水平衡時的地下水位，其值為劃分土壤水與地下水的閾值，體現(xiàn)了土壤水與地下水間的轉(zhuǎn)換關(guān)系。禹城站SMCMAX 更為敏感，說明禹城站的地下水位更淺，地下水對土壤水分狀態(tài)的影響較大。站點資料顯示禹城站地下水位埋深一般為1.5～4.0 m，欒城站地下水位埋深大部分在15 m以下。SMCWLT 用來刻畫水分脅迫對植物生長、器官凋亡等生理過程的影響，SMCREF則決定了降水與灌溉水的下滲與土壤的疏干過程，BEXP 為土壤水分特征方程的重要參數(shù)，與土壤含水狀態(tài)密切相關(guān)。欒城站的淺層土壤水（SM1、SM2） 對SMCWLT、SMCREF、BEXP 更為敏感，說明欒城站的氣候更為干旱（禹城站該輪作年干旱指數(shù)為1.929，欒城站為2.540），土壤含水率受降水與灌溉的影響大，作物生長更可能受到水分脅迫，由氣象數(shù)據(jù)根據(jù)FAO 推薦的Penman-Monteith 方法計算的日參考作物需水量與降水見圖5。

圖4 禹城站與欒城站各輸出變量參數(shù)敏感參數(shù)分析（btr = 2 : CLM）Fig.4 Parameter sensitivity analysis of each output variable for two stations

圖5 禹城站與欒城站的日參考作物需水量與降水Fig.5 Daily reference crop evapotranspiration and precipitation for two stations

2.3 氣孔阻抗的土壤水分限制因子參數(shù)化方案的比較

在植被覆蓋情況下，氣孔阻抗不僅決定了植物蒸騰過程中水汽進(jìn)入大氣的能力，而且也決定了植物的光合作用中二氧化碳進(jìn)入葉肉組織的能力。植物的氣孔開閉受到土壤水分狀態(tài)的調(diào)控，在土壤含水率較低時，植物能感應(yīng)根區(qū)土壤水

分狀態(tài)并對氣孔進(jìn)行調(diào)節(jié)，防止自身過分失水而死亡。在Noah-MP 模型中，土壤水分對氣孔開閉影響的描述存在于氣孔阻抗的土壤水分限制因子參數(shù)化方案之中。2.2 節(jié)中選用的是基于土壤水勢的CLM 方案，本節(jié)將結(jié)合2.2 節(jié)的計算結(jié)果，比較CLM 方案與基于土壤含水率的Noah 方案在參數(shù)敏感性上的異同。與2.2節(jié)的劃分原則一致，禹城站與欒城站的中/高敏感參數(shù)見表5?？紤]參數(shù)抽樣引起的變異性，在比較不同氣孔阻力參數(shù)化方案時，與比較站點間差異一樣，重點關(guān)注全局敏感度的相對大小變化，分析那些相對變化明顯的參數(shù)。

表5 禹城站與欒城站敏感性參數(shù)匯總（btr = 1 : Noah）Tab.5 Summary of sensitivity parameters for two stations（btr = 1 : Noah）

結(jié)合2 個站點來看（圖6、圖7），使用不同氣孔阻力參數(shù)化方案對SM1、SM2、NPP 的參數(shù)敏感性幾乎無影響。而將參數(shù)化方案從CLM 更換為Noah 2 個站點都出現(xiàn)了以下現(xiàn)象：①對LE 來 說LTOVRC、QE25的敏感性降低，SMCREF、SMCWLT 的敏感度升高；②對H 來說QE25 的敏感性降低。由于控制氣孔阻力的土壤水分因子參數(shù)化方案的本質(zhì)區(qū)別在于Noah 方案基于土壤含水率而CLM 方案基于水勢，在Noah 方案中氣孔阻抗由SMCREF 和SMCWLT 直接決定，因此Noah 方案中LE 對這2 個參數(shù)更敏感。模型中的水碳耦合的關(guān)鍵樞紐在于葉片的氣孔，氣孔開放程度不僅影響植物蒸騰還影響光合過程，進(jìn)而間接影響作物的生長。LTOVRC、QE25 的敏感性降低能在一定程度上反映出LE 對土壤水勢的響應(yīng)靈敏度高于土壤含水率。H 與水通量的聯(lián)系沒有LE 那么緊密，因此與土壤水直接相關(guān)的SMCREF、SMCWLT的敏感度變化不明顯。而H 對QE25 的敏感性降低但LTOVRC 的敏感度變化不明顯，這反映出氣孔開閉對H 的影響主要是通過影響植物的光合過程實現(xiàn)。對比圖6 與圖7 可以看出，在氣候更為干旱的欒城站，不同參數(shù)化方案下BEXP 的敏感度發(fā)生了非常明顯的變化，原因在于BEXP直接參與CLM方案的水勢計算過程，因此在氣候干旱地區(qū)使用該方案時需要重視對BEXP的校準(zhǔn)。

圖6 禹城站不同參數(shù)化方案的各輸出變量參數(shù)敏感參數(shù)分析Fig.6 Parameter sensitivity analysis of each output variable for different parameterization schemes at Yucheng station

圖7 欒城站不同參數(shù)化方案的各輸出變量參數(shù)敏感參數(shù)分析Fig.7 Parameter sensitivity analysis of each output variable for different parameterization schemes at Luancheng station

3 討 論

華北平原主要的土地利用類型為耕地，是一個大型的農(nóng)田生態(tài)系統(tǒng)，水熱碳通量的變化受到種植、管理等人類活動的影響。與其他植被覆蓋類型相比，陸面模式在模擬農(nóng)田時存在著巨大的挑戰(zhàn)。為了準(zhǔn)確刻畫華北平原農(nóng)田生態(tài)系統(tǒng)中的水、能量及CO2通量變化過程，降低參數(shù)不確定性對陸面模式的影響，需要開展模型參數(shù)敏感性分析。

本研究在ARSENAULT 等[6]的基礎(chǔ)上對用于敏感性分析的單站點模型進(jìn)行修改，完善預(yù)熱前后淺層土壤中的速效碳庫FASTCP、深層土壤中的遲效碳庫STBLCP 等模型狀態(tài)變量的繼承關(guān)系，代碼修改前后欒城站的敏感性分析結(jié)果見圖8。結(jié)果顯示，模型代碼修改前后，輸出變量LE、H、SM1 與SM2對不同參數(shù)的敏感程度基本一致，NEE 的敏感性參數(shù)發(fā)生巨大變化。模型代碼修改后，LE、H 對土壤參數(shù)的敏感性得以凸顯，揭示了熱通量與土壤水的潛在聯(lián)系，參數(shù)BEXP 成為最敏感參數(shù)，與HUO 等[7]的研究結(jié)果一致。BEXP 是土壤的形狀參數(shù)，與土壤的孔隙結(jié)構(gòu)有關(guān)，同時其作為土壤水分特征方程的重要參數(shù)，能夠反映土壤的持水性，與陸氣間水分交換息息相關(guān)。對于NEE 來說，模型代碼修改后，敏感性參數(shù)不再是植物生理相關(guān)參數(shù)，微生物呼吸參數(shù)MRP 的全局敏感性指標(biāo)顯著高于其他參數(shù)，這恰能體現(xiàn)農(nóng)田生態(tài)系統(tǒng)碳匯功能主要依賴土壤固碳而不是農(nóng)作物生物量固碳。同時，本研究發(fā)現(xiàn)為使凈生態(tài)系統(tǒng)碳交換量（NEE）達(dá)到穩(wěn)定值需要經(jīng)歷100 a 以上的預(yù)熱過程。修改代碼并對模型進(jìn)行充分預(yù)熱后，為土壤微生物呼吸提供能量的淺層土壤中的速效碳庫FASTCP不斷積累，這也導(dǎo)致MRP 的敏感性增強。值得注意的是，NEE 預(yù)熱穩(wěn)定后FASTCP 的值為102量(t/hm2)，遠(yuǎn)高于第2 次全國土壤普查數(shù)據(jù)中農(nóng)田表層土壤有機碳密度（26.6～32.5 t/hm2）[17]，這是Noah-MP 陸面模式對NEE 等碳循環(huán)物理量模擬能力欠缺導(dǎo)致的，開發(fā)者應(yīng)重視對碳循環(huán)模塊的優(yōu)化。在對比2種氣孔阻抗的土壤水分限制因子參數(shù)化方案對模型敏感性的影響時發(fā)現(xiàn)，基于土壤水勢的CLM 方案能更好地將干旱地區(qū)植物根區(qū)的缺水狀態(tài)反映到LE 的變化上；而YANG 等[18]也得出了類似的結(jié)論，認(rèn)為基于土壤含水率的Noah 方案在降水充沛的西雙版納站能獲得最接近觀測值的模擬結(jié)果。ARSENAULT[6]等同樣發(fā)現(xiàn)在半干旱的草原與農(nóng)田生態(tài)站CLM方案對BEXP 的敏感度高于Noah 方案，并認(rèn)為PSISAT 在不同方案間的敏感度差異也源自同樣的機理。本研究并未發(fā)現(xiàn)PSISAT 在不同方案間的敏感度差異，這可能是因為選擇的下墊面單一，只聚焦于華北平原生態(tài)系統(tǒng)。

圖8 代碼修改前后欒城站參數(shù)敏感性分析結(jié)果Fig.8 Parameter sensitivity analysis results at Luancheng station before and after code modification

本研究揭示了Noah-MP 模型模擬華北平原農(nóng)田生態(tài)系統(tǒng)水熱碳通量及土壤含水率時的敏感參數(shù)，并對比了不同氣孔阻抗的土壤水分限制因子參數(shù)化方案下參數(shù)敏感度，對參數(shù)優(yōu)化與參數(shù)化方案選擇有著指導(dǎo)意義。但本研究仍存在一些不足，例如，只考慮用戶自定義參數(shù)的敏感性，而未考慮模型嵌入?yún)?shù)，CUNTZ 等[19]發(fā)現(xiàn)對潛熱通量、總徑流及各徑流分量而言，許多模型嵌入?yún)?shù)甚至比用戶自定義參數(shù)更敏感，因此未來可以將參數(shù)敏感性分析的范圍拓展至模型嵌入?yún)?shù)。

4 結(jié) 論

（1）若不需要對參數(shù)按敏感性強弱進(jìn)行嚴(yán)格排序，在使用拉丁超立方抽樣法（LHS）篩選模型的敏感性參數(shù)時，為節(jié)約計算成本，抽樣次數(shù)可以取為參數(shù)的5倍左右。

（2）用Noah-MP 模型模擬華北平原農(nóng)田生態(tài)系統(tǒng)典型站點時，LE、H 的敏感參數(shù)有：LTOVRC、QE25、BEXP、SMCMAX、SMCREF、SMCWLT。其中，BEXP、SMCMAX、SMCREF、SMCWLT 也是SM1 和SM2 的敏感參數(shù)。NEE 對MRP、LTOVRC 的敏感性顯著高于VCMX25、QE25 及SLA 等與作物光合作用相關(guān)的參數(shù)，體現(xiàn)農(nóng)田生態(tài)系統(tǒng)碳匯功能主要依賴農(nóng)田土壤固碳。

（3）在氣候干旱的地區(qū)需要重視對土壤參數(shù)SMCMAX、SMCREF、SMCWLT、BEXP 的校準(zhǔn)，尤其在使用CLM 氣孔阻抗的土壤水分限制因子參數(shù)化方案時，BEXP 的校準(zhǔn)需要特別關(guān)注。

致謝：感謝禹城農(nóng)業(yè)綜合試驗站與欒城農(nóng)業(yè)生態(tài)系統(tǒng)試驗站為本研究提供的數(shù)據(jù)支持，感謝Arsenault Kristi R.在GitHub 平臺貢獻(xiàn)的開源代碼，使本研究能在此基礎(chǔ)上提出改進(jìn)。