汪 躍，高晶晶，譚 霄，譚 博，李龍國，劉 超

（1. 四川省水利廳規(guī)劃計(jì)劃處，四川 成都 610017； 2. 水力學(xué)與山區(qū)河流開發(fā)保護(hù)國家重點(diǎn)實(shí)驗(yàn)室，四川 成都 610065；3. 四川大學(xué)水利水電學(xué)院，四川 成都 610065）

0 引 言

鎘Cd的生物毒性強(qiáng)且分布廣泛，可通過食物鏈進(jìn)入到人體內(nèi)并積累，進(jìn)而危害人體健康［1］。施用外源硅Si 是一種有效的土壤重金屬污染治理方法，可以促進(jìn)土壤中的鎘向無效態(tài)轉(zhuǎn)化［2］，抑制植物對(duì)鎘的吸收轉(zhuǎn)運(yùn)［3-5］。此外，有研究發(fā)現(xiàn)外源硅還能夠顯著增加土壤有機(jī)碳含量，促進(jìn)土壤團(tuán)聚體的形成和穩(wěn)定，增加土壤通氣性，進(jìn)一步改善土壤結(jié)構(gòu)［6-8］。土壤水分運(yùn)動(dòng)受到團(tuán)聚體及有機(jī)質(zhì)含量的顯著影響［9］，例如土壤團(tuán)聚體可以增強(qiáng)土壤水分的入滲［10］，而有機(jī)質(zhì)含量的增加可以增強(qiáng)團(tuán)聚體穩(wěn)定性和土壤保水性［11］。由于土壤溶質(zhì)遷移，如氮、磷等營(yíng)養(yǎng)元素的流失受土壤水分運(yùn)動(dòng)影響強(qiáng)烈，施硅對(duì)土壤水分運(yùn)動(dòng)的影響可進(jìn)一步延伸到對(duì)面源污染物遷移的影響，施硅存在增加面源污染的風(fēng)險(xiǎn)，因此明確施硅對(duì)土壤水分運(yùn)動(dòng)的影響十分有必要。雖然上述研究定性描述了施硅對(duì)土壤結(jié)構(gòu)和組分的影響，但其對(duì)土壤水分運(yùn)動(dòng)的直接影響還不清楚，缺少定量研究與機(jī)理分析。

土壤水分特征曲線（簡(jiǎn)稱土-水曲線）是表示土壤水吸力與土壤含水率關(guān)系的曲線，同時(shí)也是一種重要的土壤水動(dòng)力參數(shù)，是進(jìn)行土壤水分和溶質(zhì)遷移模擬的重要基礎(chǔ)［12］。土壤水分特征曲線受多種因素影響，如土壤質(zhì)地、結(jié)構(gòu)、容重、溫度及土壤添加物等［13］。王志超等［14］研究發(fā)現(xiàn)土壤中的殘膜會(huì)影響土壤水分特征曲線，殘膜含量越高，土壤的保水能力越低。程?hào)|娟等［15］研究發(fā)現(xiàn)土壤中添加秸稈在增加土壤含水量的同時(shí)降低了土壤的比水容重。王忠江等［16］研究發(fā)現(xiàn)土壤添加生物炭時(shí)，不同質(zhì)地土壤的持水能力均得到提高。上述研究表明各種添加物均會(huì)對(duì)土壤水分特征曲線有不同程度的影響，而施加外源硅治理土壤鎘污染時(shí)，硅作為添加物對(duì)土壤水分特征曲線影響的研究尚未見報(bào)道，與此同時(shí)硅鎘反應(yīng)導(dǎo)致的有效硅含量減少以及生成的硅鎘配合物對(duì)含鎘土壤水分特征曲線是否有影響以及影響機(jī)理均不清楚。因此本文通過測(cè)量不同硅、鎘處理土壤的水分特征曲線，且進(jìn)行模型擬合，以期定量描述施硅對(duì)含鎘土壤水分運(yùn)動(dòng)的影響，并從土壤孔隙結(jié)構(gòu)、表面電位等揭示其影響機(jī)理。研究結(jié)果對(duì)于指導(dǎo)硅肥的科學(xué)施用具有重要的理論意義。

1 材料與方法

1.1 實(shí)驗(yàn)樣品

