程 強(qiáng) 徐 嬙 陳 超 薛緒掌 王忠義 孫宇瑞

(1.中國農(nóng)業(yè)大學(xué)信息與電氣工程學(xué)院， 北京 100083； 2.北京農(nóng)業(yè)智能裝備技術(shù)研究中心， 北京 100097)

0 引言

越冬期作物根區(qū)充足的蓄水量能夠減小作物對(duì)于春季降水的依賴程度，有利于春季冬小麥的生長并有可能增加產(chǎn)量[1]。土壤中水分平衡的精確預(yù)測不僅是決定水資源可用性的關(guān)鍵，也對(duì)其優(yōu)化管理起到了重要作用[2]。

在季節(jié)性凍土區(qū)，越冬期冬小麥根區(qū)水分的補(bǔ)給渠道除降水和田間灌溉外，因上層土壤凍結(jié)驅(qū)動(dòng)下層土壤水分向上運(yùn)移也是冬小麥根區(qū)水分補(bǔ)給的重要來源。隨著溫度降至0℃以下，土壤上層液態(tài)水變成冰，降低了該層的水勢。在水勢梯度作用下，下層土壤的水分開始向上運(yùn)移，形成一個(gè)水熱耦合運(yùn)移機(jī)制。以往許多學(xué)者用復(fù)雜的數(shù)學(xué)模型對(duì)該機(jī)制進(jìn)行了研究[3-8]。

在越冬期土壤的水熱耦合運(yùn)移過程中，土壤未凍水含量、含冰量和土壤溫度是必不可少的3個(gè)關(guān)鍵變量[8]。土壤凍融特征曲線(Soil freezing characteristic curve, SFC)用于表示凍土中土壤未凍水含量和溫度之間的關(guān)系，在模擬凍土的水熱耦合運(yùn)移過程中起到了重要作用。由于SFC和土壤水分特征曲線(Soil moisture characteristic curve, SMC)之間的相似性[9]，一些學(xué)者利用SMC來估算SFC進(jìn)行凍融過程中的水熱耦合運(yùn)移模擬[10]。為提高模型的預(yù)測精度，一些學(xué)者[9,11]將針式時(shí)域反射計(jì)(TDR)和溫度傳感器結(jié)合獲取SFC，但是在傳感器安裝過程中容易破壞土壤結(jié)構(gòu)。而介電管式傳感器能夠原位測定土壤未凍水含量，提高了測量精度[12-13]?；谠粶y量的土壤未凍水含量和溫度數(shù)據(jù)，得到的水熱參數(shù)有具體的物理意義，比參數(shù)辨識(shí)(反演)法所使用的模型結(jié)構(gòu)簡單，并取代田間取樣后在實(shí)驗(yàn)室測定土壤參數(shù)的傳統(tǒng)做法。

此外，由于冬季作物進(jìn)入休眠期，基本不考慮蒸騰作用，而土壤表面蒸發(fā)量作為邊界通量也對(duì)根區(qū)的蓄水和熱量平衡起著重要作用。土壤上邊界的蒸發(fā)量不僅影響地表熱通量的變化[14]，還降低了表層土壤的未凍水含量。然而，以往的研究往往認(rèn)為冬季地表蒸發(fā)量可能很小甚至可能為零[15-17]。與這些研究不同的是，F(xiàn)LERCHINGER等[4]在其模型中引入了一個(gè)地表蒸發(fā)方程，以完善越冬期水熱耦合運(yùn)移(Simutaneous heat and water, SHAW)機(jī)制，并在之后的研究中用氣象數(shù)據(jù)對(duì)模型進(jìn)行了測試[14,18-19]。 SCOTT等[1]利用HYDRUS軟件[7]來模擬越冬期的非飽和流，通過考慮由最大向上流速或彭曼潛在蒸發(fā)計(jì)算得到地表蒸發(fā)通量(以較小者為準(zhǔn))。SAITO等[20]和ZENG等[21-22]提出了水-熱-蒸汽耦合的數(shù)值模型，并證明地表蒸發(fā)量對(duì)于土壤中水分運(yùn)移有重要影響。ZHANG等[23]提出的耦合模型中計(jì)算了凍融過程中液態(tài)水和蒸發(fā)量對(duì)于水熱耦合運(yùn)移的貢獻(xiàn)率，結(jié)果證明水蒸氣在所有的土壤深度處貢獻(xiàn)率都超過了15%，驗(yàn)證了蒸發(fā)量作為土壤凍融模型的上邊界條件，在模擬過程尤其是在淺層土壤中是不可忽略的。

