何釗全 ，尚雪，張銅會，云建英

（1.陜西省紅棗重點實驗室延安大學(xué)，陜西 延安 716000；2.延安大學(xué)生命科學(xué)學(xué)院，陜西 延安 716000；3.中國科學(xué)院西北生態(tài)環(huán)境資源研究院，蘭州 730030）

0 引言

【研究意義】棗樹（Jujube）原產(chǎn)于我國，是我國特有經(jīng)濟樹種之一，具有防風(fēng)、固沙和減輕干熱風(fēng)等作用。陜北黃土高原黃河沿岸是我國優(yōu)良制干紅棗中陽木棗的主要栽培區(qū)域和適生區(qū)之一。陜北半干旱區(qū)夏季炎熱，日照較強、秋季有短促暴雨、冬季干冷，所產(chǎn)紅棗果大、營養(yǎng)豐富、口味甘甜、富含維生素，是一種保健果品[1]。近年來，紅棗種植面積持續(xù)擴大。截至2018年，陜北紅棗種植面積19×104hm2左右，產(chǎn)量74.87×104t，成為當?shù)氐闹еa(chǎn)業(yè)，經(jīng)濟、社會和生態(tài)效益顯著[2]。進而也出現(xiàn)了佳縣、延川和清澗等非常著名的全國紅棗重點縣。比如，1988年，陜西佳縣被命名為“中國紅棗名鄉(xiāng)”，是我國唯一認證的有機紅棗產(chǎn)地，具有“天然極品、百果之王”的美譽[3]。因此，探究陜北紅棗的生態(tài)適宜性區(qū)劃特征及其驅(qū)動機理，對發(fā)展區(qū)域紅棗特色產(chǎn)業(yè)具有現(xiàn)實意義。

【研究進展】近年來，該區(qū)域氣候因子、植被狀況等變化顯著，加之地貌形態(tài)、地形結(jié)構(gòu)復(fù)雜多變，會影響光、熱、水、氣等條件，對紅棗的生長發(fā)育、產(chǎn)量及其品質(zhì)產(chǎn)生不同程度的影響，最終嚴重阻礙紅棗產(chǎn)業(yè)的可持續(xù)發(fā)展[4]。國內(nèi)很多學(xué)者采用不同的統(tǒng)計分析方法對該地區(qū)氣候資源變化、氣溫降水時空變異、植被土壤特征等進行了研究，較全面地揭示了該地區(qū)的氣候地理變化趨勢。比如，一方面，1960—2010年陜北地區(qū)年均氣溫和月均氣溫呈上升趨勢，變化傾向率達0.30 ℃/10a，年降水量呈下降趨勢，變化傾向率為-12.46 mm/10a[5]，導(dǎo)致陜北產(chǎn)生輕度干旱的頻率高達68.20%，因此，氣候總體上呈暖干化趨勢，但是增長速率逐步降低[6]。另一方面，自1999年我國啟動退耕還林（草）的工程，至2015年，黃土高原造林面積已達8 954 km2[7]。在此背景下，2005—2019年，陜北地區(qū)增強型植被指數(shù)（EVI）空間分布差異明顯，總體上表現(xiàn)為，東南丘陵溝壑區(qū)植被覆蓋度較高，西北風(fēng)沙區(qū)及城鎮(zhèn)建設(shè)用地區(qū)域植被覆蓋度較低，陜北地區(qū)的植被覆蓋整體呈顯著上升趨勢，植被增加的區(qū)域面積占區(qū)域面積的64.30%，植被不變的區(qū)域占31.10%，而植被減少的區(qū)域僅占4.60%[8-9]。此外，陜北黃土地貌正負地形存在明顯差異，地面坡度的個體和空間差異反映了黃土地貌形態(tài)的空間分異[10]。研究[6,8,10]表明，黃土高原不同地貌類型在結(jié)構(gòu)和分布上各有特點，尤其是以塬、梁、峁為典型代表的黃土地貌在地貌形態(tài)、地形結(jié)構(gòu)和水系發(fā)育程度等方面表現(xiàn)出明顯的地區(qū)性特征。然而，黃土地貌形態(tài)極其復(fù)雜，具有不規(guī)則性，許多研究指標尺度依賴性大，不同尺度研究結(jié)果具有明顯差異性，是典型的土壤退化與水土流失區(qū)域。目前，通過對該區(qū)土地利用/覆被變化的氣候驅(qū)動因素、土壤侵蝕、生態(tài)系統(tǒng)生產(chǎn)力變化及脆弱生態(tài)的綜合整治，生態(tài)脆弱區(qū)條件得到顯著改善[11]。適宜的種植區(qū)域是保證結(jié)構(gòu)優(yōu)化和促進產(chǎn)業(yè)可持續(xù)發(fā)展的基礎(chǔ)。加之植被、地形生態(tài)類別多樣，氣候因子變化頻繁，紅棗的生長狀況、品質(zhì)產(chǎn)量會在不同區(qū)域有顯著異同。