邛崍市位于成都市以西75 km，屬于亞熱帶季風(fēng)濕潤(rùn)氣候，當(dāng)?shù)刂饕捎盟?小麥、水稻-玉米、水稻-油菜等輪作制度。以邛崍市前進(jìn)鎮(zhèn)為研究區(qū)，通過前期野外調(diào)查發(fā)現(xiàn)，前進(jìn)鎮(zhèn)農(nóng)田土壤總鎘含量為0.24～6.52 mg/kg，根據(jù)《土壤環(huán)境質(zhì)量農(nóng)用地土壤風(fēng)險(xiǎn)管控標(biāo)準(zhǔn)（試行）》（GB15618-2018），部分土壤已超過風(fēng)險(xiǎn)管制值。實(shí)驗(yàn)所用土壤取自前進(jìn)鎮(zhèn)水稻田（北緯30°24′，東經(jīng)103°32′），風(fēng)干后研磨并過10 目篩，其基本理化性質(zhì)如表1所示。

表1 土壤初始理化性質(zhì)Tab.1 Initial physicochemical properties of soil

1.2 實(shí)驗(yàn)方法

基于調(diào)查的邛崍水稻土鎘濃度，試驗(yàn)中設(shè)置0，2.5 和5 mg/kg 三個(gè)鎘添加水平，為確保實(shí)驗(yàn)的準(zhǔn)確性，使有效態(tài)鎘含量趨于平衡，本實(shí)驗(yàn)對(duì)添加的外源鎘進(jìn)行了30天的老化處理［17］。并參考Vasudevan 等的研究［18］，設(shè)置0，200， 400 kg SiO2/ hm2三個(gè)硅添加水平共9 種處理，每種處理3 個(gè)重復(fù)。對(duì)照組縮寫為CK，處理組縮寫為 CdiSij，i=2.5 和5，表示鎘添加水平；j=2 和4，表示SiO2添加水平。

實(shí)驗(yàn)步驟如下，稱取溴化鎘（CdBr2）試劑粉末2.27 和4.97 mg 分別與兩份1 000 g 前述初始土樣充分混合，以配制鎘添加濃度為2.5 和5 mg/kg 的土樣。取前述初始土樣以及配制好的鎘添加土樣，以1.05 g/cm3的容重裝填保持環(huán)（retaining ring，內(nèi)徑：5.5 cm， 高：1 cm），用相應(yīng)濃度的硅酸鈉（NaSiO3·9H2O）溶液逆向飽和24 h 以實(shí)現(xiàn)不同硅添加水平。飽和樣品稱重后放入壓力薄膜儀（1500F2， Soilmoisture Equipment Crop.， USA）中，分別以0，75，140，230，380，400，480，610，710，810，800，900，1 000 和1 200 kPa 進(jìn)行加壓，壓力值的設(shè)置參考相關(guān)研究確定［19］。每個(gè)壓力值下，樣品在壓力膜儀中至少加壓48 h直至不再排水，加壓結(jié)束后迅速拿出稱重記數(shù)，并使用烘干法測(cè)量最終土壤含水率。

烘干的土壤樣品研磨過10 目篩后，按照土水比1∶10 加入超純水，震蕩3 min 后過濾，通過pH 計(jì)（pHS-3C，雷磁，上海儀電，中國）測(cè)量濾液pH 值。上述烘干土樣研磨過200 目篩，取0.01 g 溶于1.5 mL 超純水中，進(jìn)行超聲處理5 min，之后經(jīng)0.45 μm 濾膜過濾后上機(jī)測(cè)試zeta 電位（Zetasizer Nano ZS， Malvern Instruments， Britain），每個(gè)樣品測(cè)試3次。

1.3 土壤當(dāng)量孔徑體積占比

土壤孔隙分布是影響水分特征曲線的重要因素之一，因此本文將通過土壤當(dāng)量孔徑占比的計(jì)算來解釋施硅對(duì)土壤水分特征曲線影響的機(jī)理。在非飽和的土-水系統(tǒng)中，根據(jù)土壤毛管理論可知，土壤中某一范圍孔徑的圓形毛管的毛管力作用會(huì)產(chǎn)生土壤水吸力，可表示為：

式中：S為土壤水吸力，hPa；σ為水的表面張力系數(shù)，常溫下為7.5×10-4N/cm；r0為毛管半徑，mm。

