環(huán)境因子變化對砒砂巖區(qū)平茬沙棘細(xì)根現(xiàn)存量的影響1)

王鑫 郭月峰 姚云峰 徐雅潔 劉曉宇

(內(nèi)蒙古農(nóng)業(yè)大學(xué)，呼和浩特，010018)

祁偉

(內(nèi)蒙古水利水電勘測設(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)象。其中很重要的原因就是對該區(qū)人工植被與其周圍環(huán)境因子之間的相互關(guān)系缺乏深入探討。所以，對平茬后植物生長與環(huán)境因子的耦合效應(yīng)逐漸引起我們的關(guān)注。

細(xì)根是植物對環(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)對環(huán)境變化的適應(yīng)機(jī)制有重要現(xiàn)實(shí)意義[16]。其在土壤中的空間分布直接影響林木的生長發(fā)育和生產(chǎn)力，細(xì)根的生長在生長季受各生態(tài)因子的顯著影響，因此反映細(xì)根生理生態(tài)功能的各指標(biāo)也會(huì)發(fā)生變化。特別對于根系的水分和養(yǎng)分吸收能力而言，細(xì)根根長密度無疑比其他指標(biāo)更有意義[17]。

目前，對于平茬后沙棘的特性與地上部分水土保持功能的研究已很多，但對平茬后沙棘根系在其生境條件下的適應(yīng)策略研究卻很少。沙棘細(xì)根儲(chǔ)存量(DRL)除了受自身生理特性的影響，還受到外部因素大氣狀況、土壤含水率、太陽輻射、水汽壓等的影響。本文運(yùn)用統(tǒng)計(jì)分析方法對砒砂巖生境條件下平茬后沙棘細(xì)根總儲(chǔ)存量與該地區(qū)環(huán)境因子進(jìn)行研究，明確平茬后沙棘細(xì)根DRL對其生長環(huán)境的響應(yīng)機(jī)制，為砒砂巖區(qū)生態(tài)建設(shè)提供新的思路，加快該區(qū)域水土流失治理，對實(shí)現(xiàn)砒砂巖的資源利用具有重要意義。

1 研究區(qū)概況

研究區(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)等。

2 研究方法

2.1 調(diào)查與取樣

于2018年選擇在內(nèi)蒙古鄂爾多斯市準(zhǔn)格爾旗暖水鄉(xiāng)13年生平茬沙棘人工林地作為試驗(yàn)地。以平茬后的沙棘為研究對象，在試驗(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月份中旬對3塊試驗(yàn)樣地的沙棘進(jìn)行平茬，于2018年、2019年4—10月每月中旬對沙棘進(jìn)行測定。微根管的安裝與平茬同時(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ì)

2.2 沙棘細(xì)根長密度的測定

根長密度表示單位土體內(nèi)植物根系的長度之和，根長密度越大表明根系占據(jù)的土壤空間越大，更有利于植物根系吸收土壤中的水分和養(yǎng)分。為更好的研究平茬后沙棘細(xì)根生長狀況，此次研究于2018年、2019年4—10月份每月中旬測定，每月數(shù)據(jù)為2 a數(shù)據(jù)均值，采用根系分析軟件RootSnap測定<2 mm的沙棘細(xì)根，計(jì)算細(xì)根根長密度。本文以單位土體的根長密度作為基本參數(shù)，具體計(jì)算如下[18]：

DRL=(LR/A)×FDO。

式中：DRL為根長密度(mm·cm-3)；A為觀測窗面積(cm2)；LR為檢測到的細(xì)根根長(mm)；FDO為田間觀測深度(cm)。各根樣長度、表面積、體積分別除以取樣土體體積計(jì)算出根長密度。

2.3 環(huán)境因子的測定

在林地附近的空曠地布設(shè)美國研發(fā)的HOBO小型自動(dòng)氣象站，設(shè)置數(shù)據(jù)采集時(shí)間間隔為10 min，定期下載觀測數(shù)據(jù)，連續(xù)測定周圍2018年、2019年4—10月份的環(huán)境因子，數(shù)據(jù)為每月數(shù)據(jù)平均值。觀測的環(huán)境因子主要包括X1為空氣溫度(Ta)；X2為空氣相對濕度(HR)；X3為水汽壓(DVP)；X4為光合有效輻射(RPA)；X5為土壤含水率(%)。土壤含水率利用中子水分測定儀測定。

2.4 數(shù)據(jù)分析

數(shù)據(jù)整理及制圖在Microsoft Excel 2016中進(jìn)行，數(shù)據(jù)分析利用SPSS 26、DPS進(jìn)行顯著性檢驗(yàn)(顯著性水平設(shè)為0.05)、單因素方差分析(ANOVA)、共線性診斷及嶺回歸分析。

3 結(jié)果與分析

3.1 沙棘細(xì)根現(xiàn)存量變化特征

4—10月份沙棘細(xì)根DRL在垂直深度上表現(xiàn)一定的分布特征(表2)，DRL在垂直深度上的季節(jié)變化特征為隨土層深度加深表現(xiàn)出逐漸降低的趨勢，各土層平均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ū)，采用方差分析對不同土層以及不同月份之間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ì)根根長密度分布特征