生態(tài)適宜性概念最早由美國景觀規(guī)劃師宗黎黎[4]提出，并隨著遙感和地理信息系統(tǒng)（GIS）技術(shù)的持續(xù)應(yīng)用和數(shù)據(jù)統(tǒng)計方法深入探究，使生態(tài)適宜性指標類別日趨系統(tǒng)和全面[12]。GIS 技術(shù)在農(nóng)業(yè)、林業(yè)、藥業(yè)等多領(lǐng)域廣泛應(yīng)用，評價對象愈發(fā)細化，且提高了區(qū)劃和評估結(jié)果的精度。整體上看，近年來，學(xué)者們針對各類作物生態(tài)適宜性評價展開了大量研究，研究方向主要集中在方法、模型及應(yīng)用性研究方面。比如，王珊等[13]利用耕地地力評價結(jié)果，對高縣玉米種植區(qū)土地進行了生態(tài)適宜性評價，明確了不同生態(tài)適宜性等級的區(qū)域分布面積，確定了生育期間≥10℃積溫和降水量為影響該區(qū)域玉米種植適宜性的主要指標，為該縣今后的玉米種植生態(tài)區(qū)劃提供科學(xué)依據(jù)。賀征兵等[14]基于陜西省黃龍縣氣象資料，數(shù)字高程模型（DEM）以及土壤數(shù)據(jù)，利用GIS 逆距離權(quán)重法的綜合插值方法對核桃種植生態(tài)氣候適宜性進行評價，確定了種植適宜區(qū)，以期在核桃栽培過程中充分發(fā)揮生態(tài)氣候資源合理利用的優(yōu)勢。莫建國等[15]基于GIS技術(shù)，利用氣候因素（平均氣溫、降水量及日照時間），確立了獼猴桃種植的氣候適宜性區(qū)劃，為河源地區(qū)進一步擴大種植面積和提高品質(zhì)提供理論支撐。孫園園等[16]以四川省稻作區(qū)為分析樣區(qū)，利用GIS 空間差值法，研究了優(yōu)質(zhì)稻生產(chǎn)的氣候生態(tài)條件適宜性及其空間分布格局，為科學(xué)規(guī)劃優(yōu)化優(yōu)質(zhì)稻生產(chǎn)布局和結(jié)構(gòu)奠定案例基礎(chǔ)。王迅等[17]通過分析四川晚熟柑橘適宜性區(qū)劃，結(jié)合對代表性產(chǎn)區(qū)的實地調(diào)研，分析出四川發(fā)展晚熟柑橘的氣候優(yōu)勢點，并得到了晚熟柑橘氣候適宜性精細區(qū)劃圖。上述分析結(jié)論更加充分驗證了作物生態(tài)適宜性評價的必要性和現(xiàn)實意義。但是，評價因子側(cè)重點的差異性會影響評價方法的不同，最終導(dǎo)致評價方向和適用性范疇的不同。

【切入點】鑒于氣候、植被和立地條件是影響紅棗生長的關(guān)鍵因素，而人們對陜北紅棗生長在不同區(qū)域上的氣候、地理等限制性因子認識有限，對具體的立地環(huán)境缺少研究，勢必影響了紅棗產(chǎn)業(yè)的持續(xù)發(fā)展。同時，氣候、地形、土壤等與紅棗生長的關(guān)系極力促進紅棗生態(tài)適宜性精細區(qū)劃和評價?！緮M解決的關(guān)鍵問題】因此，本文采用GIS 技術(shù)和層次分析法，探討陜北地區(qū)紅棗分布與生態(tài)環(huán)境之間的關(guān)系，對紅棗生態(tài)適宜區(qū)進行等級劃分，確定紅棗生態(tài)適宜性區(qū)劃圖；且通過地理探測器模型和生態(tài)學(xué)統(tǒng)計方法篩選影響紅棗生態(tài)適宜性的主要生態(tài)因子及其影響機理。研究結(jié)果可為陜北紅棗合理引種栽培提供理論依據(jù)。

1 材料與方法

1.1 研究區(qū)概況

陜北地區(qū)位于黃土高原中部，地處35°02′—39°35′N，107°15′—110°15′E，總面積92 521.40 km2（圖1）。年均降水量約350～600 mm，年平均氣溫8～12 ℃，為典型的半干旱地區(qū)。土壤以栗鈣土、黃綿土為主。該區(qū)域平均海拔高，氣候干燥，且光照充足、晝夜溫差大，是我國優(yōu)質(zhì)紅棗的重要產(chǎn)區(qū)。

圖1 研究區(qū)位置Fig.1 Location of the study area

1.2 數(shù)據(jù)來源