本文基于原位測量的土壤未凍水含量和溫度(2011—2012年)，提出一種改進(jìn)型的凍融過程中水熱耦合運(yùn)移模型參數(shù)估算方法，并將估算參數(shù)用于模型精度的驗(yàn)證過程(2012—2013年)?？紤]到地表蒸發(fā)量對(duì)土壤水分運(yùn)移與再分布的重要性，本文擬構(gòu)建土壤凍融過程水熱耦合運(yùn)移模型用于模擬越冬期小麥根區(qū)蓄水與能量交換，并探討冬季地表蒸發(fā)量對(duì)土壤水熱耦合模型預(yù)測精度的影響。

1 材料與方法

1.1 模型介紹

1.1.1水流通量傳輸過程

FLERCHINGER等[4]在1989年提出使用一階偏微分方程來模擬凍土和非凍土中的水熱耦合運(yùn)移，其中水流通量的傳輸方程[24](向下為正)為

(1)

(2)

式中Ks——飽和導(dǎo)水率，m/h

θs——飽和含水量，m3/m3

b——土壤孔隙分布參數(shù)

θi——土壤冰含量，m3/m3

θl——未凍水含量，m3/m3

ρi——冰密度，kg/m3

ρl——液態(tài)水密度，kg/m3

qv——蒸汽通量，m/h

K——非飽和導(dǎo)水率，m/h

ψ——土壤基質(zhì)勢，m

U——水通量的匯項(xiàng)，m3/(m3·h)

z——土壤深度，mt——時(shí)間，h

式(1)各項(xiàng)依次代表：體積含水率的變化量；體積含冰量的變化量；流入土壤層的凈液態(tài)水通量；進(jìn)入土壤層的水蒸氣含量；由于根系吸收而形成的一個(gè)匯項(xiàng)。

土壤基質(zhì)勢ψ可通過Campbell模型[25]計(jì)算，即

(3)

式中ψe——空氣進(jìn)入勢，kPa

由Richards方程可以得到

(4)

式中D——土壤中水分的擴(kuò)散率，m/h

地表蒸發(fā)量作為土壤表層的水流通量，影響了凍融過程中的水熱耦合運(yùn)移，則由式(1)、(4)可得

(5)

由達(dá)西定律可得

(6)

式中ql——水流通量，m/h

在上邊界處ql為地表蒸發(fā)量。該模型以上邊界上的水流通量(即蒸發(fā)量)作為驅(qū)動(dòng)力，計(jì)算下一個(gè)迭代時(shí)間上的液態(tài)水含量；另一方面，凍土中水分的流動(dòng)同時(shí)會(huì)產(chǎn)生熱對(duì)流，進(jìn)而影響土壤中的熱量傳輸。

1.1.2熱量傳輸過程

土壤基質(zhì)中熱量傳輸(包括液態(tài)水產(chǎn)生的熱對(duì)流和蒸發(fā)吸收的潛熱)方程為

(7)

式中Cs——土壤體積熱容，J/(m3·℃)

T——土壤溫度，℃

ρv——水蒸氣密度，kg/m3

Lf——融化潛熱，kJ/kg

Lv——蒸發(fā)潛熱，kJ/kg

λs——土壤熱傳導(dǎo)率，W/(m·℃)

S——土壤層的源項(xiàng)，m3/(m3·h)

式(7)中各項(xiàng)依次代表：因溫度升高儲(chǔ)存的比熱；水凍結(jié)成冰所需要的潛熱；土壤層中的凈蒸發(fā)潛熱；進(jìn)入土壤層的凈熱傳導(dǎo)；由于液態(tài)水通量產(chǎn)生的熱對(duì)流；土壤層的源項(xiàng)(可能包括太陽輻射和長波輻射)。當(dāng)不考慮地表蒸發(fā)的影響時(shí)，忽略凈蒸發(fā)潛熱項(xiàng)。

土壤體積熱容Cs是土壤各組分的體積熱容之和，其計(jì)算公式為

Cs=∑ρjcjθj

(8)

式中ρj、cj、θj——土壤的第j個(gè)組成成分的密度、比熱容和體積分?jǐn)?shù)

