王鑫 郭月峰 姚云峰 徐雅潔 劉曉宇
(內(nèi)蒙古農(nóng)業(yè)大學(xué),呼和浩特,010018)
祁偉
(內(nèi)蒙古水利水電勘測(cè)設(shè)計(jì)院)
內(nèi)蒙古鄂爾多斯市準(zhǔn)格爾旗暖水鄉(xiāng)圪秋溝流域內(nèi)廣泛分布著砒砂巖[1],砒砂巖是一種多分布于晉陜蒙交界地區(qū),形成于晚古生代二疊紀(jì)、中生代三疊紀(jì)、侏羅紀(jì)和白堊紀(jì)的陸相碎屑巖[2]。該地區(qū)是黃河流域生態(tài)環(huán)境最脆弱、水土流失最嚴(yán)重的區(qū)域,屬于我國典型的風(fēng)水復(fù)合侵蝕區(qū),給當(dāng)?shù)剞r(nóng)業(yè)生產(chǎn)和生態(tài)環(huán)境帶來災(zāi)難性的影響[3]。砒砂巖抗蝕能力很弱,干旱條件下堅(jiān)硬如石;遇到降雨,砒砂巖浸水膨脹崩解、潰散成沙;遇風(fēng)則又極易剝蝕風(fēng)化[4-5]。加之該區(qū)域自然條件惡劣,降雨量較少,水分嚴(yán)重匱乏,這使植物的扎根十分困難,致使部分樹木生長不良,有枯枝現(xiàn)象,降低了植物的成活率,甚至導(dǎo)致植物大面積枯死,進(jìn)而影響了林分的可持續(xù)發(fā)展[6]。
沙棘(Hippophaerhamnoides)屬于胡頹子科(Elaeanaeeae)沙棘屬(Hippophae)的灌木或小喬木,喜陽光,不耐陰,耐干旱和貧瘠,其根系十分發(fā)達(dá),分蘗萌生能力強(qiáng),繁殖快,生物量大,因此,沙棘是黃土高原丘陵地區(qū)水土保持、防風(fēng)固沙、改善生態(tài)環(huán)境的先鋒樹種,它已成為一種兼具生態(tài)效益、經(jīng)濟(jì)效益和社會(huì)效益的樹種[7]。在砒砂巖地種植沙棘,沙棘人工林可有效降低土壤密度,增加土壤孔隙率和土壤持水率,因此,在該區(qū)域得到了大面積的種植推廣[8-9]。然而,由于內(nèi)蒙古砒砂巖地區(qū)特殊的土壤性狀和特定的地質(zhì)條件加之干旱缺水等原因,導(dǎo)致該區(qū)大面積種植的沙棘林在生長到10 a左右時(shí),開始出現(xiàn)生長衰退、生產(chǎn)力下降等現(xiàn)象。其中很重要的原因就是對(duì)該區(qū)人工植被與其周圍環(huán)境因子之間的相互關(guān)系缺乏深入探討。所以,對(duì)平茬后植物生長與環(huán)境因子的耦合效應(yīng)逐漸引起我們的關(guān)注。
細(xì)根是植物對(duì)環(huán)境變化響應(yīng)最敏感和最活躍的部分[10],其形態(tài)參數(shù)和養(yǎng)分含量不僅影響細(xì)根生理生態(tài)特性[11-12]、反映林地養(yǎng)分有效性[13],還表達(dá)了植物生存環(huán)境的變化信息[14],因此,了解細(xì)根的特性有助于更好地理解植物的功能特性和環(huán)境適應(yīng)性[15]。在全球氣候變化日益加劇的時(shí)代背景下,研究其性狀指標(biāo)對(duì)環(huán)境變化的適應(yīng)機(jī)制有重要現(xiàn)實(shí)意義[16]。其在土壤中的空間分布直接影響林木的生長發(fā)育和生產(chǎn)力,細(xì)根的生長在生長季受各生態(tài)因子的顯著影響,因此反映細(xì)根生理生態(tài)功能的各指標(biāo)也會(huì)發(fā)生變化。特別對(duì)于根系的水分和養(yǎng)分吸收能力而言,細(xì)根根長密度無疑比其他指標(biāo)更有意義[17]。
目前,對(duì)于平茬后沙棘的特性與地上部分水土保持功能的研究已很多,但對(duì)平茬后沙棘根系在其生境條件下的適應(yīng)策略研究卻很少。沙棘細(xì)根儲(chǔ)存量(DRL)除了受自身生理特性的影響,還受到外部因素大氣狀況、土壤含水率、太陽輻射、水汽壓等的影響。本文運(yùn)用統(tǒng)計(jì)分析方法對(duì)砒砂巖生境條件下平茬后沙棘細(xì)根總儲(chǔ)存量與該地區(qū)環(huán)境因子進(jìn)行研究,明確平茬后沙棘細(xì)根DRL對(duì)其生長環(huán)境的響應(yīng)機(jī)制,為砒砂巖區(qū)生態(tài)建設(shè)提供新的思路,加快該區(qū)域水土流失治理,對(duì)實(shí)現(xiàn)砒砂巖的資源利用具有重要意義。
研究區(qū)位于內(nèi)蒙古鄂爾多斯市準(zhǔn)格爾旗暖水鄉(xiāng)圪秋溝流域內(nèi),地理坐標(biāo)為39°42′~39°50′N,110°25′~110°48′E,該流域內(nèi)地形起伏較大,溝壑較多,梁卯起伏;其土壤侵蝕強(qiáng)烈,水土流失比較嚴(yán)重。平均海拔達(dá)到1 044 m,屬于典型的中溫帶半干旱大陸性季風(fēng)氣候,平均日照時(shí)間3 000 h,無霜期148 d。年平均降水量約400 mm,密集分布在7、8月份;年均蒸發(fā)量2 093 mm,年均氣溫為7.5 ℃,≥10 ℃積溫為2 900~3 500 ℃。土壤類型大部分以黃綿土為主,伴隨著砒砂巖景觀,主要為栗鈣土和風(fēng)沙土,該流域主要以具有水土保持和防風(fēng)固沙效用的人工植被為主,主要造林樹種有沙棘、油松(Pinustabuliformis)、檸條錦雞兒(Caraganakorshinskii)、紫花苜蓿(Medicagosativa)和山杏(Prunussibirica)等,人工沙棘林林下植被主要有羊草(Leymuschinensis)、克氏針茅(Stipakrylovii)、糙隱子草(Cleistogenessquarrosa)等。
于2018年選擇在內(nèi)蒙古鄂爾多斯市準(zhǔn)格爾旗暖水鄉(xiāng)13年生平茬沙棘人工林地作為試驗(yàn)地。以平茬后的沙棘為研究對(duì)象,在試驗(yàn)地上選擇立地條件、管護(hù)措施、生長狀況等基本條件一致的3塊平茬沙棘人工林地作為試驗(yàn)樣地,3塊試驗(yàn)樣地面積均為50 m×50 m,坡度小于5°,試驗(yàn)樣地內(nèi)沙棘株距×行距為2 m×4 m。為了更好地研究平茬后沙棘細(xì)根與其生境環(huán)境因子的關(guān)系,選擇全面平茬,留茬高度為10 cm。2018年3月份中旬對(duì)3塊試驗(yàn)樣地的沙棘進(jìn)行平茬,于2018年、2019年4—10月每月中旬對(duì)沙棘進(jìn)行測(cè)定。微根管的安裝與平茬同時(shí)進(jìn)行,為避免邊緣效應(yīng),將微根管集中安裝在試驗(yàn)樣地中心。所以,本試驗(yàn)在每個(gè)平茬樣地內(nèi)選擇3叢沙棘標(biāo)準(zhǔn)叢(株高83 cm、南北冠幅108 cm、東西冠幅130 cm、平均分支數(shù)為7);在距沙棘0 cm處安裝微根管,3個(gè)平茬樣地共計(jì)9根,在深度為1 m、半徑為2 m范圍內(nèi)采集根系信息。樣地概況見表1。
表1 沙棘樣地統(tǒng)計(jì)
根長密度表示單位土體內(nèi)植物根系的長度之和,根長密度越大表明根系占據(jù)的土壤空間越大,更有利于植物根系吸收土壤中的水分和養(yǎng)分。為更好的研究平茬后沙棘細(xì)根生長狀況,此次研究于2018年、2019年4—10月份每月中旬測(cè)定,每月數(shù)據(jù)為2 a數(shù)據(jù)均值,采用根系分析軟件RootSnap測(cè)定<2 mm的沙棘細(xì)根,計(jì)算細(xì)根根長密度。本文以單位土體的根長密度作為基本參數(shù),具體計(jì)算如下[18]:
DRL=(LR/A)×FDO。
式中:DRL為根長密度(mm·cm-3);A為觀測(cè)窗面積(cm2);LR為檢測(cè)到的細(xì)根根長(mm);FDO為田間觀測(cè)深度(cm)。各根樣長度、表面積、體積分別除以取樣土體體積計(jì)算出根長密度。
在林地附近的空曠地布設(shè)美國研發(fā)的HOBO小型自動(dòng)氣象站,設(shè)置數(shù)據(jù)采集時(shí)間間隔為10 min,定期下載觀測(cè)數(shù)據(jù),連續(xù)測(cè)定周圍2018年、2019年4—10月份的環(huán)境因子,數(shù)據(jù)為每月數(shù)據(jù)平均值。觀測(cè)的環(huán)境因子主要包括X1為空氣溫度(Ta);X2為空氣相對(duì)濕度(HR);X3為水汽壓(DVP);X4為光合有效輻射(RPA);X5為土壤含水率(%)。土壤含水率利用中子水分測(cè)定儀測(cè)定。
數(shù)據(jù)整理及制圖在Microsoft Excel 2016中進(jìn)行,數(shù)據(jù)分析利用SPSS 26、DPS進(jìn)行顯著性檢驗(yàn)(顯著性水平設(shè)為0.05)、單因素方差分析(ANOVA)、共線性診斷及嶺回歸分析。
4—10月份沙棘細(xì)根DRL在垂直深度上表現(xiàn)一定的分布特征(表2),DRL在垂直深度上的季節(jié)變化特征為隨土層深度加深表現(xiàn)出逐漸降低的趨勢(shì),各土層平均DRL由大到小順序?yàn)镈RL(0~15 cm)、DRL(15~30 cm)、DRL(30~45 cm)、DRL(45~60 cm)、DRL(60~100 cm),細(xì)根大部分分布在0~30 cm范圍內(nèi),0~30 cm土層是沙棘細(xì)根的主要分布區(qū)和生長活躍區(qū),采用方差分析對(duì)不同土層以及不同月份之間DRL的差異進(jìn)行顯著性檢驗(yàn);經(jīng)過方差分析表明,各土層間DRL差異顯著(P<0.05),DRL季節(jié)動(dòng)態(tài)變化隨土層深度變化表現(xiàn)出顯著性差異(P<0.05)。通過該區(qū)4—10月份根長密度分布圖我們可以看到在0~100 cm土層范圍內(nèi),月平均DRL的波動(dòng)范圍為0.59~7.23 mm·cm-2,各土層最小值均出現(xiàn)在4月份,最低為0.59 mm·cm-2;各土層DRL均在生長季8月份達(dá)到峰值,最高為7.23 mm·cm-2。
表2 沙棘細(xì)根根長密度分布特征
為直觀的看出沙棘細(xì)根生長與其周圍環(huán)境因子的變化特征,選擇空氣溫度(Ta)、空氣相對(duì)濕度(HR)、水汽壓(DVP)、光合有效輻射(RPA)、土壤含水率等指標(biāo)分別與4—10月份沙棘細(xì)根現(xiàn)存量繪制變化曲線。
由圖1a可知,沙棘細(xì)根DRL與Ta的季節(jié)變化趨勢(shì)大致相同,其中平均氣溫極小值出現(xiàn)在4月份(5.32 ℃),4月份之后氣溫開始回升,受土壤有機(jī)質(zhì)活性影響,此時(shí)細(xì)根生長緩慢,而在5—8月份生長季中光和有效輻射充足,隨著雨季的到來,土壤有機(jī)質(zhì)活性增加,因此平均氣溫極大值在8月份達(dá)到最高(23.67 ℃),最高氣溫較DRL更早達(dá)到峰值,到了9月份,隨著氣溫的下降,細(xì)根生長緩慢。因此,DRL最大值一般會(huì)出現(xiàn)在秋季,最小值出現(xiàn)在春季。
由圖1b可知,HR的變化趨勢(shì)與沙棘細(xì)根DRL的動(dòng)態(tài)變化趨勢(shì)呈負(fù)相關(guān),即沙棘細(xì)根DRL隨著HR的減小而呈逐漸上升的趨勢(shì)。進(jìn)入4月份后隨著太陽輻射強(qiáng)度與氣溫的升高,HR逐漸減小,沙棘細(xì)根DRL持續(xù)上升,HR在8月份左右降為最低值,約為47.69%,較5月份降低了15.18%,此時(shí)沙棘細(xì)根DRL達(dá)到峰值。
由圖1c可知,隨著生長季內(nèi)物候變化,月平均DVP總體表現(xiàn)為前期(4—6月份)較低、中后期(7—8月份)較高、末期(9—10月份)迅速下降的變化趨勢(shì)。隨著氣溫和空氣相對(duì)濕度的變化,作為蒸騰作用主要驅(qū)動(dòng)力的DVP在初夏季節(jié)(6月份)出現(xiàn)迅速上升,月均上升速度最快時(shí)月提高了0.61 kPa,進(jìn)入7月份上升趨勢(shì)有所減緩,在8月份平均DVP達(dá)到峰值(1.81 kPa),9月份月平均DVP下降速度最快,較8月份降低了45%;因此,表現(xiàn)出月平均DVP與沙棘細(xì)根DRL呈正相關(guān)。
如圖1d所示是沙棘細(xì)根DRL與RPA的連日變化,可以看出沙棘細(xì)根DRL與RPA的連月變化呈正相關(guān),因此,RPA的最小值出現(xiàn)在旱季(4月份),最大值則出現(xiàn)在雨季(8月份),該區(qū)域7—9月份為雨季,降雨次數(shù)較多,陰雨天較其他月份多,因此RPA增長量沒有明顯大幅度變化趨勢(shì)。
如圖1e所示,沙棘細(xì)根DRL的季節(jié)變化趨勢(shì)大致與土壤含水率變化趨勢(shì)相近,均呈現(xiàn)單峰變化趨勢(shì)。沙棘細(xì)根土壤含水率與沙棘細(xì)根DRL呈正相關(guān)關(guān)系,其中4—8月份為土壤含水率增長階段,從5.2%增長至9.8%,8月份之后土壤含水率開始下降,從圖中我們可以看到細(xì)根DRL在5、6月份增長速率最大,此時(shí)土壤含水率維持在8%左右。
圖1 不同環(huán)境因子對(duì)細(xì)根根系現(xiàn)存量影響的季節(jié)動(dòng)態(tài)
3.3.1 建立多元線性回歸模型
沙棘DRL不僅受植物自身的生理特性的影響,與其周圍環(huán)境因子也有密切的聯(lián)系,所以理論上沙棘細(xì)根生長有很大部分取決于植物周圍環(huán)境因子的影響,由此建立沙棘DRL回歸模型:以沙棘DRL(Y)為因變量,選取Ta(X1)、HR(X2)、DVP(X3)、RPA(X4)、土壤含水量(X5)為自變量。對(duì)因變量Y與各自變量之間進(jìn)行相關(guān)性分析,結(jié)果表明Y與X1、X2、X3、X4、X5都存在顯著相關(guān)性,相關(guān)系數(shù)分別為0.755、-0.859、0.790、0.761、0.923,其中Y與X5的相關(guān)系數(shù)最大,為0.923,總之,用Y與這些自變量做多元線性回歸是合理的。對(duì)Y與5個(gè)自變量進(jìn)行線性回歸分析,結(jié)果見表3、表4,由此得到回歸方程。
表3 變量系數(shù)
表4 單因素方差分析
Y=0.546-0.544X1-30.530X2+2.152X3+0.037X4+
217.377X5。
根據(jù)擬合出來的模型結(jié)果,可決系數(shù)(R2)為0.999,調(diào)整后的可決系數(shù)(R2)為0.995,模型的擬合效果較好;模型的整體顯著性F檢驗(yàn)統(tǒng)計(jì)量為258.907,概率P值為0.047(小于0.05),但用所有自變量做回歸效果不好。
3.3.2 多重共線性診斷
通過表3我們可以看到5個(gè)自變量的方差膨脹系數(shù)分別為50.487、6.032、25.174、17.566、4.254。除X2和X5外,其他變量的方差擴(kuò)大因子都大于10,即變量之間存在嚴(yán)重的多重共線性。表中最大的條件數(shù)為217.377,自變量間存在多重共線性。當(dāng)自變量之間存在較強(qiáng)的多重共線性時(shí),求得的多重線性回歸模型很不穩(wěn)定。此時(shí),嶺回歸分析有可能較好地解決前述提及的問題。
3.3.3 沙棘細(xì)根現(xiàn)存量與環(huán)境因子嶺回歸分析
首先對(duì)原始數(shù)據(jù)進(jìn)行中心化和標(biāo)準(zhǔn)化,首次回歸分析,將嶺參數(shù)k默認(rèn)從0到1,步長為0.01,共有101個(gè)值。
嶺回歸的概念:
(1)對(duì)全部的5個(gè)自變量做嶺回歸分析,并得到嶺跡圖,見圖2。
X1為Ta,X2為HR,X3為DVP、X4為RPA、X5為土壤含水率。
(2)經(jīng)過這樣的分析,保留4個(gè)變量(X1、X2、X4、X5),再做嶺跡圖(圖3)。
X1為Ta,X2為HR,X3為RPA、X4為土壤含水率。
(3)對(duì)剩余的4個(gè)變量重新進(jìn)行嶺跡分析,由于X3的絕對(duì)值很小,并且它與X1的相關(guān)系數(shù)也比較高,因此考慮剔除掉X3;經(jīng)過這樣的分析,最后保留3個(gè)變量(X1、X4、X5),再做嶺跡圖(圖4)。
在嶺跡圖上,k=0.21以后,各參數(shù)逐漸穩(wěn)定。因而選取k=0.21,計(jì)算此時(shí)的嶺回歸估計(jì)結(jié)果為:復(fù)決定系數(shù)R2=0.941,調(diào)整相關(guān)R2=0.940,F(xiàn)=15.935,P=0.024<0.05呈顯著相關(guān),可以說回歸方程擬合效果不錯(cuò),各解釋變量的系數(shù)也較為合理了。
綜上得到嶺回歸方程為:
Y=19.825+0.004X1-31.539X2+120.560X5。
運(yùn)用嶺回歸分析方法成功解決了沙棘細(xì)根現(xiàn)存量影響因素分析的回歸建模問題。采用它可以通過允許小的偏差而換取高于無偏估計(jì)量的精度,因?yàn)樗咏鎸?shí)值的可能性較大。
X1為Ta,X2為HR,X3為土壤含水率。
王亞楠[19]提出在0~30 cm土壤潮濕,并且溫度較高,這樣的環(huán)境適合土壤微生物生存,因此土壤養(yǎng)分的有效性較高。本研究中我們得到平茬后沙棘細(xì)根DRL生長季內(nèi)在0~30 cm分布最為集中,這可能因?yàn)?~30 cm土壤養(yǎng)分的有效性高,所以細(xì)根DRL較大,同樣證實(shí)了王亞楠提出的觀點(diǎn)。
馬明呈[20]通過研究季節(jié)變化對(duì)沙棘根系生長的影響發(fā)現(xiàn),在濕潤季節(jié),表土層(0~15 cm)中細(xì)根分布比例顯著高于干旱季節(jié),這一現(xiàn)象可能與細(xì)根生長與含水量升高(雨季開始)相關(guān)。本研究表明平茬后沙棘細(xì)根DRL的季節(jié)變化呈秋季單峰型,平茬后沙棘細(xì)根DRL隨季節(jié)變化與土層深度存在一定負(fù)相關(guān)性,隨著季節(jié)變化,沙棘細(xì)根DRL均隨著土層深度的增加而逐漸減少;本研究區(qū)屬于典型的砒砂巖區(qū),因此氣候干燥,年均降水量少,降雨密集分布在7、8月份;春季為旱季,溫度較低,并且4月份氣溫開始回升,受到土壤有機(jī)質(zhì)活性影響,細(xì)根生長較為緩慢,因此各土層最小值均出現(xiàn)在4月份,夏季溫度升高,隨著雨季的到來,土壤有機(jī)質(zhì)活性增加,DRL也相應(yīng)增加;7、8月份沙棘細(xì)根DRL較高,也可能是因?yàn)樵谠摃r(shí)期土壤水分、空氣溫度、濕度、光照等生長條件良好,植物為了滿足地上部分的生長需要做出積極響應(yīng),使得細(xì)根大量生長,吸收更多的水分和養(yǎng)分。
沙棘細(xì)根生長不僅受自身生理特征的影響,周圍環(huán)境因子對(duì)沙棘生長發(fā)育也有很大的影響,所以理論上平茬后沙棘細(xì)根DRL有很大部分取決于植物周圍環(huán)境因子的影響;DRL的季節(jié)動(dòng)態(tài)受樹木生長的生理需求、根系本身的物候特征、環(huán)境條件(氣候特征)的綜合影響,發(fā)現(xiàn)細(xì)根生物量的季節(jié)變化與枝葉物候節(jié)律一致[21]。本研究表明,平茬后沙棘細(xì)根DRL與Ta、HR、DVP、RPA、土壤含水率等指標(biāo)均具有顯著相關(guān)性,Ta、DVP、RPA、土壤含水率與沙棘細(xì)根DRL呈正相關(guān),而HR與沙棘細(xì)根DRL呈負(fù)相關(guān),8月份前隨著HR降低,DRL逐漸升高,8月份之后則呈相反趨勢(shì)。RPA的強(qiáng)度可誘導(dǎo)氣孔啟閉,同時(shí)又直接影響著Ta、DVP等環(huán)境因子的變化,因此Ta、DVP、RPA均與沙棘細(xì)根DRL呈正相關(guān)關(guān)系,說明沙棘細(xì)根DRL受其周圍環(huán)境因子影響,但各環(huán)境因子的影響程度不同,想要對(duì)其深入了解,需作進(jìn)一步研究。
對(duì)沙棘細(xì)根DRL與環(huán)境因子建模我們發(fā)現(xiàn)回歸模型與實(shí)際情況不相符和,其回歸系數(shù)不穩(wěn)定,通過方差擴(kuò)大因子分析我們發(fā)現(xiàn)環(huán)境因子之間相互影響,使得回歸模型不準(zhǔn)確,因此我們采用嶺回歸對(duì)各自變量進(jìn)行分析。嶺回歸分析結(jié)果剔除了共線性嚴(yán)重的變量X3、X4,得到的嶺回歸方程符合實(shí)際情況,回歸方程也比較穩(wěn)定,根據(jù)分析結(jié)果最終選取Ta、HR、土壤含水率作為DRL的解釋變量,其對(duì)平茬后沙棘DRL影響程度由大到小為土壤含水率、HR、Ta。因此,7、8月份細(xì)根DRL較高可能主要是因?yàn)榇藭r(shí)土壤含水率較高,為滿足自身的生長需求,細(xì)根大量生長以便吸收水分和養(yǎng)分。在其他變量不變的情況下,空氣溫度每增加1%,沙棘細(xì)根長密度增加0.004%;空氣相對(duì)濕度每增加1%,沙棘細(xì)根長密度降低31.539%;土壤含水率每增加1%,沙棘細(xì)根長密度增加120.560%。嶺回歸法解決多重共線性問題有其獨(dú)到之處,與其他方法不盡相同。無論采取什么方法,都應(yīng)從實(shí)際情況出發(fā),選擇對(duì)解決實(shí)際問題有利而簡(jiǎn)單的方法,不僅可以對(duì)分析各變量之間的作用和聯(lián)系帶來意想不到的幫助,而且可以達(dá)到事半功倍的效果[22]。
通過以上分析得出如下結(jié)論:(1)平茬后沙棘細(xì)根DRL生長季內(nèi)在0~30 cm分布最為集中,細(xì)根DRL的季節(jié)變化呈秋季單峰型,先增加后減少,在7、8月份平茬后沙棘細(xì)根DRL達(dá)到最大值。(2)平茬后沙棘細(xì)根DRL與Ta、HR、DVP、RPA、土壤含水率等指標(biāo)均具有顯著相關(guān)性,Ta、DVP、RPA、土壤含水率與沙棘細(xì)根DRL呈正相關(guān),而HR與沙棘細(xì)根DRL呈負(fù)相關(guān)。(3)平茬后沙棘細(xì)根DRL與各環(huán)境因子的嶺回歸方程為:Y=19.825+0.004X1-31.539X2+120.560X5,其中Y為平茬后沙棘細(xì)根長密度(DRL);X1為空氣溫度(Ta);X2為空氣相對(duì)濕度(HR);X5為土壤含水率;根據(jù)分析結(jié)果最終選取Ta、HR、土壤含水率作為DRL的解釋變量,其對(duì)平茬后沙棘DRL影響程度由大到小為土壤含水率、HR、Ta。