本文采用的數(shù)據(jù)主要有數(shù)字高程數(shù)據(jù)（DEM）、研究區(qū)縣級行政矢量邊界、年平均氣溫、年降水量、年積溫（≥10 ℃）、土地利用數(shù)據(jù)。其中，DEM收集于地理空間數(shù)據(jù)云（http://www.gscloud.cn/）；研究區(qū)縣級行政矢量邊界來源于地球大數(shù)據(jù)科學(xué)工程數(shù)據(jù)共享服務(wù)系統(tǒng)；年平均氣溫、年降水量、年積溫（≥10 ℃）及土壤類型等柵格數(shù)據(jù)收集于中國科學(xué)院資源環(huán)境科學(xué)數(shù)據(jù)中心（http://www.resdc.cn/）。上述矢量和柵格數(shù)據(jù)均解譯完成。

1.3 研究方法

1.3.1 適宜性環(huán)境因子選取

本研究根據(jù)差異性、定量性原則，通過分析研究區(qū)地理和氣候資料，參考其他學(xué)者研究方法，選取地形、氣候、土壤和植被4 個尺度的生態(tài)環(huán)境因子作為陜北紅棗生態(tài)適宜性評價指標。首先，730～1 100 m被確定為適宜陜北紅棗生長的海拔高程，東南、南、西南方向被選定為陽坡，0～25°被選定為適宜坡度[18]。其次，9.6～11.10 ℃、3 000～3 200 ℃、430～580 mm依次被選定為適宜陜北紅棗生長的年平均氣溫、年積溫（≥10 ℃）、年降水量[18]。此外，黃綿土、黑壚土和灰鈣土被選定為適宜土壤類型[19]，溫帶落葉闊葉林被選定為適合陜北紅棗生長的地理類別[20]。其中，選擇溫帶落葉闊葉林作為植被因子，主要是因為，相比灌木和草本植物，它們具有樹木種類一致性；保護灌木和草本植物所表征的生態(tài)景觀格局多樣性的同時，并能夠減弱不同樹種間的生長競爭性[20]。

1.3.2 環(huán)境因子提取

利用ArcGIS10.6軟件將獲取的環(huán)境數(shù)據(jù)統(tǒng)一重采樣為30 m 分辨率，并獲得陜北DEM、氣象、土壤及土地利用柵格圖。通過該軟件的空間分析工具（Spatial Analyst）按照適宜紅棗生長的環(huán)境因子的閾值范圍，對紅棗的10 個適宜性環(huán)境因子進行空間化處理。在所有適宜性環(huán)境因子提取并柵格可視化后，依次通過“創(chuàng)建隨機點”和“多值提取至點”2 個步驟實現(xiàn)它們的點數(shù)據(jù)的獲取，以便后期統(tǒng)計分析使用。

選取極差法分別對數(shù)據(jù)進行標準化處理，將各個環(huán)境因子值域確定為[0-1]，具體計算式為：

正相關(guān)指標：

負相關(guān)指標：

式中：aij代表經(jīng)過標準化處理后第i個單元的第j個指標的值；xij代表第i個單元的第j個指標的原始值；xmax代表第i個單元的第j個指標的最大值；xmin代表第i個單元的第j個指標的最小值。

1.3.3 適宜性環(huán)境因子權(quán)重確定

采用層次分析法將適宜性環(huán)境因子按順序分解為不同層次結(jié)構(gòu)，采用1-9 標度評分標準賦值建立判斷矩陣，確定各個指標的權(quán)重（表1）。通過對其結(jié)果進行一致性檢驗（一致性比率≤0.1），證實判斷矩陣一致性好，指標權(quán)重結(jié)果合理。

表1 陜北紅棗生態(tài)適宜性的影響因子權(quán)重Table 1 Weight of influencing factors of ecological suitability of Jujube in Northern Shaanxi

1.3.4 空間疊加分析

基于各環(huán)境因子的權(quán)重大小，利用ArcGIS10.6軟件對上述提取的10 個環(huán)境因子進行空間綜合疊加分析，最終得出紅棗生態(tài)適宜性指數(shù)（S）并進行等級劃分，即優(yōu)適宜區(qū)（S≥90%）、次適宜區(qū)（90%＞S≥80%）、一般適宜區(qū)（80%＞S≥70%）及非適宜區(qū)（S＜70%）。紅棗生態(tài)適宜性指數(shù)的柵格圖可視化后，同樣通過“創(chuàng)建隨機點”和“多值提取至點”2 個步驟實現(xiàn)數(shù)值的獲取，用于地理探測和數(shù)據(jù)統(tǒng)計分析。

1.3.5 地理探測分析

地理探測器（GeoDetector）是一種能夠探測空間分異性的新型統(tǒng)計學(xué)方法[21]，主要從以下幾方面進行空間探測。

因子探測方法是探測各因子的解釋力大小，模型如式（3）：

式中：q值代表各環(huán)境因子對研究區(qū)增強型植被指數(shù)（EVI）的解釋力，值域為[0-1]，該值越大，說明該環(huán)境因子對EVI的解釋力越強；h=1，2，…，L為影響因子的分層數(shù)；Nh和N分別代表各環(huán)境因子的層h和整個研究區(qū)的樣本數(shù)；σh和σ分別代表層h和整個研究區(qū)的EVI值的方差。