DE VRIES[26]提出了用于計(jì)算土壤的熱傳導(dǎo)率λs的理論，該理論認(rèn)為非常濕潤的土壤可以看作是混合了液態(tài)水、土壤顆粒、冰晶以及分散的空氣的連續(xù)介質(zhì)。這種理想模型下的土壤熱導(dǎo)率的計(jì)算公式為

(9)

式中mj、λj為土壤的第j個(gè)組成成分的權(quán)重因子、熱傳導(dǎo)率，例如砂粒含量、粉粒含量、黏粒含量、有機(jī)質(zhì)含量、冰、液態(tài)水、空氣含量等。DE VRIES[26]討論了關(guān)于權(quán)重因子mj的計(jì)算方法。

1.2 參數(shù)估算

土壤凍融特征線(SFC)是指凍土中的溫度(0℃之下)和未凍水含量之間的關(guān)系，在構(gòu)建凍融過程中土壤中水熱耦合運(yùn)移，尤其是優(yōu)化數(shù)值模型的參數(shù)時(shí)發(fā)揮了重要作用[12]。當(dāng)土壤中出現(xiàn)冰時(shí)，液態(tài)水和固態(tài)冰的壓強(qiáng)與溫度有關(guān)，它們之間的關(guān)系為

(10)

式中pl——液相的壓強(qiáng)，Pa

pi——冰的壓強(qiáng)，Pa

π——土壤溶質(zhì)的滲透壓，Pa

TK——絕對(duì)溫度，K

假設(shè)土壤在低鹽度條件下，滲透壓為零；且冰中為零壓力，則方程(10)可簡化為[27]

(11)

式中g(shù)——重力加速度，m/s2

假設(shè)未凍土中水流通量和土壤水分特征線在凍土中同樣適用，則由方程(3)、(11)可得

(12)

凍土中同時(shí)測量得到的溫度和液態(tài)水含量確定了SFC，根據(jù)方程(12)可以得到T<0時(shí)，θl與T之間為冪函數(shù)關(guān)系(圖1)，通過曲線擬合可得到相應(yīng)的統(tǒng)計(jì)學(xué)參數(shù)，進(jìn)而反演出ψe和b的值，分別為-127 kPa和2.5。

圖1 土壤凍融特征曲線 Fig.1 Soil freezing characteristic curves

圖2 2011—2012年和2012—2013年兩個(gè)越冬期土壤凍融上邊界條件 Fig.2 Upper boundary conditions of air temperature and soil surface evaporation for soil freezing-thawing process over two winters of 2011—2012 and 2012—2013

1.3 邊界條件

模型中的邊界條件分為第1類邊界條件和第2類邊界條件。其中第1類邊界條件是指在上邊界給出了空氣溫度值；第2類邊界條件是指在上邊界給出了向上的水流通量(即蒸發(fā)量)的值(圖2)。從2011—2012年和2012—2013年兩個(gè)冬季的空氣溫度和蒸發(fā)量的值可知，2011—2012年冬季空氣最低溫度為-18.16℃，在11月24日左右溫度降至0℃左右，之后溫度一直下降至0℃以下，進(jìn)入凍結(jié)期，直至3月中旬時(shí)溫度上升至0℃以上，且存在明顯的融化潛熱階段；2012—2013年空氣最低溫度為-22.68℃，凍融期約為2012年11月29日至2013年3月1日，空氣溫度波動(dòng)較大，融化潛熱階段不明顯。

1.4 研究區(qū)概況

試驗(yàn)區(qū)位于北京市昌平區(qū)小湯山精準(zhǔn)農(nóng)業(yè)基地(116°26′39″E，40°10′43″N)，屬于大陸性季風(fēng)氣候，年平均降雨量為587 mm，年平均溫度為11.7℃。該基地屬于季節(jié)性凍土區(qū)，每年11—12月土壤開始凍結(jié)，2—3月開始融化，越冬期(11—3月)平均降雨量約為42 mm，土壤平均溫度約為-1.1℃。按照國際制土壤質(zhì)地分級(jí)標(biāo)準(zhǔn)，試驗(yàn)區(qū)土壤被劃分為粉砂質(zhì)壤土，其砂粒、粉粒和粘粒所占比例分別為39.9%、46.6%和13.5%，土壤飽和含水量為0.41 m3/m3，飽和導(dǎo)水率為0.01 m/h。