以便于區(qū)別土壤真實(shí)的孔徑，此時(shí)把毛管直徑稱為當(dāng)量孔徑，那么當(dāng)量孔徑與土壤水吸力的關(guān)系為［21］：

式中：D為當(dāng)量孔徑，mm。

θ1、θ2（θ1＞θ2）時(shí)的土壤當(dāng)量孔徑可計(jì)算得到為D1、D2，則可計(jì)算得出土壤中當(dāng)量孔徑范圍在D2～D1之間土壤孔隙所占的體積與孔隙總體積占比是θ1-θ2。

1.4 土壤水分特征曲線模型

土壤水分特征曲線是在一定溫度條件下，土壤含水率與土壤水吸力的關(guān)系曲線，是土壤最重要的水力特性之一［13］。由于土壤水分特征曲線的影響因素復(fù)雜，至今尚沒有從理論上建立土壤含水率和土壤水吸力之間的關(guān)系，通常用經(jīng)驗(yàn)公式來描述，使用比較廣泛的模型有Brooks-Corey（BC）模型［20］、van Genuchten（VG）模型［21］和Dual-Porosity（DP）［22］等。本研究應(yīng)用美國國家鹽改中心（US Salinity Laboratory）開發(fā)的RETC 軟件［23］，基于最小二乘法對(duì)不同處理土壤的水分特征曲線數(shù)據(jù)進(jìn)行不同模型的擬合，確定模型參數(shù)。

（1）Brooks-Core模型：

式中：Se為土壤的飽和度；θ為土壤體積含水率；θr為土壤殘余體積含水率；θs為土壤飽和體積含水率；α為進(jìn)氣值的倒數(shù)；h為土壤吸力；λ為與土壤孔隙分布有關(guān)的參數(shù)。

（2）Van- Genuchten模型：

式中：m和n為影響水分特征曲線形態(tài)的經(jīng)驗(yàn)參數(shù)，m= 1 - 1/n(n＞1)，其余符號(hào)含義與式（4）相同。

（3）Dual-Porosity模型：

式中：w1和w2為兩個(gè)區(qū)域的權(quán)重因子，不同的區(qū)域代表不同的孔隙系統(tǒng)，其中0＜w1，w2＜1，∑W= 1；α1和α2為兩個(gè)區(qū)域進(jìn)氣值的倒數(shù)；m1和m2、n1和n2為兩個(gè)區(qū)域的土壤孔隙尺寸的分布參數(shù)；其余符號(hào)含義與式（4）相同。

1.5 數(shù)據(jù)統(tǒng)計(jì)與分析

試驗(yàn)數(shù)據(jù)用Excel 和origin 軟件進(jìn)行分析和作圖，顯著性分析使用SPSS軟件進(jìn)行分析，顯著水平為P＜0.05。不同模型對(duì)不同處理的數(shù)據(jù)的擬合效果不同，根據(jù)擬合模型的預(yù)測(cè)值與實(shí)測(cè)值，可計(jì)算決定系數(shù)R2，均方根誤差（Root Mean Square Error，RMSE）和幾何平均數(shù)（Geometric Mean，GMER）這3 個(gè)統(tǒng)計(jì)指標(biāo)來綜合評(píng)價(jià)不同模型在不同處理下的擬合效果。

2 結(jié)果和分析

2.1 硅、鎘對(duì)土壤持水特性的影響

2.1.1 硅、鎘對(duì)土壤水分特征曲線的影響

圖1反映了不同處理下土壤水吸力隨體積含水率的變化趨勢(shì)。由圖可見，Cd2.5和Cd5處理與CK 相比，幾乎重合，說明鎘的添加對(duì)土壤水分特征曲線的影響不大；Si2和Si4處理下的水分特征曲線相比于CK 整體向右偏移明顯，即同一吸力下加硅處理土壤具有更高的含水率，說明加硅后土壤的持水能力有所增強(qiáng)；在0～710 kPa 吸力段，Cd2.5Si2、Cd2.5Si4、Cd5Si2和Cd5Si4處理下的水分特征曲線相比于其他處理略向右偏移，但是在710～1 200 kPa吸力段卻向左劇烈偏移，說明加硅會(huì)使含鎘土壤高吸力段的持水能力相對(duì)于其他處理大幅減弱。

圖1 不同處理下土壤水分特征曲線Fig.1 Water characteristic curves of different treatments