生態(tài)探測方法是探究任意2 個環(huán)境因子x1和x2對紅棗生態(tài)適宜性指數(shù)的影響是否有顯著差異。2 種因子之間，如果行因子和列因子有顯著性差異，則為Y，沒有為N。用F統(tǒng)計量進行衡量：

式中：Nx1和Nx2分別代表環(huán)境因子x1和x2的樣本數(shù)量；SSWx1和SSWx2分別代表由2 個環(huán)境因子形成分層的層內(nèi)方差之和；L1和L2分別代表變量x1和x2分層數(shù)目。

交互作用探測方法用于識別各因子共同作用（增強或減弱）或相互獨立作用對紅棗生態(tài)適宜性的解釋力大小。在計算2 個環(huán)境因子交互的q值后比較q（x1）、q（x2）和q（x1∩x2）大小。

1.3.6 數(shù)據(jù)處理

采用ArcGIS 10.6 繪制紅棗適宜性區(qū)劃等級圖，采用地理探測器模型（Geographical detector model）和典型相關(guān)分析（CCA）探析環(huán)境因子對紅棗生態(tài)適宜性指數(shù)的影響特征。利用Origin 2022 繪制主要環(huán)境因子和紅棗生態(tài)適宜性指數(shù)之間的相關(guān)性，利用韋恩圖（Venn diagram）確定主要環(huán)境因子對紅棗生態(tài)適宜性指數(shù)的解釋率大小。

2 結(jié)果與分析

2.1 陜北紅棗生態(tài)適宜性區(qū)劃等級空間格局

圖2為陜北紅棗生態(tài)適宜性區(qū)劃等級的空間格局。由圖2可知，紅棗最優(yōu)適宜區(qū)和次適宜區(qū)面積小，均分布在陜北南部的富縣及鄰近區(qū)域；一般適宜區(qū)集中分布在陜北中南部沿黃河西岸區(qū)域，以綏德縣、清澗縣、延川縣及延長縣為代表。最優(yōu)適宜區(qū)、次適宜區(qū)、一般適宜區(qū)的面積分別達到8 548.02、50 732.19、887 126.76 hm2?？梢钥闯?，一般適宜區(qū)面積依次是最優(yōu)適宜區(qū)和次適宜區(qū)面積的103.80 倍和17.50 倍。充分說明，陜北紅棗生長主要表現(xiàn)為一般適宜區(qū)，而有利于紅棗生長的最優(yōu)生態(tài)適宜區(qū)較少。

圖2 陜北紅棗生態(tài)適宜性區(qū)劃等級的空間格局Fig.2 Spatial pattern of zoning levels of ecological suitability of Jujube in Northern Shaanxi

2.2 陜北紅棗生態(tài)適宜性影響因子的地理探測

2.2.1 因子探測分析

通過計算各環(huán)境因子的q值，提取各環(huán)境因子對紅棗生態(tài)適宜性現(xiàn)狀的影響力，探測結(jié)果如表2所示。各環(huán)境因子對紅棗生態(tài)適宜性的影響程度排序為：年積溫＞年平均氣溫＞年降水量＞高程＞坡度＞坡向＞其他（土壤因子和植被因子）。由此可見，年積溫、年平均氣溫、年降水量及高程的q值均大于0.50，解釋力均大于45%，是影響紅棗生態(tài)適宜性的主要環(huán)境因子；坡度和坡向的q值在0.20～0.40 區(qū)間，是紅棗生態(tài)適宜性的次要影響因子；土壤類型和植被因子的q值均小于0.10，對紅棗生態(tài)適宜性的影響較小。

表2 各影響因子探測q 值Table 2 q values of factor detection

2.2.2 生態(tài)探測分析

通過生態(tài)探測分析（表3）得出，氣候因子（年降水量、年積溫、年平均氣溫）與地形因子（高程、坡度、坡向）對紅棗生態(tài)適宜性變化的影響存在顯著性差異，與土壤類別（黃綿土、黑壚土、灰鈣土）、植被類型（溫帶落葉闊葉林帶）對紅棗適宜性的影響無顯著性差異。這充分證實了因子探測結(jié)果的正確性，即可以表明，在影響陜北紅棗的主要生態(tài)環(huán)境因子方面，氣候因子（年降水量、年積溫、年平均氣溫）對陜北紅棗生態(tài)適宜性的影響最為突出，且其影響程度明顯高于地形因子中的坡度和坡向。此外，土壤類別（黃綿土、黑壚土、灰鈣土）與植被類型（溫帶落葉闊葉林帶）對紅棗生態(tài)適宜性的影響也無顯著性差異。結(jié)合因子分析結(jié)果表明，土壤類別和植被類型對紅棗生態(tài)適宜性影響小且影響程度相近，因此，在后期的紅棗區(qū)域規(guī)模培育過程中不作為重點關(guān)注因素，而著重科學(xué)規(guī)劃氣候因子和地形因子中的高程。