進(jìn)行田間試驗(yàn)時(shí)，在每個(gè)測量點(diǎn)處，將6根PVC管(長2 m、直徑50 mm)以2 m的間隔垂直插入土壤中，使用介電管式傳感器(Diviner 2000型, Sentek Sensor Technologies)進(jìn)行不同深度土壤水分的測量，測量間隔為10 cm。土壤未凍水含量可通過介電常數(shù)計(jì)算得到[22]，測量點(diǎn)的液態(tài)水含量取6個(gè)點(diǎn)的平均值，其水平和垂直方向上的其他位置數(shù)據(jù)通過插值計(jì)算獲得。土壤溫度使用一組數(shù)字化傳感器DS18B20型(DS18B20+， Maxim Integrated)進(jìn)行測量，其精度為0.5℃(在-55～125℃之間)。在距介電管式傳感器1 m遠(yuǎn)的位置挖洞，直徑為5 cm，將溫度傳感器以10 cm的間隔放入洞中。在放置溫度傳感器時(shí)，不斷地將該位置的原土壤填回至洞中[8]。

土壤蒸發(fā)量作為模型的邊界條件，使用基地的稱重式蒸滲儀系統(tǒng)得到。試驗(yàn)所用蒸滲儀數(shù)據(jù)為24套中型蒸滲儀(長1 m、寬0.75 m、高2.3 m)測量的平均值。每套蒸滲儀采用杠桿式稱重系統(tǒng)，在利用平衡塊抵消土箱和土體質(zhì)量后，使用質(zhì)量傳感器測量土壤中水分質(zhì)量，以反映土壤的蒸發(fā)量的變化[28]。傳感器的測量頻率為每15 min一次，靈敏度為0.05～0.1 mm。

1.5 誤差分析

(1)均方根誤差(Root mean square error, RMSE)用于反映估算值和實(shí)測值之間的總體差異，對(duì)特大或特小誤差反映敏感。當(dāng)RMSE越接近于0時(shí)，表明估算誤差越小，模擬精度越高。其計(jì)算公式為

(13)

式中Mi——實(shí)測值Si——模擬值

N——實(shí)測樣本數(shù)

(2)平均偏差(Mean bias error, MBE)用于反映估算值和實(shí)測值之間的平均偏差，正值表示蒸散量被高估，負(fù)值表示被低估。MBE越接近于0時(shí)，精度越高。其計(jì)算公式為

(14)

(3)d值表示模擬值和實(shí)測值之間的一致性和精確度，取值范圍0≤d≤1，當(dāng)d越接近1時(shí)，代表模型的模擬值精確度越接近100%。其計(jì)算公式為

(15)

2 結(jié)果與討論

2.1 2011—2013年兩年冬季凍融過程的實(shí)測和模擬結(jié)果比較

春季冬小麥返青期之前其根系分布范圍約為0～20 cm[29-30]，因此根據(jù)土壤10 cm和20 cm處未凍水含量和溫度的模擬值和實(shí)測值的對(duì)比情況來分析冬小麥根區(qū)的水熱運(yùn)移(圖3、4)。在2011—2012年冬季，基于原位測量的土壤未凍水含量和土壤溫度(負(fù)溫條件時(shí))得到的SFC參數(shù)和考慮了地表蒸發(fā)量的水熱耦合運(yùn)移模型進(jìn)行數(shù)值模擬，土壤溫度模擬值與實(shí)測值擬合度較高，模擬效果較好(d值均在0.95左右，表1)。由于凍土中未凍水含量主要受溫度影響，未凍水含量的模擬精度也達(dá)到較高水平，其中，10 cm處d值為0.930，RMSE為0.036 m3/m3，20 cm處d值為0.870，RMSE為0.055 m3/m3。在溫度逐漸下降至0℃，土壤20 cm處的未凍水逐漸凍結(jié)的過程中，模擬值出現(xiàn)一定的低估，但是土壤溫度的模擬值和實(shí)測值之間該現(xiàn)象并不明顯。出現(xiàn)這種情況的原因可能是在SFC的標(biāo)定過程中(圖1)，使用整個(gè)凍融過程的凍土溫度和未凍水含量數(shù)據(jù)進(jìn)行曲線擬合，但SFC曲線在凍結(jié)-融化過程中存在滯后效應(yīng)，土壤凍結(jié)時(shí)速度較快，難以滿足穩(wěn)態(tài)或準(zhǔn)穩(wěn)態(tài)的假設(shè)條件[10,31]，在凍結(jié)階段，相同的溫度條件下擬合曲線的未凍水含量值相對(duì)于實(shí)測值偏低。在融化階段，由于不同深度處土壤溫度的高估導(dǎo)致未凍水含量出現(xiàn)不同程度的高估，如2011—2012年冬季20 cm處融化階段的溫度高估了0.3℃，但是未凍水含量卻上升了約0.1 m3/m3。0℃附近較小的溫度變化可能引起未凍水含量的較大差異(圖1)，這主要是受融化潛熱的影響，而擬合的SFC曲線在融化階段的相對(duì)高估也是影響未凍水含量預(yù)測精度的原因之一。