2.1.2 硅、鎘對(duì)不同形態(tài)土壤水的影響

在不同的吸力范圍內(nèi)，土壤水的形態(tài)及其作用力有所不同。在8～630 kPa 吸力段，土壤水以毛管水為主，主要受土壤孔隙分布影響；在630～3 100 kPa吸力段，土壤水以膜狀水為主，主要受土壤顆粒表面的范德華力、水分子內(nèi)聚力等影響［13］。圖1顯示含鎘土壤加硅之后在710 kPa 吸力前后持水性有顯著變化，說明毛管水和膜狀水受到的影響并不相同。為分析土壤對(duì)毛管水吸持性的影響，根據(jù)蔣文君等的研究［24］，我們以140 kPa作為低吸力段和高吸力段的分界點(diǎn)進(jìn)行0～710 kPa 吸力范圍內(nèi)的當(dāng)量孔徑孔隙體積占比計(jì)算，結(jié)果如表2所示。結(jié)果顯示，低吸力段大孔隙當(dāng)量孔徑介于0.002 14～0.002 50 mm，高吸力段小孔隙當(dāng)量孔徑介于0.000 49～0.002 14 mm。CK 處理低吸力段孔隙體積占比為18.65%，高吸力段孔隙體積占比為11.35%。與CK 相比，單獨(dú)加鎘對(duì)土壤孔隙沒有顯著影響；單獨(dú)加硅之后土壤大孔隙體積占比隨著添加濃度的增加而增大，但小孔徑體積占比沒有顯著變化。含鎘土壤加硅之后，土壤大孔隙體積占比排序?yàn)镃d2.5Si4＞Cd5Si4＞Cd2.5Si2＞Cd5Si2，相同含鎘量下，施硅量越高則大孔隙體積占比越高，而小孔隙體積占比無明顯變化，與單獨(dú)加硅處理展現(xiàn)的規(guī)律一致。但是在相同施硅量下，鎘含量越高則大孔隙體積占比越低，雖然該趨勢(shì)沒有統(tǒng)計(jì)上的顯著差異，但從某種程度上說明鎘與硅的反應(yīng)對(duì)硅增加土壤大孔隙數(shù)量的效應(yīng)產(chǎn)生影響。有研究表明施硅可以顯著降低土壤中交換態(tài)和鐵錳結(jié)合態(tài)鎘含量，顯著提高碳酸鹽結(jié)合態(tài)和殘?jiān)鼞B(tài)鎘的含量，同時(shí)伴隨著有效硅含量的降低［25］。因此，施硅增加土壤大當(dāng)量孔徑孔隙主要是通過有效硅實(shí)現(xiàn)，硅-鎘反應(yīng)導(dǎo)致的形態(tài)變化會(huì)削弱硅對(duì)土壤孔隙的影響。

在710～1 200 kPa 吸力段，Cd2.5Si2、Cd2.5Si4、Cd5Si2和Cd5Si4處理的水分特征曲線向左劇烈偏移，說明含鎘土壤加硅之后對(duì)膜狀水的吸持減弱而單獨(dú)加鎘/硅均無影響。

表2 不同處理不同吸力段當(dāng)量孔徑體積占比Tab.2 Equivalent diameter distributions of soil pore for different treatments

2.2 硅、鎘對(duì)土壤水分特征曲線影響的機(jī)理

為進(jìn)一步分析硅、鎘對(duì)土壤水分特征曲線影響的機(jī)理，本研究測(cè)定了土壤的zeta 電位及pH 值，結(jié)果如表3所示。結(jié)果顯示，所有處理土壤的zeta電位值處于-25.22～-45.17 mV 之間，這與其他研究者測(cè)定的土壤zeta 電位值區(qū)間（-15～-60 mV）大體一致［26］。與CK 相比，單獨(dú)加鎘土壤的zeta 電位無明顯變化，單獨(dú)施硅土壤的zeta電位絕對(duì)值顯著增大，且與施硅量成正比，例如Si2處理土壤的zeta 電位絕對(duì)值較CK 提高了38.67%，Si4提高了74.33%。施硅后含鎘土壤的zeta 電位絕對(duì)值較不加硅處理（CK 和單獨(dú)加鎘）土壤顯著升高，較單獨(dú)加硅土壤略有降低但不顯著。同時(shí)，施硅會(huì)顯著增加土壤的pH 值，并且zeta 電位絕對(duì)值與pH 值完全正相關(guān)。施硅增加zeta 電位絕對(duì)值一方面是因?yàn)楣杷岱肿咏Y(jié)合至土壤表面后，由于其含有大量電負(fù)性較高的含氧基團(tuán)（Si-OH），使得土壤zeta 電位絕對(duì)值增大［24］；另一方面硅酸鈉水解反應(yīng)使得土壤pH值升高，也會(huì)導(dǎo)致zeta電位絕對(duì)值增大［27］。