表3 各環(huán)境因子影響紅棗生態(tài)適宜性的顯著性差異Table 3 Significant differences of environmental factors influencing ecological suitability of Jujube

2.2.3 交互探測分析

本文利用交互作用探測器對各環(huán)境因子進行探測，獲取各因子之間的相互關(guān)系（表4）。可以發(fā)現(xiàn)，各環(huán)境因子對紅棗生態(tài)適宜性的影響存在交互作用，不存在相互獨立起作用的因子（表5）。

由表4可發(fā)現(xiàn)，影響紅棗生態(tài)適宜性的主要和次要環(huán)境因子，即年降水量、高程、年平均氣溫、年積溫、坡度及坡向之間表現(xiàn)出極顯著的交互作用，而與其他環(huán)境因子間交互作用不顯著。由表5可知，絕大部分的環(huán)境因子交互q值大于其單因子的q值，因子互動效應(yīng)呈相互增強和非線性增強關(guān)系。其中，年降水量、年平均氣溫、高程、坡度、坡向與土壤類型的交互作用，均大于其因子間的簡單加和作用，呈現(xiàn)出對紅棗生態(tài)適宜性的非線性增強的影響效應(yīng)。而與其他環(huán)境因子的交互作用均小于因子加和作用，表現(xiàn)出相互增強的交互效應(yīng)。此外，年積溫、土壤類型、植被類型與其他環(huán)境因子的交互作用均小于其因子加和作用，說明它們對陜北紅棗生態(tài)適宜性的影響形式表現(xiàn)為相互增強的交互效應(yīng)。總之，各環(huán)境因子對紅棗生態(tài)適宜性的影響具有交互效應(yīng)，并非獨立存在的，也不是簡單的疊加，而是呈相互增強或非線性增強效應(yīng)。

表4 紅棗生態(tài)適宜性影響因子的交互作用Table 4 Interaction of factors influencing ecological suitability of Jujube

表5 影響紅棗生態(tài)適宜性的環(huán)境因子的交互效應(yīng)Table 5 Interaction effect of environmental factors influencing ecological suitability of Jujube

2.3 影響陜北紅棗生態(tài)適宜性的因子關(guān)系

本文在確定陜北地區(qū)紅棗生態(tài)適宜性區(qū)劃特征的基礎(chǔ)上，采用因子探測明晰了各環(huán)境因子對紅棗生態(tài)適宜性指數(shù)的影響差異特征。其中，年積溫、年平均氣溫、年降水量及高程的影響解釋率均大于45%，是調(diào)控紅棗生態(tài)適宜性的主要環(huán)境因子?；诘湫拖嚓P(guān)分析得出，年降水量、高程和年積溫極顯著影響紅棗生態(tài)適宜性指數(shù)，而其他環(huán)境因子的影響程度較?。▓D3）。這充分證實了因子探測結(jié)果的正確性，確定了降水量和溫度對紅棗種植的重要性。

圖3 環(huán)境因子對紅棗生態(tài)適宜性的影響Fig.3 Influence of environmental factors on ecological suitability of Jujube

通過相關(guān)性分析得出，年降水量、高程、年積溫及紅棗生態(tài)適宜性指數(shù)存在極顯著的相關(guān)性。其中，年降水量和年積溫與紅棗生態(tài)適宜性指數(shù)呈極顯著的非線性負相關(guān)關(guān)系，而高程和生態(tài)適宜性指標呈現(xiàn)極顯著的非線性正相關(guān)關(guān)系（圖4）。說明主要環(huán)境因子在不同數(shù)值范圍內(nèi)與紅棗的生態(tài)適宜性指數(shù)表現(xiàn)出不同的協(xié)同/權(quán)衡的動態(tài)變化關(guān)系。其中，在0～340 mm 和大于475 mm 的年降水量閾值內(nèi)，生態(tài)適宜性指數(shù)會隨年降水量的增加而增大，表現(xiàn)出正向協(xié)同關(guān)系，在340～475 mm 閾值內(nèi)，生態(tài)適宜性分布指數(shù)會隨年降水量的增加而減小，權(quán)衡關(guān)系顯著；年積溫在0～2 600 ℃和大于3 750 ℃閾值內(nèi)，生態(tài)適宜性指數(shù)隨著年積溫的增加而增大，呈正向協(xié)同關(guān)系，2 600～3 750 ℃閾值內(nèi)，生態(tài)適宜性指數(shù)隨著年積溫的增加而減小，權(quán)衡關(guān)系亦顯著；高程在0～900 m 和大于1 300 m 閾值內(nèi)，適宜性分布指數(shù)隨著高程的增加而減小，二者為權(quán)衡關(guān)系，900～1 300 m 閾值內(nèi)，生態(tài)適宜性指數(shù)隨著高程的增加而增大，因此表現(xiàn)為明顯的正向協(xié)同關(guān)系。整體上看，紅棗生態(tài)適宜性與年降水量和年積溫均表現(xiàn)為顯著的權(quán)衡關(guān)系，而與高程呈顯著的正向協(xié)同關(guān)系。此外，可以看出，利于紅棗生態(tài)適宜性分布和利于紅棗生長的主要環(huán)境因子（年降水量、高程及年積溫）的適宜閾值范圍不同。這主要是由于生態(tài)適宜性指數(shù)是由文中選定的10 個環(huán)境因子（地形、氣候、土壤及植被因子）綜合疊加得到的表征紅棗生態(tài)適宜性分布環(huán)境狀況和面積的指數(shù)，更具全面性和復(fù)雜性；同時，生態(tài)適宜性指數(shù)與紅棗生長是表征陜北紅棗種植培育狀況的2 類不同的影響途徑。因此，二者對主要環(huán)境因子的響應(yīng)節(jié)律表現(xiàn)出一定的關(guān)聯(lián)性和差異性。