將2012—2013年的土壤表層蒸發(fā)量、空氣溫度、土壤初始溫度和未凍水含量作為輸入數(shù)據(jù)，利用數(shù)值模型和已估算的SFC參數(shù)可得到未凍水含量的預(yù)測值，將該預(yù)測值和實(shí)測值進(jìn)行對(duì)比分析來驗(yàn)證模型的適用性。整個(gè)凍融過程中，預(yù)測值能夠較好地追蹤到土壤中溫度和未凍水含量的動(dòng)態(tài)變化，土壤溫度的模擬值整體偏低，在土壤凍結(jié)后以及融化過程中精度相對(duì)提高(圖3、4)。未凍水含量在10 cm和20 cm處模擬值的d值分別為0.924和0.774，RMSE分別為0.046 m3/m3和0.071 m3/m3(表2)。在凍結(jié)階段，土壤溫度的模擬值出現(xiàn)了明顯的低估，從而導(dǎo)致未凍水含量在凍結(jié)過程中的低估(圖3)。

圖3 2011—2012和2012—2013兩個(gè)越冬期土壤10 cm和20 cm處未凍水含量模擬值與實(shí)測值比較 Fig.3 Comparisons of simulated and measured soil unfrozen water contents at 10 cm and 20 cm over two winters of 2011—2012 and 2012—2013

圖4 2011—2012年和2012—2013年兩個(gè)越冬期土壤10 cm和20 cm處溫度模擬值與實(shí)測值比較 Fig.4 Comparisons of simulated and measured soil temperatures at 10 cm and 20 cm across two winters of 2011—2012 and 2012—2013

2012—2013年土壤凍融過程水熱耦合運(yùn)移模型的數(shù)值模擬結(jié)果具有較高精度(表2)，驗(yàn)證了這種基于原位測量數(shù)據(jù)的改進(jìn)型模型參數(shù)估算方法的可行性與適用性。傳統(tǒng)的田間取樣方法是指在觀測點(diǎn)周圍取土樣，然后在實(shí)驗(yàn)室環(huán)境下測定水熱參數(shù)。然而，土壤的空間變異性會(huì)引起較大的誤差，取樣過程還會(huì)對(duì)土壤的結(jié)構(gòu)造成破壞，影響測量精度。此外，從田間取樣到實(shí)驗(yàn)室測定參數(shù)還耗費(fèi)了大量的人力和時(shí)間成本。CHENG 等[12]比較了由原位測量的SFC推導(dǎo)出的土壤水分特征曲線(0.806≤R2≤0.994)和田間取樣后實(shí)驗(yàn)室壓力盤測定的土壤水分特征曲線(0.980≤R2≤0.997)，結(jié)果表明該方法的估算精度與傳統(tǒng)方法相比在重疊區(qū)差異很小，證明了該方法的有效性。

2.2 蒸發(fā)量對(duì)模型預(yù)測結(jié)果的影響