表3 不同處理土壤的zeta電位與pH值Tab.3 zeta potential and pH value of soils in different treatments

已有研究表明土壤zeta電位的絕對(duì)值決定了土壤膠體的抗分散能力和穩(wěn)定性，可以用來表示土壤顆粒之間彼此相互作用的強(qiáng)度。根據(jù)DLVO 理論［28，29］，zeta 電位絕對(duì)值越大，在靜電斥力的作用下，土壤膠體就越分散，分散的膠體更容易隨著水流運(yùn)動(dòng)而堵塞大孔隙，導(dǎo)致土壤的小孔隙體積占比增大［30］。但是由于壓力膜儀脫水過程中保持環(huán)高度僅為1 cm，難以形成連續(xù)水流帶動(dòng)細(xì)小顆粒運(yùn)動(dòng)，便沒有發(fā)生土壤孔隙堵塞，因此加硅對(duì)小孔隙占比沒有影響。土壤團(tuán)聚體的形成包括土壤膠體顆粒絮凝形成微團(tuán)粒以及微團(tuán)粒之間的膠結(jié)。土壤膠體顆粒能否絮凝取決于顆粒之間的交互作用能［31］。加硅后土壤zeta 電位絕對(duì)值增加，土壤膠體顆粒之間的靜電斥力增大，顆粒之間相互排斥不易聚集，似乎與郭磊等［6］的研究結(jié)論相反。實(shí)際上可溶性硅酸鹽在水中發(fā)生水解反應(yīng)可生成硅酸，硅酸進(jìn)一步形成無定形的二氧化硅沉淀，并以膠態(tài)粒子或凝膠出現(xiàn)。雖然施硅增加zeta 電位使得土壤膠體無法絮凝形成新的土壤微團(tuán)粒，但硅酸凝膠可使已有微團(tuán)粒膠結(jié)在一起，微團(tuán)粒之間的孔隙是施硅增加土壤大孔隙的原因。

對(duì)土壤膜狀水而言，土壤顆粒表面的總位能決定著對(duì)膜狀水的吸持，主要受靜電斥力和范德華引力支配。范德華引力只與作用距離和哈梅克常數(shù)（Hamaker constant）有關(guān)。靜電斥力受電解質(zhì)類型、濃度、pH值等影響，其本質(zhì)是因?yàn)檫@些環(huán)境條件改變使顆粒表面電位值發(fā)生改變，從而影響顆粒間的靜電斥力［30］。大多數(shù)研究以zeta 電位近似為表面電位，含鎘與不含鎘土壤加硅之后zeta電位絕對(duì)值都顯著增加，靜電斥力增大，但是只有含鎘土壤加硅之后才表現(xiàn)出對(duì)膜狀水吸持的減弱，單獨(dú)加硅并沒有影響，說明靜電斥力不起主要作用。郭磊［6］通過X 射線衍射儀、傅立葉變換紅外光譜儀以及拉曼光譜儀等表征手段，確認(rèn)了水溶性硅-鎘配合物的生成。可以合理推測(cè)是硅-鎘反應(yīng)生成水溶性硅-鎘配合物溶解在水中，可能引起水分子與土壤顆粒之間范德華力的減小，從而減弱土壤對(duì)膜狀水的吸持。

2.3 施硅對(duì)含鎘土壤水分特征曲線模型的影響

2.3.1 不同水分特征曲線模型的比較