圖4 年降水量、年積溫及高程與紅棗生態(tài)適宜性指數(shù)的關(guān)系Fig.4 Relationships between annual precipitation,annual accumulated temperature,elevation,and ecological suitability index of Jujube

進一步利用韋恩圖發(fā)現(xiàn)，年降水量、高程及年積溫對紅棗生態(tài)適宜性指數(shù)具有明顯的單一和交互解釋率。三者單一的解釋率分別達22%、15%、11%，且年降水量與高程，年降水量與年積溫，及高程與年積溫對紅棗生態(tài)適宜性指數(shù)的交互解釋率依次為11%、9%、7%。三者總共能夠解釋紅棗生態(tài)適宜性指數(shù)的78%（圖5），更加驗證了它們作為主要環(huán)境因子對紅棗生態(tài)適宜性的作用。

圖5 年降水量、年積溫和高程對紅棗生態(tài)適宜性指數(shù)的解釋率Fig.5 Interpretation rate of annual precipitation,annual accumulated temperature and elevation to ecological suitability index of Jujube

3 討論

3.1 陜北紅棗生態(tài)適宜性區(qū)劃特征

研究區(qū)紅棗生態(tài)適宜性區(qū)劃等級在空間上從東南部到西北部由高到低遞減。本研究表明，陜北紅棗最優(yōu)適宜區(qū)和次適宜區(qū)為陜北南部區(qū)域，主要以富縣為代表；富縣鄰邊區(qū)域分布較少，或只有小植株零星分布。二者分布面積分別為8 548.02、50 732.19 hm2。一般適宜區(qū)集中分布在陜北中南部沿黃河西岸區(qū)域，主要以綏德縣、清澗縣、延川縣及延長縣為代表，并向四周逐漸減少。一般適宜區(qū)面積顯著大于最優(yōu)適宜區(qū)和次適宜區(qū)，高達887 126.76 hm2，是陜北地區(qū)紅棗適宜性區(qū)劃的代表性特征。而陜北地區(qū)的北部（橫山區(qū)、榆陽區(qū)、神木市、府谷縣）和西部區(qū)縣（定邊縣、吳起縣、志丹縣、靖邊縣）是紅棗非適宜區(qū)。

以往紅棗種植區(qū)劃指標的選擇大多僅考慮了氣象因素，導(dǎo)致區(qū)劃結(jié)果粗略[22-23]。在分析紅棗與氣象要素關(guān)系的基礎(chǔ)上，綜合地形和土壤環(huán)境要素，以及植被的影響，實現(xiàn)了陜北地區(qū)紅棗種植生態(tài)適宜性的精細化區(qū)劃，有利于區(qū)域紅棗合理區(qū)劃和集約化經(jīng)營管理。通過實地調(diào)研查證，確定該研究的評價結(jié)果基本符合該地區(qū)紅棗種植的實際情況。優(yōu)適宜區(qū)和次適宜區(qū)，棗林分布分散，種植面積小，需要發(fā)揮地理優(yōu)勢，大力培育，增加規(guī)模；一般適宜區(qū)棗林種植面積大，但是紅棗品質(zhì)仍需要著重提升。

3.2 陜北紅棗生態(tài)適宜性的驅(qū)動機理