3.2 沙棘細(xì)根現(xiàn)存量與環(huán)境因子的關(guān)系

為直觀的看出沙棘細(xì)根生長與其周圍環(huán)境因子的變化特征，選擇空氣溫度(Ta)、空氣相對濕度(HR)、水汽壓(DVP)、光合有效輻射(RPA)、土壤含水率等指標(biāo)分別與4—10月份沙棘細(xì)根現(xiàn)存量繪制變化曲線。

由圖1a可知，沙棘細(xì)根DRL與Ta的季節(jié)變化趨勢大致相同，其中平均氣溫極小值出現(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的變化趨勢與沙棘細(xì)根DRL的動(dòng)態(tài)變化趨勢呈負(fù)相關(guān)，即沙棘細(xì)根DRL隨著HR的減小而呈逐漸上升的趨勢。進(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月份)迅速下降的變化趨勢。隨著氣溫和空氣相對濕度的變化，作為蒸騰作用主要驅(qū)動(dòng)力的DVP在初夏季節(jié)(6月份)出現(xiàn)迅速上升，月均上升速度最快時(shí)月提高了0.61 kPa，進(jìn)入7月份上升趨勢有所減緩，在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增長量沒有明顯大幅度變化趨勢。

如圖1e所示，沙棘細(xì)根DRL的季節(jié)變化趨勢大致與土壤含水率變化趨勢相近，均呈現(xiàn)單峰變化趨勢。沙棘細(xì)根土壤含水率與沙棘細(xì)根DRL呈正相關(guān)關(guān)系，其中4—8月份為土壤含水率增長階段，從5.2%增長至9.8%，8月份之后土壤含水率開始下降，從圖中我們可以看到細(xì)根DRL在5、6月份增長速率最大，此時(shí)土壤含水率維持在8%左右。

圖1 不同環(huán)境因子對細(xì)根根系現(xiàn)存量影響的季節(jié)動(dòng)態(tài)

3.3 沙棘細(xì)根現(xiàn)存量與環(huán)境因子回歸分析

3.3.1 建立多元線性回歸模型

沙棘DRL不僅受植物自身的生理特性的影響，與其周圍環(huán)境因子也有密切的聯(lián)系，所以理論上沙棘細(xì)根生長有很大部分取決于植物周圍環(huán)境因子的影響，由此建立沙棘DRL回歸模型：以沙棘DRL(Y)為因變量，選取Ta(X1)、HR(X2)、DVP(X3)、RPA(X4)、土壤含水量(X5)為自變量。對因變量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與這些自變量做多元線性回歸是合理的。對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)境因子嶺回歸分析

首先對原始數(shù)據(jù)進(jìn)行中心化和標(biāo)準(zhǔn)化，首次回歸分析，將嶺參數(shù)k默認(rèn)從0到1，步長為0.01，共有101個(gè)值。

嶺回歸的概念：

(1)對全部的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)對剩余的4個(gè)變量重新進(jìn)行嶺跡分析，由于X3的絕對值很小，并且它與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為土壤含水率。

4 結(jié)論與討論

王亞楠[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é)變化對沙棘根系生長的影響發(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)境因子對沙棘生長發(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月份之后則呈相反趨勢。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)境因子的影響程度不同，想要對其深入了解，需作進(jìn)一步研究。

對沙棘細(xì)根DRL與環(huán)境因子建模我們發(fā)現(xiàn)回歸模型與實(shí)際情況不相符和，其回歸系數(shù)不穩(wěn)定，通過方差擴(kuò)大因子分析我們發(fā)現(xiàn)環(huán)境因子之間相互影響，使得回歸模型不準(zhǔn)確，因此我們采用嶺回歸對各自變量進(jìn)行分析。嶺回歸分析結(jié)果剔除了共線性嚴(yán)重的變量X3、X4，得到的嶺回歸方程符合實(shí)際情況，回歸方程也比較穩(wěn)定，根據(jù)分析結(jié)果最終選取Ta、HR、土壤含水率作為DRL的解釋變量，其對平茬后沙棘DRL影響程度由大到小為土壤含水率、HR、Ta。因此，7、8月份細(xì)根DRL較高可能主要是因?yàn)榇藭r(shí)土壤含水率較高，為滿足自身的生長需求，細(xì)根大量生長以便吸收水分和養(yǎng)分。在其他變量不變的情況下，空氣溫度每增加1%，沙棘細(xì)根長密度增加0.004%；空氣相對濕度每增加1%，沙棘細(xì)根長密度降低31.539%；土壤含水率每增加1%，沙棘細(xì)根長密度增加120.560%。嶺回歸法解決多重共線性問題有其獨(dú)到之處，與其他方法不盡相同。無論采取什么方法，都應(yīng)從實(shí)際情況出發(fā)，選擇對解決實(shí)際問題有利而簡單的方法，不僅可以對分析各變量之間的作用和聯(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為空氣相對濕度(HR)；X5為土壤含水率；根據(jù)分析結(jié)果最終選取Ta、HR、土壤含水率作為DRL的解釋變量，其對平茬后沙棘DRL影響程度由大到小為土壤含水率、HR、Ta。