使用RETC 軟件分別對(duì)各處理中3 個(gè)重復(fù)的平均值進(jìn)行VG 模型、BC 模型和DP 模型擬合，擬合結(jié)果如表4 所示。不同處理土壤在各模型擬合的表現(xiàn)并不相同，各模型擬合的不同處理下R2的平均值排序?yàn)椋篤G＞BC＞DP，不同處理下RMSE的平均值排序?yàn)椋築C＞VG＞DP，不同處理下GMER的平均值排序?yàn)椋篤G＞BC＞DP，其中R2越接近1，RMSE越小，GMER越大，模型擬合效果越好。綜合來看，VG 模型擬合效果最好，更適用于定量描述施硅對(duì)含鎘土壤水分特征曲線的影響。

表4 不同水分特征曲線模型擬合效果評(píng)價(jià)Tab.4 Evaluation of fitting performance for different models

2.3.2 Van-Genuchten模型擬合參數(shù)變化特征

利用RETC軟件對(duì)各處理的每個(gè)重復(fù)進(jìn)行VG模型擬合，并對(duì)擬合參數(shù)做方差分析，結(jié)果如表5 所示。與CK 相比，Cd2.5和Cd5處理對(duì)θs和α無顯著影響，Si2和Si4處理使θs和α顯著增大，且θs與施硅量成正比。θs是土壤飽和含水率，可以反映土壤孔隙度。α代表進(jìn)氣值的倒數(shù)，其值越大，土壤排出的水就越少［16］。這兩個(gè)水力參數(shù)對(duì)加硅的響應(yīng)與生物炭添加類似，都可以提高土壤的保水性，但是兩者原理并不相同。添加生物炭對(duì)土壤水分特征曲線的影響以物理變化為主，生物炭屬于多孔介質(zhì)，其本身可以吸持更多水分，而施硅對(duì)土壤水分特征曲線的影響既有物理變化也有化學(xué)變化。施硅對(duì)pH、zeta 電位的影響以及可溶性硅酸鹽的水解反應(yīng)均可影響土壤團(tuán)聚體的形成，改變土壤結(jié)構(gòu)，從而影響土壤水分特征曲線。與CK 相比，單獨(dú)施硅還可以顯著提升土壤殘余含水率θr，這是因?yàn)榉稚⒌耐寥滥z體顆粒比表面積大，能夠吸附更多的水分。但含鎘土壤施硅后θr顯著下降，這與小節(jié)2.1.2 土壤膜狀水吸持規(guī)律一致，進(jìn)一步說明硅-鎘配合物減小范德華力從而減少了土壤對(duì)水分的吸持。

3 結(jié) 論

本文通過壓力膜儀測(cè)量了不同硅、鎘處理土壤的水分特征曲線并進(jìn)行模型模擬及機(jī)理分析，主要結(jié)論如下。

（1）與對(duì)照組相比，鎘添加對(duì)土壤水分特征曲線沒有影響，單獨(dú)施硅以及含鎘土壤施硅會(huì)使曲線偏移，分別導(dǎo)致整體土壤含水率增加以及710～1 200 kPa吸力段的土壤含水率下降。

（2）相比于對(duì)照組，鎘添加對(duì)土壤孔隙分布沒有影響，施硅能夠增加含鎘土壤與不含鎘土壤大當(dāng)量孔徑孔隙占比，提高土壤低吸力段持水能力，同時(shí)施硅對(duì)小當(dāng)量孔徑孔隙沒有影響。

（3）施硅能夠顯著增加土壤zeta 電位絕對(duì)值，分散土壤顆粒，不利于土壤微團(tuán)粒的形成，但施硅后形成的硅酸凝膠可促進(jìn)微團(tuán)粒的膠結(jié)從而增加土壤大孔隙占比，提高低吸力范圍內(nèi)的土壤持水性；含鎘土壤施硅后形成的硅-鎘配合物可能導(dǎo)致范德華力減弱，從而降低了土壤高吸力范圍內(nèi)的持水性。

（4）相比于DP 模型與BC 模型，VG 模型能更好地描述硅、鎘對(duì)土壤水分特征曲線的影響，且施硅對(duì)含鎘土壤水力特性的影響可通過影響水分特征曲線的飽和含水率和殘余含水率進(jìn)行定量表達(dá)。

因此，通過外源硅治理水稻土鎘污染時(shí)應(yīng)當(dāng)注意施硅導(dǎo)致的土壤大孔隙占比增加從而引起潛在的面源污染物遷移問題以及硅-鎘反應(yīng)導(dǎo)致的土壤膜狀水吸持下降對(duì)作物水分利用的影響。