2000年以來，陜北地區(qū)生態(tài)系統(tǒng)格局變化顯著。主要表現(xiàn)為林木、草地的面積占比分別增加19%、7%。生態(tài)系統(tǒng)格局已由退耕還林（還草）生態(tài)工程實施前的耕地主導(dǎo)型轉(zhuǎn)化為林木主導(dǎo)型[11]。陜北棗樹作為該區(qū)域優(yōu)勢特色經(jīng)濟林木樹種之一，在區(qū)域生態(tài)系統(tǒng)調(diào)控方面發(fā)揮著重要的現(xiàn)實意義。本文通過因子探測分析和典型相關(guān)分析明晰了年積溫、年平均氣溫、年降水量及高程的影響解釋率均大于45%，對陜北紅棗生態(tài)適宜性指數(shù)具有極顯著的影響，是調(diào)控紅棗生態(tài)適宜性的主要環(huán)境因子。陜北地區(qū)的降水量和溫度在生態(tài)系統(tǒng)格局變化過程中波動顯著，相互制約，主導(dǎo)著區(qū)域氣象節(jié)律。1960—2017年，陜北地區(qū)年際地溫、氣溫顯著上升（P＜0.01），同時氣溫通過地-氣的熱量交換調(diào)控地溫變化[24]。這很大程度上因為退耕還林（還草）生態(tài)工程的實施。該地區(qū)造林引起微弱的增溫（（0.03±0.22）℃）和地表溫度季節(jié)性變化。春季及夏季的降溫程度較大，冬季降溫作用微弱。春季白天地表降溫最強，其次是夏季，秋冬季節(jié)降溫較弱[25]。地表溫度變化與降水梯度密切相關(guān)。年降水量、年平均氣溫增加趨勢的氣候傾向率分別達到 14.73 mm/（10 a）、0.25 ℃/（10 a）（P＞0.05）[2]。季節(jié)上，陜北地區(qū)夏季和秋季氣象干旱頻率高，春季和冬季氣象干旱頻率相對較低[26]。夏季極端降水約占夏季總降水量的54%，且極端降水占比在將近70%的站點表現(xiàn)出增長趨勢[27]。由此可見，陜北地區(qū)雖然具有植被和季節(jié)性降水量增加所導(dǎo)致的濕潤化氣候特征的發(fā)展趨勢，且有利于地表林木生長和恢復(fù)。但是由于氣溫仍然處于年際持續(xù)升高趨勢，所以該區(qū)域的氣候節(jié)律整體上仍然處于暖干化狀態(tài)[5]。因此，在陜北半干旱地區(qū)，主要以旱生植被和樹種為主。紅棗耐旱性較強，在熱量充足的條件下，降水對其生長的影響顯得更為重要。降水可以緩解水分脅迫，使得生長季末期延遲[28]，且降水和溫度對棗樹生長的影響占主導(dǎo)位置，這與劉靜等[29]和吉珍霞等[30]的觀點一致。故氣象的年際和季節(jié)性變化更替，對棗林的生長產(chǎn)生定期、不定量的影響。

溫度和降水的變化能夠顯著解釋林木覆蓋的變化。其中，溫度與陜北地區(qū)林木覆蓋度之間呈極顯著正相關(guān)的區(qū)域占84%，極顯著負相關(guān)區(qū)域占14%；降水與林木覆蓋度之間呈極顯著正相關(guān)的區(qū)域占54%，極顯著負相關(guān)區(qū)域占37%[31]。這一點與本研究結(jié)果相似，即年降水量和年積溫與生態(tài)適宜性指數(shù)極顯著負相關(guān)（權(quán)衡關(guān)系），而高程和生態(tài)適宜性指數(shù)極顯著正相關(guān)（正向協(xié)同關(guān)系）。此外，陜北地區(qū)林木對氣候變化的響應(yīng)存在明顯的空間差異性，不同氣候因子對不同林木覆蓋類型的影響程度不同[32—33]。這主要是因為陜北地區(qū)氣候特征本身具有顯著的區(qū)域差異性。陜北北部和中部黃河沿岸的府谷縣、神木縣、吳堡縣及綏德縣等區(qū)域，近年來降水量增多，尤其在紅棗成熟采摘期的集中降水現(xiàn)象的增加；加之黃河支流水域的區(qū)域微氣候特征的影響，使得紅棗出現(xiàn)裂果和棗果含水量增加的現(xiàn)象，直接減弱了紅棗的甜度和芳香風(fēng)味[20]。因此，紅棗的生態(tài)適宜性指數(shù)呈從陜北南部和東部區(qū)縣向西部和北部區(qū)縣遞減的空間布局差異性。結(jié)合各環(huán)境因子，進行地理探測中的生態(tài)探測分析，進一步得出，年降水量、年積溫、年平均氣溫與高程、坡度、坡向?qū)t棗生態(tài)適宜性的影響存在顯著性差異，與土壤類別（黃綿土、黑壚土、灰鈣土）、植被類型對紅棗生態(tài)適宜性的影響無顯著性差異。這與梁軼等[20]的研究結(jié)果一致，即基于陜北紅棗的喜光和耐旱的氣候特性，坡向和坡度也會顯著影響紅棗的產(chǎn)量和品質(zhì)，所以需要對坡向類型（陽坡、陰坡）進行優(yōu)化，對坡度適宜性進行人工修定。由此可知，在栽培種植過程中，需要穩(wěn)定造林面積，同時根據(jù)環(huán)境因子對紅棗的影響差異性和相關(guān)性，協(xié)同穩(wěn)定調(diào)控區(qū)域氣象變化幅度，以期實現(xiàn)棗林的平穩(wěn)生長。