為分析越冬期地表蒸發(fā)量對(duì)土壤水分運(yùn)移的影響程度，在忽略蒸發(fā)量影響的假設(shè)下運(yùn)用模型模擬出一組未凍水含量的變化情況。根據(jù)2012—2013年冬季凍融過程中土壤未凍水含量有無蒸發(fā)時(shí)的對(duì)比(圖3c、3d)，可以發(fā)現(xiàn)在不考慮蒸發(fā)時(shí)，未凍水含量和溫度的模擬值均高估，且隨著凍融過程的進(jìn)行，地表蒸發(fā)量對(duì)未凍水含量的累計(jì)影響越來越明顯。在融化階段不考慮蒸發(fā)時(shí)，比考慮蒸發(fā)影響高估約0.024 m3/m3，約占土壤未凍水含量的16%。對(duì)比不同土壤深度處蒸發(fā)量的影響，發(fā)現(xiàn)蒸發(fā)量對(duì)于土壤凍融過程中未凍水含量的影響隨著土壤深度的增加而逐漸減小。在模型中加入地表蒸發(fā)量作為邊界條件后，模型的模擬精度普遍提高，尤其是在土壤表層。對(duì)有無蒸發(fā)條件下模擬值與實(shí)測值之間進(jìn)行誤差分析(表2)可以發(fā)現(xiàn)：考慮蒸發(fā)量的影響后，10 cm處的模擬值的d值由0.892增加至0.924，RMSE由0.059 m3/m3降為0.046 m3/m3，而20 cm處的影響較小。因此，對(duì)越冬期小麥根區(qū)的未凍水含量的模擬應(yīng)該考慮蒸發(fā)量的影響，尤其是接近土壤表層[17]。

表1 2011—2012年冬季土壤未凍水含量和溫度誤差分析 Tab.1 Analysis of soil unfrozen water content and temperature in winter of 2011—2012

越冬期土壤凍融過程中，2012—2013年的蒸發(fā)量較前一年較高(圖2)，對(duì)模型的預(yù)測精度影響較大。出現(xiàn)該現(xiàn)象的原因是該年冬季降水量相對(duì)較大，使土壤表層含水量增加，從而增加了地表蒸發(fā)量。因此，模型的預(yù)測精度也會(huì)受到氣候條件等的影響。

表2 2012—2013年冬季有無蒸發(fā)時(shí)土壤未凍水含量和溫度模擬值誤差對(duì)比分析 Tab.2 Analysis of soil unfrozen water content and temperature with or without soil evaporation in winter of 2012—2013

3 結(jié)論

(1)提出了一種土壤凍融過程水熱耦合運(yùn)移參數(shù)原位估算方法，借助2011—2012年原位測定的試驗(yàn)數(shù)據(jù)擬合SFC曲線，對(duì)土壤水熱參數(shù)進(jìn)行最優(yōu)估算，然后檢驗(yàn)了該方法的適用性。結(jié)果表明：2012—2013年冬季土壤凍融過程中模型的模擬值與實(shí)測值擬合度較高，能夠準(zhǔn)確預(yù)測越冬期小麥根區(qū)未凍水含量與溫度的變化規(guī)律，10 cm處土壤未凍水含量和溫度的d值分別達(dá)到0.924和0.918，20 cm處達(dá)到0.774和0.833?？梢姡摲椒ü浪愕乃疅釁?shù)能夠保證模型的預(yù)測精度，從而取代田間取樣并在實(shí)驗(yàn)室測定土壤參數(shù)的傳統(tǒng)做法，比參數(shù)辨識(shí)(反演)法所使用的模型結(jié)構(gòu)簡單，不破壞土壤結(jié)構(gòu)，節(jié)省了人力和時(shí)間成本。

(2)為檢驗(yàn)地表蒸發(fā)量對(duì)模型精度的影響，對(duì)2012—2013年越冬期的模擬值和實(shí)測值進(jìn)行誤差分析，10 cm處的未凍水含量預(yù)測值的d值由0.892增加為0.924，RMSE由0.059 m3/m3降為0.046 m3/m3。結(jié)果表明，考慮地表蒸發(fā)量對(duì)土壤凍融過程中水熱運(yùn)移的影響后，土壤表層未凍水含量預(yù)測精度明顯提高。此外，隨著土壤深度的增加，蒸發(fā)量的影響逐漸減小。

(3)除土壤邊界條件和氣候因素外，SFC曲線在凍融過程的滯后效應(yīng)可能是造成模擬過程中凍結(jié)階段低估和融化階段高估的重要原因。且土壤在經(jīng)歷多次凍融循環(huán)后，紋理結(jié)構(gòu)和孔隙分布等物理性質(zhì)可能會(huì)產(chǎn)生改變，進(jìn)而影響SFC。因此，SFC曲線及其參數(shù)的標(biāo)定和修正是土壤凍融過程中進(jìn)行水熱耦合運(yùn)移數(shù)值模擬的重要步驟，也需要在未來的研究中找出更加適當(dāng)和精確的模型標(biāo)定與參數(shù)估算方法。