進一步交互探測分析發(fā)現(xiàn)，年降水量、年平均氣溫、高程、坡度、坡向與土壤類型呈對紅棗適宜性的非線性增強的影響效應(yīng)。而與其他環(huán)境因子呈相互增強的交互效應(yīng)。這主要是因為地形、植被、氣候及土壤因子之間的關(guān)系本身具有時空復(fù)雜性和關(guān)聯(lián)性[5-6,8,10]。地形通過影響植被和氣候分布格局，來影響土壤類型，進而綜合性調(diào)控紅棗生態(tài)適宜性生長和分布區(qū)域及態(tài)勢。此外，年降水量、高程和年積溫對紅棗適宜性具有明顯的單一和交互解釋率。三者能夠整體性解釋紅棗生態(tài)適宜性指數(shù)的78%，足以證實它們對紅棗生態(tài)適宜性的顯著作用。因此，各環(huán)境因子交互作用對紅生態(tài)棗適宜性的影響是相互增強或非線性增強效應(yīng)，并非獨立存在和簡單疊加。基于環(huán)境因子和林木、植被的交互關(guān)系，生態(tài)環(huán)境的持續(xù)規(guī)?；G化，也會促進環(huán)境因子的良性循環(huán)。比如，2000年以來，植被恢復(fù)對陜北地區(qū)局地降水量的增加顯現(xiàn)了積極效應(yīng)。觀測數(shù)據(jù)顯示，2000—2015年該區(qū)域降水量以7.84 mm/a2的速率顯著增加，大規(guī)模植被恢復(fù)對研究區(qū)降水量增加速率的貢獻高達37.40%[34]。

陜北地區(qū)退耕還林還草工程實施效果顯著，工程實施后共有297 066.15 hm2耕地轉(zhuǎn)化為林地和草地，林草覆蓋率由57.30%增長至60.50%[35]。該區(qū)域生態(tài)系統(tǒng)格局已從耕地型轉(zhuǎn)化為林木和草地型[11,36]。充分說明，2000年以來，人為對植被、林地覆蓋的影響也向正作用方向發(fā)展。張翀等[37]利用地理探測器分析得出，土地利用變化就是影響林地、植被覆蓋變化的重要人為主導(dǎo)因子。據(jù)統(tǒng)計，1997—2018年，陜北地區(qū)未利用地面積占比減少5%[11]。在后續(xù)的土地利用規(guī)劃中，可以進一步調(diào)整區(qū)域土地利用模式和結(jié)構(gòu)，在適宜區(qū)開發(fā)未利用地，增加紅棗的規(guī)模栽培，促進土地利用結(jié)構(gòu)的生態(tài)轉(zhuǎn)變進程。這對于區(qū)域生態(tài)環(huán)境的持續(xù)改善具有現(xiàn)實意義。

綜上分析認為，對于陜北優(yōu)適宜區(qū)和次適宜區(qū)，由于實地考察證實，該2 類區(qū)域紅棗種植密度和面積小，需要提升生產(chǎn)培育面積和規(guī)模，引領(lǐng)優(yōu)質(zhì)紅棗發(fā)展形成優(yōu)勢突出和特色鮮明的紅棗產(chǎn)業(yè)帶；一般適宜區(qū)域，應(yīng)以培肥地力為重點，通過提升耕地肥力和質(zhì)量促進紅棗生長。非適宜區(qū)，應(yīng)及時調(diào)整農(nóng)業(yè)產(chǎn)業(yè)結(jié)構(gòu)，改種其他作物或退耕還林；亦或根據(jù)影響因子作用大小，針對性地進行溫室紅棗種植培育。通過主要影響因子的人工合理調(diào)控，并配合水肥等農(nóng)藝措施，達到陜北紅棗生長和品質(zhì)的改善。

4 結(jié)論

1）陜北紅棗生態(tài)適宜性區(qū)劃等級從東南部到西北部由高到低遞減。最優(yōu)適宜區(qū)和次適宜區(qū)集中分布在富縣，分布面積達8 548.02 hm2；一般適宜區(qū)分布在綏德縣、清澗縣及延川縣等陜北沿黃河西岸區(qū)域，面積顯著高于最優(yōu)適宜區(qū)和次適宜區(qū)。

2）氣候因子（年降水量、年積溫、年平均氣溫）和地形因子（高程、坡度、坡向）對紅棗生態(tài)適宜性影響較大，且存在顯著性差異。其中，年降水量、高程和年積溫（≥10 ℃）為主要驅(qū)動因子，三者對紅棗生態(tài)適宜性指數(shù)的解釋率分別達到22%、15%、11%，且總和解釋率高達78%。

3）各環(huán)境因子對紅棗生態(tài)適宜性的影響并非獨立存在的，也不是簡單的疊加，而是呈相互增強或非線性增強的交互作用效應(yīng)。其中，年降水量和年積溫與紅棗生態(tài)適宜性指數(shù)極顯著負相關(guān)（權(quán)衡關(guān)系），高程與生態(tài)適宜性指數(shù)極顯著正相關(guān)（正向協(xié)同關(guān)系）。